CN112539852B - 基于生物传热模型的磁共振温度成像方法及系统 - Google Patents
基于生物传热模型的磁共振温度成像方法及系统 Download PDFInfo
- Publication number
- CN112539852B CN112539852B CN202011211171.0A CN202011211171A CN112539852B CN 112539852 B CN112539852 B CN 112539852B CN 202011211171 A CN202011211171 A CN 202011211171A CN 112539852 B CN112539852 B CN 112539852B
- Authority
- CN
- China
- Prior art keywords
- heat transfer
- magnetic resonance
- temperature
- biological heat
- transfer model
- 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.)
- Active
Links
- 238000012546 transfer Methods 0.000 title claims abstract description 87
- 238000003384 imaging method Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 33
- 230000008569 process Effects 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000010438 heat treatment Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000012544 monitoring process Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 230000009191 jumping Effects 0.000 claims description 3
- 230000003902 lesion Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000005764 inhibitory process Effects 0.000 abstract description 6
- 210000001519 tissue Anatomy 0.000 description 13
- 238000010586 diagram Methods 0.000 description 8
- 206010028980 Neoplasm Diseases 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000004800 polyvinyl chloride Substances 0.000 description 2
- 238000000015 thermotherapy Methods 0.000 description 2
- 206010020843 Hyperthermia Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036031 hyperthermia Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000415 inactivating effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229920000915 polyvinyl chloride Polymers 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/36—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using magnetic elements, e.g. magnets, coils
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于生物传热模型的磁共振温度成像方法及系统，该方法包括以下步骤：S101，在hybrid算法中引入生物传热模型，修改正则约束项；S102，在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；S103，利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像。根据本发明的基于生物传热模型的磁共振温度成像方法，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
Description
技术领域
本发明涉及生物医学工程技术领域，尤其是涉及一种基于生物传热模型的磁共振温度成像方法及系统。
背景技术
热疗技术作为一种微创甚至无创的新型肿瘤治疗手段，具有死亡率低、复发率低、术后恢复迅速等优点，肿瘤热疗的核心目的为在对肿瘤细胞进行灭活的同时，保障健康组织细胞的安全性。为了实现这一目的，肿瘤热疗过程中需要密切关注被加热区域及周边正常组织的温度变化情况，从而根据不同组织的热耐受度评估手术进度，进行合理手术规划。
磁共振温度成像技术(magnetic resonance temperature imaging,MRTI)由于能够对含水组织进行近实时、全局、无损的温度成像等有点，而且磁共振成像具有高的软组织对比度和成像分辨率，这一技术同时结合了测温及影响引导技术，因而在肿瘤热疗中具有最广泛的运用前景。
在临床治疗的磁共振温度监控过程中，传统的磁共振成像序列需要较长的数据获取时间，在采集过程中温度是变化的，而重建的温度分布仅仅是采集过程中实际温度的拟合值，即时间平均效应，这会大大地降低该技术在临床中的运用前景。为了提升数据获取速度与时间分辨率，磁共振快速成像序列成为了一种选择，但是较高的时间分辨率与成像速度会导致空间分辨率的下降，扫描获得的图像质量较差，信噪比较低，为了解决这一问题，引入了生物传热模型进行辅助。
生物传热(Bio-Heat Transfer,BHT)模型能够利用组织的初始条件、边界条件和热物理学参数，在已知外界热源输入的情况下，对生物组织的温度场变化情况进行模拟，从而预测未来某一时刻组织的温度分布。但是当硬件设施和选用序列一定时，时间分辨率、空间分辨率和图像质量三者在成像过程中不可兼得。因此，上述技术存在改进空间。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此，本发明的一个目的在于提出一种基于生物传热模型的磁共振温度成像方法，该基于生物传热模型的磁共振温度成像方法，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
本发明第二个目的提出了一种采用上述基于生物传热模型的磁共振温度成像方法的系统。
本发明第三个目的提出了一种非临时性可读存储介质。
根据本发明实施例的基于生物传热模型的磁共振温度成像方法，包括以下步骤：
S101，在hybrid算法中引入生物传热模型，修改正则约束项；
S102，在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；
S103，利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像。
根据本发明的基于生物传热模型的磁共振温度成像方法，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，在步骤S101中，修改正则约束项的过程包括：
S201，对病灶组织加热过程进行生物传热仿真模拟，得到测温时刻的模拟温度场；
S202，通过磁共振中温度与PRF之间的强线性关系预测对应的相位变化量；
S203，利用相位变化量对hybrid算法中的正则约束项进行修正，使相位场变化量的范数约束修改为相位场的测量值与预测值之差的范数约束。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，对磁共振温度成像模型的参数采用拟牛顿法，通过求解Jacobbi矩阵实现对生物传热模型中参数的优化处理。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，步骤S102包括：
S301，根据生物传热模型初始参数确定判据条件和阈值；
S302，根据判据条件，计算终止条件值和磁共振测得到温度分布值是否吻合，满足，终止运算，不满足，进入S303；
S303，根据当前生物传热模型的参数值计算出此时的温度分布值，此时的温度分布值和磁共振测得到温度分布值作差，通过差分法计算Jacobbi矩阵；
S304，根据Jacobbi矩阵，利用拟牛顿法更新参数向量，跳转至S302，通过若干次迭代，收敛得到准确参数值。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，步骤S103包括：
S401，用加热探针对病灶组织进行温度变化检测，筛选出温度大于阈值的部分作为监控区域；
S402，使用Kalman滤波器对监控区域进行滤波处理；
S403，将Kalman滤波器对监控区域的滤波过程放置在GPU上进行运算处理。
根据本发明的第二方面的基于生物传热模型的磁共振温度成像系统，采用如第一方面任一种所述的基于生物传热模型的磁共振温度成像方法，包括：
第一模块，用于执行在hybrid算法中引入生物传热模型、修改正则约束项；
第二模块，用于执行在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；
第三模块，用于执行利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像。
进一步地，第一模块包括：
第一模块A部，用于对病灶组织加热过程进行生物传热仿真模拟，得到测温时刻的模拟温度场；
第一模块B部，用于通过磁共振中温度与PRF之间的强线性关系预测对应的相位变化量；
第一模块C部，用于利用相位变化量对hybrid算法中的正则约束项进行修正，使相位场变化量的范数约束修改为相位场的测量值与预测值之差的范数约束。
进一步地，第二模块包括：
第二模块A部，用于根据生物传热模型初始参数确定判据条件和阈值；
第二模块B部，用于根据判据条件，计算终止条件值和磁共振测得到温度分布值是否吻合；
第二模块C部，用于根据当前生物传热模型的参数值计算出此时的温度分布值；
第二模块D部，用于根据Jacobbi矩阵，利用拟牛顿法更新参数向量。
综上，根据本发明第二方面的基于生物传热模型的磁共振温度成像系统，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
根据本发明的第三方面的非临时性可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现第一方面任一种所述的基于生物传热模型的磁共振温度成像方法。所述非临时性可读存储介质与上述的基于生物传热模型的磁共振温度成像方法相对于现有技术所具有的优势相同，在此不再赘述。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是根据本发明实施例的基于生物传热模型的磁共振温度成像方法的示意图；
图2是根据本发明实施例的修改正则约束项的过程示意图；
图3是根据本发明实施例的磁共振温度成像模型对生物传热模型的参数进行优化处理的示意图；
图4是根据本发明实施例的利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合并迭代更新图像的示意图；
图5是根据本发明实施例的基于生物传热模型的磁共振温度成像系统的结构示意图；
图6是根据本发明实施例的第一模块的结构示意图；
图7是根据本发明实施例的第二模块的结构示意图；
图8是根据本发明实施例的多次迭代后的更新参数向量的结果示意图。
附图标记：
10-基于生物传热模型的磁共振温度成像系统，1-第一模块，11-第一模块A部，12-第一模块B部，13-第一模块C部，2-第二模块，21-第二模块A部，22-第二模块B部，23-第二模块C部，24-第二模块D部，3-第三模块。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
下面参考图1-图8描述根据本发明实施例的基于生物传热模型的磁共振温度成像方法。如图1所示，根据本发明实施例的基于生物传热模型的磁共振温度成像方法，包括以下步骤：
S101，在hybrid(混合)算法中引入生物传热模型，修改正则约束项；进一步地，采用hybrid(混合)算法不仅能够有效抑制成像过程中的运动伪影问题，还能够解决主磁场随时间的变化带来的磁场漂移积累问题。在现有技术中，hybrid(混合)算法受温度的影响很大，即，hybrid(混合)算法在重建的温度分布中，如某区域温度升高值较小，则该区域的温度升高值会被抑制，从而导致对取值的误判，而在本发明中，hybrid(混合)算法通过引入生物传热模型，可以有效地修改正则约束项，从而提升，hybrid(混合)算法取值的准确性。
S102，在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；
S103，利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像。进一步地，在临床温度监控中往往同时需要较高的时间分辨率以及空间分辨率，本发明通过引入Kalman滤波器，利用生物传热模型模型对磁共振测量得到的温度图进行滤波，最后重建出信噪比和测温精度更高的温度图。
根据本发明的基于生物传热模型的磁共振温度成像方法，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，如图2所示，在步骤S101中，修改正则约束项的过程包括：
S201，对病灶组织加热过程进行生物传热仿真模拟，得到测温时刻的模拟温度场；
S202，通过磁共振中温度与PRF(pulse repetition frequency脉冲重复频率)之间的强线性关系预测对应的相位变化量；
S203，利用相位变化量对hybrid算法中的正则约束项进行修正，使相位场变化量的范数约束修改为相位场的测量值与预测值之差的范数约束。
进一步地，在现有技术中，在非加热区域，生物传热模型模拟的温度为0，此时hybrid目标方程保持不变，hybrid(混合)算法判断是准确的。但是当在加热温度上升较少的区域，当hybrid(混合)算法误判该区域的温度为0时，这样就会导致实际测量的温度值发生偏差。在本发明中，通过生物传热模型的修正，可以使得hybrid(混合)算法对加热温度上升较少的区域作出正确的温度判断，从而使得实际的测量温度值更加准确，有利于降低误差率。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，对磁共振温度成像模型的参数采用拟牛顿法，并通过求解Jacobbi矩阵实现对生物传热模型中参数的优化处理。这样有利于提升预测的准确度。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，如图3所示，步骤S102包括：
S301，根据生物传热模型初始参数确定判据条件和阈值；
S302，根据判据条件，计算终止条件值和磁共振测得到温度分布值是否吻合，满足，终止运算，不满足，进入S303；
S303，根据当前生物传热模型的参数值计算出此时的温度分布值，此时的温度分布值和磁共振测得到温度分布值作差，通过差分法计算Jacobbi矩阵；
S304，根据Jacobbi矩阵，利用拟牛顿法更新参数向量，跳转至S302，通过若干次迭代，收敛得到准确参数值。例如，在一个具体的实施例中，如图8所示，通过三到四次的迭代即可收敛到准确参数值。需要说明的是，PVC即聚氯乙烯(Polyvinyl chloride)。
根据本发明一个实施例的基于生物传热模型的磁共振温度成像方法，如图4所示，步骤S103包括：
S401，用加热探针对病灶组织进行温度变化检测，筛选出温度大于阈值的部分作为监控区域；
S402，使用Kalman滤波器对监控区域进行滤波处理；
S403，将Kalman滤波器对监控区域的滤波过程放置在GPU(Graphics ProcessingUnit图形处理器)上进行运算处理。例如，在一个具体的实施例中，滤波过程通过放置在GPU上进行运算处理，成像时长可降低到每一帧需要五秒的时间，从而有利于满足临床手术需求。
综上所述，根据本发明的基于生物传热模型的磁共振温度成像方法，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
根据本发明的第二方面的基于生物传热模型的磁共振温度成像系统10，采用如第一方面任一种的基于生物传热模型的磁共振温度成像方法，如图5所示，包括：
第一模块1，用于执行在hybrid算法中引入生物传热模型、修改正则约束项；
第二模块2，用于执行在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；
第三模块3，用于执行利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像。
进一步地，如图6所示，第一模块1包括：
第一模块A部11，用于对病灶组织加热过程进行生物传热仿真模拟，得到测温时刻的模拟温度场；
第一模块B部12，用于通过磁共振中温度与PRF之间的强线性关系预测对应的相位变化量；
第一模块C部13，用于利用相位变化量对hybrid算法中的正则约束项进行修正，使相位场变化量的范数约束修改为相位场的测量值与预测值之差的范数约束。
进一步地，如图7所示，第二模块2包括：
第二模块A部21，用于根据生物传热模型初始参数确定判据条件和阈值；
第二模块B部22，用于根据判据条件，计算终止条件值和磁共振测得到温度分布值是否吻合；
第二模块C部23，用于根据当前生物传热模型的参数值计算出此时的温度分布值；
第二模块D部24，用于根据Jacobbi矩阵，利用拟牛顿法更新参数向量。
综上，根据本发明第二方面的基于生物传热模型的磁共振温度成像系统10，根据本发明的基于生物传热模型的磁共振温度成像方法，一方面能够保留hybrid算法对于运动以及温漂的抑制能力；另一方面还能够大大提升温度变化较小区域的重建误差，增加hybrid算法的适用性，提升预测的准确度。
本发明第三方面还提供了一种非临时性可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本发明第一方面的基于生物传热模型的磁共振温度成像方法，从而具有抑制运动以及温漂的能力强、误差率更低和预测的准确度高等优点。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
Claims (5)
1.一种基于生物传热模型的磁共振温度成像方法，其特征在于，包括以下步骤：
S101，在hybrid算法中引入生物传热模型，修改正则约束项；
S102，在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；
S103，利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像；
在步骤S101中，修改正则约束项的过程包括：
S201，对病灶组织加热过程进行生物传热仿真模拟，得到测温时刻的模拟温度场；
S202，通过磁共振中温度与脉冲重复频率(PRF)之间的强线性关系预测对应的相位变化量；
S203，利用相位变化量对hybrid算法中的正则约束项进行修正，使相位场变化量的范数约束修改为相位场的测量值与预测值之差的范数约束；
对磁共振温度成像模型的参数采用拟牛顿法，通过求解Jacobbi矩阵实现对生物传热模型中参数的优化处理；
步骤S102包括：
S301，根据生物传热模型初始参数确定判据条件和阈值；
S302，根据判据条件，计算终止条件值和磁共振测得的温度分布值是否吻合，满足，终止运算，不满足，进入S303；
S303，根据当前生物传热模型的参数值计算出此时的温度分布值，此时的温度分布值和磁共振测得的温度分布值作差，通过差分法计算Jacobbi矩阵；
S304，根据Jacobbi矩阵，利用拟牛顿法更新参数向量，跳转至S302，通过若干次迭代，收敛得到准确参数值；
步骤S103包括：
S401，用加热探针对病灶组织进行温度变化检测，筛选出温度大于阈值的部分作为监控区域；
S402，使用Kalman滤波器对监控区域进行滤波处理；
S403，将Kalman滤波器对监控区域的滤波过程放置在图形处理器(GPU)上进行运算处理。
2.一种基于生物传热模型的磁共振温度成像系统，其特征在于，该系统采用根据权利要求1所述的基于生物传热模型的磁共振温度成像方法，包括：
第一模块，用于执行在hybrid算法中引入生物传热模型、修改正则约束项；
第二模块，用于执行在温度分布计算过程中利用磁共振温度成像模型对生物传热模型的参数进行优化处理与实时更新；
第三模块，用于执行利用Kalman滤波器将生物传热模型与磁共振温度成像模型进行结合，迭代更新图像。
3.根据权利要求2所述的基于生物传热模型的磁共振温度成像系统，其特征在于，第一模块包括：
第一模块A部，用于对病灶组织加热过程进行生物传热仿真模拟，得到测温时刻的模拟温度场；
第一模块B部，用于通过磁共振中温度与脉冲重复频率(PRF)之间的强线性关系预测对应的相位变化量；
第一模块C部，用于利用相位变化量对hybrid算法中的正则约束项进行修正，使相位场变化量的范数约束修改为相位场的测量值与预测值之差的范数约束。
4.根据权利要求3所述的基于生物传热模型的磁共振温度成像系统，其特征在于，第二模块包括：
第二模块A部，用于根据生物传热模型初始参数确定判据条件和阈值；
第二模块B部，用于根据判据条件，计算终止条件值和磁共振测得的温度分布值是否吻合；
第二模块C部，用于根据当前生物传热模型的参数值计算出此时的温度分布值；
第二模块D部，用于根据Jacobbi矩阵，利用拟牛顿法更新参数向量。
5.一种非临时性可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现根据权利要求1所述的基于生物传热模型的磁共振温度成像方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011211171.0A CN112539852B (zh) | 2020-11-03 | 2020-11-03 | 基于生物传热模型的磁共振温度成像方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011211171.0A CN112539852B (zh) | 2020-11-03 | 2020-11-03 | 基于生物传热模型的磁共振温度成像方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112539852A CN112539852A (zh) | 2021-03-23 |
| CN112539852B true CN112539852B (zh) | 2022-09-27 |
Family
ID=75013694
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011211171.0A Active CN112539852B (zh) | 2020-11-03 | 2020-11-03 | 基于生物传热模型的磁共振温度成像方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112539852B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101507603A (zh) * | 2008-10-14 | 2009-08-19 | 清华大学 | 一种磁共振测温的方法和装置 |
| DE102010039737A1 (de) * | 2010-08-25 | 2012-03-01 | Siemens Aktiengesellschaft | Verfahren zur Ermittlung eines die Temperaturverteilung innerhalb eines zu vermessenden Objekts beschreibenden Untersuchungsbilddatensatzes und Magnetresonanzeinrichtung |
| CN106618571A (zh) * | 2016-11-16 | 2017-05-10 | 深圳先进技术研究院 | 一种磁共振成像方法和系统 |
| CN110140042A (zh) * | 2016-05-31 | 2019-08-16 | Q生物公司 | 张量场映射 |
-
2020
- 2020-11-03 CN CN202011211171.0A patent/CN112539852B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101507603A (zh) * | 2008-10-14 | 2009-08-19 | 清华大学 | 一种磁共振测温的方法和装置 |
| DE102010039737A1 (de) * | 2010-08-25 | 2012-03-01 | Siemens Aktiengesellschaft | Verfahren zur Ermittlung eines die Temperaturverteilung innerhalb eines zu vermessenden Objekts beschreibenden Untersuchungsbilddatensatzes und Magnetresonanzeinrichtung |
| CN110140042A (zh) * | 2016-05-31 | 2019-08-16 | Q生物公司 | 张量场映射 |
| CN106618571A (zh) * | 2016-11-16 | 2017-05-10 | 深圳先进技术研究院 | 一种磁共振成像方法和系统 |
Non-Patent Citations (1)
| Title |
|---|
| 基于Kalman滤波和生物传热模型的实时磁共振温度成像精度提升;吴锦超等;《清华大学学报(自然科学版)》;20191025;第60卷(第4期);第334-340页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112539852A (zh) | 2021-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110364238B (zh) | 基于机器学习的对比剂管理 | |
| US10610302B2 (en) | Liver disease assessment in medical imaging | |
| CN104244818B (zh) | 在非侵入式治疗期间的基于参考的运动跟踪 | |
| Roney et al. | Variability in pulmonary vein electrophysiology and fibrosis determines arrhythmia susceptibility and dynamics | |
| JP6027546B2 (ja) | 医用画像診断装置及び医用画像診断装置を用いた位相決定方法 | |
| JP5865522B2 (ja) | 冷却期間中に取得された磁気共鳴データを用いた処置計画の修正 | |
| US11747424B2 (en) | Magnetic resonance imaging apparatus, image processing apparatus, and image processing method | |
| JP2006198407A (ja) | 磁気共鳴（ｍｒ）イメージングにおける動き補償方法および動き補償システム | |
| US9715726B2 (en) | Method and system for B0 drift and respiratory motion compensation in echo-planar based magnetic resonance imaging | |
| Corrado et al. | Quantifying atrial anatomy uncertainty from clinical data and its impact on electro-physiology simulation predictions | |
| Mesin et al. | Semi-automated tracking and continuous monitoring of inferior vena cava diameter in simulated and experimental ultrasound imaging | |
| US9320488B2 (en) | Method and device for correction of movement artifacts in a computed tomography image | |
| CN110751187A (zh) | 异常区域图像生成网络的训练方法和相关产品 | |
| CN112539852B (zh) | 基于生物传热模型的磁共振温度成像方法及系统 | |
| US20140185898A1 (en) | Image generation method and apparatus | |
| CN116802695A (zh) | 生成解剖腔体的解剖模型 | |
| JP2006026076A (ja) | 磁気共鳴イメージング装置 | |
| Kim et al. | Motion-robust, multi-slice, real-time MR thermometry for MR-guided thermal therapy in abdominal organs | |
| US11906608B2 (en) | Automated detection of water-fat swaps in Dixon magnetic resonance imaging | |
| WO2022144263A1 (en) | Quality measure for a mapping function | |
| Noorda et al. | Subject-specific liver motion modeling in MRI: a feasibility study on spatiotemporal prediction | |
| JP2022534031A (ja) | 運動アーチファクト検出の方法 | |
| JP6873134B2 (ja) | Ｓｅｎｓｅイメージングにおける画像アーチファクトの除去 | |
| Gifford et al. | Tests of scanning model observers for myocardial SPECT imaging | |
| US20240177270A1 (en) | Image processing apparatus and method |
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 |