incorporates Gaussian smoothing of the observed anomalies with a frequency response of exp(-ν2/ω2) where ν is the frequency of the observed anomaly, and ω a frequency parameter which must be chosen with regard to the spacing of the observation points. The linear summation formula is modified to correct for irregularities in the distribution of observation points. The method has been developed for the interpolation of local variations in the deflcxtion components between astrogeodetic points.