CN114879107A - 一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 - Google Patents
一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 Download PDFInfo
- Publication number
- CN114879107A CN114879107A CN202210594734.1A CN202210594734A CN114879107A CN 114879107 A CN114879107 A CN 114879107A CN 202210594734 A CN202210594734 A CN 202210594734A CN 114879107 A CN114879107 A CN 114879107A
- Authority
- CN
- China
- Prior art keywords
- intrinsic
- time domain
- coil
- receiving coil
- frequency receiving
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000002599 functional magnetic resonance imaging Methods 0.000 title claims abstract description 44
- 238000003384 imaging method Methods 0.000 claims abstract description 66
- 238000004088 simulation Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 16
- 230000035945 sensitivity Effects 0.000 claims description 33
- 230000002123 temporal effect Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 6
- 239000004033 plastic Substances 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 4
- -1 polytetrafluoroethylene Polymers 0.000 claims description 4
- 229920001343 polytetrafluoroethylene Polymers 0.000 claims description 4
- 239000004810 polytetrafluoroethylene Substances 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 24
- 230000008859 change Effects 0.000 abstract description 3
- 230000015654 memory Effects 0.000 description 12
- 230000008878 coupling Effects 0.000 description 7
- 238000010168 coupling process Methods 0.000 description 7
- 238000005859 coupling reaction Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 230000005520 electrodynamics Effects 0.000 description 6
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000005672 electromagnetic field Effects 0.000 description 4
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 241000080590 Niso Species 0.000 description 1
- 239000004809 Teflon Substances 0.000 description 1
- 229920006362 Teflon® Polymers 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 210000005013 brain tissue Anatomy 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000008367 deionised water Substances 0.000 description 1
- 229910021641 deionized water Inorganic materials 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
- G01R33/3671—Electrical details, e.g. matching or coupling of the coil to the receiver involving modulation of the quality factor of the RF coil
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种fMRI中射频接收线圈本征时域稳定性参数的测量方法，包括：通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈‑模体的间距，以模拟线圈与模体之间的相对运动；采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，并进行处理重建，得到本征成像数据；根据本征成像数据计算得到线圈的本征时域稳定性参数。本发明可排除所有可能存在于功能磁共振采集中复杂的生理噪声，只考虑由于射频接收线圈与被试模体之间的相对运动带来的时域噪声，获取的本征时域稳定性参数可用于指导功能磁共振成像应用中的专用射频接收线圈的设计与使用，提高功能磁共振成像的图像性能。
Description
技术领域
本发明属于功能磁共振成像技术领域，具体涉及一种fMRI中射频接收线圈本征时域稳定性参数的测量方法。
背景技术
在7T(tesla，特斯拉)超高场磁共振成像系统中，由于磁场强度的提高，图像信噪比以及功能磁共振成像(fMRI，functional Magnetic Resonance Imaging))中检测神经元活动时的对比度也随之提高，因此超高场磁共振被广泛应用于亚毫米功能成像中。然而，在超高场磁共振成像中，随着磁场强度的提高，被试运动利用的时域噪声也随之增强，从而极大削弱了超高场下可能达到的成像潜能。迄今为止，功能磁共振成像在缓解时域噪声方面主要致力于改进数据获取方法，如图像后处理和运动校正算法，但很少能从根本上解决磁场强度相关的时域噪声问题。
其中，与磁场强度相关的时域噪声可以归因于随着场强提升的射频工作频率，该射频工作频率使得电磁场与被试之间产生更复杂的相互作用。具体的，成像被试与射频线圈之间由于电动力学耦合而作为介质负载，当射频接收线圈的空间位置相对于成像被试发生变化时，电动力学耦合会受到干扰；当线圈与被成像被试的距离保持不变时，电动力学耦合才会达到稳态。
现有技术中，MRI(Magnetic Resonance Imaging，磁共振成像)获取序列和图像后处理算法中常假设射频接收线圈和被试之间具有恒定的耦合水平，但这在fMRI中显然不适用。原因在于：功能磁共振扫描中人体大脑内部不可避免地发生变化，即使受试者的头部是静止的，也会有脑组织的非刚性运动，以及血液和脑脊液等的流动，从而干扰成像被试与射频接收线圈之间的相互作用。尽管有学者研究了射频接收线圈参与影响功能磁共振成像的时域信噪比(tSNR，time SIGNAL NOISE RATIO)，但这些研究均无法将射频接收线圈与其他噪声，如生理噪声相互剥离，从而无法获知射频接收线圈本身带来怎样的时域噪声，进而无法获知该射频接收线圈是否适用于功能磁共振成像中。
因此，亟需提出一种在功能磁共振成像中，能够有效测量射频接收线圈本身带来的时域噪声的方法。
发明内容
本发明的目的是提供一种fMRI中射频接收线圈本征时域稳定性参数的测量方法，用于解决现有技术的方法无法将射频接收线圈与其他噪声相互剥离，进而无法获知射频接收线圈本身带来怎样的时域噪声的技术问题。
为了实现上述目的，本发明采用以下技术方案：
第一方面，本发明提供一种fMRI中射频接收线圈本征时域稳定性参数的测量方法，包括：
通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，其中，所述介质件需要满足电导率小于1×10-15S/m，相对介电常数介于1-3之间；
采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，并对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据；
根据所述本征成像数据计算得到线圈的本征时域稳定性参数，所述本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性。
在一种可能的设计中，所述射频接收线圈至少包括单通道接收线圈和/或多通道接收阵列线圈。
在一种可能的设计中，所述被试仿真模体包括水模体，所述水模体内部设有模拟人体电特性的填充液，所述电特性至少包括介电常数和电导率。
在一种可能的设计中，所述介质件包括填充有聚四氟乙烯材料的塑料垫。
在一种可能的设计中，采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，包括：
在插入每一厚度的介质件时，利用短TE序列对所述被试仿真模体进行质子密度加权成像，设置N个时域采样点/次，采集得到不同体素分辨率下的多组回波EPI图像。
在一种可能的设计中，当采用单通道接收线圈时，对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据，包括：
利用傅立叶变换对原始K空间的回波数据进行处理重建，得到本征成像数据，其中，所述本征成像数据的矩阵大小与k空间矩阵大小一致。
在一种可能的设计中，当采用多通道接收阵列线圈时，对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据，包括：
利用平方和方法对原始K空间的三维多通道回波数据进行加权组合重建，得到多通道信息压缩后的本征成像数据。
在一种可能的设计中，根据所述本征成像数据计算得到线圈本征时域稳定性参数，包括：
根据使用不同体素分辨率采集得到的本征成像数据的时域数据，分别计算对应不同体素分辨率下的本征时域信噪比tSNR*，包括：
计算每一像素点的时域均值和时域标准差，将所述时域均值作为线圈的本征时域灵敏度均值
并将所述时域标准差作为线圈的本征时域噪声标准差σt *′；
根据所述本征时域灵敏度均值
和所述本征时域噪声标准差σt *′，计算得到本征时域信噪比tSNR*，其中，
在一种可能的设计中，在根据所述本征成像数据计算得到线圈的本征时域稳定性参数之前，所述方法还包括：
将在射频电源开启时采集得到的本征成像数据作为线圈本征灵敏度数据S*；
将在射频电源关闭时采集得到的本征成像数据作为热噪声数据，计算得到热噪声标准差
根据所述线圈本征灵敏度数据S*和所述热噪声标准差
计算得到本征信噪比SNR*，
在一种可能的设计中，根据所述本征成像数据计算得到线圈本征时域稳定性参数，包括：
对所述线圈本征灵敏度数据S*的时域标准差进行归一化计算，得到本征时域灵敏度稳定性λ*，
其中，σs *表示线圈本征灵敏度的时域标准差；
对所述热噪声数据的时域标准差进行归一化计算，得到本时域热噪声稳定性α*，
其中，
表示被试仿真模体噪声的时域标准差。
第二方面，本发明提供一种fMRI中射频接收线圈本征时域稳定性参数的测量装置，包括：
相对运动模拟模块，用于通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，其中，所述介质件需要满足电导率小于1×10-15S/m，相对介电常数介于1-3之间；
成像数据采集模块，用于采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，并对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据；
时域参数计算模块，用于根据所述本征成像数据计算得到线圈的本征时域稳定性参数，所述本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性。
第三方面，本发明提供一种计算机设备，包括依次通信相连的存储器、处理器和收发器，其中，所述存储器用于存储计算机程序，所述收发器用于收发消息，所述处理器用于读取所述计算机程序，执行如第一方面任意一种可能的设计中所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法。
第四方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，执行如第一方面任意一种可能的设计中所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法。
第五方面，本发明提供一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行如第一方面任意一种可能的设计中所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法。
有益效果：
本发明通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，从而可以排除所有可能存在于功能磁共振采集中复杂的生理噪声，只考虑由于射频接收线圈与被试模体之间的相对运动带来的时域噪声；通过采集在插入不同厚度的介质件时被试仿真模体的多个回波EPI图像，并对多个回波图像进行处理重建，得到本征成像数据；根据本征成像数据计算得到线圈的本征时域稳定性参数，本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性，该本征时域稳定性参数可用于指导功能磁共振成像应用中的专用射频接收线圈的设计与使用，提高功能磁共振成像的图像性能。
附图说明
图1为本实施例中的fMRI中射频接收线圈本征时域稳定性参数的测量方法的流程图。
具体实施方式
为使本说明书实施例的目的、技术方案和优点更加清楚，下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
由于射频接收线圈的电磁场决定了磁共振信号本征灵敏度和本征热噪声水平，基于此已有本征信噪比SNR*的概念被提出，通过本征信噪比SNR*可以根据数值电磁仿真计算来确定实际磁共振扫描时图像信噪比的上限；然而，本征信噪比无法评价时域信号带来的动态噪声，因此，射频接收线圈的本征时域性能应当被量化，并且可以用来评估射频接收线圈对时域噪声的贡献。
此外，由于在MRI扫描中由射频接收线圈的电动力学所决定的磁共振信号本征灵敏度与本征热噪声水平会有波动，本实施例将其定义为本征时域噪声。由于UHF(特高频)的射频波长较短，由此产生的本征时域噪声可能是与磁场强度相关的时域噪声的主要组成部分，从而恶化了亚毫米空间尺度下的成像性能，因此，本征时域噪声的确定对于提高成像性能具有重要意义。
此外，现代MRI系统使用射频接收线圈的相控阵列设计，以实现在更大的视野范围上获得比较高的图像信噪比，同时可实现并行成像技术以提高图像编码效率。但线圈不同单元之间与被试载荷之间存在复杂的电动力学耦合，不可避免地会影响并行成像的灵敏度、热噪声和噪声放大因子(g因子)。这种电动力学耦合的扰动会导致相控阵列单元间相互耦合发生变化，从而使本征时域噪声更加复杂。目前还没有方法能够将射频接收线圈与其他噪声，如生理噪声相互剥离，从而无法获知射频接收线圈本身带来怎样的时域噪声，进而无法获知该射频接收线圈是否适用于功能磁共振成像中。基于此，本实施例提出了一种fMRI中射频接收线圈本征时域稳定性参数的测量方法，并通过该方法测量得到射频接收线圈本征时域稳定性参数，用于指导功能磁共振成像应用中的专用射频接收线圈的设计与使用，提高功能磁共振成像的图像性能。以下通过实施例对该方法进行说明，具体如下：
如图1所示，第一方面，本实施例提供一种fMRI中射频接收线圈本征时域稳定性参数的测量方法，包括但不限于由步骤S101～S103实现：
步骤S101.通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，其中，所述介质件需要满足电导率小于1×10-15S/m，介电常数介于1-3之间，即介电性能对电磁场造成的影响非常小；
在步骤S101中，所述射频接收线圈至少包括单通道接收线圈和/或多通道接收阵列线圈，优选的，在本实施例中，选用三个直径为2cm、3.5cm和5cm的单通道接收线圈以及一个覆盖面为5cm直径圆的16通道接收阵列线圈，当然，可以理解的是，本实施例不限于采用上述直径或上述通道数的线圈，任意能够实现本实施例发明目的的线圈结构，均属于本发明的保护范围之内，此处不做限定。
在步骤S101中，所述被试仿真模体包括水模体，所述水模体内部设有模拟人体人体电特性的填充液，所述电特性至少包括介电常数和电导率，优选的，所述水模体采用底面直径为11cm，长为20cm的均匀圆柱水模体，且该水模体内的填充液为含有37.5％的NiSO4·5H2O和0.5％的NaCl的去离子水；当然，可以理解的，本实施例的被试仿真模体不限于采用水模体，任意能够模拟被试人体电特性的模体，均适用于本实施例的方法中，其中，该仿真模体的电导率和介电常数与人体一致，其体积或形状等参数可根据模拟不同人体而有所不同。
在步骤S101中，所述介质件包括填充有聚四氟乙烯材料的塑料垫，优选的，本实施例分别选用厚度为3mm、6mm和8mm的塑料垫，并在每次对被试仿真模体进行扫描之前插入其中一种厚度的塑料垫，则通过插入不同厚度的塑料垫来改变线圈与被试仿真模体之间的间距，从而模拟线圈和被试仿真模体的相对运动；此外，为了保证介质件对射频接收线圈的电磁场特性的影响尽量降到最低，本实施例选用介电性能低于阈值的介质件，优选的，所述介质件采用聚四氟乙烯材料，即特氟龙，介电常数εr＝2.1，当然，可以理解的是，本实施例还能够选用其它低介电常数的材料，此处不做限定。
步骤S102.采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，并对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据；
其中，需要说明的是，优选的，所述回波图像为单次回波EPI(Echo PlanarImaging，平面回波成像)图像。
在步骤S102一种具体的实施方式中，采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，包括：
在插入每一厚度的介质件时，利用短TE序列对所述被试仿真模体进行质子密度加权成像，设置N个时域采样点/次，采集得到不同体素分辨率下的多组回波图像。
其中，优选的，本实施例设置30时域采样点/次，并分别采集体素大小为0.85×0.85×1mm3,1.15×1.15×1mm3,1.45×1.45×1mm3,1.75×1.75×1mm3、2×2×1mm3的回波图像，从而可以获得不同的本征信噪比。需要说明的是，由于在采集每组数据时，仅改变了体素大小这一扫描参数，而尽可能不影响其他参数，例如成像范围、对比度等，从而可以通过仅改变采集图像的体素大小来获取不同的本征信噪比。那么，本实施例每一体素大小下能够采集得到一组30个在时间上离散的数据点，且该30个数据点的本征信噪比一致，从而可以对信噪比数据进行校准。
在步骤S102一种具体的实施方式中，当采用单通道接收线圈时，对多个所述回波EP I图像的原始K空间数据进行处理重建，得到本征成像数据，包括：
利用傅立叶变换对原始K空间的回波数据进行处理重建，得到本征成像数据，其中，所述本征成像数据的矩阵大小与K空间矩阵大小一致。
具体的，原始K空间的回波数据Nx*Ny通过傅立叶变换进行重建得到图像数据Nx’*Ny’，其中，Nx，Ny分别为回波采样数和相位编码数，Nx’*Ny’分别为对应图像列数和行数，且该图像的矩阵大小与k空间矩阵大小一致。
在步骤S102一种具体的实施方式中，当采用多通道接收阵列线圈时，对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据，包括：
利用平方和方法对K空间的三维多通道回波数据进行加权组合重建，得到多通道信息压缩后的本征成像数据。
具体的，对于K空间的三维多通道回波数据Nx*Ny*Nch，Nch为接收线圈通道数目，需要将其进行压维处理，即压缩至二维数据，即通过平方和方法进行加权组合重建，得到本征成像数据。
步骤S103.根据所述本征成像数据计算得到线圈的本征时域稳定性参数，所述本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性。
其中，需要说明的是，所述本征时域信噪比是指仅考虑射频接收线圈和被试样本相对运动，由射频接收线圈灵敏度、热噪声以及二者的时域波动等因素所决定的时域信噪比性能参数，用于评价特定射频接收线圈是否适用于功能磁共振成像这种时域获取应用中，即同时考虑射频接收线圈和被试样本相对运动带来的时域噪声，以及射频接收线圈的灵敏度和热噪声。
在步骤S103一种具体的实施方式中，根据所述本征成像数据计算得到线圈本征时域稳定性参数，包括：
根据使用不同体素分辨率采集得到的本征成像数据的时域数据，分别计算对应不同体素分辨率下的本征时域信噪比tSNR*，包括：
计算每一像素点的时域均值和时域标准差，将所述时域均值作为线圈的本征时域灵敏度均值
并将所述时域标准差作为线圈的本征时域噪声标准差σt *′；
根据所述本征时域灵敏度均值
和所述本征时域噪声标准差σt *′，计算得*到本征时域信噪比tSNR*，其中，
在步骤S103一种具体的实施方式中，在根据所述本征成像数据计算得到线圈的本征时域稳定性参数之前，所述方法还包括：
将在射频电源开启时采集得到的本征成像数据作为线圈本征灵敏度数据S*；
将在射频电源关闭时采集得到的本征成像数据作为热噪声数据，计算得到热噪声标准差
根据所述线圈本征灵敏度数据S*和所述热噪声标准差
计算得到本征信噪比SNR*，
在一种可能的设计中，根据所述本征成像数据计算得到线圈本征时域稳定性参数，包括：
对所述线圈本征灵敏度数据S*的时域标准差进行归一化计算，得到本征时域灵敏度稳定性λ*，
其中，σs *表示线圈本征灵敏度的时域标准差；
对所述热噪声数据的时域标准差进行归一化计算，得到本时域热噪声稳定性α*，
其中，
表示被试仿真模体噪声的时域标准差。
其中，需要说明的是，所述本时域热噪声稳定性α*还可以通过拟合本征信噪比SNR*和tSNR*的关系来得到，为验证模型参数α*是否可以表征热噪声的时域稳定性，对每个接收线圈的仅含噪声数据计算热噪声的归一化时域方差Var，然后通过线性回归分析方法分析其与α*之间的关系，若二者相同或相近，则说明模型参数α*可以表征热噪声的时域稳定性。
基于上述公开的内容，本实施例通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，从而可以排除所有可能存在于功能磁共振采集中复杂的生理噪声，只考虑由于射频接收线圈与被试模体之间的相对运动带来的时域噪声；通过采集在插入不同厚度的介质件时被试仿真模体的多个回波EPI图像，并对多个回波图像进行处理重建，得到本征成像数据；根据本征成像数据计算得到线圈的本征时域稳定性参数，本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性，该本征时域稳定性参数可用于指导功能磁共振成像应用中的专用射频接收线圈的设计与使用，提高功能磁共振成像的图像性能。
第二方面，本发明提供一种fMRI中射频接收线圈本征时域稳定性参数的测量装置，包括：
相对运动模拟模块，用于通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，其中，所述介质件需要满足电导率小于1×10-15S/m，介电常数介于1-3之间；
成像数据采集模块，用于采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，并对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据；
时域参数计算模块，用于根据所述本征成像数据计算得到线圈的本征时域稳定性参数，所述本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性。
本实施例第二方面提供的前述装置的工作过程、工作细节和技术效果，可以参见如上第一方面或第一方面中任意一种可能设计所述的方法，于此不再赘述。
第三方面，本发明提供一种计算机设备，包括依次通信相连的存储器、处理器和收发器，其中，所述存储器用于存储计算机程序，所述收发器用于收发消息，所述处理器用于读取所述计算机程序，执行如第一方面任意一种可能的设计中所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法。
具体举例的，所述存储器可以但不限于包括随机存取存储器(Random-AccessMemory，RAM)、只读存储器(Read-Only Memory，ROM)、闪存(Flash Memory)、先进先出存储器(First Input First Output，FIFO)和/或先进后出存储器(First Input Last Output，FILO)等等；所述处理器可以不限于采用型号为STM32F105系列的微处理器；所述收发器可以但不限于为WiFi(无线保真)无线收发器、蓝牙无线收发器、GPRS(General Packet RadioService，通用分组无线服务技术)无线收发器和/或ZigBee(紫蜂协议，基于IEEE802.15.4标准的低功耗局域网协议)无线收发器等。此外，所述计算机设备还可以但不限于包括有电源模块、显示屏和其它必要的部件。
本实施例第三方面提供的前述计算机设备的工作过程、工作细节和技术效果，可以参见如上第一方面或第一方面中任意一种可能设计所述的方法，于此不再赘述。
第四方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，执行如第一方面任意一种可能的设计中所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法。
其中，所述计算机可读存储介质是指存储数据的载体，可以但不限于包括软盘、光盘、硬盘、闪存、优盘和/或记忆棒(Memory Stick)等，所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。
本实施例第四方面提供的前述计算机设备的工作过程、工作细节和技术效果，可以参见如上第一方面或第一方面中任意一种可能设计所述的方法，于此不再赘述。
第五方面，本发明提供一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行如第一方面任意一种可能的设计中所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法。
本实施例第五方面提供的前述计算机可读存储介质的工作过程、工作细节和技术效果，可以参见如上第一方面或第一方面中任意一种可能设计所述的方法，于此不再赘述。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，包括：
通过在射频接收线圈与被试仿真模体之间分别插入不同厚度的介质件来改变线圈-模体间距，以模拟线圈与模体之间的相对运动，其中，所述介质件需要满足电导率小于1×10-15S/m，相对介电常数介于1-3之间；
采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，并对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据；
根据所述本征成像数据计算得到线圈的本征时域稳定性参数，所述本征时域稳定性参数至少包括本征时域信噪比、本征时域灵敏度稳定性和本征时域热噪声稳定性。
2.根据权利要求1所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，所述射频接收线圈至少包括单通道接收线圈和/或多通道接收阵列线圈。
3.根据权利要求1所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，所述被试仿真模体包括水模体，所述水模体内部设有模拟人体电特性的填充液，所述电特性至少包括介电常数和电导率。
4.根据权利要求1所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，所述介质件包括填充有聚四氟乙烯材料的塑料垫。
5.根据权利要求1所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，采集在插入不同厚度的介质件时所述被试仿真模体的多个回波EPI图像，包括：
在插入每一厚度的介质件时，利用短TE序列对所述被试仿真模体进行质子密度加权成像，设置N个时域采样点/次，采集得到不同体素分辨率下的多组回波EPI图像。
6.根据权利要求2所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，当采用单通道接收线圈时，对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据，包括：
利用傅立叶变换对原始K空间的回波数据进行处理重建，得到本征成像数据，其中，所述本征成像数据的矩阵大小与K空间矩阵大小一致。
7.根据权利要求2所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，当采用多通道接收阵列线圈时，对多个所述回波EPI图像的原始K空间数据进行处理重建，得到本征成像数据，包括：
利用平方和方法对原始K空间的三维多通道回波数据进行加权组合重建，得到多通道信息压缩后的本征成像数据。
8.根据权利要求5所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，根据所述本征成像数据计算得到线圈本征时域稳定性参数，包括：
根据使用不同体素分辨率采集得到的本征成像数据的时域数据，分别计算对应不同体素分辨率下的本征时域信噪比tSNR*，包括：
计算每一像素点的时域均值和时域标准差，将所述时域均值作为线圈的本征时域灵敏度均值
并将所述时域标准差作为线圈的本征时域噪声标准差σt *′；
根据所述本征时域灵敏度均值
和所述本征时域噪声标准差σt *′，计算得到本征时域信噪比tSNR*，其中，
9.根据权利要求8所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，在根据所述本征成像数据计算得到线圈的本征时域稳定性参数之前，所述方法还包括：
将在射频电源开启时采集得到的本征成像数据作为线圈本征灵敏度数据S*；
将在射频电源关闭时采集得到的本征成像数据作为热噪声数据，计算得到热噪声标准差
根据所述线圈本征灵敏度数据S*和所述热噪声标准差
计算得到本征信噪比SNR*，
10.根据权利要求9所述的fMRI中射频接收线圈本征时域稳定性参数的测量方法，其特征在于，根据所述本征成像数据计算得到线圈本征时域稳定性参数，包括：
对所述线圈本征灵敏度数据S*的时域标准差进行归一化计算，得到本征时域灵敏度稳定性λ*，
其中，σs *表示线圈本征灵敏度的时域标准差；
对所述热噪声数据的时域标准差进行归一化计算，得到本时域热噪声稳定性α*，
其中，
表示被试仿真模体噪声的时域标准差。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210594734.1A CN114879107B (zh) | 2022-05-27 | 2022-05-27 | 一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 |
| PCT/CN2022/099606 WO2023226115A1 (zh) | 2022-05-27 | 2022-06-17 | 一种fmri中射频接收线圈本征时域稳定性参数的测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210594734.1A CN114879107B (zh) | 2022-05-27 | 2022-05-27 | 一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114879107A true CN114879107A (zh) | 2022-08-09 |
| CN114879107B CN114879107B (zh) | 2023-01-03 |
Family
ID=82679123
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210594734.1A Active CN114879107B (zh) | 2022-05-27 | 2022-05-27 | 一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN114879107B (zh) |
| WO (1) | WO2023226115A1 (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5552707A (en) * | 1993-08-30 | 1996-09-03 | Hitachi Medical Corporation | RF probe which can vary distance between receiving coils which face each other |
| CN1238463A (zh) * | 1998-06-05 | 1999-12-15 | 通用电器横河医药系统株式会社 | 磁共振成象线圈和磁共振成象装置 |
| US20110080167A1 (en) * | 2009-10-02 | 2011-04-07 | Siemens Aktiengesellschaft | Method for determining motion parameters of an object in a magnetic field |
| CN102028470A (zh) * | 2010-12-23 | 2011-04-27 | 上海辰光医疗科技有限公司 | 用于磁共振成像系统的眼球射频线圈装置 |
| CN102908193A (zh) * | 2011-08-03 | 2013-02-06 | 上海联影医疗科技有限公司 | 一种提高信噪比的四肢定位装置及其四肢接收装置 |
| CN106125029A (zh) * | 2016-06-21 | 2016-11-16 | 中国科学院深圳先进技术研究院 | 多通道磁共振射频线圈性能评估方法及其装置 |
| CN110940945A (zh) * | 2019-12-02 | 2020-03-31 | 浙江大学 | 具有高时域信号稳定性的磁共振成像射频线圈组件 |
| CN113203895A (zh) * | 2021-03-31 | 2021-08-03 | 北京品驰医疗设备有限公司 | 植入式医疗仪器用体外电磁感应线圈的测试系统 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8798719B2 (en) * | 2009-12-02 | 2014-08-05 | Qing X. Yang | Method of utilization of high dielectric constant (HDC) materials for reducing SAR and enhancing SNR in MRI |
| CN101794329B (zh) * | 2009-12-07 | 2011-09-14 | 清华大学 | 一种计算磁共振成像rf线圈信噪比的方法 |
| US20130099786A1 (en) * | 2010-07-02 | 2013-04-25 | Koninklijke Philips Electronics N.V. | Parallel magnetic resonance imaging using undersampled coil data for coil sensitivity estimation |
| CN103536288A (zh) * | 2012-07-12 | 2014-01-29 | 上海联影医疗科技有限公司 | 磁共振系统及其射频线圈性能检测方法 |
| CN102928796B (zh) * | 2012-09-28 | 2014-12-24 | 清华大学 | 快速扩散磁共振成像和重建方法 |
| US10429475B2 (en) * | 2013-03-12 | 2019-10-01 | The General Hospital Corporation | Method for increasing signal-to-noise ratio in magnetic resonance imaging using per-voxel noise covariance regularization |
| CN103698732B (zh) * | 2013-12-12 | 2016-03-09 | 深圳先进技术研究院 | 磁共振射频线圈性能评测方法和系统 |
-
2022
- 2022-05-27 CN CN202210594734.1A patent/CN114879107B/zh active Active
- 2022-06-17 WO PCT/CN2022/099606 patent/WO2023226115A1/zh unknown
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5552707A (en) * | 1993-08-30 | 1996-09-03 | Hitachi Medical Corporation | RF probe which can vary distance between receiving coils which face each other |
| CN1238463A (zh) * | 1998-06-05 | 1999-12-15 | 通用电器横河医药系统株式会社 | 磁共振成象线圈和磁共振成象装置 |
| US20110080167A1 (en) * | 2009-10-02 | 2011-04-07 | Siemens Aktiengesellschaft | Method for determining motion parameters of an object in a magnetic field |
| CN102028470A (zh) * | 2010-12-23 | 2011-04-27 | 上海辰光医疗科技有限公司 | 用于磁共振成像系统的眼球射频线圈装置 |
| CN102908193A (zh) * | 2011-08-03 | 2013-02-06 | 上海联影医疗科技有限公司 | 一种提高信噪比的四肢定位装置及其四肢接收装置 |
| CN106125029A (zh) * | 2016-06-21 | 2016-11-16 | 中国科学院深圳先进技术研究院 | 多通道磁共振射频线圈性能评估方法及其装置 |
| CN110940945A (zh) * | 2019-12-02 | 2020-03-31 | 浙江大学 | 具有高时域信号稳定性的磁共振成像射频线圈组件 |
| CN113203895A (zh) * | 2021-03-31 | 2021-08-03 | 北京品驰医疗设备有限公司 | 植入式医疗仪器用体外电磁感应线圈的测试系统 |
Non-Patent Citations (2)
| Title |
|---|
| YANG GAO ET AL.: "《Intrinsic Temporal Performance of the RF Receive Coil in Magnetic Resonance Imaging》", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
| 徐斌 等: "《非人灵长类动物超高场磁共振脑成像技术进展》", 《中国生物医学工程学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114879107B (zh) | 2023-01-03 |
| WO2023226115A1 (zh) | 2023-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6557710B2 (ja) | 核磁気共鳴（ｎｍｒ）フィンガープリンティング | |
| Jaspan et al. | Compressed sensing MRI: a review of the clinical literature | |
| JP6232043B2 (ja) | Ｋ空間領域の中央および周辺のサンプリング最中での呼吸運動の検出を用いた磁気共鳴画像の再構成 | |
| US10627468B2 (en) | Nuclear magnetic resonance (NMR) fingerprinting | |
| CN100591269C (zh) | 数据校正装置、数据校正方法、磁共振成像装置和x射线ct装置 | |
| CN107656224B (zh) | 一种磁共振成像方法、装置及系统 | |
| CN110346743B (zh) | 一种磁共振弥散加权成像方法和装置 | |
| US20180231626A1 (en) | Systems and methods for magnetic resonance fingerprinting for quantitative breast imaging | |
| Nogueira et al. | Breast DWI at 3 T: influence of the fat-suppression technique on image quality and diagnostic performance | |
| US11047942B2 (en) | Magnetic resonance imaging method and system with optimal variable flip angles | |
| Wang et al. | T1 estimation for aqueous iron oxide nanoparticle suspensions using a variable flip angle SWIFT sequence | |
| Cristobal‐Huerta et al. | Compressed Sensing 3D‐GRASE for faster High‐Resolution MRI | |
| CN110992440A (zh) | 弱监督磁共振快速成像方法和装置 | |
| Utzschneider et al. | Towards accelerated quantitative sodium MRI at 7 T in the skeletal muscle: Comparison of anisotropic acquisition-and compressed sensing techniques | |
| WO2013086951A1 (zh) | 动态对比度增强磁共振成像方法和系统 | |
| CN112370040A (zh) | 磁共振成像方法、装置、存储介质及电子设备 | |
| Kratzer et al. | Sodium relaxometry using 23Na MR fingerprinting: A proof of concept | |
| CN114966507B (zh) | 一种fMRI中射频接收线圈本征时域稳定性评价方法 | |
| CN114879107B (zh) | 一种fMRI中射频接收线圈本征时域稳定性参数的测量方法 | |
| CN113129401B (zh) | 一种参数化磁共振成像的图像重建方法 | |
| Cabini et al. | Fast deep learning reconstruction techniques for preclinical magnetic resonance fingerprinting | |
| CN111681184B (zh) | 一种神经黑色素图像重建方法、装置、设备和存储介质 | |
| US8717023B2 (en) | Relaxometry quantification self-justification fitting | |
| Liu et al. | A three‐dimensional Magnetic Resonance Spin Tomography in Time‐domain protocol for high‐resolution multiparametric quantitative magnetic resonance imaging | |
| CN115137347A (zh) | 三维超短回波时间磁共振指纹成像的髓鞘定量成像方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |