Frequency filtering of spatially modulated or &ldquotagged&rdquo MRI data in the spatial frequency k-space domain with subsequent 2DFT to the spatial domain and pixel-by-pixel arithmetic calculations provide robust ratio values that can be subjected to inverse trigonometric functions to derive B1 maps for an MRI system.