本發明提出一種核磁共振t1噪聲的壓制方法，包括如下步驟：載入核磁共振多維譜數據S(m,n)，并計算直接維投影Psubgt;1/subgt;；將直接維投影Psubgt;1/subgt;采用連續小波變換求得一階導數譜Psubgt;2/subgt;，并對Psubgt;2/subgt;采用滑動窗口法得到譜峰區域，并將該區域作為t1噪聲區域；取直接維投影Psubgt;1/subgt;的非t1噪聲區域的連續d列數據，分別求取各列的噪聲水平，并將各列噪聲水平的平均值σsubgt;1/subgt;作為t1噪聲壓制的目標噪聲水平；取t1噪聲區域對應的矩陣S中的某列，采用連續小波變換得到一階導數Psubgt;3/subgt;，并用滑動窗口算法得到該列的噪聲區域Psubgt;4/subgt;；對t1噪聲區域Psubgt;4/subgt;采用Whittaker算法進行平滑，使該區域的噪聲水平被壓制到目標值σsubgt;1/subgt;；繼續取t1噪聲區域對應的矩陣S的下一列數據，計算噪聲區域Psubgt;4/subgt;并進行Whittaker算法進行平滑，直到t1噪聲區域的平滑全部完成。

【技術領域】

本發明涉及核磁共振數據處理技術領域，尤其涉及一種核磁共振t1噪聲的壓制方法。

【背景技術】

核磁共振在生物大分子、分子動力學、藥物研發等領域取得了重要應用。蛋白質等復雜結構的大分子使核磁共振譜圖重疊嚴重，分辨率降低，對譜圖的辨認造成困難。核磁共振多維譜是解決譜峰重疊的重要手段，但由于儀器的不穩定性使間接維的噪聲水平增大，形成t1噪聲，對譜圖的辨認形成干擾。為提高核磁共振多維譜圖的數據質量，如何壓制t1噪聲是核磁共振數據處理的重要課題。

現有的核磁共振t1噪聲壓制方法主要分為三種類型：1.在脈沖序列執行期間使用線性梯度場進行相干路徑選擇(Horne?T?J,Morris?G?A.P-Type?gradient-enhanced?COSYexperiments?show?lowert1?noise?than?N-type[J].Magnetic?Resonance?inChemistry,1997,35(10):680-686.)；2.通過抑制對角峰，以減小較強的對角峰形成的t1噪聲(Denk?W,Wagner?G,Rance?M,et?al.Combined?suppression?of?diagonal?peaks?andt1?ridges?in?two-dimensional?nuclear?overhauser?enhancement?spectra[J].Journal?of?Magnetic?Resonance,1985,62(2):350-355.)；3.在數據處理階段采用參考去卷積(Gibbs?A,Morris?G?A,Swanson?AG,et?al.Suppression?of?t1Noise?in?2D?NMRSpectroscopy?by?Reference?Deconvolution[J].Journal?of?Magnetic?Resonance,Series?A,1993,101(3):351-356)、數據平滑(Glaser?S?J,Kalbitzer?H?R.Improvementof?two-dimensional?NMR?spectra?by?weighted?mean?t1-ridge?subtraction?andantidiagonal?reduction[J].Journal?of?Magnetic?Resonance,1986,68(2):350-354.)、數據相關性去噪(Simon?Poulding,Adrian?J.Charlton,James?Donarski.Removal?of?t1noise?from?metabolomic?2D?1H–13C?HSQC?NMR?spectra?by?Correlated?TraceDenoising[J].Journal?of?Magnetic?Resonance,189(2):190-199.)等。其中，基于數據處理的t1噪聲壓制方法應用較為廣泛，但是，目前的方法普遍執行效率較低，且對譜峰的正負、寬窄等不同特征的普適性不強。

鑒于此，實有必要提供一種新型的核磁共振t1噪聲的壓制方法以克服上述缺陷。

【發明內容】