CN109541511A - 温升预测方法和装置 - Google Patents
温升预测方法和装置 Download PDFInfo
- Publication number
- CN109541511A CN109541511A CN201811431311.8A CN201811431311A CN109541511A CN 109541511 A CN109541511 A CN 109541511A CN 201811431311 A CN201811431311 A CN 201811431311A CN 109541511 A CN109541511 A CN 109541511A
- Authority
- CN
- China
- Prior art keywords
- electrode
- image
- sequence
- particular range
- rms
- 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 60
- 238000007689 inspection Methods 0.000 claims abstract description 45
- 238000003384 imaging method Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000006698 induction Effects 0.000 abstract description 29
- 206010037660 Pyrexia Diseases 0.000 abstract description 27
- 238000002595 magnetic resonance imaging Methods 0.000 description 45
- 238000010586 diagram Methods 0.000 description 8
- 210000000746 body region Anatomy 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 229920001817 Agar Polymers 0.000 description 5
- 239000008272 agar Substances 0.000 description 5
- 238000004070 electrodeposition Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000005611 electricity Effects 0.000 description 4
- 230000005311 nuclear magnetism Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 210000004705 lumbosacral region Anatomy 0.000 description 3
- 239000013307 optical fiber Substances 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 208000018737 Parkinson disease Diseases 0.000 description 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 239000000499 gel Substances 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 229920003023 plastic Polymers 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 201000006517 essential tremor Diseases 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005404 monopole Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 229920002635 polyurethane Polymers 0.000 description 1
- 239000004814 polyurethane Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本公开提供一种温升预测方法和装置,用于预测电极在磁共振成像检查、即MRI检查下的温升,所述温升预测方法包括:获取在所述MRI检查下所述电极周围的特定范围的B1图像;将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以计算所述电极的接收系数;确定在所述MRI检查下要扫描的序列,并确定所述序列的射频场的大小、即B1+rms;根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场;以及根据所述特定范围内的B1场来预测所述电极的温升。根据本公开的实施例提供的温升预测方法和装置,可以在考虑到电极的形状、角度和放置位置等的情况下更加精确地表征RF感应发热,从而确定具备电极的医疗设备在MRI检查下的RF感应发热。
Description
技术领域
本公开涉及医疗设备领域,尤其涉及一种温升(温度升高)预测方法和装置。
背景技术
诸如DBS(Deep Brain Stimulation,深部脑刺激)等的IMD(Implantable MedicalDevice,植入式医疗设备)是一种新兴的运动障碍治疗方法。目前,IMD广泛用于PD(Parkinson’s disease,帕金森病)、肌张力障碍和特发性震颤等疾病。随着MRI(MagneticResonance Imaging,磁共振成像)和越来越多的植入式医疗设备的发展,使用可植入设备的MRI检查的数量正在迅速增加。然而,DBS在MRI中应用存在一定的风险,其中RF(RadioFrequency,射频)感应发热是主要潜在风险之一。
目前,SAR(Specific Absorption Rate,比吸收率)被用作表征RF感应发热的主要指标,但是不存在用于估计SAR的通用标准,并且不同MRI厂商的SAR估计方法也有所不同。另一方面,SAR的计算非常保守,对采用的序列和参数限制较大,从而导致植入诸如起搏器等的植入式医疗设备的患者无法进行常规的核磁扫描。此外,B1+rms是用于在MRI中激发质子共振的平均RF磁场,并表征RF场的大小,并且与患者状态无关,因此目前被视为RF感应发热的新指标。然而,该指标在使用中仍然受到很多限制,这是因为B1+rms是基于理想信号源的假设而忽略了患者对B1场的干扰。此外,植入式医疗设备中的电极的形状、角度和放置位置等因素都可能影响RF感应发热,然而,由这些因素所产生的影响不能用SAR和B1+rms表征。因此,利用目前的SAR和B1+rms指标有时不能准确地表征RF感应发热。
发明内容
有鉴于此,在一个方面,本公开提出了一种温升预测方法,用于预测电极在磁共振成像检查、即MRI检查下的温升,所述温升预测方法包括:获取在所述MRI检查下所述电极周围的特定范围的B1图像;将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以确定所述电极的接收系数;确定在所述MRI检查下要扫描的序列,并确定所述序列的射频场的大小、即B1+rms;根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场;以及根据所述特定范围内的B1场来预测所述电极的温升。
对于上述温升预测方法,在一种可能的实现方式中,所述背景射频场的B1图像是以下B1图像中任一:在无所述电极的情况下的所述特定范围的B1图像;以及与所述电极的距离在特定范围以外的范围的B1图像。
对于上述温升预测方法,在一种可能的实现方式中,获取在所述MRI检查下所述电极周围的特定范围的B1图像,包括:获取所述电极的尖端部分的特定范围内的B1图像。
对于上述温升预测方法,在一种可能的实现方式中,确定所述序列的B1+rms包括:在所述序列为一个的情况下,确定该序列的B1+rms;以及在所述序列为多个的情况下,确定多个序列的B1+rms的最大值或平均值。
对于上述温升预测方法,在一种可能的实现方式中,将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以计算所述电极的接收系数包括以下方式中任一:将所述B1图像中的B1场与所述背景射频场进行对比以得到所述B1场与所述背景射频场的比值,并将所得到的比值中大于第一阈值的比值的均方值作为所述电极的接收系数,或者将所得到的比值中小于第二阈值的比值的均方值的倒数作为所述电极的接收系数,其中所述第二阈值小于所述第一阈值;以及将所述B1图像中的B1场与所述背景射频场进行对比,以将所述B1场的均方值与所述背景射频场的均方值的比值计算为所述电极的接收系数。
对于上述温升预测方法,在一种可能的实现方式中,根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场,包括:将所述接收系数与所述序列的B1+rms相乘,以确定所述特定范围内的B1场。
对于上述温升预测方法,在一种可能的实现方式中,根据所述特定范围内的B1场来预测所述电极的温升,包括:对所述特定范围内的B1场求平方;以及根据对所述特定范围内的B1场求平方所得到的值来预测所述电极的温升。
对于上述温升预测方法,在一种可能的实现方式中,所述电极位于植入式医疗设备内。
对于上述温升预测方法,在一种可能的实现方式中,在获取在所述MRI检查下所述电极周围的特定范围的B1图像之前,还包括:获取所述电极的位置和方向。
对于上述温升预测方法,在一种可能的实现方式中,在根据所述特定范围内的B1场来预测所述电极的温升之后,还包括:根据所预测出的所述电极的温升来判断是否发出警报提醒。
在另一方面,本公开提出了一种温升预测装置,用于预测电极在磁共振成像检查、即MRI检查下的温升,所述温升预测装置包括:B1图像获取单元,用于获取在所述MRI检查下所述电极周围的特定范围的B1图像;接收系数确定单元,用于将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以确定所述电极的接收系数;B1+rms确定单元,用于确定在所述MRI检查下要扫描的序列,并确定所述序列的射频场的大小、即B1+rms;B1场确定单元,用于根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场;以及预测单元,用于根据所述特定范围内的B1场来预测所述电极的温升。
对于上述温升预测装置,在一种可能的实现方式中,所述背景射频场的B1图像是以下B1图像中任一:在无所述电极的情况下的所述特定范围的B1图像;以及与所述电极的距离在特定范围以外的范围的B1图像。
对于上述温升预测装置,在一种可能的实现方式中,所述B1图像获取单元用于:获取所述电极的尖端部分的特定范围内的B1图像。
对于上述温升预测装置,在一种可能的实现方式中,所述B1+rms确定单元用于:在所述序列为一个的情况下,确定该序列的B1+rms;以及在所述序列为多个的情况下,确定多个序列的B1+rms的最大值或平均值。
对于上述温升预测装置,在一种可能的实现方式中,所述接收系数计算单元用于进行以下操作中任一:
将所述B1图像中的B1场与所述背景射频场进行对比以得到所述B1场与所述背景射频场的比值,并将所得到的比值中大于第一阈值的比值的均方值作为所述电极的接收系数,或者将所得到的比值中小于第二阈值的比值的均方值的倒数作为所述电极的接收系数,其中所述第二阈值小于所述第一阈值;以及将所述B1图像中的B1场与所述背景射频场进行对比,以将所述B1场的均方值与所述背景射频场的均方值的比值计算为所述电极的接收系数。
对于上述温升预测装置,在一种可能的实现方式中,所述B1场确定单元用于:将所述接收系数与所述序列的B1+rms相乘,以确定所述特定范围内的B1场。
对于上述温升预测装置,在一种可能的实现方式中,所述预测单元用于:对所述特定范围内的B1场求平方;以及根据对所述特定范围内的B1场求平方所得到的值来预测所述电极的温升。
对于上述温升预测装置,在一种可能的实现方式中,所述电极位于植入式医疗设备内。
对于上述温升预测装置,在一种可能的实现方式中,还包括:电极位置和方向获取单元,用于获取所述电极的位置和方向。
对于上述温升预测装置,在一种可能的实现方式中,还包括:警报单元,用于根据所预测出的所述电极的温升来判断是否发出警报提醒。
根据本公开的温升预测方法和装置,通过根据接收系数K和序列的B1+rms来计算电极的特定范围的B1场、然后根据该B1场来预测RF感应发热情况,由此可以在考虑到电极的形状、角度和放置位置等的情况下预测RF感应发热。因此,无论电极的形状、角度和放置位置等如何,应用本公开的上述实施例都可以更加精确地表征RF感应发热,从而确定具备电极的医疗设备在MRI检查下的RF感应发热。
根据下面参考附图对示例性实施例的详细说明,本公开的其它特征及方面将变得清楚。
附图说明
包含在说明书中并且构成说明书的一部分的附图与说明书一起示出了本公开的示例性实施例、特征和方面,并且用于解释本公开的原理。
图1示出根据本公开一实施例的温升预测方法的流程图。
图2示出根据本公开的另一实施例的温升预测方法的流程图。
图3示出将电极与测量温度升高的光纤并排捆绑的示意图。
图4示出将对局部B1场求平方得到的值与温度升高进行拟合而得到的直线的示意图。
图5示出根据本公开一实施例的针对温升预测方法的有限元数值模拟的示意图。
图6示出根据本公开一实施例的针对温升预测方法的有限元数值模拟的结果的示意图。
图7示出根据本公开一实施例的温升预测装置的结构框图。
图8示出根据本公开的另一实施例的温升预测装置的结构框图。
具体实施方式
以下将参考附图详细说明本公开的各种示例性实施例、特征和方面。附图中相同的附图标记表示功能相同或相似的元件。尽管在附图中示出了实施例的各种方面,但是除非特别指出,不必按比例绘制附图。
在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。
另外,为了更好的说明本公开,在下文的具体实施方式中给出了众多的具体细节。本领域技术人员应当理解,没有某些具体细节,本公开同样可以实施。在一些实例中,对于本领域技术人员熟知的方法、手段、元件和电路未作详细描述,以便于凸显本公开的主旨。
图1示出根据本公开一实施例的温升预测方法的流程图。其中该温升预测方法主要用于预测电极在MRI检查下的温升。如图1所示,该温升预测方法包括:
步骤S101、获取在MRI检查下电极周围的特定范围的B1图像;
步骤S102、将该B1图像与在MRI检查下的背景射频场的B1图像进行对比,以计算电极的接收系数K;
步骤S103、确定在MRI检查下要扫描的序列,并确定序列的B1+rms;
步骤S104、根据接收系数K和序列的B1+rms来确定特定范围内的B1场;以及
步骤S105、根据特定范围内的B1场来预测电极的温升。
以下将具体阐述上述各步骤。
对于上述步骤S101,首先,例如对患者进行MRI检查,以获取电极周围的B1图像。在一种可能的实现方式中,电极位于植入式医疗设备中。该植入式医疗设备是诸如为心脏起搏器等的医疗设备,并且可被植入患者体内的相应位置处,以对患者进行相应的治疗或者辅助患者器官功能的正常运行。此外,还可以直接获取之前对该患者的扫描测试结果,由此可以从该扫描测试结果中获取到电极周围的B1图像。
在一种可能的实现方式中,由于电极的RF感应发热主要集中在电极的尖端部分,因此在步骤S101中主要获取电极的尖端部分周围的B1图像。
在获取电极周围的B1图像之后,从该B1图像中提取特定范围的B1图像(换句话说,局部B1图像)。该特定范围可以是电极周围的长方体区域,该长方体区域的长度、宽度、高度例如分别为50mm、20mm和20mm。优选地,该长方体区域的长度的范围例如为10mm~80mm,该长方体区域的宽度的范围例如为6mm~40mm,以及该长方体区域的高度的范围例如为6mm~40mm。此外,该特定范围还可以是上述长方体的某一个截面。以上以示例方式说明了特定范围,但是本公开的特定范围不限于此,并且该特定范围具体可以根据应用场景来适当设置。
由此,可以获取到电极周围的特定范围内的B1图像。
对于上述步骤S102,可以将所获取到的B1图像与标准无干扰的图像进行对比,以得到电极的接收系数K。该标准无干扰的图像可以是受电极的电影较小的B1图像,例如是在MRI检查下的背景射频场的B1图像。在一种可能的实现方式中,背景射频场的B1图像是以下B1图像中任一:在无电极的情况下的特定范围的B1图像;以及与电极的距离在特定范围以外的范围的B1图像。其中,在该特定范围以外的范围的B1图像受到电极的影响较小。具体而言,在步骤S102中,预先获取受电极影响较小的标准无干扰的图像。然后,将上述步骤S101中所获取到的特定范围的B1图像与该标准无干扰的图像进行对比,以得到电极的接收系数K。
在一种可能的实现方式中,可以利用如下方式来计算电极的接收系数K。可以将电极周围B1场信噪比大于第一阈值或小于第二阈值的位置的B1场均方根作为电极的接收系数。具体而言,可以将B1图像中的B1场与背景射频场进行对比,以得到B1场与背景射频场的比值(即,信噪比),并将所得到的比值中大于第一阈值的比值的均方值作为电极的接收系数K,或者将所得到的比值中小于第二阈值的比值的均方值的倒数作为电极的接收系数K,其中第二阈值小于第一阈值。例如,第一阈值可以为1.2,以及第二阈值可以为0.5。然而,第一阈值和第二阈值的具体数值不限于此,其可以根据具体的应用场景而设置成其它适当的值。
在一种可能的实现方式中,还可以利用如下方式来计算电极的接收系数K。具体而言,将B1图像中的B1场与背景射频场进行对比,以将B1场的均方值与背景射频场的均方值的比值计算为电极的接收系数K。
需要说明的是,以上列举了两种计算接收系数K的具体方式,但是本领域技术人员能够知道,根据电极周围的B1图像与在MRI检查下的背景射频场的B1图像,可以采用任何其它适当的方式来计算电极的接收系数K,只要所计算出的接收系数能够适当地反映电极的形状、角度和放置位置对RF感应发热的影响即可。
由于电极的形状、角度和放置位置等因素均会影响具有该电极的医疗设备在MRI检查下的RF感应发热,因此通过将电极周围的B1图像与标准无干扰的B1图像相比较、由此计算出接收系数K,可以考虑到这些因素对电极的RF感应发热的影响。这样,可以在后续步骤中应用接收系数K来精确地预测RF感应发热。
对于上述步骤S103,确定患者在此次MRI检查中要扫描的序列,例如T1、T2、BOLD序列等。然后,对MRI检查设备进行适当设置,以根据该设置来对患者进行MRI检查,并获取(读出)所要扫描的序列的射频场的大小,即B1+rms。
在要扫描的序列是一个的情况下,可以获取该序列的B1+rms作为本步骤中需要的B1+rms。例如,在要扫描T1序列的情况下,可以获取T1序列的B1+rms。例如,在要扫描T2序列的情况下,可以获取T2序列的B1+rms。又如,在要扫描BOLD序列的情况下的情况下,可以获取BOLD序列的B1+rms。
另一方面,在要扫描的序列是多个的情况下,首先可以分别获取这些序列的B1+rms,然后取这些序列中的B1+rms的平均值或最大值作为本步骤中需要的B1+rms。例如,在要扫描T1序列和T2序列的情况下,可以分别获取T1序列和T2序列的B1+rms,然后取所得到的B1+rms的平均值或最大值作为本步骤中需要的B1+rms。优选地,取这些序列中的B1+rms的最大值、即max(B1+rms)作为本步骤中需要的B1+rms,这样可以最大限度地降低患者在MRI检查中的风险。
此外,还可以预先读取各序列的B1+rms。然后,将各序列的B1+rms的平均值和/或最大值存储在存储器中。这样,在需要获取序列的B1+rms时,可以直接从存储器中读取序列的平均值和/或最大值。
需要说明的是,以上通过示例的方式列举了T1序列、T2序列和BOLD序列等几种序列以及它们的组合,但是本公开不限于此,序列可以是患者体内的任意其它序列,并且所要扫描的序列可以是任意序列的任意组合。
由此,通过上述步骤S103,可以得到序列的B1+rms。
对于上述步骤S104,可以根据步骤S102中所计算出的电极的接收系数K和步骤S103中所得到的序列的B1+rms来计算特定范围的B1场。
在一种可能的实现方式中,可以将电极的接收系数K与序列的B1+rms相乘来计算出特定范围的B1场(local B1),即local B1=K*B1+rms。
对于上述步骤S105中,可以根据步骤S104中所计算出的特定范围的B1场来预测电极的温度升高,从而预测电极的RF感应发热情况。
本发明人经过多次实验和研究发现,电极的特定范围的B1场与电极的温度升高存在一定的关系。具体而言,电极的特定范围的B1场的平方与电极的温度升高存在线性关系。
因此,在一种可能的实现方式中,在步骤S104中得到电极的特定范围的B1场之后,可以对该B1场求平方,然后根据求平方所得到的值和该线性关系来预测出温度升高,由此预测出电极的RF感应发热情况。
这样,根据上述实施例,通过根据接收系数K和序列的B1+rms来计算电极的特定范围的B1场、然后根据该B1场来预测RF感应发热情况,由此可以在考虑到电极的形状、角度和放置位置等的情况下预测RF感应发热。因此,无论电极的形状、角度和放置位置等如何,应用本公开的上述实施例都可以精确地表征RF感应发热,从而确定具备电极的医疗设备在MRI检查下的RF感应发热。
图2是示出根据本公开的另一实施例的温升预测方法的流程图。图2中与图1相同的附图标记表示相同的含义,并且在本实施例中将省略其说明。
本实施例的温升预测方法与上述实施例的温升预测方法主要区别在于,在上述步骤S101之前,还可以包括以下步骤:
步骤S100、获取电极的位置和方向。
具体而言,在获取电极周围的B1图像之前,需要先获取电极的位置和方向。例如,通过MRI检查来对患者进行扫描,由此扫描定位出位于植入式医疗设备内的电极的路径位置,并进一步确定出该电极的路径方向。
这样,在获取到电极的位置和方向之后,可以便于获取到在MRI检查下该电极(或者该电极的尖端部分)周围部分的B1图像。
此外,在一种可能的实现方式中,本实施例的温升预测方法在步骤S105之后,还可以包括以下步骤:
步骤S106、根据所预测出的电极的温升来判断是否发出警报提醒。
具体而言,可以针对电极的RF感应发热情况设置安全阈值,并将步骤S105中所预测出的电极的RF感应发热与该安全阈值进行比较。例如,可以用温度升高来表征RF感应发热情况,相应地可以设置相应的温度升高作为安全阈值。在所预测出的RF感应发热、即温度升高大于或等于安全阈值的情况下,可以发出警报,从而提醒医护人员该患者进行MRI检查存在一定的风险,这样医护人员例如可以终止MRI检查。另一方面,在温度升高小于安全阈值的情况下,表明可以安全地对患者进行MRI检查,这时可以不发出警告,或者也可以发出表明MRI检查安全的提醒。
需要说明的是,可以针对不同的患者设置不同的安全阈值。然而,也可以针对不同的患者设置相同的安全阈值。
这样,通过本实施例的温升预测方法,可以直观地提示患者在MRI检查中的风险,由此尽可能的避免患者在MRI检查中遭受到的危害。
需要说明的是,以上以预测位于植入式医疗设备内的电极在MRI检查、即核磁下的温升的示例来说明了本公开,但是本公开不限于此,并且本公开的温升预测方法同样一般的医疗设备内的电极、或者单独应用于单极。换句话说,通过本公开的温升预测方法,同样可以预测一般的医疗设备内的电极在核磁下的温升,甚至可以预测单独的电极在核磁下的温升。
以下将通过实验的方式来进一步阐述本发明的温升预测方法。
为了确定本公开的温升预测方法的可行性,采用体模的方法进行了验证。首先将电极与测量温度升高的光纤并排捆绑,使得光纤测温仪与电极尖端触点紧密接触。图3示出将电极与测量温度升高的光纤并排捆绑的示意图。
然后,配置琼脂溶液,配置方法为:琼脂为15g/L,CuSO4为5mmol/L,以及NaCl为1g/L。其中,加入CuSO4的目的是缩短T1弛豫时间,从而模拟人体组织。
然后,将并排捆绑的电极与光纤垂直悬挂,以使得其大概与塑料圆柱筒的轴线重合,并调节电极位置的高度,使电极顶端与塑料圆柱筒的底面距离约为5cm。为了防止高温使电极外层的聚氨酯管变形,在琼脂溶液的温度低于70℃之后,再将琼脂溶液缓慢注入塑料圆柱筒,并静置以待其冷却凝结。将凝胶模在室温条件下放置12小时,从而确保其在扫描时温度达到平衡状态。
之后,进行核磁实验,具体阐述如下。
将上述的琼脂模型沿着核磁共振仪的Z轴放置在体模中。然后在体模中注入凝胶,使液面达到9cm。针对核磁共振仪设置以下两种序列。序列一为B1场测量序列,参数如下:
TR(Repetition Time,重复时间)=2000ms,TE(Echo Time,回波时间)=shortest(最短),flip angle(翻转角)=60°,FOV(Field of View,视野区域)=180mm×180mm,slice thickness(层厚)=5mm,slices(层数)=10,以及dynamic scans(重复扫描次数)=1,voxel size(体素大小)=2.5x2.5,B1+rms=0.72uT(9%)。
序列二选取T2TSE序列,该序列的B1+rms比较大,因此能够在相同的条件下对电极进行更有效的加热,从而提高结果的信噪比。该序列参数为:
TE=90ms,flip angle=90deg,TR=shortest,FOV=200×200×60,voxel size=0.7×0.7×2,slices=30,fast imaging mode(快速成像模式)=TSE。
为了探究电极在不同路径形态(形状)、摆放位置、摆放角度的情况下,电极尖端的B1场与温度升高之间的关系,设置了不同的摆放方式,如下表所示。
编号 | 电极位置 | 形态 | 角度 |
1 | 体模左腰部 | 拉直 | 0 |
2 | 体模左腰部 | 团圈 | 0 |
3 | 体模左腰部 | 直角 | 0 |
4 | 体模中间部分 | 拉直 | 0 |
5 | 体模中间部分 | 团圈 | 0 |
6 | 体模中间部分 | 拉直 | 30 |
7 | 体模中间部分 | 拉直 | 60 |
8 | 体模中间部分 | S形 | 0 |
接着,进行数据拟合。
首先,将前7种情况局部B1场求平方,然后与温度升高进行拟合,可以得到如图4所示的直线。图4示出将对局部B1场求平方得到的值与温度升高进行拟合而得到的直线的示意图。然后,根据测量得到的拟合解析式为:
y=0.2938x
其中x为局部B1场均方根,y为实际的温升。根据第8种情况下采集得到的B1场计算预测的温度升高约为6.66℃,这与实际测量得到的温度6.35℃接近度相对误差小于5%。也就是说,根据本实施例的温升预测方法预测出的电极的温度升高与实际测量到的温度升高非常接近。
由此可知,根据本公开的实施例的温升预测方法,可以精确地预测电极的RF感应发热情况。
以下将通过数值计算的方式来研究本公开的温升预测方法。
本研究通过有限元数值模拟。图5示出根据本公开一实施例的针对温升预测方法的有限元数值模拟的示意图。首先,构建一个八通道鸟笼线圈以产生RF场,并将RF场调整为以128MHz工作(图5的a)。构建ASTM体模以模拟人体的头部和躯干,并且位于鸟笼线圈的等的中心内部,并且直径为1.3mm的电极位于其中。4个铂铱环与实际DBS电极相同。绝缘厚度为0.65mm,电极内部的螺旋线简化为直线。体模介质的电导率和介电常数分别设定为0.48s/m和78。绝缘的电导率和介电常数分别设定为0和3。坐标原点位于体模躯干的中心。将未植入体模的全身SAR调整为1W/kg。为了覆盖复杂的现实条件,如图5的b所示,模拟了不同的电极位置和角度。还模拟了如下的不同形状的电极(图5的c):(i)直线放置;(ii)具有直角弯曲的电极;(iii)具有螺旋弯曲的电极;(iv)与金属外壳连接的电极;(v)具有“S”弯曲的电极。所有弯曲距离电极尖端至少100mm,以确保尖端附近的长方体范围相同。将距离电极的尖端接触面0.3mm的点设定为温度探头,以对RF感应发热进行评价。
在数据后处理中,在该研究中,如图5的d所示,电极附近的长方体区域中的局部B1场的平均值被用作RF感应发热的表征指标。长方体区域的长度、宽度和高度分别为50mm、20mm和20mm。为了确定局部B1作为表征指标的可行性和稳健性,模拟了电极的不同形状、角度和放置位置的9种不同电极条件。接着,计算局部B1场平均值的平方,并将所得到的值与电极的温度升高进行拟合。通过拟合所得到的相关系数用于估计局部B1场作为RF发热的表征的实用性。
图6示出根据本公开一实施例的针对温升预测方法的有限元数值模拟的结果的示意图。图6的a示出局部B1场分布。结果表明,B1场的分布受电极的形状、角度和放置位置的显着影响,而电极尖端周围的分布模式相似。接着,计算并绘制局部B1场平均值的平方与温度升高之间的相关性,结果如图6的b所示。如图6的b所示,示出了强线性相关性,并且相关系数R约为0.9992,这意味着即使对于复杂的放置,局部B1场平均值也是RF感应发热的可行表征指标。
由此,通过数值计算的方式,也验证了根据本公开的实施例的温升预测方法,可以精确地预测电极的RF感应发热。
图7示出根据本公开一实施例的温升预测装置的结构框图。温升预测装置70用于预测电极在磁共振成像检查、即MRI检查下的温升。如图7所示,所述温升预测装置70主要包括:B1图像获取单元71,用于获取在所述MRI检查下所述电极周围的特定范围的B1图像;接收系数计算单元72,用于将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以计算所述电极的接收系数;B1+rms确定单元73,用于确定在所述MRI检查下要扫描的序列,并确定所述序列的射频场的大小、即B1+rms;B1场确定单元74,用于根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场;以及预测单元75,用于根据所述特定范围内的B1场来预测所述电极的温升。
在一种可能的实现方式中,所述背景射频场的B1图像是以下B1图像中任一:在无所述电极的情况下的所述特定范围的B1图像;以及与所述电极的距离在特定范围以外的范围的B1图像。
在一种可能的实现方式中,所述B1图像获取单元71可以用于获取所述电极的尖端部分的特定范围内的B1图像。
在一种可能的实现方式中,所述B1+rms确定单元73可以用于:在所述序列为一个的情况下,确定该序列的B1+rms;以及在所述序列为多个的情况下,确定多个序列的B1+rms的最大值或平均值。
在一种可能的实现方式中,所述接收系数计算单元72用于进行以下操作中任一:将所述B1图像中的B1场与所述背景射频场进行对比以得到所述B1场与所述背景射频场的比值,并将所得到的比值中大于第一阈值的比值的均方值作为所述电极的接收系数,或者将所得到的比值中小于第二阈值的比值的均方值的倒数作为所述电极的接收系数,其中所述第二阈值小于所述第一阈值;以及将所述B1图像中的B1场与所述背景射频场进行对比,以将所述B1场的均方值与所述背景射频场的均方值的比值计算为所述电极的接收系数。
在一种可能的实现方式中,所述B1场确定单元74用于将所述接收系数与所述序列的B1+rms相乘,以确定所述特定范围内的B1场。
在一种可能的实现方式中,所述预测单元75用于:对所述特定范围内的B1场求平方;以及根据对所述特定范围内的B1场求平方所得到的值来预测所述电极的温升。
在一种可能的实现方式中,所述电极位于植入式医疗设备内。
根据本公开的实施例的温升预测装置70的具体实现方法详见上述实施例中所阐述的温升预测方法,这里不再赘述。
这样,根据上述实施例,通过根据接收系数K和序列的B1+rms来计算电极的特定范围的B1场、然后根据该B1场来预测RF感应发热情况,由此可以在考虑到电极的形状、角度和放置位置等的情况下预测RF感应发热。因此,无论电极的形状、角度和放置位置等如何,应用本公开的上述实施例都可以精确地表征RF感应发热,从而确定具备电极的医疗设备在MRI检查下的RF感应发热。
图8示出根据本公开的另一实施例的温升预测装置的结构框图。图8中与图7所示的组件相同的组件具有相同的含义,这里将省略其说明。
如图8所示,图8中的温升预测装置80与图7中的温升预测装置70的主要区别在于,除了图7中的温升预测装置70的组件以外,所述温升预测装置80还可以包括:电极位置和方向获取单元76,用于获取所述电极的位置和方向。
在一种可能的实现方式中,所述温升预测装置80还可以包括:警报单元76,用于根据所预测出的所述电极的温升来判断是否发出警报提醒。
根据本公开的实施例的温升预测装置80的具体实现方法详见上述实施例中所阐述的温升预测方法,这里不再赘述。
这样,通过本实施例的温升预测装置,可以直观地提示患者在MRI检查中的风险,由此尽可能的避免患者在MRI检查中遭受到的危害。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种温升预测方法,用于预测电极在磁共振成像检查、即MRI检查下的温升,其特征在于,所述温升预测方法包括:
获取在所述MRI检查下所述电极周围的特定范围的B1图像;
将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以计算所述电极的接收系数;
确定在所述MRI检查下要扫描的序列,并确定所述序列的射频场的大小、即B1+rms;
根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场;以及
根据所述特定范围内的B1场来预测所述电极的温升。
2.根据权利要求1所述的温升预测方法,其特征在于,所述背景射频场的B1图像是以下B1图像中任一:在无所述电极的情况下的所述特定范围的B1图像;以及与所述电极的距离在特定范围以外的范围的B1图像。
3.根据权利要求1所述的温升预测方法,其特征在于,获取在所述MRI检查下所述电极周围的特定范围的B1图像,包括:
获取所述电极的尖端部分的特定范围内的B1图像。
4.根据权利要求1所述的温升预测方法,其特征在于,确定所述序列的B1+rms包括:
在所述序列为一个的情况下,确定该序列的B1+rms;以及
在所述序列为多个的情况下,确定多个序列的B1+rms的最大值或平均值。
5.根据权利要求1所述的温升预测方法,其特征在于,将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以计算所述电极的接收系数包括以下方式中任一:
将所述B1图像中的B1场与所述背景射频场进行对比以得到所述B1场与所述背景射频场的比值,并将所得到的比值中大于第一阈值的比值的均方值作为所述电极的接收系数,或者将所得到的比值中小于第二阈值的比值的均方值的倒数作为所述电极的接收系数,其中所述第二阈值小于所述第一阈值;以及
将所述B1图像中的B1场与所述背景射频场进行对比,以将所述B1场的均方值与所述背景射频场的均方值的比值计算为所述电极的接收系数。
6.根据权利要求1所述的温升预测方法,其特征在于,根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场,包括:
将所述接收系数与所述序列的B1+rms相乘,以确定所述特定范围内的B1场。
7.根据权利要求1所述的温升预测方法,其特征在于,根据所述特定范围内的B1场来预测所述电极的温升,包括:
对所述特定范围内的B1场求平方;以及
根据对所述特定范围内的B1场求平方所得到的值来预测所述电极的温升。
8.根据权利要求1所述的温升预测方法,其特征在于,所述电极位于植入式医疗设备内。
9.根据权利要求1至8中任一项所述的温升预测方法,其特征在于,在根据所述特定范围内的B1场来预测所述电极的温升之后,还包括:
根据所预测出的所述电极的温升来判断是否发出警报提醒。
10.一种温升预测装置,用于预测电极在磁共振成像检查、即MRI检查下的温升,其特征在于,所述温升预测装置包括:
B1图像获取单元,用于获取在所述MRI检查下所述电极周围的特定范围的B1图像;
接收系数计算单元,用于将所述B1图像与在所述MRI检查下的背景射频场的B1图像进行对比,以计算所述电极的接收系数;
B1+rms确定单元,用于确定在所述MRI检查下要扫描的序列,并确定所述序列的射频场的大小、即B1+rms;
B1场确定单元,用于根据所述接收系数和所述序列的B1+rms来确定所述特定范围内的B1场;以及
预测单元,用于根据所述特定范围内的B1场来预测所述电极的温升。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811431311.8A CN109541511B (zh) | 2018-11-26 | 2018-11-26 | 温升预测方法和装置 |
PCT/CN2019/100406 WO2020107944A1 (zh) | 2018-11-26 | 2019-08-13 | 温升预测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811431311.8A CN109541511B (zh) | 2018-11-26 | 2018-11-26 | 温升预测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109541511A true CN109541511A (zh) | 2019-03-29 |
CN109541511B CN109541511B (zh) | 2020-02-14 |
Family
ID=65850530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811431311.8A Active CN109541511B (zh) | 2018-11-26 | 2018-11-26 | 温升预测方法和装置 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109541511B (zh) |
WO (1) | WO2020107944A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020107944A1 (zh) * | 2018-11-26 | 2020-06-04 | 清华大学 | 温升预测方法和装置 |
CN115458143A (zh) * | 2022-09-16 | 2022-12-09 | 兰州大学 | 基于神经网络的无源植入式医疗器械射频发热评估方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11691005B2 (en) * | 2018-12-20 | 2023-07-04 | Medtronic, Inc. | Medical device and MRI systems |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101109720A (zh) * | 2006-07-19 | 2008-01-23 | 西门子(中国)有限公司 | 测量磁性材料磁感应强度对温度变化特性的方法及装置 |
CN101268943A (zh) * | 2007-03-19 | 2008-09-24 | 株式会社东芝 | 磁共振摄像装置及方法、磁共振摄像维护装置及*** |
CN101415459A (zh) * | 2004-08-25 | 2009-04-22 | 麦克罗文逊股份有限公司 | 可植入物装置的热脱离*** |
WO2011014465A2 (en) * | 2009-07-31 | 2011-02-03 | Proteus Biomedical, Inc. | Lead for use in rf field |
CN102866373A (zh) * | 2011-12-12 | 2013-01-09 | 中国科学院深圳先进技术研究院 | 磁共振温度成像中温度测量的矫正方法和*** |
CN103513204A (zh) * | 2012-06-29 | 2014-01-15 | 西门子(深圳)磁共振有限公司 | 一种磁共振成像中k空间数据的轨迹校正方法和装置 |
CN104606780A (zh) * | 2015-01-19 | 2015-05-13 | 清华大学 | 一种mri相容的植入式医疗器械及其连接方法和连接机构 |
EP2944972A1 (en) * | 2014-05-12 | 2015-11-18 | Max-Delbrück-Centrum für Molekulare Medizin (MDC) | Method for the assessment of radiofrequency induced power absorption of conducting implants during MRI |
CN106667487A (zh) * | 2016-12-23 | 2017-05-17 | 清华大学 | 一种监测有源植入物周围组织温度的方法和磁共振成像*** |
CN106896334A (zh) * | 2017-04-18 | 2017-06-27 | 清华大学 | 一种预先评估mr下有源植入物周围组织温度的方法和磁共振成像*** |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7123012B2 (en) * | 2002-11-29 | 2006-10-17 | Advanced Imaging Research, Inc. | Multiple tuned radio frequency coil for resonance imaging and spectroscopic analysis |
CN100425199C (zh) * | 2006-08-25 | 2008-10-15 | 北京源德生物医学工程有限公司 | 永磁磁共振图像导引体外高能聚焦超声***和方法 |
CN104083823B (zh) * | 2014-06-27 | 2017-09-26 | 清华大学 | 一种mri兼容的植入式电极 |
CN205181369U (zh) * | 2015-11-30 | 2016-04-27 | 上海联影医疗科技有限公司 | 一种pet-mr扫描装置 |
CN109541511B (zh) * | 2018-11-26 | 2020-02-14 | 清华大学 | 温升预测方法和装置 |
-
2018
- 2018-11-26 CN CN201811431311.8A patent/CN109541511B/zh active Active
-
2019
- 2019-08-13 WO PCT/CN2019/100406 patent/WO2020107944A1/zh active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101415459A (zh) * | 2004-08-25 | 2009-04-22 | 麦克罗文逊股份有限公司 | 可植入物装置的热脱离*** |
CN101109720A (zh) * | 2006-07-19 | 2008-01-23 | 西门子(中国)有限公司 | 测量磁性材料磁感应强度对温度变化特性的方法及装置 |
CN101268943A (zh) * | 2007-03-19 | 2008-09-24 | 株式会社东芝 | 磁共振摄像装置及方法、磁共振摄像维护装置及*** |
WO2011014465A2 (en) * | 2009-07-31 | 2011-02-03 | Proteus Biomedical, Inc. | Lead for use in rf field |
CN102866373A (zh) * | 2011-12-12 | 2013-01-09 | 中国科学院深圳先进技术研究院 | 磁共振温度成像中温度测量的矫正方法和*** |
CN103513204A (zh) * | 2012-06-29 | 2014-01-15 | 西门子(深圳)磁共振有限公司 | 一种磁共振成像中k空间数据的轨迹校正方法和装置 |
EP2944972A1 (en) * | 2014-05-12 | 2015-11-18 | Max-Delbrück-Centrum für Molekulare Medizin (MDC) | Method for the assessment of radiofrequency induced power absorption of conducting implants during MRI |
CN104606780A (zh) * | 2015-01-19 | 2015-05-13 | 清华大学 | 一种mri相容的植入式医疗器械及其连接方法和连接机构 |
CN106667487A (zh) * | 2016-12-23 | 2017-05-17 | 清华大学 | 一种监测有源植入物周围组织温度的方法和磁共振成像*** |
CN106896334A (zh) * | 2017-04-18 | 2017-06-27 | 清华大学 | 一种预先评估mr下有源植入物周围组织温度的方法和磁共振成像*** |
Non-Patent Citations (4)
Title |
---|
CHRISTINA L. SAMMET 等: "RF-related heating assessment of extracranial neurosurgical implants at 7T", 《MAGNETIC RESONANCE IMAGING》 * |
MO, XIAOLONG 等: "Study of deep brain stimulation lead resonant length in 3.0 T MRI RF magnetic field", 《ELECTRONICS LETTERS》 * |
任新颖: "肿瘤热疗中组织温度场测量方法及其关键技术研究", 《中国博士学位论文全文数据库 医药卫生科技辑》 * |
李波,姜长青,张峰,李路明,马羽: "基于伪影模板的术后电极定位方法", 《清华大学学报(自然科学版)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020107944A1 (zh) * | 2018-11-26 | 2020-06-04 | 清华大学 | 温升预测方法和装置 |
CN115458143A (zh) * | 2022-09-16 | 2022-12-09 | 兰州大学 | 基于神经网络的无源植入式医疗器械射频发热评估方法 |
CN115458143B (zh) * | 2022-09-16 | 2023-05-23 | 兰州大学 | 基于神经网络的无源植入式医疗器械射频发热评估方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2020107944A1 (zh) | 2020-06-04 |
CN109541511B (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nordbeck et al. | Spatial distribution of RF‐induced E‐fields and implant heating in MRI | |
CN109541511A (zh) | 温升预测方法和装置 | |
Neufeld et al. | Measurement, simulation and uncertainty assessment of implant heating during MRI | |
CN106896334B (zh) | 预先评估mr下有源植入物周围组织温度的方法和mri*** | |
Calcagnini et al. | In vitro investigation of pacemaker lead heating induced by magnetic resonance imaging: role of implant geometry | |
Griffin et al. | Safely assessing radiofrequency heating potential of conductive devices using image‐based current measurements | |
Bhusal et al. | Measurements and simulation of RF heating of implanted stereo‐electroencephalography electrodes during MR scans | |
CN106667487B (zh) | 一种监测有源植入物周围组织温度的方法和磁共振成像*** | |
Shrivastava et al. | Effect of the extracranial deep brain stimulation lead on radiofrequency heating at 9.4 Tesla (400.2 MHz) | |
US10912482B2 (en) | Method for determining at least one type and/or condition of cells and system | |
Lucano et al. | A numerical investigation on the effect of RF coil feed variability on global and local electromagnetic field exposure in human body models at 64 MH z | |
Yang et al. | Impact of electrode structure on RF-induced heating for an AIMD implanted lead in a 1.5-Tesla MRI system | |
US11428767B2 (en) | Computer-implemented method for providing an actuation sequence and method for acquiring measurement data by a magnetic resonance device | |
McCabe et al. | Technique to assess the compatibility of medical implants to the RF field in MRI | |
Dharmadhikari et al. | Evaluation of radiofrequency safety by high temperature resolution MR thermometry using a paramagnetic lanthanide complex | |
Karanasiou et al. | Single-frequency and multiband microwave radiometry for feasible brain conductivity variation imaging during reactions to external stimuli | |
Dixit et al. | Thermo‐acoustic ultrasound for noninvasive temperature monitoring at lead tips during MRI | |
Brown | Radiofrequency heating near medical devices in magnetic resonance imaging | |
Chen et al. | Planar MIT sensor array with gradiometers for local hemorrhage detection | |
Talcoth et al. | Electromagnetic modeling of pacemaker lead heating during MRI | |
McCabe | Implantable electrode structures and their RF effects in MRI machines | |
CN108469543A (zh) | 感生电压测量方法及*** | |
Brown et al. | RF-induced Heating Near Active Implanted Medical Devices in MRI: Impact of Tissue Simulating Medium | |
Li et al. | Using scaling approach to estimate MRI RF field induced heating for small medical implant | |
Bassen et al. | Evaluation of unintended electrical stimulation from MR gradient fields |
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 |