主权项 |
一种基于统计的核磁共振一维谱信噪比提高方法,其特征在于,包括以下步骤:步骤1、核磁共振仪器载入设定的脉冲序列和采样参数,对核磁共振FID信号进行多次扫描采集;步骤2、将步骤1中各次扫描采集到的核磁共振FID信号分别保存;步骤3、将步骤2中各次扫描对应的核磁共振FID信号进行傅立叶变换,得到各次扫描对应的谱数据并组成二维数组Spec(n,m),其中n表示扫描总次数,m表示各次扫描对应的谱数据的点数;步骤4、根据步骤3中二维数组Spec(n,m)计算各次扫描的谱数据的相同数据点索引位置的数据点的变异系数CV<sub>j</sub>,将所有变异系数CV<sub>j</sub>组成变异系数数组CV,其中j表示各次扫描对应的数据点索引位置,j∈[1,m];步骤5、将步骤3得到二维数组Spec(n,m)对应各次扫描对应的谱数据的相同数据点索引位置的数据值进行累加,得到累加后的各次扫描对应的谱数据的相同数据点索引位置的数据值AddedSpec<sub>j</sub>,各AddedSpec<sub>j</sub>组成累加谱数据组AddedSpec;步骤6、计算步骤4中得到的变异系数数组CV的噪声水平ε,步骤7、根据步骤4中得到的变异系数数组CV和步骤6得到的噪声水平ε通过以下公式获得权重数组w:<maths num="0001"><math><![CDATA[<mrow><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><mi>α</mi><mo>+</mo><mfrac><mrow><mn>1</mn><mo>-</mo><mi>α</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mi>β</mi><mo>/</mo><mi>ϵ</mi><mo>*</mo><mo>(</mo><mrow><msub><mi>CV</mi><mi>i</mi></msub><mo>-</mo><mi>γ</mi><mo>*</mo><mi>ϵ</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001164296130000011.GIF" wi="713" he="127" /></maths>其中,α、β、γ均为预设的系数,α表示累加谱数据组AddedSpec的噪声系数,β表示抑制噪声的程度,CV<sub>i</sub>为步骤4得到的变异系数数组CV的第i个数值,ε为步骤6得到的噪声水平,γ表示噪声与峰进行区分的系数,w<sub>i</sub>为权重数组w的第i个数值,步骤8、将步骤7中得到的权重数组w,与步骤5中得到的累加谱数据组AddedSpec逐点相乘,得到新的谱数据weightedSpec。 |