Ultra-low frequency (ULF) wave observations across the heliosphere often rely on the sign of correlations between plasma (density/pressure) and magnetic field perturbations to distinguish between fast and slow magnetosonic modes. However, the assumptions behind this magnetohydrodynamic result are not always valid, particularly within the magnetosphere which is inhomogeneous and supports standing waves along the geomagnetic field. Through theory and a global simulation, we find both effects can result in anomalous plasma–magnetic field correlations. The interference pattern in standing waves can lead both body and surface magnetosonic waves to have different cross-phases than their constituent propagating waves. Furthermore, if the scale of gradients in the background are shorter than the wavelength or the waves are near-incompressible, then advection by the wave of inhomogeneities can overcome the wave's inherent sense of compression. These effects need to be allowed for and taken into account when applying the typical diagnostic to observations.