[發明專利]一種核磁共振譜采樣截斷偽峰的抑制方法

有效

• 申請號: 201610310469.4
• 申請日: 2016-05-11
• 公開（公告）號: CN105928965B
• 公開（公告）日: 2017-09-15
• 發明（設計）人: 蔣濱;孫鵬;羅凡;張許;楊運煌;李從剛;周欣;劉買利
• 申請（專利權）人: 中國科學院武漢物理與數學研究所
• 主分類號: G01N24/08
• 分類號: G01N24/08
• 代理公司: 武漢宇晨專利事務所42001
• 代理人: 李鵬,王敏鋒

1.一種核磁共振譜采樣截斷偽峰的抑制方法，其特征在于，包括以下步驟：

步驟1：設定第一百分比閾值p和第二百分比閾值q，第一百分比閾值p的取值范圍為95％～98％，第二百分比閾值q的取值范圍為0.005％～0.015％，

步驟2：獲得的時域信號a，時域信號a的采樣數據點點數為N₁，對時域信號a進行填零處理，填零后的時域信號a的采樣數據點點數為N₂，然后對進填零后的時域信號a進行傅立葉變換后獲得NMR譜圖A，

步驟3：求得NMR譜圖A中最高信號峰的強度值I，以最高信號峰的強度值I與第一百分比閾值p的乘積pI為第一比較閾值，將NMR譜圖A中所有強度值大于第一比較閾值的采樣數據點對應的強度值分別減去第一比較閾值，并將NMR譜圖A中所有強度值小于等于第一比較閾值的采樣數據點對應的強度值設定為0，得到NMR譜圖B；將NMR譜圖A中所有強度值大于第一比較閾值的采樣數據點對應的強度值設定為第一比較閾值，且NMR譜圖A中所有強度值小于等于第一比較閾值的采樣數據點對應的強度值不變，得到NMR譜圖C，

步驟4、將步驟2獲得的NMR譜圖A中最高信號峰的強度值I與第二百分比閾值q的乘積qI為第二比較閾值，

若NMR譜圖B中最高信號峰的強度值小于第二比較閾值，則進入步驟7；

若NMR譜圖B中最高信號峰的強度值大于等于第二比較閾值，則進入步驟5，

步驟5：將NMR譜圖B進行逆傅立葉變換得到時域信號b，并對時域信號b作截斷處理，即將N₂個采樣數據點的時域信號b只保留前面的N₁個采樣數據點，

步驟6：將時域信號a減去截斷處理后的時域信號b獲得新的時域信號a，將新的時域信號a作為步驟2中的時域信號a，循環迭代進入步驟2，

步驟7：將各個循環迭代中獲得的NMR譜圖B，以及最后一次循環迭代中獲得的NMR譜圖C相疊加，即得到NMR譜圖D，

步驟8：對NMR譜圖D作窗口平滑處理修正，得到最終NMR譜圖E。