Method for removing effects of shallow resistivity structures in electromagnetic survey data to produce a low frequency resistivity anomaly map, or alternatively imaging resistivity structures at their correct depth levels. The method involves solving Maxwell's electromagnetic field equations by either forward modeling or inversion, and requires at least two survey data sets, one taken at the source frequency selected to penetrate to a target depth, the other a higher frequency able to penetrate only shallow depths.