Skip to content

Commit

Permalink
Version submitted to G3
Browse files Browse the repository at this point in the history
  • Loading branch information
leouieda committed Jun 12, 2023
1 parent a8b0cee commit 629ce77
Show file tree
Hide file tree
Showing 4 changed files with 21 additions and 20 deletions.
4 changes: 2 additions & 2 deletions paper/abstract-plain.tex
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
Very small magnetic particles in rocks and other materials can store information about what the Earth’s magnetic field was like in the past.
\noindent{Very small magnetic particles in rocks and other materials can store information about what the Earth’s magnetic field was like in the past.
But not all particles are good recorders of this magnetic information, and some may have recorded different overlapping directions and strengths.
So it is important to measure each particle separately in order to identify and separate the good recorders from the bad ones.
A device called a ``quantum diamond microscope'' is able to measure the magnetic field near the surface of a rock sample at microscopic scale.
We propose a new method for processing data from this microscope that is able to find out the individual magnetizations of large amounts of small magnetic particles automatically.
We created a computer program to execute the method, which calculates the 3D position and magnetization of each particle using the simple model of a magnetic dipole.
We tested the method on simulated data, using fake magnetic particles for which we know the correct magnetization and position, and real data, both of which showed good results in most cases.
The method we created has the potential to enable the widespread study of the magnetism of natural materials with more detail than before.
The method we created has the potential to enable the widespread study of the magnetism of natural materials with more detail than before.}
27 changes: 14 additions & 13 deletions paper/abstract.tex
Original file line number Diff line number Diff line change
@@ -1,13 +1,14 @@
Paleomagnetic data are usually obtained from whole cylindric samples, where the signal results from the sum of magnetic moments from hundreds of thousands to millions of magnetic particles within the sample volume.
This usually includes both stable and unstable remanence carriers.
Recently, magnetic microscopy techniques allowed the investigation of individual grains by directly imaging their magnetic field.
However, the determination of the magnetic moments of individual grains is hindered by the intrinsic ambiguity in the inversion of potential field data, as well as by the large number of grains found in any one microscopy image.
We present a fast, semi-automated algorithm capable of estimating the position and magnetization of each ferromagnetic (l.s) source using only the magnetic microscopy data.
Our algorithm works in three steps: (i) we first apply image processing techniques to identify and isolate data window boundaries for each source; (ii) with these window boundaries, the position of the sources is estimated using Euler deconvolution; and finally (iii) using the position information, the algorithm is able to estimate the magnetic dipole moment direction and intensity for each source through an overdetermined linear inverse problem using a dipolar approximation.
The method does not require any type of additional information about the sample or the sources.
Sensitivity tests were run to estimate the stability of our routine to the depth of particles, signal-to-noise ratio, and non-dipolarity of the sources.
Tests with simple synthetic data show the high effectiveness of the methodology for recovering the position and magnetic information for both dipolar and non-dipolar sources.
More complex synthetic data including over 100 different magnetic particles were devised to emulate real rock data.
Results obtained on these data also show the feasibility and robustness of the algorithm to semi-automatically estimate the position and magnetic moment of a large number of particles.
This is further confirmed through an application to real data in which we are able to retrieve the expected bimodal isothermal remanent directions that were induced in the sample.
Given its semi-automatic nature, its low processing cost, and the possibility of simultaneous inversion of the magnetic moment of a great number of magnetic particles, the methodology here proposed is a step forward in enabling paleomagnetic applications of magnetic microscopy.
% Paleomagnetic data are usually obtained from whole cylindric samples, where the signal results from the sum of magnetic moments from hundreds of thousands to millions of magnetic particles within the sample volume.
% This usually includes both stable and unstable remanence carriers.
% Recently, magnetic microscopy techniques allowed the investigation of individual grains by directly imaging their magnetic field.
% However, the determination of the magnetic moments of individual grains is hindered by the intrinsic ambiguity in the inversion of potential field data, as well as by the large number of grains found in any one microscopy image.
% We present a fast, semi-automated algorithm capable of estimating the position and magnetization of each ferromagnetic (l.s) source using only the magnetic microscopy data.
% Our algorithm works in three steps: (i) we first apply image processing techniques to identify and isolate data window boundaries for each source; (ii) with these window boundaries, the position of the sources is estimated using Euler deconvolution; and finally (iii) using the position information, the algorithm is able to estimate the magnetic dipole moment direction and intensity for each source through an overdetermined linear inverse problem using a dipolar approximation.
% The method does not require any type of additional information about the sample or the sources.
% Sensitivity tests were run to estimate the stability of our routine to the depth of particles, signal-to-noise ratio, and non-dipolarity of the sources.
% Tests with simple synthetic data show the high effectiveness of the methodology for recovering the position and magnetic information for both dipolar and non-dipolar sources.
% More complex synthetic data including over 100 different magnetic particles were devised to emulate real rock data.
% Results obtained on these data also show the feasibility and robustness of the algorithm to semi-automatically estimate the position and magnetic moment of a large number of particles.
% This is further confirmed through an application to real data in which we are able to retrieve the expected bimodal isothermal remanent directions that were induced in the sample.
% Given its semi-automatic nature, its low processing cost, and the possibility of simultaneous inversion of the magnetic moment of a great number of magnetic particles, the methodology here proposed is a step forward in enabling paleomagnetic applications of magnetic microscopy.
Paleomagnetic data is collected from bulk samples, containing a mixture of stable and unstable magnetic particles. Recently, magnetic microscopy techniques have allowed the examination of individual magnetic grains. However, accurately determining the magnetic moments of these grains is difficult and time-consuming due to the inherent ambiguity of potential field data and the large number of grains in each image. Here we introduce a fast and semi-automated algorithm that estimates the position and magnetization of each magnetic source solely based on magnetic microscopy data, without requiring any additional information. The algorithm follows a three-step process. Firstly, it employs image processing techniques to identify and isolate data window boundaries for each magnetic source. Secondly, the algorithm uses Euler deconvolution to estimate their positioning. Lastly, the algorithm solves an overdetermined linear inverse problem, using a dipolar approximation, to determine the direction and intensity of the magnetic dipole moment for each source. To validate the algorithm, we conducted tests analyzing particle depth, intensity, and non-dipolarity. Additionally, synthetic data experiments were performed to demonstrate the effectiveness of the methodology in recovering the position and magnetic information of both dipolar and non-dipolar sources. These experiments included complex samples designed to simulate real rock data. Then, we applied it to real data, where it accurately retrieved the expected directions induced in the sample. The semi-automated nature of our algorithm, combined with its low processing cost and ability to simultaneously analyze the magnetic moments of numerous particles, represents a significant advancement in facilitating paleomagnetic applications of magnetic microscopy.
2 changes: 1 addition & 1 deletion paper/content.tex
Original file line number Diff line number Diff line change
Expand Up @@ -711,7 +711,7 @@ \subsection{Applicability to real-world scenarios}
\centering
\includegraphics[width=1\linewidth]{figures/complex-synthetic-stereograms.png}
\caption{ Comparison of true and estimated dipole magnetic moments and directions for the complex synthetic sample.
a) Simulation of a thin section of rock with 103 particles uniformly magnetized (dipolar sources) by two different induced fields. The average directions of the induced fields are $D=\ang{0}$ and $I=\ang{0}$ and $D=\ang{180}$ and $I=\ang{0}$, yielding two stable directions. b) All estimated vector directions of the identified sources ($M=99$). c) The Filtering criteria used to select the magnetic directions of the sources with the best fitting ($M=96$), determined by the coefficient of determination ($R^2 \geq 0.85$) and signal-to-noise ratio ($\text{SNR} \geq 5$).
a) Simulation of a thin section of rock with 103 particles uniformly magnetized (dipolar sources) by two different induced fields. The average directions of the induced fields are $D=\ang{0}$ / $I=\ang{0}$ and $D=\ang{180}$ / $I=\ang{0}$, yielding two stable directions. b) All estimated vector directions of the identified sources ($M=99$). c) The Filtering criteria used to select the magnetic directions of the sources with the best fitting ($M=96$), determined by the coefficient of determination ($R^2 \geq 0.85$) and signal-to-noise ratio ($\text{SNR} \geq 5$).
}
\label{complex-synthetic-stereograms}
\end{figure}
Expand Down
8 changes: 4 additions & 4 deletions paper/manuscript-agu.tex
Original file line number Diff line number Diff line change
Expand Up @@ -60,7 +60,7 @@


\begin{document}

\justifying

% Import the title, etc.
\input{information.tex}
Expand Down Expand Up @@ -145,9 +145,9 @@
% \end{keypoints}

\begin{keypoints}
\item enter point 1 here
\item enter point 2 here
\item enter point 3 here
\item We present a fast algorithm to estimate the position and magnetization of individual sources using only magnetic microscopy data.
\item The algorithm was validated using synthetic samples that included both dipolar and non-dipolar sources with varying depths and intensities.
\item The successful application of the algorithm on a speleothem sample represents an advancement of this methodology for paleomagnetic studies.
\end{keypoints}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Expand Down

0 comments on commit 629ce77

Please sign in to comment.