1. A method for terrain correction for potential field geophysical data measured from at least one potential field sensor at at least one measurement position based on representation of the terrain response as a 3D Cauchy-type surface integral over the topography of the earth surface, the method comprising: a. measuring at least one component of gravity and/or magnetic scalar, vector, and/or tensor data with at least one potential field sensor in at least one measurement position along at least one survey line; b. measuring the terrain as a digital elevation model; c. selecting an appropriate physical property function to describe the physical property distribution within the terrain; d. forward model the terrain response from the digital elevation model with a physical property function for the at least one component of gravity and/or magnetic vector and/or tensor data with at least one potential field sensor in at least one measurement position along at least one survey line, using 3D Cauchy-type surface integral over the terrain surface; e. calculating the system terrain response for the at least one component of gravity and/or magnetic vector and/or tensor data with at least one potential field sensor in at least one measurement position along at least one survey line; and f. calculating the terrain corrected potential field data.