一种增强磁共振谱仪动态范围的方法

resonance imaging:quantization noise reduction and data memory saving,j.ieee aes magazine,june 1990将k空间相位编码行低频分量进行压缩，高频分量进行扩张。通过仿真表明，非线性放大将量化噪声降低了10到20倍，折合动态范围增加了4位左右。根据论文j.bollenbeck,m.vester,r.oppelt,h.kroeckel,w.schnell.a high performance multi
‑
channel rf receiver for magnet resonance imaging systems,j.proc.intl.soc.mag.reson.med.13(2005)可知，2005年西门子公司提出了一种非线性压缩的具体实现方案，在mri接收机adc之前设计了一种信号幅度压缩器/扩展器(压缩扩展器)，该技术规避了对接收机内部自动增益控制的需求，当adc转换速率为10ms/s时，mri dr能够做到164db/hz，对应有效位数(enob)为16位。如果adc前没有压缩扩展器则需要18位分辨率才能达到164db/hz的dr。由于西门子公司提供的方案必须在模拟电路中进行信号幅值的动态压缩，因此这种设计在工程实现上是一项极具挑战性的任务，设计指标除了满足基本功能(dr和压缩特性)外，还必须综合考虑非线性压缩信号可恢复性、电压偏移、温度稳定性等因素的影响。


技术实现要素：

7.为解决上述技术问题，本发明提出一种增强磁共振谱仪动态范围的方法，在相位编码和频率编码方向都用线性梯度脉冲进行编码，并且本发明不使用非线性压缩的模拟电路对一个相位编码行的不同频率编码所采集信号进行非线性放大。
8.本发明采用的技术方案为：一种增强磁共振谱仪动态范围的方法，包括：
9.s1、通过预扫描得到不同增益设置下k空间的回波信号的幅值分布情况；
10.s2、对不同增益对应的增益差和相位进行迭代校正；
11.s3、正式扫描启动后，将k空间相位编码行分成n部分，每部分设置一个接收增益，得到完整的k空间数据。
12.步骤s1所述与扫描的频率编码步数和相位编码步数设置为正式扫描的1/4。
13.步骤s1具体为：按照6db增益步进从大到小分别设置n个增益控制参数g1、g2
…
gn；其中相位编码有m步，在每个增益参数控制下，通过扫描得到一组k空间数据；依次从最大增益g1开始，分析对应增益参数控制下扫描采集到的k空间数据，找出所有回波幅值未溢出的相位编码行；后续增益控制参数依次在剩余的若干步相位编码里找出所有回波幅值未溢出的相位编码行；
14.记g1到g
n
增益参数控制下剩余的未溢出的相位编码步数分别是m1、m2…
m
n
，其中m1+m2…
+m
n
＝m。
15.步骤s2具体为：从最大增益g1开始，对相邻的两个增益控制参数控制下各自对应的相位编码步数所采集回波对应行做增益和相位比较，得到相邻的两个增益控制参数控制下的统计增益差因子与回波相位差。
16.步骤s3包括：将k空间每次相位编码行的频率编码步数和总相位编码步数扩充为预扫描时的4倍。
17.步骤s3还包括对对不同增益步进控制下相位编码行数据进行归一化，具体的：增益g1控制下的相位编码为4m1步数据保持不变；后续增益控制参数控制下对应相位编码步数的数据根据该增益控制参数之前的所有统计增益差因子与回波相位差进行归一化。
18.归一化计算过程为：g
i
控制下相位编码为4m
i
步数据各个点的复信号乘以步数据各个点的复信号乘以
19.本发明的有益效果：本发明基于常规的adc(16bit位宽&80mhz采样率)进行磁共振回波信号的采集，在相位编码和频率编码方向都用线性梯度脉冲进行编码，并且本发明不使用非线性压缩的模拟电路对一个相位编码行的不同频率编码所采集信号进行非线性放大。本发明通过预扫描得到不同增益设置下k空间的回波信号的幅值分布情况，正式扫描启动后，将k空间相位编码行分成n部分，每部分设置一个接收增益，使得相应部分的相位编码行对应的回波信号到达adc进行模数转换时能够充分利用adc的动态范围；本发明具备以下优点：
20.1、降低了对磁共振谱仪接收机adc的性能要求，磁共振接收机在固定增益控制下adc的最佳采样位至少需要19位才能覆盖mr信号的全动态范围，接收机在可变增益控制下adc的最佳采样位宽为15位以上即可覆盖mr信号的全动态范围；
21.2、本发明在不修改磁共振接收机电路的情况下，通过优化序列的扫描参数提高了图像的动态范围；
22.3、整个序列扫描的时间只增加5％左右，而重建图像的信噪比提高了10～20db(adc位宽相当于提高4位)；
23.4、本发明的迭代信号增益和相位校正方法，规避了噪声基底对幅值较小回波的影响，减少了校正误差。
附图说明
24.图1为本发明的方案流程图；
25.图2为本发明实施例提供的预扫描(prescan)的k空间示意图；
26.图3为本发明实施例提供的进行不同增益步进下增益差和相位的迭代校正示意图；
27.其中，图3(a)为通过直方图得到不同增益步进下增益差的示意图，图3(b)为通过另外一幅直方图得到不同增益步进下相位的示意图；
28.图4为校正后的预扫描(prescan)的k空间示意图；
29.图5为本发明实施例提供的正式扫描的k空间示意图；
30.图6为本发明方法的仿真结果；
31.其中，图6(a)为固定增益下不同adc位宽与采集信噪比的关系，图6(b)为动态调整增益下不同adc位宽与采集信噪比的关系。
具体实施方式
32.为便于本领域技术人员理解本发明的技术内容，下面结合附图对本发明内容进一步阐释。
33.本发明基于常规的adc(16bit位宽&80mhz采样率)进行磁共振回波信号的采集，在相位编码和频率编码方向都用线性梯度脉冲进行编码，并且本发明不使用非线性压缩的模拟电路对一个相位编码行的不同频率编码所采集信号进行非线性放大。
34.本发明通过预扫描得到不同增益设置下k空间的回波信号的幅值分布情况，正式
扫描启动后，将k空间相位编码行分成n部分，每部分设置一个接收增益，使得相应部分的相位编码行对应的回波信号到达adc进行模数转换时能够充分利用adc的动态范围。
35.本发明的实现过程如图1所示，包括以下三个步骤：
36.1、在预扫描时，如图2所示，横坐标为频率步数，纵坐标为相位编码步数，设置扫描矩阵大小为正式扫描的1/4，即将频率编码步数和相位编码步数分别设置为正式扫描的1/4。按照6db增益步进从大到小分别设置n个增益控制参数g1、g2…
g
n
，即在这n个接收增益控制参数下各扫描一次；这样设计的好处可以大幅度减少预扫描的时间，其中相位编码有m步，在每个增益参数控制下，通过扫描得到一组k空间数据。依次从最大增益g1开始，分析该增益参数控制下扫描采集到的k空间数据，找出所有回波幅值未溢出的相位编码行，合计m1步相位编码；然后分析g2控制下扫描采集到的k空间数据，从剩余的m
‑
m1步相位编码里找出所有回波幅值未溢出的相位编码行，合计m2步相位编码；按照以上规律分别得到g3到g
n
增益参数控制下剩余的未溢出的相位编码步数分别是m3、m4…
m
n
，其中m1+m2…
+m
n
＝m。
37.2、进行不同增益步进下增益差和相位的迭代校正。将数据格式设置为32位浮点。g1控制下m1步相位编码所采集回波与g2控制下同样这些m1步相位编码所采集回波对应行做增益和相位比较。如图3所示，通过直方图得到g1与g2控制下的统计增益差因子为δg1(如图3(a)所示)，同样通过另外一幅直方图得到g2相对于g1控制下的回波相位差为φ1(如图3(b)所示)，依次规律，分别得到g2与g3的统计增益差因子为δg2，g3相对于g2控制下的回波相位差为φ2，
……
，以及g
n
‑1与g
n
的统计增益差因子为δg
n
‑1，g
n
相对于g
n
‑1控制下的回波相位差为φ
n
‑1。
38.图4为根据图3得到相邻接收增益间的增益差和相位差，对k空间各部分在不同接收增益控制下k空间各部分的增益和相位补偿实施方法，即k空间各部分在不同接收增益控制下得到的数据归一化到同一接收增益下；图4中右侧是增益差和相位差补偿值，用于补偿图4中对应括号区域的数据。
39.3、启动正式扫描，如图5所示，将k空间每次相位编码行的频率编码步数和总相位编码步数由预扫描的m扩充为4m，则增益g1控制下的相位编码为4m1步，g2控制下相位编码为4m2步，以此类推，g
n
控制下的相位编码步数为4m
n
步。扫描采集得到完整的k空间数据，分别对不同增益步进控制下相位编码行数据进行归一化。增益g1控制下的相位编码为4m1步数据保持不变；g2控制下相位编码为4m2步数据各个点的复信号乘以g3控制下相位编码为4m3步数据各个点的复信号乘以以此类推，g
n
控制下相位编码为4m
n
步数据各个点的复信号乘以图5中g1,
…
,g
n
为对应k空间区域的接收增益参数。
40.本实施例中k空间横坐标方向是频率编码，纵坐标方向是相位编码，每一行就是相位编码的一步。
41.如图6所示的仿真结果，表达了固定增益和动态调整增益下的adc位宽与采集信号信噪比的对比关系，图6(a)为固定增益下不同adc位宽与采集信噪比的关系，图6(b)为动态调整增益下不同adc位宽与采集信噪比的关系，其中gc为动态增益(gain control的缩写)。通过图6仿真可知谱仪在扫描过程中通过增益动态调整，15位位宽的adc就能实现相当于扫描过程中增益固定且adc位宽为19位的采集信号动态范围，且该方法不需要对谱仪进行任
何硬件修改。
42.本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。