Regularizations on the shape of the prisms. To obtain stabilized solutions, we impose the zeroth- and first-order Tikhonov The lowest values of the goal function form the set of candidate solutions. We run successive inversions for a range of tentative total-magnetization intensitiesĪnd depths to the top of the shallowest prism. The position of these vertices, the horizontal location of each prism, and the depth extent of all prisms are the parameters to be estimated by solving a constrained nonlinear inverse problem of minimizing a goal function. The horizontal cross-section of each prism is a polygon defined by a given number ofĮqui-angularly spaced vertices from 0º to 360º, whose polygon verticesĪre described by polar coordinates with an origin defined by a horizontal locationīecause our method estimates the radii of each polygon vertex we refer to it as We approximate the source by an ensemble of vertically juxtaposed right prisms, all of them with the same total-magnetization vector and depth extent. The total-magnetization direction is known. We present a method for inverting total-field anomaly data to estimate the geometry ofĪ uniformly magnetized 3-D geological source in the subsurface. The blue prisms represent the true model and the red prisms represent the estimated model. This method estimates the geometry of a 3D magnetic source from total field anomaly data.įigure 1: Result for complex model simulation.
#Pingas rangers code#
This repository contains the manuscript and supplementary code and data for the article "Magnetic radial inversion for 3-D source geometry estimation" submitted for publication in Geophysical Journal International.
Magnetic radial inversion for 3-D source geometry estimation