CN113655534A - Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition - Google Patents

Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition Download PDF

Info

Publication number
CN113655534A
CN113655534A CN202110794109.7A CN202110794109A CN113655534A CN 113655534 A CN113655534 A CN 113655534A CN 202110794109 A CN202110794109 A CN 202110794109A CN 113655534 A CN113655534 A CN 113655534A
Authority
CN
China
Prior art keywords
tensor
signal
decomposition
order
magnetic resonance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110794109.7A
Other languages
Chinese (zh)
Other versions
CN113655534B (en
Inventor
刘欢
王泽华
董浩斌
赵昌峰
王晓斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences
Original Assignee
China University of Geosciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences filed Critical China University of Geosciences
Priority to CN202110794109.7A priority Critical patent/CN113655534B/en
Publication of CN113655534A publication Critical patent/CN113655534A/en
Application granted granted Critical
Publication of CN113655534B publication Critical patent/CN113655534B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

The invention provides a nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition, which comprises the following steps: acquiring a multichannel signal: enabling one channel to start sampling at intervals of delta t through time delay sampling to obtain a multi-channel signal; converting each channel signal into Hankel matrix to form third-order tensor
Figure DDA0003162185570000011
To third order tensor
Figure DDA0003162185570000012
Performing Tucker decomposition and processing to obtain third order tensor
Figure DDA0003162185570000013
Tensor of third order
Figure DDA0003162185570000014
Recovering the signal into a multi-channel signal, and performing CP tensor decomposition processing on a second-order tensor X formed by the multi-channel signal to finally obtain a high signal-to-noise ratio signal X after multi-channel signal fusionnew. The invention has the beneficial effects that: the method effectively overcomes the limitation of the current algorithm under strong noise interference, simultaneously ensures the real-time property and universality of the algorithm, is applicable to instruments such as a proton magnetometer and a nuclear magnetic resonance water finder, improves the precision of frequency domain measurement and the accuracy of hydrological parameters, and effectively improves the precision and the accuracy of related geoscience instruments.

Description

Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition
Technical Field
The invention relates to the field of signal processing, in particular to a nuclear magnetic resonance (FID) signal noise suppression method and method based on multi-linear singular value tensor decomposition.
Background
Free Induction Decay (FID) signals are signals generated in nuclear magnetic resonance technology, and in related instruments (such as proton precession magnetometers, nuclear magnetic resonance water detectors and the like) in the field of geophysical exploration by using the nuclear magnetic resonance technology, required geomagnetic field values, hydrological parameters and the like are obtained by means of extracting parameters of the nuclear magnetic resonance FID signals, measuring frequencies, inverting hydrological parameters and the like, but the obtained FID signals are very weak, and the FID signals are easily submerged in environmental noise and artificial noise. How to suppress various types of noise becomes a major limitation of nmr techniques in related geosciences. Therefore, in order to improve the measurement accuracy and accuracy of the geoscience instrument, it is necessary to perform noise suppression on the acquired FID signal to obtain more accurate information such as frequency, hydrological parameters, and the like.
The FID signal noise suppression means in the current stage mainly performs noise suppression through algorithms, and commonly used algorithms include SVD, wavelet transform, PCA, EMD, and the like, but the algorithms still have the following limitations in practical application: 1) the requirement on the data volume is high, and the noise suppression effect is greatly reduced under the condition of few sampling points; 2) the analysis processing is carried out in a time-frequency domain, the calculation amount is large, and the real-time performance cannot be guaranteed;
3) the universality is weak, and the algorithms perform noise suppression on specific types of noise, but in a practical situation, a plurality of types of noise exist simultaneously. Therefore, when the obtained FID signal is noisy and has more noise types, these algorithms have not been able to perform effective noise suppression, so that the instrument cannot work normally.
Disclosure of Invention
In view of the above technical defects, the present invention provides a method for suppressing noise of nuclear magnetic resonance FID signals based on nonlinear singular value tensor decomposition, including the following steps:
s101: acquiring a multichannel signal: enabling one channel to start sampling at intervals of delta t through time delay sampling to obtain a multi-channel signal;
s102: noise suppression based on multilinear singular value tensor decomposition:
converting each channel signal into Hankel matrix to form third-order tensor
Figure BDA0003162185550000021
To third order tensor
Figure BDA0003162185550000022
Performing Tucker decomposition and processing to obtain third order tensor
Figure BDA0003162185550000023
Tensor of third order
Figure BDA0003162185550000024
Recovering the signal into a multi-channel signal, and performing CP tensor decomposition processing on a second-order tensor X formed by the multi-channel signal to finally obtain a high signal-to-noise ratio signal X after multi-channel signal fusionnew。
Further, the third order tensor obtained in step S102
Figure BDA0003162185550000025
The specific process of performing the Tucker decomposition treatment comprises the following steps:
to tensor
Figure BDA0003162185550000026
Performing Tucker decomposition to obtain a core tensor
Figure BDA0003162185550000027
With a plurality of adjoint matrices U(n)；
To core tensor
Figure BDA0003162185550000028
Carrying out zero setting processing on singular values with relatively small median values to obtain a processed core tensor
Figure BDA0003162185550000029
The processed core tensor
Figure BDA00031621855500000210
A plurality of adjoint matrixes U obtained in the last step(n)Carrying out tensor modular multiplication operation to obtain third-order tensor
Figure BDA00031621855500000211
I.e. a part of the noise suppression is done.
Further, for the third order tensor
Figure BDA00031621855500000212
The restoration is a multichannel signal X, and the specific process of CP tensor decomposition is as follows:
tensor of third order
Figure BDA00031621855500000213
Transforming into a multi-channel signal matrix X (namely a second-order tensor) through a Hankel inverter; and decomposing the second-order tensor X into a form of two groups of one-dimensional vector outer products by using CP tensor decomposition to obtain a vector a and a vector b, and calculating the vector a and the vector b to obtain a fused denoising FID signal so as to complete all noise suppression.
Further, the operation process of the a, b vector is as follows:
Figure BDA00031621855500000214
wherein mean (a) x b represents that the mean value of the vector a is multiplied by the vector b to obtain a fused signal xnew。
The invention has the beneficial effects that: the method effectively overcomes the limitation of the current algorithm under strong noise interference, simultaneously ensures the real-time property and universality of the algorithm, is applicable to instruments such as a proton magnetometer and a nuclear magnetic resonance water finder, improves the precision of frequency domain measurement and the accuracy of hydrological parameters, and effectively improves the precision and the accuracy of related geoscience instruments.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of a multi-channel equal delay sampling method;
fig. 3 is a comparison of the effect of the FID signal measured in the extreme environment after noise suppression by SVD, PCA, and multi-linear singular value tensor decomposition (MLSVD).
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be further described with reference to the accompanying drawings.
Referring to fig. 1, a method for suppressing noise of a nuclear magnetic resonance FID signal based on nonlinear singular value tensor decomposition includes the following steps:
s101: acquiring a multichannel signal: enabling one channel to start sampling at intervals of delta t through time delay sampling to obtain a multi-channel signal;
in the embodiment of the invention, a multi-channel equal delay sampling method is adopted to obtain a plurality of related FID signals, one channel is enabled to start sampling every time delta t through delay sampling, and the sampling time is equal, so that a plurality of signals related to the original signals can be obtained.
For example, when the number of channels is 5, 5 FID signals x are obtained1,x2,x3,x4,x5A schematic diagram of the multi-channel equal delay sampling method is shown in fig. 2.
The method actually yields 5 vectors x1,x2,x3,x4,x5There are different signal and noise components:
Figure BDA0003162185550000041
wherein xsRepresenting a signal component, xnRepresenting the noise component. The correlation among signal components of multi-channel data acquired by a multi-channel equal delay sampling method is strong, the correlation among noise components is low, and compared with multi-channel sampling, the correlation of noise data is reduced.
S102: noise suppression based on multilinear singular value tensor decomposition:
converting each channel signal into Hankel matrix to form third-order tensor
Figure BDA00031621855500000411
To third order tensor
Figure BDA00031621855500000412
And performing Tucker decomposition processing and Hankel inverse transformation to obtain a multichannel signal X, performing CP tensor decomposition on the X to complete signal fusion, and completing noise suppression.
In the invention, after acquiring multichannel data, the multichannel data is constructed into a third-order tensor, and the noise processing is carried out by utilizing the decomposition of a multilinear singular value tensor:
1) and transforming the signals of each channel into a Hankel matrix, wherein a three-dimensional array formed by a plurality of Hankel matrices is a third-order tensor. After the third-order tensor is constructed, the tensor needs to be subjected to multi-linear singular value tensor decomposition. Assuming a third order tensor is derived:
Figure BDA0003162185550000042
2) to tensor
Figure BDA0003162185550000043
Performing Tucker decomposition: tensor
Figure BDA0003162185550000044
Decomposed into a core tensor
Figure BDA0003162185550000045
With N adjoint matrices
Figure BDA0003162185550000046
The result of successive modulo n multiplications.
3) Tensor
Figure BDA0003162185550000047
Obtaining a core tensor after Tucker decomposition
Figure BDA0003162185550000048
The values also represent the signal and noise components of the data. The smaller value representing noise is set to zero and restored to new tensor
Figure BDA0003162185550000049
The effect of noise processing can be obtained.
4) Will tensor
Figure BDA00031621855500000410
Obtaining a second-order tensor (matrix) X through Hankel inverse transformation, decomposing the second-order tensor (matrix) X through a CP tensor to obtain two vectors a and b, and calculating the vectors a and b:
Figure BDA0003162185550000051
mean (a) x b denotes multiplying the mean of vector a by vector b to obtain the fused signal xnewNamely the fused FID signal.
The present invention provides an embodiment as follows:
fig. 3 is a comparison of the effect of the FID signal measured in the extreme environment after noise suppression by SVD, PCA and multi-linear singular value tensor decomposition (MLSVD), and it can be seen from the time domain diagram (SVD is the uppermost and lowermost part on the left side of fig. 3, MLSVD is the middle trapezoidal part on the left side of fig. 3, and PCA is the part between the two) and the frequency domain diagram that the noise reduction effect of the FID signal noise suppression algorithm based on the multi-linear singular value decomposition is obvious, and the useful information of the FID signal can still be recovered under the condition that other algorithms cannot work. The signal envelope is clearer and the attenuation characteristic is more obvious in the time domain; the frequency spectrum curve is smoother and more approximate to the real FID signal frequency spectrum in the frequency domain.
The noise suppression effects of the SVD, PCA and the multi-linear singular value tensor decomposition algorithms at different sampling points are sorted, the comparison results are shown in Table 1, and it can be seen that the multi-linear singular value tensor decomposition algorithms have obvious noise suppression effects compared with other algorithms at different sampling points, and when the noise is large, the SVD and the PCA cannot work, and the multi-linear singular value tensor decomposition algorithms can still effectively perform noise suppression.
Table one: and decomposing the noise reduction effect of the SVD, the PCA and the multi-linear singular value tensor under different conditions.
Figure BDA0003162185550000052
Figure BDA0003162185550000061
Compared with the existing FID signal noise suppression algorithm, the method has the following characteristics:
(1) a multi-channel signal acquisition method based on time delay is adopted to acquire more signal noise characteristic information;
(2) the noise singular value is calculated and processed by using the multi-linear singular value tensor decomposition, and the real-time performance of the algorithm is ensured while the noise suppression effect is ensured;
(3) information fusion of multi-channel signals is realized through CP tensor decomposition, and the noise suppression effect of the invention is further improved.
The invention provides a nuclear magnetic resonance signal noise suppression method based on multi-linear singular value tensor decomposition, which effectively overcomes the limitation of the current algorithm under strong noise interference, ensures the real-time property and universality of the algorithm, is suitable for instruments such as a proton magnetometer and a nuclear magnetic resonance water finder, improves the precision of frequency domain measurement and the accuracy of hydrological parameters, and effectively improves the precision and the accuracy of related geoscience instruments.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (4)

1. A nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition is characterized by comprising the following steps: the method comprises the following steps:
s101: acquiring a multichannel signal: enabling one channel to start sampling at intervals of delta t through time delay sampling to obtain a multi-channel signal;
s102: noise suppression based on multilinear singular value tensor decomposition:
converting each channel signal into Hankel matrix to form third-order tensor
Figure FDA0003162185540000011
To third order tensor
Figure FDA0003162185540000012
Performing Tucker decomposition and processing to obtain third order tensor
Figure FDA0003162185540000013
Tensor of third order
Figure FDA0003162185540000014
Recovering to multi-channel signals, and performing CP tensor decomposition processing on a second-order tensor X formed by the multi-channel signals to finally obtain the height of the multi-channel signals after fusionSignal to noise ratio signal xnew。
2. The method for suppressing noise of nuclear magnetic resonance (FID) signals based on the nonlinear singular value tensor decomposition as claimed in claim 1, wherein the method comprises the following steps: the third order tensor obtained in step S102
Figure FDA0003162185540000015
The specific process of performing the Tucker decomposition treatment comprises the following steps:
to tensor
Figure FDA0003162185540000016
Performing Tucker decomposition to obtain a core tensor
Figure FDA0003162185540000017
With a plurality of adjoint matrices U(n)；
To core tensor
Figure FDA0003162185540000018
Carrying out zero setting processing on singular values with relatively small median values to obtain a processed core tensor
Figure FDA0003162185540000019
The processed core tensor
Figure FDA00031621855400000110
A plurality of adjoint matrixes U obtained in the last step(n)Carrying out tensor modular multiplication operation to obtain third-order tensor
Figure FDA00031621855400000111
I.e. a part of the noise suppression is done.
3. The method for suppressing noise of nuclear magnetic resonance (FID) signals based on the nonlinear singular value tensor decomposition as claimed in claim 1, wherein the method comprises the following steps: to third order tensor
Figure FDA00031621855400000112
The restoration is a multichannel signal X, and the specific process of CP tensor decomposition is as follows:
tensor of third order
Figure FDA00031621855400000113
Transforming into a multi-channel signal matrix X (namely a second-order tensor) through a Hankel inverter; and decomposing the second-order tensor X into a form of two groups of one-dimensional vector outer products by using CP tensor decomposition to obtain a vector a and a vector b, and calculating the vector a and the vector b to obtain a fused denoising FID signal so as to complete all noise suppression.
4. The method for suppressing noise of nuclear magnetic resonance (FID) signals based on the nonlinear singular value tensor decomposition as claimed in claim 3, wherein: the operation process of the a and b vectors is as follows:
Figure FDA0003162185540000021
wherein mean (a) x b represents that the mean value of the vector a is multiplied by the vector b to obtain a fused signal xnew。
CN202110794109.7A 2021-07-14 2021-07-14 Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition Active CN113655534B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110794109.7A CN113655534B (en) 2021-07-14 2021-07-14 Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110794109.7A CN113655534B (en) 2021-07-14 2021-07-14 Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition

Publications (2)

Publication Number Publication Date
CN113655534A true CN113655534A (en) 2021-11-16
CN113655534B CN113655534B (en) 2022-05-17

Family

ID=78477383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110794109.7A Active CN113655534B (en) 2021-07-14 2021-07-14 Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition

Country Status (1)

Country Link
CN (1) CN113655534B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114235411A (en) * 2021-12-28 2022-03-25 频率探索智能科技江苏有限公司 Bearing outer ring defect positioning method

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070018644A1 (en) * 2003-01-30 2007-01-25 Flexman John H Signal processing for detection of nqr signals
US7253627B1 (en) * 2006-07-19 2007-08-07 Univ King Fahd Pet & Minerals Method for removing noise from nuclear magnetic resonance signals and images
CN106772648A (en) * 2016-12-05 2017-05-31 中国地质大学（武汉） Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio
CN110333545A (en) * 2019-06-28 2019-10-15 中国地质大学（武汉） Improve the system and method for rotary proton class sensor tuning precision and signal-to-noise ratio
US20200049782A1 (en) * 2017-02-14 2020-02-13 The United States of America, as represented by the Secretary, Department of Health and Human Ser.. Denoising of dynamic magnetic resonance spectroscopic imaging using low rank approximations in the kinetic domain
CN110826017A (en) * 2019-09-25 2020-02-21 中国地质大学（武汉） Signal denoising method based on parameter optimization Hankel matrix and singular value decomposition
CN111257929A (en) * 2020-02-17 2020-06-09 成都理工大学 Singular value attenuation rank reduction denoising method
CN111929733A (en) * 2020-08-21 2020-11-13 电子科技大学 Seismic signal regularization processing method based on slice sampling
CN111985501A (en) * 2020-07-30 2020-11-24 河海大学 Hyperspectral image feature extraction method based on self-adaptive high-order tensor decomposition

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070018644A1 (en) * 2003-01-30 2007-01-25 Flexman John H Signal processing for detection of nqr signals
US7253627B1 (en) * 2006-07-19 2007-08-07 Univ King Fahd Pet & Minerals Method for removing noise from nuclear magnetic resonance signals and images
CN106772648A (en) * 2016-12-05 2017-05-31 中国地质大学（武汉） Improve the method and system of rotary proton class sensor tuning precision and signal to noise ratio
US20200049782A1 (en) * 2017-02-14 2020-02-13 The United States of America, as represented by the Secretary, Department of Health and Human Ser.. Denoising of dynamic magnetic resonance spectroscopic imaging using low rank approximations in the kinetic domain
CN110333545A (en) * 2019-06-28 2019-10-15 中国地质大学（武汉） Improve the system and method for rotary proton class sensor tuning precision and signal-to-noise ratio
CN110826017A (en) * 2019-09-25 2020-02-21 中国地质大学（武汉） Signal denoising method based on parameter optimization Hankel matrix and singular value decomposition
CN111257929A (en) * 2020-02-17 2020-06-09 成都理工大学 Singular value attenuation rank reduction denoising method
CN111985501A (en) * 2020-07-30 2020-11-24 河海大学 Hyperspectral image feature extraction method based on self-adaptive high-order tensor decomposition
CN111929733A (en) * 2020-08-21 2020-11-13 电子科技大学 Seismic signal regularization processing method based on slice sampling

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HUAN LIU,等: "Improving the Signal-to-Noise-Ratio of Free Induction Decay Signals Using a New Multilinear Singular Value Decomposition-Based Filter", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 *
TIANYU QIU，等: "Review and Prospect: NMR Spectroscopy Denoising & Reconstruction with Low Rank Hankel Matrices and Tensors", 《MAGNETIC RESONANCE IN CHEMISTRY》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114235411A (en) * 2021-12-28 2022-03-25 频率探索智能科技江苏有限公司 Bearing outer ring defect positioning method
CN114235411B (en) * 2021-12-28 2023-06-30 频率探索智能科技江苏有限公司 Bearing outer ring defect positioning method

Also Published As

Publication number Publication date
CN113655534B (en) 2022-05-17

Similar Documents

Publication Publication Date Title
JP2735064B2 (en) Waveform analyzer
US20100157731A1 (en) Automatic dispersion extraction of multiple time overlapped acoustic signals
CN109885903B (en) Model-based ground nuclear magnetic resonance signal spike noise removing method
CN112083508B (en) Artificial source electromagnetic exploration signal denoising method and system with noise reference channel
CN112487604B (en) Long-time nonlinear drift compensation method for output data of marine gravimeter
CN113655534B (en) Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition
CN114091538B (en) Intelligent noise reduction method for discrimination loss convolutional neural network based on signal characteristics
Chen et al. Improved methods for exponential parameter estimation in the presence of known poles and noise
CN109165432B (en) Magnetic resonance spectrum reconstruction method based on partial singular value sum
Zhou et al. Geomagnetic sensor noise reduction for improving calibration compensation accuracy based on improved HHT algorithm
LIN et al. Segmented time-frequency peak filtering for random noise reduction of MRS oscillating signal
CN113640891B (en) Singular spectrum analysis-based transient electromagnetic detection data noise filtering method
Zala et al. High-resolution signal and noise field estimation using the L 1 (least absolute values) norm
CN112285793B (en) Magnetotelluric denoising method and system
Baziw Incorporation of iterative forward modeling into the principle phase decomposition algorithm for accurate source wave and reflection series estimation
Xu et al. Multichannel reflectivity inversion with sparse group regularization based on HPPSG algorithm
Cai et al. De-noising for NMR oil well logging signals based on empirical mode decomposition and independent component analysis
CN112817056B (en) Magnetotelluric signal denoising method and system
CN111008356B (en) Gamma energy spectrum set analysis method for deducting background based on WTS VD algorithm
US5873050A (en) Method of analyzing capabilities of multiple-suppression computer seismic data processing software
CN117078785B (en) Rapid non-Cartesian magnetic resonance intelligent imaging method
CN114972562B (en) Fast magnetic resonance imaging method combining coil sensitivity estimation and image reconstruction
CN112505779B (en) Method for removing collected footprints based on feature decomposition
JP2022190994A (en) Magnetic resonance imaging apparatus, noise removal method and image processing apparatus
CN116821593A (en) Nuclear magnetic resonance data noise suppression method based on iterative Kalman filtering

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant