CN104898172A - 一种基于互相关的核磁共振全波信号噪声滤除方法 - Google Patents
一种基于互相关的核磁共振全波信号噪声滤除方法 Download PDFInfo
- Publication number
- CN104898172A CN104898172A CN201510256496.3A CN201510256496A CN104898172A CN 104898172 A CN104898172 A CN 104898172A CN 201510256496 A CN201510256496 A CN 201510256496A CN 104898172 A CN104898172 A CN 104898172A
- Authority
- CN
- China
- Prior art keywords
- noise
- signal
- correlation
- magnetic resonance
- nuclear magnetic
- 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
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种基于互相关的核磁共振全波信号噪声滤除方法,是利用噪声与拉莫尔频率的正弦信号不相关,而FID幅度衰减正弦信号与拉莫尔频率的正弦信号具有相关性的特点,通过互相关运算滤除噪声,然后拟合互相关波形的包络,并重构不含噪声的互相关波形,最后利用解卷积算法提取核磁共振全波数据中的FID信号。该方法运算数据计算量小,可以同时压制工频及其谐波噪声、随机噪声和尖峰噪声,明显提高核磁共振全波数据信噪比,有利于扩展核磁共振找水仪的应用范围和探测深度,而且该发明不需要额外的参考消噪装置,节约了成本,一次采集的全波数据通过互相关即可达到理想的滤波效果,节省测量时间。
Description
技术领域
本发明涉及一种核磁共振测深(Magnetic Resonance Sounding,MRS)信号噪声滤除方法,具体是基于互相关噪声研制原理的核磁共振全波信号噪声的滤除方法。
背景技术
核磁共振地下水探测利用人工激发的电磁场使地下水中的氢核形成宏观磁矩,这一宏观磁矩在地磁场中产生旋进运动,使用线圈接收宏观磁矩进动产生的电磁信号,该电磁信号的幅度随时间按指数形式衰减,成为自由感应衰减(FreeInduction Decay,FID)信号,通过对FID进行分析可以探测地下是否存在地下水。核磁共振找水是一种快速、直接的地下水探测方法。FID信号的表达式为:
上式中,ω0表示拉莫尔角频率,由当地的地磁场强度决定,一般为1000-3000Hz;FID信号的三个特征参数E0(初始振幅)、T2 *(平均弛豫时间)和(初始相位)分别与含水层的分布、厚度、平均含水率、渗透率和导电性等信息有着密切的关系,准确测量FID信号是进行核磁共振地下水探测的关键,测量FID信号时需要在野外勘探过程中将线圈接收信号的时域波形记录下来,该时域波形成为核磁共振全波信号。进行地下水探测时,核磁共振信号非常微弱,一般为nV级,而周围电磁场噪声却非常强烈,如工频及其谐波噪声、随机噪声和尖峰噪声,即使在信号检测电路中使用了带通滤波器,核磁共振全波信号的信噪比依然很低,一般低于0dB,这严重限制了核磁共振地下水探测的应用。
目前,核磁共振全波信号一般采用数据叠加方法,该方法能抑制大部分的随机噪声,但对工频噪声和尖峰噪声抑制效果差,且需要多次采集核磁共振信号,测量效率低。CN103823244A公开了一种磁共振三分量消噪装置,该发明基于参考通道噪声与主通道含噪信号中噪声的相关性来实现信噪分离,但由于某些噪声的无规律性和不稳定性以及混合机制的不确定性,限制了算法的应用,而且增加参考线圈和参考通道使仪器系统庞大复杂。CN104459809A公开了一种基于独立成分分析的全波核磁共振信号噪声滤除方法,首先利用核磁共振测深探水仪采集FID信号,通过频谱分析获得采集信号中含有的工频谐波干扰频率,采用数字正交法构造输入通道信号解决欠定盲源分离问题;然后,将构造的输入通道信号与采集的全波信号一并作为输入信号进行独立成分分析,得到分离的FID信号;最后采用频谱校正法解决分离后MRS信号幅值不定问题,进而提取去噪后的FID信号。该发明主要针对核磁共振全波信号中工频谐波干扰或某一单频干扰,不能滤除随机噪声和尖峰噪声,该方法算法过程复杂,计算量大,而且输出结果不稳定。
发明内容
本发明的目的就在于针对上述现有技术的不足,提供一种基于互相关的核磁共振全波信号噪声滤除方法。
本发明的目的是通过以下方式实现的:
一种基于互相关的核磁共振全波信号噪声滤除方法,包括以下步骤:
A、读取核磁共振接收机野外采集的全波数据s(k)及其采样率fs;
B、使用快速傅里叶变换求取全波数据的频谱S(ω);
C、读取频谱S(ω)中FID信号的频率f0所对应的相位θ0;
D、读取频谱S(ω)中,拉莫尔频率f0左右两侧的两个工频谐波频率f1和f2,及其所对应的幅度a1,a2,相位θ1和θ2;
E、重构两个工频谐波噪声的波形n1(k)和n2(k);
F、在全波数据s(k)中减去n1(k)和n2(k),得到x(k);
G、将x(k)与参考信号cos(2πkf0/fs)进行互相关运算,得到R(m);
H、利用希尔伯特变换获取R(m)的包络信号,En(m);
I、将包络信号En(m)分为长度相等的两段,并分别使用双指数函数拟合包络曲线;
J、使用双指数函数的拟合结果重构包络信号,Env(m)=a·exp(bmf0/fs)+c·exp(dmf0/fs),a、b、c、d为拟合参数;
K、构造不含噪声的互相关波形rf(m)=Env(m)·cos(2πmf0/fs+θ0);
使用参考信号cos(2πkf0/fs)对rf(m)进行解卷积运算,得到的解卷积结果ff(k)即为滤除噪声的核磁共振全波信号。
有益效果:本发明提供的基于互相关核磁共振全波信号噪声滤除方法,是利用噪声与拉莫尔频率的正弦信号不相关,而FID幅度衰减正弦信号与拉莫尔频率的正弦信号具有相关性的特点,通过互相关运算滤除噪声,然后拟合互相关波形的包络,并重构不含噪声的互相关波形,最后利用解卷积算法提取核磁共振全波数据中的FID信号。该方法运算数据计算量小,可以同时压制工频及其谐波噪声、随机噪声和尖峰噪声,明显提高核磁共振全波数据信噪比,有利于扩展核磁共振找水仪的应用范围和探测深度,而且该发明不需要额外的参考消噪装置,节约了成本,一次采集的全波数据通过互相关即可达到理想的滤波效果,节省测量时间。
附图说明
附图1是基于互相关的核磁共振全波信号噪声滤除方法的流程图。
附图2是野外采集的不含噪声的核磁共振全波数据的波形图。
附图3是野外采集的含有噪声的核磁共振全波数据的波形图。
附图4是含有噪声的核磁共振全波数据的频谱图。
附图5是核磁共振全波数据与参考信号进行互相关后的波形图。
附图6是互相关波形的包络拟合曲线。
附图7是重构后得到的不含噪声的互相关波形图。
附图8是解卷积得到的滤除噪声的核磁共振全波数据波形图。
具体实施方式
下面结合附图和实施例作进一步详细说明:
如图1所示,一种基于互相关的核磁共振全波信号噪声滤除方法,包括以下步骤:
A、读取核磁共振接收机野外采集的全波数据,核磁共振找水仪接收机野外工作时以恒定的采样率fs采集数据,读取接收机记录的全波数据s(k),k=1,2,3…N;
B、使用快速傅里叶变换求取全波数据的频谱S(ω),在Matlab中使用fft(s,N)命令可得到全波的频谱S(ω),S(ω)包含N个频率所对应的复数,频率间隔为fs/N;
C、读取频谱S(ω)中FID信号的频率f0所对应的相位θ0,S(ω)包含N个频率所对应的复数代表了该频率的幅度和相位信息,即
上式中,abs表示取模运算符,arctan为反正切运算符,im为取虚部运算符,re为取实部运算符;
D、读取频谱S(ω)中,拉莫尔频率f0左右两侧的两个工频谐波频率f1和f2,全波数据与参考信号作互相关运算时,如果在频谱有上接近拉莫尔频率的周期信号,该周期信号会影响互相关运算的结果,这是由于当周期噪声的频率接近信号频率时,周期噪声与FID信号的相关性增强,互相关噪声抑制效果变差。因此,需要将频率最接近拉莫尔频率的两个工频谐波滤除掉。读取这两个工频谐波噪声的幅度a1,a2,相位θ1和θ2;
E、重构两个工频谐波噪声的波形n1(k)和n2(k),由于这两个工频谐波为频率已知的周期信号,当其幅度和相位已知时,其波形就可唯一确定,可使用如下公式重构这两个工频谐波噪声的波形:
F、在全波数据s(k)中减去n1(k)和n2(k),得到x(k)
x(k)=s(k)-n1(k)-n2(k)
得到的x(k)不含有频率最接近拉莫尔频率的两个工频谐波噪声,但依然含有其他频率的工频噪声、随机噪声和尖峰噪声;
G、将x(k)与参考信号cos(2πkf0/fs)进行互相关运算,得到R(k),计算公式为:
相关检测是一种噪声抑制的有效方法,它基于信号与噪声的统计特性进行检测,互相关函数是两个时域信号相似性的一种度量,被检测信号与参考信号具有相关性,而与工频及其谐波噪声、随机噪声和尖峰噪声的相关性较差或没有相关性,利用这一差异可把被检测信号与噪声区分开了。虽然互相关检测一般应用于周期的正弦信号,但是核磁共振的FID信号是幅度按指数规律衰减正弦信号,与频率为拉莫尔频率的参考正弦信号cos(2πkf0/fs)具有较强相关性,而与随机噪声和尖峰噪声没有相关性,与频率不等于拉莫尔频率的周期噪声相关性较差,因此,互相关检测适用于核磁共振全波数据噪声滤除。两个长度均为N的数据进行互相关运算后得到的互相关数组长度为2N-1;
H、利用希尔伯特变换获取R(m)的包络信号En(m),互相关数据R(m)的希尔伯特变换为:
信号经希尔伯特变换后,在频域各频率分量的幅度保持不变,但相位将出现90°相移。得到解解析信号为:z(m)=R(m)+jR′(m)
解析信号z(k)的瞬时振幅即为互相关R(m)波形的包
络En(m);
I、将En(m)分为两段,每段长度均为N,即En1(k)=En(1:N),En2(k)=En(N:2N-1),分别使用双指数函数拟合En1(k)和En2(k),采用方法为最小二乘拟合,即
a1、b1、c1、d1;a1、b1、c1、d1为拟合得到的参数;
J、使用双指数函数的拟合结果重构包络信号,Env(1:N)=a1·exp(b1kf0/fs)+c1·exp(d1kf0/fs),Env(N:2N-1)=a2·exp(b2kf0/fs)+c2·exp(d2kf0/fs),拟合后的包络信号将不含有噪声;
K、构造不含噪声的互相关波形rf(m)=Env(m)·cos(2πmf0/fs+θ0),Env(m)为拟合得到的包络,cos(2πmf0/fs+θ0)为FID信号相位相同的正弦信号,其频率为拉莫尔频率;
使用参考信号cos(2πkf0/fs)对rf(k)进行解卷积运算,互相关函数可表示为卷积运算:
R12=f1(t)*f2(t),
“*”表示卷积运算,可见将f2(t)反褶(变量取符号)与f1(t)做卷积即得到f1(t)与f2的互相关函数,同时,由于参考信号cos(2πkf0/fs)的反褶为其本身,R(m)可以看作x(k)与参考信号cos(2πkf0/fs)的卷积运算,重构的互相关波形rf(m)是不含噪声的核磁全波数据与参考信号的卷积。因此,通过反卷积可以得到不含噪声的核磁全波数据。在Matlab是使用解卷积命令deconv即可得到ff(k):ff=deconv(rf,cos(2πf0t)),得到的解卷积结果ff(k)即为滤除噪声的核磁共振全波信号。
野外采集的核磁共振全波数据含有FID的波形为:
上式中,FID信号的初始幅度为100nV,平均弛豫时间T2 *为300ms,拉莫尔频率f0为2326Hz,初始相位为30°,接收机采样率为100kHz,信号采集时间为500ms,采样点数为50000。不含噪声的FID信号波形见图2所示。
采集到的核磁共振全波数据中除了微弱的FID信号外,还包含工频及其谐波噪声、随机噪声和尖峰噪声,噪声的表达式为:
n(k)=10-6square(t,50)+wag(1,N,-140)+10-6pulstran(t,0.05:0.05:0.5,'gauspuls')
上式中,Matlab中square用于产生50Hz方波,来模拟工频工频及其谐波噪声,方波的幅度为1μV;wgn(1,N,-140)用于产生标准差为100nV的零均值随机噪声;pulstran函数用于产生尖峰噪声,尖峰噪声的幅度为1μV;这些噪声为加性噪声叠加在FID信号上,核磁共振全波数据s(k)的表达式为:
s(k)=f(k)+n(k)
该核磁共振全波数据的信噪比为-30dB。包含噪声的核磁共振全波数据波形如图3所示。
在Matlab中使用函数fft函数对全波数据x(k)进行快速傅里叶变换,得到的频谱图如图4所示,频谱数据S(ω)包含50000个频率所对应的复数,频率间隔为2Hz。
从频谱数组S(ω)中找到第1164个元素为,对此复数取相角得到θ0为31°。
根据拉莫尔频率2326Hz,得到最接近拉莫尔频率的两个工频谐波频率分别为2250Hz和2350Hz,从频谱数组S(ω)中找到第1126个元素为1.3216e-09-1.3472e-08i,对此复数取模得到a1为27.73nV,θ1为-1.5°;从频谱数组S(ω)中找到第1176个元素为6.5767e-10-1.3849e-08i,对此复数取模得到a2为27.07nV,θ2为-1.4°。重构得到的这两个工频谐波波形为:
从核磁全波数据s(k)中减去这两个工频谐波噪声得到x(k),计算式如下:
x(k)=s(k)-n1(k)-n2(k)
由于减去的这两个谐波噪声相对总体噪声而言很小,核磁共振全波数据的波形几乎不会发生明显变化。将x(k)与参考信号cos(2πkf0/fs)按照如下公式进行互相关运算,在Matlab中使用xcorr命令即可得到两个数组的互相关函数,得到的R(m)波形如图5所示,可见互相关运算后的波形与原始的FID信号存在差异,波形可分为两段:前一段呈现幅度按指数规律增加的正弦信号,后一段呈现幅度按指数规律衰减的正弦信号。
利用希尔伯特变换获取R(m)的包络信号En(m),将En(m)分为两段,每段长度均为50000,即En1(k)=En(1:50000),En2(k)=En(50000:99999),分别使用双指数函数拟合En1(k)和En2(k),采用方法为最小二乘拟合,得到的拟合结果为:
a1=0.001364,b1=0.1022,c1=-0.001367,d1=-3.625;a2=-0.0001759,b2=0.4319,c2=0.001402,d2=-3.628。
使用双指数函数的拟合结果重构包络信号,Env(1:N)=a1·exp(b1kf0/fs)+c1·exp(d1kf0/fs),Env(N:2N-1)=a2·exp(b2kf0/fs)+c2·exp(d2kf0/fs),拟合后的包络信号将不含有噪声,图6是互相关波形的包络拟合曲线,包络拟合曲线分为两段。
构造不含噪声的互相关波形rf(m)=Env(m)·cos(2πmf0/fs+θ0),Env(m)为拟合得到的包络,图7是重构后得到的不含噪声的互相关波形图,重构后的互相关波形比原互相关波形光滑。
使用参考信号cos(2πkf0/fs)对rf(k)进行解卷积运算,在Matlab是使用解卷积命令deconv即可得到ff(k):ff=deconv(rf,cos(2πf0t)),得到的解卷积结果ff(k)即为滤除噪声的核磁共振全波信号,图8是解卷积得到的滤除噪声的核磁共振全波数据波形图,对比图2和图8发现二者波形一致,说明基于互相关的核磁共振全波信号噪声滤除方法能够恢复淹没在噪声中的FID信号。
Claims (1)
1.一种基于互相关的核磁共振全波信号噪声滤除方法,其特征在于,包括以下步骤:
A、读取核磁共振接收机野外采集的全波数据s(k)及其采样率fs;
B、使用快速傅里叶变换求取全波数据的频谱S(ω);
C、读取频谱S(ω)中FID信号的频率f0所对应的相位θ0;
D、读取频谱S(ω)中,拉莫尔频率f0左右两侧的两个工频谐波频率f1和f2,及其所对应的幅度a1,a2,相位θ1和θ2;
E、重构两个工频谐波噪声的波形n1(k)和n2(k);
F、在全波数据s(k)中减去n1(k)和n2(k),得到x(k);
G、将x(k)与参考信号cos(2πkf0/fs)进行互相关运算,得到R(m);
H、利用希尔伯特变换获取R(m)的包络信号,En(m);
I、将包络信号En(m)分为长度相等的两段,并分别使用双指数函数拟合包络曲线;
J、使用双指数函数的拟合结果重构包络信号,Env(m)=a·exp(bmf0/fs)+c·exp(dmf0/fs),a、b、c、d为拟合参数;
K、构造不含噪声的互相关波形rf(m)=Env(m)·cos(2πmf0/fs+θ0);使用参考信号cos(2πkf0/fs)对rf(m)进行解卷积运算,得到的解卷积结果ff(k)即为滤除噪声的核磁共振全波信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510256496.3A CN104898172B (zh) | 2015-05-19 | 2015-05-19 | 一种基于互相关的核磁共振全波信号噪声滤除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510256496.3A CN104898172B (zh) | 2015-05-19 | 2015-05-19 | 一种基于互相关的核磁共振全波信号噪声滤除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104898172A true CN104898172A (zh) | 2015-09-09 |
CN104898172B CN104898172B (zh) | 2017-05-10 |
Family
ID=54030932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510256496.3A Active CN104898172B (zh) | 2015-05-19 | 2015-05-19 | 一种基于互相关的核磁共振全波信号噪声滤除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104898172B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597557A (zh) * | 2016-11-03 | 2017-04-26 | 长安大学 | 一种基频可控的宽频带核磁共振及其最佳响应提取方法 |
CN106646637A (zh) * | 2016-12-27 | 2017-05-10 | 吉林大学 | 一种去除核磁信号中尖峰噪声的方法 |
CN106772646A (zh) * | 2016-12-19 | 2017-05-31 | 吉林大学 | 一种地面核磁共振信号提取方法 |
CN107045149A (zh) * | 2017-04-17 | 2017-08-15 | 吉林大学 | 一种基于双奇异值分解的全波核磁共振信号噪声滤除方法 |
CN108345039A (zh) * | 2018-01-12 | 2018-07-31 | 吉林大学 | 一种消除地面核磁共振数据中邻频谐波干扰的方法 |
CN108363109A (zh) * | 2018-01-12 | 2018-08-03 | 吉林大学 | 地面核磁共振数据的谐波干扰基频快速搜索方法 |
CN109100813A (zh) * | 2018-08-14 | 2018-12-28 | 吉林大学 | 一种基于协同滤波消除地面核磁共振数据中尖峰噪声的方法 |
CN109557015A (zh) * | 2016-11-01 | 2019-04-02 | 北京信息科技大学 | 流式细胞仪光谱重叠信号数字处理方法 |
CN109782219A (zh) * | 2018-12-26 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 一种基于相关性的信号累积测相方法 |
CN109885903A (zh) * | 2019-01-29 | 2019-06-14 | 吉林大学 | 一种基于模型的地面核磁共振信号尖峰噪声去除方法 |
CN110133559A (zh) * | 2019-05-21 | 2019-08-16 | 辽宁开普医疗系统有限公司 | 一种磁共振b0场扰动补偿系统及方法 |
CN111239837A (zh) * | 2020-02-20 | 2020-06-05 | 吉林大学 | 基于mcmc的地面核磁共振信号参数提取方法 |
CN111308561A (zh) * | 2020-03-11 | 2020-06-19 | 中国科学院地质与地球物理研究所 | 一种电磁信号强噪声去除方法 |
CN112180454A (zh) * | 2020-10-29 | 2021-01-05 | 吉林大学 | 一种基于ldmm的磁共振地下水探测随机噪声抑制方法 |
US11317822B2 (en) | 2020-07-31 | 2022-05-03 | Synex Medical Inc. | Weak signal detection system and method |
CN115017933A (zh) * | 2022-03-11 | 2022-09-06 | 吉林大学 | 基于类周期性小波系数恢复的核磁信号尖峰噪声抑制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1324306A (en) * | 1970-09-18 | 1973-07-25 | Commissariat Energie Atomique | Detection of submerged magnetic bodies for example underwater wrecks and submarines |
CN101251606A (zh) * | 2008-04-09 | 2008-08-27 | 吉林大学 | 弱信号检测仪器中有用信号频带内工频谐波干扰抑制电路 |
CN104502984A (zh) * | 2014-12-15 | 2015-04-08 | 吉林大学 | 特定频率噪声对消地下核磁共振探测装置及探测方法 |
CN104614778A (zh) * | 2015-01-27 | 2015-05-13 | 吉林大学 | 基于ica的核磁共振地下水探测信号噪声消除方法 |
-
2015
- 2015-05-19 CN CN201510256496.3A patent/CN104898172B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1324306A (en) * | 1970-09-18 | 1973-07-25 | Commissariat Energie Atomique | Detection of submerged magnetic bodies for example underwater wrecks and submarines |
CN101251606A (zh) * | 2008-04-09 | 2008-08-27 | 吉林大学 | 弱信号检测仪器中有用信号频带内工频谐波干扰抑制电路 |
CN104502984A (zh) * | 2014-12-15 | 2015-04-08 | 吉林大学 | 特定频率噪声对消地下核磁共振探测装置及探测方法 |
CN104614778A (zh) * | 2015-01-27 | 2015-05-13 | 吉林大学 | 基于ica的核磁共振地下水探测信号噪声消除方法 |
Non-Patent Citations (3)
Title |
---|
林君,等: "延时叠加工频谐波消除系统设计与实验", 《实验室研究与探索》 * |
田宝凤,等: "基于参考线圈和变步长自适应的磁共振信号噪声压制方法", 《地球物理学报》 * |
田宝凤,等: "核磁共振信号工频谐波的自适应滤除方法", 《吉林大学学报(信息科学版)》 * |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109557015A (zh) * | 2016-11-01 | 2019-04-02 | 北京信息科技大学 | 流式细胞仪光谱重叠信号数字处理方法 |
CN109557015B (zh) * | 2016-11-01 | 2021-05-25 | 北京信息科技大学 | 流式细胞仪光谱重叠信号数字处理方法 |
CN106597557A (zh) * | 2016-11-03 | 2017-04-26 | 长安大学 | 一种基频可控的宽频带核磁共振及其最佳响应提取方法 |
CN106597557B (zh) * | 2016-11-03 | 2018-11-23 | 长安大学 | 一种基频可控的宽频带核磁共振及其最佳响应提取方法 |
CN106772646B (zh) * | 2016-12-19 | 2018-11-06 | 吉林大学 | 一种地面核磁共振信号提取方法 |
CN106772646A (zh) * | 2016-12-19 | 2017-05-31 | 吉林大学 | 一种地面核磁共振信号提取方法 |
CN106646637A (zh) * | 2016-12-27 | 2017-05-10 | 吉林大学 | 一种去除核磁信号中尖峰噪声的方法 |
CN107045149A (zh) * | 2017-04-17 | 2017-08-15 | 吉林大学 | 一种基于双奇异值分解的全波核磁共振信号噪声滤除方法 |
CN107045149B (zh) * | 2017-04-17 | 2018-10-16 | 吉林大学 | 一种基于双奇异值分解的全波核磁共振信号噪声滤除方法 |
CN108363109A (zh) * | 2018-01-12 | 2018-08-03 | 吉林大学 | 地面核磁共振数据的谐波干扰基频快速搜索方法 |
CN108345039A (zh) * | 2018-01-12 | 2018-07-31 | 吉林大学 | 一种消除地面核磁共振数据中邻频谐波干扰的方法 |
CN108345039B (zh) * | 2018-01-12 | 2019-07-23 | 吉林大学 | 一种消除地面核磁共振数据中邻频谐波干扰的方法 |
CN108363109B (zh) * | 2018-01-12 | 2019-06-21 | 吉林大学 | 地面核磁共振数据的谐波干扰基频快速搜索方法 |
CN109100813B (zh) * | 2018-08-14 | 2019-07-12 | 吉林大学 | 一种基于协同滤波消除地面核磁共振数据中尖峰噪声的方法 |
CN109100813A (zh) * | 2018-08-14 | 2018-12-28 | 吉林大学 | 一种基于协同滤波消除地面核磁共振数据中尖峰噪声的方法 |
CN109782219A (zh) * | 2018-12-26 | 2019-05-21 | 中国电子科技集团公司第二十研究所 | 一种基于相关性的信号累积测相方法 |
CN109885903A (zh) * | 2019-01-29 | 2019-06-14 | 吉林大学 | 一种基于模型的地面核磁共振信号尖峰噪声去除方法 |
CN110133559A (zh) * | 2019-05-21 | 2019-08-16 | 辽宁开普医疗系统有限公司 | 一种磁共振b0场扰动补偿系统及方法 |
CN110133559B (zh) * | 2019-05-21 | 2021-05-11 | 辽宁开普医疗系统有限公司 | 一种磁共振b0场扰动补偿系统及方法 |
CN111239837B (zh) * | 2020-02-20 | 2021-01-05 | 吉林大学 | 基于mcmc的地面磁共振信号参数提取方法 |
CN111239837A (zh) * | 2020-02-20 | 2020-06-05 | 吉林大学 | 基于mcmc的地面核磁共振信号参数提取方法 |
CN111308561A (zh) * | 2020-03-11 | 2020-06-19 | 中国科学院地质与地球物理研究所 | 一种电磁信号强噪声去除方法 |
US11317822B2 (en) | 2020-07-31 | 2022-05-03 | Synex Medical Inc. | Weak signal detection system and method |
US11931137B2 (en) | 2020-07-31 | 2024-03-19 | Synex Medical Inc. | Weak signal detection system and method |
CN112180454A (zh) * | 2020-10-29 | 2021-01-05 | 吉林大学 | 一种基于ldmm的磁共振地下水探测随机噪声抑制方法 |
CN112180454B (zh) * | 2020-10-29 | 2023-03-14 | 吉林大学 | 一种基于ldmm的磁共振地下水探测随机噪声抑制方法 |
CN115017933A (zh) * | 2022-03-11 | 2022-09-06 | 吉林大学 | 基于类周期性小波系数恢复的核磁信号尖峰噪声抑制方法 |
CN115017933B (zh) * | 2022-03-11 | 2024-04-05 | 吉林大学 | 基于类周期性小波系数恢复的核磁信号尖峰噪声抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104898172B (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898172B (zh) | 一种基于互相关的核磁共振全波信号噪声滤除方法 | |
CN107783200B (zh) | 一种联合emd与tfpf算法的全波磁共振信号随机噪声消减方法 | |
CN109828318B (zh) | 一种基于变分模态分解的磁共振测深信号噪声滤除方法 | |
CN101881628A (zh) | 基于混沌系统和小波阈值去噪的微弱周期信号的检测方法 | |
CN106568589A (zh) | 一种基于经验小波变换碰摩声发射消噪方法 | |
CN103823244A (zh) | 磁共振三分量消噪装置及消噪方法 | |
Ghanati et al. | Joint application of a statistical optimization process and empirical mode decomposition to magnetic resonance sounding noise cancelation | |
CN108345039A (zh) | 一种消除地面核磁共振数据中邻频谐波干扰的方法 | |
CN110989020B (zh) | 一种音频大地电磁数据噪声干扰的滤波方法及系统 | |
CN108254794B (zh) | 一种基于建模反恢复技术的磁共振消噪方法及装置 | |
Gong et al. | Nonlinear vibration feature extraction based on power spectrum envelope adaptive empirical Fourier decomposition | |
Qin et al. | Dense framelets with two generators and their application in mechanical fault diagnosis | |
LIN et al. | Segmented time-frequency peak filtering for random noise reduction of MRS oscillating signal | |
CN108957552A (zh) | 基于ss-pca的地震数据海浪噪声压制方法 | |
CN103040459B (zh) | 一种多通道微弱生理信息记录系统中工频干扰的高保真滤除方法 | |
Qin et al. | A new method for multicomponent signal decomposition based on self-adaptive filtering | |
CN109782363B (zh) | 一种基于时域建模与频域对称的磁共振信号消噪方法 | |
CN110456420B (zh) | 基于近端参考线圈核磁共振地下水探测信号噪声消除方法 | |
Naghizadeh | A unified method for interpolation and de-noising of seismic records in the f-k domain | |
CN115291151B (zh) | 一种基于低相关分段的高精度磁共振信号频率测量方法 | |
CN108680620B (zh) | 一种用于电学层析成像系统的幅值解调方法 | |
CN112764095B (zh) | 基于广义能量比的vsp数据q值计算方法及系统 | |
CN117851756A (zh) | 一种基于二次频率耦合算法的磁共振测深信号提取方法 | |
Shahi et al. | Time-based noise removal from magnetic resonance sounding signals | |
CN111580188B (zh) | 一种大地电磁时间域标定方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |