CN117310581A - 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 - Google Patents
一种核磁共振信号衰减拟合方法、系统、设备及存储介质 Download PDFInfo
- Publication number
- CN117310581A CN117310581A CN202311316023.9A CN202311316023A CN117310581A CN 117310581 A CN117310581 A CN 117310581A CN 202311316023 A CN202311316023 A CN 202311316023A CN 117310581 A CN117310581 A CN 117310581A
- Authority
- CN
- China
- Prior art keywords
- water
- diffusion coefficient
- coefficient
- exp
- attenuation
- 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
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000009792 diffusion process Methods 0.000 claims abstract description 114
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 114
- 230000006870 function Effects 0.000 claims description 40
- 238000005315 distribution function Methods 0.000 claims description 37
- 239000004519 grease Substances 0.000 claims description 31
- 238000005259 measurement Methods 0.000 claims description 10
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 claims description 9
- 238000002592 echocardiography Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 abstract description 9
- 238000012545 processing Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 206010061218 Inflammation Diseases 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/443—Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明属于数据处理技术领域，具体公开了一种核磁共振信号衰减拟合方法、系统、设备及存储介质，通过将CPMG序列核磁共振测量得到的由若干个回波构成的时序信号代入基于先验分布函数所构建多指数衰减函数拟合公式中进行拟合计算，得到所需水T1衰减系数、水T2衰减系数和水扩散系数，可以实现更高效、稳定的核磁共振信号衰减拟合。本发明通过加入对被测量物质的先验知识，在不增加变量数量的条件下提高拟合的准确性和稳定度，可以替代现有的使用单一常数求解的衰减拟合方法，有效提升拟合值的稳定性，使得到的衰减系数和扩散系数可以很好地反映真实情况。
Description
技术领域
本发明属于数据处理技术领域，具体涉及一种核磁共振信号衰减拟合方法、系统、设备及存储介质。
背景技术
核磁共振通过拟合电磁时序信号的衰减过程得到测量对象的相关参数(包括T1、T2衰减系数和扩散系数等)，可以达到对测量对象的成分和性质进行非侵入式检测。当对测量对象同时激发的不同成分比例感兴趣时，可以利用不同成分(水和油脂)的衰减系数的差异而引起的信号衰减过程的差异，将接收的全部信号视为不同衰减信号的线性加权组合，通过拟合加权系数得到不同成分的占比。在低分辨率的核磁共振应用场景下，通常把近似性质的物质视作同一类物质(自由水和结合水视作水类，区分于油脂类)，在现有拟合模型中用单一参数代表这一类物质的平均属性。
现有核磁共振信号衰减拟合方法把每一类物质对应的系数视作单一未知常数来求解。然而对于较复杂的人体组织，结合水和自由水的衰减系数实为跨度较大的分布。且在非均匀场的情况下，磁场梯度的空间分布并不均匀，导致对梯度较敏感的水扩散系数在空间上并不相同。使用单一常数求解得到的水衰减系数和扩散系数并不能很好的反映真实值。另一方面，通过将结合水和自由水视作两种物质分别拟合的方法或者使用拓展指数拟合的方法会显著增加变量的数量，降低拟合值的稳定性(拟合值对随机噪音更敏感)。
发明内容
本发明的目的是提供一种核磁共振信号衰减拟合方法、系统、设备及存储介质，用以解决现有技术中存在的上述问题。
为了实现上述目的，本发明采用以下技术方案：
第一方面，提供一种核磁共振信号衰减拟合方法，包括：
构建基于位置参数和尺度参数的扩散系数先验分布函数，并根据扩散系数先验分布函数构建多指数衰减函数拟合公式；
获取通过CPMG序列核磁共振测量得到的由若干个回波构成的时序信号，并根据时序信号确定若干信号值；
将若干信号值分别代入多指数衰减函数拟合公式中，得到多元方程组；
对多元方程组进行求解，得到水T1衰减系数、水T2衰减系数和水扩散系数；
输出水T1衰减系数、水T2衰减系数和水扩散系数。
在一个可能的设计中，所述多指数衰减函数拟合公式为
S＝W*(1-exp(TR/R1，W))*exp(-TE/T2，W)*∫rp(x|Dw，σ)g(x)dx
+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝exp(-TE*b*x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
在一个可能的设计中，所述多指数衰减函数拟合公式为
S＝W*∫xp(x|T1，W，σ)g(x)dx*exp(-TE*(1/T2，F+b*DF))
+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝1-exp(-TR/x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
在一个可能的设计中，所述多指数衰减函数拟合公式为
S＝W*(1-exp(TR/T1，W))*∫xp(x|T2，W，σ)g(x)dx*exp(-TE*b*Dw)
+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝exp(-TE/x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
在一个可能的设计中，所述扩散系数值变量x的取值范围为(0，200]，所述对多元方程组进行求解包括：求得多元方程组在扩散系数值变量x取值范围(0，200]内的近似最优解。
在一个可能的设计中，所述尺度参数σ设定为σ＝2*10-3，所述位置参数为水扩散系数的平均值。
在一个可能的设计中，所述扩散系数先验分布函数采用高斯函数。
第二方面，提供一种核磁共振信号衰减拟合系统，构建单元、获取单元、代入单元、解算单元和输出单元，其中：
构建单元，用于构建基于位置参数和尺度参数的扩散系数先验分布函数，并根据扩散系数先验分布函数构建多指数衰减函数拟合公式；
获取单元，用于获取通过CPMG序列核磁共振测量得到的由若干个回波构成的时序信号，并根据时序信号确定若干信号值；
代入单元，用于将若干信号值分别代入多指数衰减函数拟合公式中，得到多元方程组；
解算单元，用于对多元方程组进行求解，得到水T1衰减系数、水T2衰减系数和水扩散系数；
输出单元，用于输出水T1衰减系数、水T2衰减系数和水扩散系数。
第三方面，提供一种核磁共振信号衰减拟合设备，包括：
存储器，用于存储指令；
处理器，用于读取所述存储器中存储的指令，并根据指令执行上述第一方面中任意一种所述的方法。
第四方面，提供一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，使得所述计算机执行第一方面中任意一种所述的方法。同时，还提供一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行第一方面中任意一种所述的方法。
有益效果：本发明通过将CPMG序列核磁共振测量得到的由若干个回波构成的时序信号代入基于先验分布函数所构建多指数衰减函数拟合公式中进行拟合计算，得到所需水T1衰减系数、水T2衰减系数和水扩散系数，可以实现更高效、稳定的核磁共振信号衰减拟合。本发明通过加入对被测量物质的先验知识，在不增加变量数量的条件下提高拟合的准确性和稳定度，可以替代现有的使用单一常数求解的衰减拟合方法，有效提升拟合值的稳定性，使得到的衰减系数和扩散系数可以很好地反映真实情况。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例1中方法的步骤示意图；
图2为本发明实施例1中函数p在不同设定尺度参数σ下的分布情况；
图3为本发明实施例1中样本中水的扩散系数的估计值统计结果；
图4为本发明实施例1中样本中水的变异系数统计结果；
图5为本发明实施例2中系统的构成示意图；
图6为本发明实施例3中设备的构成示意图。
具体实施方式
在此需要说明的是，对于这些实施例方式的说明用于帮助理解本发明，但并不构成对本发明的限定。本文公开的特定结构和功能细节仅用于描述本发明的示例实施例。然而，可用很多备选的形式来体现本发明，并且不应当理解为本发明限制在本文阐述的实施例中。
应当理解，除非另有明确的规定和限定，术语“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在实施例中的具体含义。
在下面的描述中提供了特定的细节，以便于对示例实施例的完全理解。然而，本领域普通技术人员应当理解可以在没有这些特定细节的情况下实现示例实施例。例如可以在框图中示出系统，以避免用不必要的细节来使得示例不清楚。在其他实施例中，可以不以非必要的细节来示出众所周知的过程、结构和技术，以避免使得实施例不清楚。
实施例1：
本实施例提供一种核磁共振信号衰减拟合方法，可应用于相应的核磁共振信号衰减拟合处理终端，如图1所示，方法包括以下步骤：
S1.构建基于位置参数和尺度参数的扩散系数先验分布函数，并根据扩散系数先验分布函数构建多指数衰减函数拟合公式。
具体实施时，可使用随机变量描述同一类物质扩散系数的分布。随机变量服从的概率函数一般可以用两个参数描述——位置参数(μ)和尺度参数(σ)。测量随机变量X得到观测值x的概率可以表达为Pr(X＝x)＝p(x|μ,σ)。基于位置参数和尺度参数的扩散系数先验分布函数p的位置参数即为水(包括自由水和结合水)的扩散系数平均值。扩散系数先验分布函数p视为水扩散系数值的分布，可采用高斯函数，形成高斯分布，如图2所示，即为扩散系数先验分布函数在不同设定尺度参数σ下的分布情况，横轴为水扩散系数的值，纵轴为具有对应值的水占所有自由水以及结合水的比率。
扩散系数先验分布函数构建多指数衰减函数拟合公式可表示为
S＝W*(1-exp(TR/R1，W))*exp(-TE/T2，W)*∫xp(x|Dw，σ)g(x)dx
+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝exp(-TE*b*x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
或者，所述多指数衰减函数拟合公式表示为
S＝W*∫xp(x|T1，W，σ)g(x)dx*exp(-TE*(1/T2，F+b*DF))
+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝1-exp(-TR/x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
或者，所述多指数衰减函数拟合公式表示为
S＝W*(1-exp(TR/T1，W))*∫xp(x|T2，W，σ)g(x)dx*exp(-TE*b*DW)
+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝exp(-TE/x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
S2.获取通过CPMG序列核磁共振测量得到的由若干个回波构成的时序信号，并根据时序信号确定若干信号值。
具体实施时，可通过CPMG序列(自旋回讯磁振脉冲序列)核磁共振测量得到的由1508个回波构成的时序信号，并根据时序信号确定1508个信号值S，时序信号服从多指数衰减函数拟合公式的指数衰减信号。
S3.将若干信号值分别代入多指数衰减函数拟合公式中，得到多元方程组。
具体实施时，将1508个信号值S分别代入多指数衰减函数拟合公式中，得到相应的多元方程组。在本实施例中油脂T1衰减系数，T1，F，油脂T2衰减系数T2，F，以及油脂扩散系数DF均为根据已有样品中油脂先验知识所设定好的对应常数。TR为设定的重复时间常数，TE为设定的回波时间常数,b为设定的磁场梯度对扩散系数的影响常数。在实际CPMG序列核磁共振测量时，给定TR和TE(16个有不同TR和TE的值的组合)以及b(4个不同值，对应TE)。x的取值范围由被描述物质的性质和描述精度的需求决定，本实施例中使用(0，200]区间内间隔均匀的2000个值，即0.1，0.2，...199.9，200。
S4.对多元方程组进行求解，得到水T1衰减系数、水T2衰减系数和水扩散系数。
具体实施时，求解要先确定水扩散系数函数p所服从的分布和尺度参数，根据已有经验统计，在同一个样品上的20次测量结果可以得到如图3和图4所示的关于水扩散系数统计性质，图3为样本中水扩散系数的估计值，图4为样本中水的变异系数(水扩散系数标准差除以均值，越小说明数值越稳定)。可以看出，在σ＝2*10-3时，得到的水扩散系数值有较好的统计性质(估计值偏移较小，且值较稳定)。因此，可采用σ＝2*10-3，对应的水扩散系数函数p如图2中对应σ的分布所示。基于此，通过解多元方程组1508个等式得到的近似最优解，即得到水T1衰减系数、水T2衰减系数和水扩散系数Dw。
S5.输出水T1衰减系数、水T2衰减系数和水扩散系数。
具体实施时，将求解得到的水T1衰减系数T1，W、水T2衰减系数T2，W和水扩散系数Dw进行输出展示及保存，以便后续利用水T1衰减系数、水T2衰减系数和水扩散系数Dw进行物质性质判断。示例性地，如较高的T1，W和T2，W意味着较多的自由水，在人体组织中可认为有炎症的存在，在岩芯样本中可视为较大的孔隙度，较低的Dw对应较弱的分子扩散，一般对应液体的流动受阻。
本实施例方法通过加入对被测量物质的先验知识，在不增加变量数量的条件下提高拟合的准确性和稳定度，可以替代现有的使用单一常数求解的衰减拟合方法，有效提升拟合值的稳定性，使得到的衰减系数和扩散系数可以很好地反映真实情况。
实施例2：
本实施例提供一种核磁共振信号衰减拟合系统，如图2所示，包括构建单元、获取单元、代入单元、解算单元和输出单元，其中：
构建单元，用于构建基于位置参数和尺度参数的扩散系数先验分布函数，并根据扩散系数先验分布函数构建多指数衰减函数拟合公式；
获取单元，用于获取通过CPMG序列核磁共振测量得到的由若干个回波构成的时序信号，并根据时序信号确定若干信号值；
代入单元，用于将若干信号值分别代入多指数衰减函数拟合公式中，得到多元方程组；
解算单元，用于对多元方程组进行求解，得到水T1衰减系数、水T2衰减系数和水扩散系数；
输出单元，用于输出水T1衰减系数、水T2衰减系数和水扩散系数。
实施例3：
本实施例提供一种核磁共振信号衰减拟合设备，如图5所示，在硬件层面，包括：
数据接口，用于建立处理器与核磁共振测量端的数据对接；
存储器，用于存储指令；
处理器，用于读取所述存储器中存储的指令，并根据指令执行实施例1中的核磁共振信号衰减拟合方法。
可选地，该设备还包括内部总线。处理器与存储器和数据接口可以通过内部总线相互连接，该内部总线可以是ISA(Industry Standard Architecture，工业标准体系结构)总线、PCI(Peripheral Component Interconnect，外设部件互连标准)总线或EISA(Extended Industry Standard Architecture，扩展工业标准结构)总线等。所述总线可以分为地址总线、数据总线、控制总线等。
所述存储器可以但不限于包括随机存取存储器(Random Access Memory，RAM)、只读存储器(Read Only Memory，ROM)、闪存(Flash Memory)、先进先出存储器(First InputFirst Output，FIFO)和/或先进后出存储器(First In Last Out，FILO)等。所述处理器可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(NetworkProcessor，NP)等；还可以是数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field－Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
实施例4：
本实施例提供一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，使得所述计算机执行实施例1中的核磁共振信号衰减拟合方法。其中，所述计算机可读存储介质是指存储数据的载体，可以但不限于包括软盘、光盘、硬盘、闪存、优盘和/或记忆棒(Memory Stick)等，所述计算机可以是通用计算机、专用计算机、计算机网络或者其他可编程系统。
本实施例还提供一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行实施例1中的核磁共振信号衰减拟合方法。其中，所述计算机可以是通用计算机、专用计算机、计算机网络或者其他可编程系统。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种核磁共振信号衰减拟合方法，其特征在于，包括：
构建基于位置参数和尺度参数的扩散系数先验分布函数，并根据扩散系数先验分布函数构建多指数衰减函数拟合公式；
获取通过CPMG序列核磁共振测量得到的由若干个回波构成的时序信号，并根据时序信号确定若干信号值；
将若干信号值分别代入多指数衰减函数拟合公式中，得到多元方程组；
对多元方程组进行求解，得到水T1衰减系数、水T2衰减系数和水扩散系数；
输出水T1衰减系数、水T2衰减系数和水扩散系数。
2.根据权利要求1所述的一种核磁共振信号衰减拟合方法，其特征在于，所述多指数衰减函数拟合公式为S＝W*(1-exp(TR/T1，W))*exp(-TE/T2，W)*∫xp(x|Dw，σ)g(x)dx+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝exp(-TE*b*x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，DW表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
3.根据权利要求1所述的一种核磁共振信号衰减拟合方法，其特征在于，所述多指数衰减函数拟合公式为S＝W*∫xp(x|T1，W，σ)g(x)dx*exp(-TE*(1/T2，F+b*DF))+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝1-exp(-TR/x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数，b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
4.根据权利要求1所述的一种核磁共振信号衰减拟合方法，其特征在于，所述多指数衰减函数拟合公式为
S＝W*(1-exp(TR/T1，W))*∫xp(x|T2，W，σ)g(x)dx*exp(-TE*b*DW)+F*(1-exp(TR/T1，F))*exp(-TE*(1/T2，F+b*DF))
g(x)＝exp(-TE/x)
其中，S表示信号值，T1，W表示水T1衰减系数，T2，W表示水T2衰减系数，Dw表示水扩散系数，p(·)表示扩散系数先验分布函数，x表示扩散系数先验分布函数的扩散系数值变量，σ表示设定的尺度参数，T1，F表示设定的油脂T1衰减系数，T2，F表示设定的油脂T2衰减系数，DF表示设定的油脂扩散系数，TR为设定的重复时间常数，TE为设定的回波时间常数,b为设定的磁场梯度对扩散系数的影响常数，exp(·)表示以自然常数e为底的指数函数，W表示水的权重，F表示油脂的权重。
5.根据权利要求2-4任意一项所述的一种核磁共振信号衰减拟合方法，其特征在于，所述扩散系数值变量x的取值范围为(0，200]，所述对多元方程组进行求解包括：求得多元方程组在扩散系数值变量x取值范围(0，200]内的近似最优解。
6.根据权利要求2-4任意一项所述的一种核磁共振信号衰减拟合方法，其特征在于，所述尺度参数σ设定为σ＝2*10-3，所述位置参数为水扩散系数的平均值。
7.根据权利要求2-4任意一项所述的一种核磁共振信号衰减拟合方法，其特征在于，所述扩散系数先验分布函数采用高斯函数。
8.一种核磁共振信号衰减拟合系统，其特征在于，包括构建单元、获取单元、代入单元、解算单元和输出单元，其中：
构建单元，用于构建基于位置参数和尺度参数的扩散系数先验分布函数，并根据扩散系数先验分布函数构建多指数衰减函数拟合公式；
获取单元，用于获取通过CPMG序列核磁共振测量得到的由若干个回波构成的时序信号，并根据时序信号确定若干信号值；
代入单元，用于将若干信号值分别代入多指数衰减函数拟合公式中，得到多元方程组；
解算单元，用于对多元方程组进行求解，得到水T1衰减系数、水T2衰减系数和水扩散系数；
输出单元，用于输出水T1衰减系数、水T2衰减系数和水扩散系数。
9.一种核磁共振信号衰减拟合设备，其特征在于，包括：
存储器，用于存储指令；
处理器，用于读取所述存储器中存储的指令，并根据指令执行权利要求1-7任意一项所述的核磁共振信号衰减拟合方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，使得所述计算机执行权利要求1-7任意一项所述的核磁共振信号衰减拟合方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311316023.9A CN117310581B (zh) | 2023-10-11 | 2023-10-11 | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311316023.9A CN117310581B (zh) | 2023-10-11 | 2023-10-11 | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117310581A true CN117310581A (zh) | 2023-12-29 |
| CN117310581B CN117310581B (zh) | 2024-05-10 |
Family
ID=89296951
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311316023.9A Active CN117310581B (zh) | 2023-10-11 | 2023-10-11 | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117310581B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5565775A (en) * | 1995-06-23 | 1996-10-15 | Exxon Research And Engineering Company | Producible fluid volumes in porous media determined by pulsed field gradient nuclear magnetic resonance |
| CN1971617A (zh) * | 2005-03-25 | 2007-05-30 | 西门子共同研究公司 | 先验约束的均值平移分析 |
| US20100271021A1 (en) * | 2009-04-27 | 2010-10-28 | Kecheng Liu | Method and apparatus for diffusion tensor magnetic resonance imaging |
| CN104471423A (zh) * | 2012-06-29 | 2015-03-25 | Cr发展公司 | 在组织微细管网络内的水的相对量的量化 |
| CN110786854A (zh) * | 2019-11-05 | 2020-02-14 | 广州互云医院管理有限公司 | 一种水脂混合体系下的反转回复序列t1测量方法 |
| CN111721795A (zh) * | 2020-06-29 | 2020-09-29 | 无锡鸣石峻致医疗科技有限公司 | 一种基于核磁共振系统的物质测量方法及系统 |
| CN111855726A (zh) * | 2019-04-24 | 2020-10-30 | 中国科学院化学研究所 | 一种溶液内部小分子与体系内组分相互作用的核磁检测方法 |
| CN112415454A (zh) * | 2020-11-09 | 2021-02-26 | 无锡鸣石峻致医疗科技有限公司 | 一种adc-t2二维图谱的测量方法、装置、计算机设备及非均匀场磁共振系统 |
| CN112462311A (zh) * | 2020-11-09 | 2021-03-09 | 无锡鸣石峻致医疗科技有限公司 | 一种横向磁化矢量衰减时间常数的校正测量方法、装置、计算机设备及非均匀场磁共振系统 |
| CN112581385A (zh) * | 2020-11-30 | 2021-03-30 | 南方医科大学 | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 |
| CN115452874A (zh) * | 2021-06-08 | 2022-12-09 | 中国石油天然气股份有限公司 | 基于核磁共振对火山岩进行测量的方法、设备及存储介质 |
-
2023
- 2023-10-11 CN CN202311316023.9A patent/CN117310581B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5565775A (en) * | 1995-06-23 | 1996-10-15 | Exxon Research And Engineering Company | Producible fluid volumes in porous media determined by pulsed field gradient nuclear magnetic resonance |
| CN1971617A (zh) * | 2005-03-25 | 2007-05-30 | 西门子共同研究公司 | 先验约束的均值平移分析 |
| US20100271021A1 (en) * | 2009-04-27 | 2010-10-28 | Kecheng Liu | Method and apparatus for diffusion tensor magnetic resonance imaging |
| CN104471423A (zh) * | 2012-06-29 | 2015-03-25 | Cr发展公司 | 在组织微细管网络内的水的相对量的量化 |
| CN111855726A (zh) * | 2019-04-24 | 2020-10-30 | 中国科学院化学研究所 | 一种溶液内部小分子与体系内组分相互作用的核磁检测方法 |
| CN110786854A (zh) * | 2019-11-05 | 2020-02-14 | 广州互云医院管理有限公司 | 一种水脂混合体系下的反转回复序列t1测量方法 |
| CN111721795A (zh) * | 2020-06-29 | 2020-09-29 | 无锡鸣石峻致医疗科技有限公司 | 一种基于核磁共振系统的物质测量方法及系统 |
| US20230157567A1 (en) * | 2020-06-29 | 2023-05-25 | Wuxi Marvel Stone Healthcare Co., Ltd. | Nuclear magnetic resonance system-based substance measurement method and system |
| CN112415454A (zh) * | 2020-11-09 | 2021-02-26 | 无锡鸣石峻致医疗科技有限公司 | 一种adc-t2二维图谱的测量方法、装置、计算机设备及非均匀场磁共振系统 |
| CN112462311A (zh) * | 2020-11-09 | 2021-03-09 | 无锡鸣石峻致医疗科技有限公司 | 一种横向磁化矢量衰减时间常数的校正测量方法、装置、计算机设备及非均匀场磁共振系统 |
| CN112581385A (zh) * | 2020-11-30 | 2021-03-30 | 南方医科大学 | 基于多先验约束扩散峰度成像张量估计方法、介质和设备 |
| CN115452874A (zh) * | 2021-06-08 | 2022-12-09 | 中国石油天然气股份有限公司 | 基于核磁共振对火山岩进行测量的方法、设备及存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| GREBENKOV, DS: "Multiexponential attenuation of the CPMG spin echoes due to a geometrical confinement", JOURNAL OF MAGNETIC RESONANCE, vol. 180, no. 1, 31 May 2006 (2006-05-31), pages 118 - 216 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117310581B (zh) | 2024-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Robson et al. | Comprehensive quantification of signal‐to‐noise ratio and g‐factor for image‐based and k‐space‐based parallel imaging reconstructions | |
| Kiselev et al. | Is the “biexponential diffusion” biexponential? | |
| Nehrke | On the steady‐state properties of actual flip angle imaging (AFI) | |
| Wang et al. | A fast and non-destructive LF-NMR and MRI method to discriminate adulterated shrimp | |
| CN112462311B (zh) | 一种横向磁化矢量衰减时间常数的校正测量方法、装置、计算机设备及非均匀场磁共振系统 | |
| CN112881959A (zh) | 一种用于磁共振成像的梯度涡流补偿方法及系统 | |
| Yokoo et al. | Quantitative R2* MRI of the liver with rician noise models for evaluation of hepatic iron overload: Simulation, phantom, and early clinical experience | |
| Bourne et al. | Biexponential diffusion decay in formalin‐fixed prostate tissue: preliminary findings | |
| WO2022095142A1 (zh) | 一种adc-t2二维图谱的测量方法、装置、计算机设备及非均匀场磁共振系统 | |
| Akhondi‐Asl et al. | Fast myelin water fraction estimation using 2D multislice CPMG | |
| Vegh et al. | Selective channel combination of MRI signal phase | |
| Giovannetti et al. | Classical and lateral skin effect contributions estimation in strip MR coils | |
| Lee et al. | Echo time‐range effects on gradient‐echo based myelin water fraction mapping at 3T | |
| Wiens et al. | Computationally rapid method of estimating signal‐to‐noise ratio for phased array image reconstructions | |
| CN117310581B (zh) | 一种核磁共振信号衰减拟合方法、系统、设备及存储介质 | |
| Kuchel et al. | Analytical solutions and simulations for spin-echo measurements of diffusion of spins in a sphere with surface and bulk relaxation | |
| CN111366600B (zh) | 一种检测水果甜度的方法、装置、系统及存储介质 | |
| Mayer et al. | Fast parallel spiral chemical shift imaging at 3T using iterative SENSE reconstruction | |
| Stinson et al. | Dixon‐type and subtraction‐type contrast‐enhanced magnetic resonance angiography: A theoretical and experimental comparison of SNR and CNR | |
| CN109736793B (zh) | 一种评价储层孔隙结构的方法及设备 | |
| Tran-Gia et al. | Consideration of slice profiles in inversion recovery Look–Locker relaxation parameter mapping | |
| Li et al. | Improved estimation of myelin water fractions with learned parameter distributions | |
| CN115435895A (zh) | 一种基于概率分布函数的本底噪声测量方法 | |
| Negreiros et al. | Ultra low cost analog BIST using spectral analysis | |
| CN115184852A (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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20240415 Address after: Room 1002, 1st Floor, Building 4 (3 # B), Hefei Innovation and Entrepreneurship Park, No. 268 Furong Road, Hefei Economic and Technological Development Zone, Anhui Province, 230000 Applicant after: Anhui Junde Medical Technology Co.,Ltd. Country or region after: China Address before: Room 705, Jinbo block, 8 Hongyi Road, Xinwu District, Wuxi City, Jiangsu Province, 214000 Applicant before: Wuxi mingshijunzhi Medical Technology Co.,Ltd. Country or region before: China |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |