Dear Martin,
This may not be exactly what you want, but you don't really say what your
intended application is, so it's hard to tell. But to go from potential to field
component you normally use a directional derivative. So you need a directional
integral to go the other way.
If you pole reduce and then do the pseudogravity transform, that's what you've
done. Of course that's given you the potential of the pole-reduced field, and
maybe you will need to fiddle with the scaling, but maybe that's good enough for
you?
There's plenty of commercial software out there to do that (GETgrid, Geosoft and
Intrepid all spring to mind and all offer academic pricing). If you want the
actual algorithms, they're all Fourier filters of course. They are nicely
decribed in enough detail to implement (I have done so) in Gunn, P J, 1975.
Linear transformations of gravity and magnetic fields. Geophys. Prosp., v. 23,
p. 300-312.
Regards
Alan Reid
mail wrote:
> Dear group,
>
> Can anyone recommend an algorithm for the transformation of gridded magnetic
> field anomalies into their magnetic potential equivalents? I've checked
> Blakely, 1996, but although the basic principle is described there doesn't
> seem to be an actual algorithm.
>
> Thanks,
>
> Martin Roseveare
> ______________________________________
>
> ArchaeoPhysica Ltd
> Reconnaissance & Geophysics for Archaeology
>
> Phone: +44 (0) 7050 369789
> E-mail: [log in to unmask]
> Web: www.archaeophysica.co.uk
--
Alan B Reid PhD
Reid Geophysics
49 Carr Bridge Drive
Leeds LS16 7LB
UK
Phone: +(44)-113-261-0252
Fax: +(44)-113-293-0062
e-mail: mailto:[log in to unmask]
website: http://www.reid-geophys.co.uk
|