CN107249690B - 高强度聚焦超声损伤判定方法及装置 - Google Patents
高强度聚焦超声损伤判定方法及装置 Download PDFInfo
- Publication number
- CN107249690B CN107249690B CN201580001252.9A CN201580001252A CN107249690B CN 107249690 B CN107249690 B CN 107249690B CN 201580001252 A CN201580001252 A CN 201580001252A CN 107249690 B CN107249690 B CN 107249690B
- Authority
- CN
- China
- Prior art keywords
- displacement
- treatment
- magnetic resonance
- tissue
- focused ultrasound
- 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
- 238000002604 ultrasonography Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title abstract description 28
- 238000006073 displacement reaction Methods 0.000 claims abstract description 107
- 230000006378 damage Effects 0.000 claims abstract description 26
- 230000005855 radiation Effects 0.000 claims abstract description 26
- 230000008859 change Effects 0.000 claims abstract description 22
- 230000005284 excitation Effects 0.000 claims abstract 2
- 230000001575 pathological effect Effects 0.000 claims description 31
- 238000003384 imaging method Methods 0.000 claims description 14
- 230000011218 segmentation Effects 0.000 claims description 4
- 230000005404 monopole Effects 0.000 claims description 3
- 238000000264 spin echo pulse sequence Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 11
- 230000003902 lesion Effects 0.000 abstract description 3
- 238000002203 pretreatment Methods 0.000 abstract description 2
- 230000001338 necrotic effect Effects 0.000 abstract 1
- 210000001519 tissue Anatomy 0.000 description 84
- 238000010586 diagram Methods 0.000 description 14
- 238000010438 heat treatment Methods 0.000 description 11
- 238000002560 therapeutic procedure Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 235000021336 beef liver Nutrition 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 6
- 235000015277 pork Nutrition 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 210000005228 liver tissue Anatomy 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000012045 magnetic resonance elastography Methods 0.000 description 4
- 230000003685 thermal hair damage Effects 0.000 description 4
- 208000027418 Wounds and injury Diseases 0.000 description 3
- 238000004925 denaturation Methods 0.000 description 3
- 230000036425 denaturation Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 206010028851 Necrosis Diseases 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 102000004169 proteins and genes Human genes 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000017074 necrotic cell death Effects 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 239000011120 plywood Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 230000000451 tissue damage Effects 0.000 description 1
- 231100000827 tissue damage Toxicity 0.000 description 1
Classifications
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgical Instruments (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
一种高强度聚焦超声损伤判定方法及装置，包括：S101、获得采用高强度聚焦超声对病变组织治疗前和治疗后的磁共振声辐射力位移检测序列；S102，S103、分别在治疗前和治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗前和治疗后的病变组织内位移；S104、当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。该方法不受温度的影响，可实现快速准确定性组织是否已发生坏死；无需复杂的外部激励源，不易受噪声影响；操作方法简单，不会对组织产生附加的热量及损伤，检测方法安全。
Description
技术领域
本发明涉及生物医学工程技术领域，尤其一种高强度聚焦超声损伤判定方法及装置。
背景技术
高强度聚焦超声(High intensity focused ultrasound，HIFU)可利用超声波的穿透性和能量沉积性将体外低能量超声波汇聚于特定深度的靶组织，使焦域内温度迅速升高至60℃，引发蛋白质变性和组织凝固性坏死，而不损伤周围正常组织，HIFU治疗已逐渐发展为一种重要的肿瘤治疗手段。高强度聚焦超声治疗过程中的图像引导、监控及治疗效果评判是实现聚焦超声安全治疗的基本保障，也是高效、精细和个性化治疗的基础。目前常用的高强度聚焦超声治疗效果评判方法包括：
(1)基于组织T1、T2、质子密度进行损伤检测：T1-W、T2-W等都是磁共振监控高强度聚焦超声治疗的常用序列，当蛋白质受高温作用变性时，组织T1、T2和质子密度值会发生变化。然而这些参数除了受组织性质改变本身影响外，也受温度的影响，比如组织T1值随温度的升高而增长。因此基于组织T1、T2及质子密度进行损伤判断时，为了准确定性组织是否已发生坏死，还是温度影响引起的图像信号改变，需要等高强度聚焦超声引起的温度扩散到热平衡。
(2)基于超声弹性成像进行损伤检测：弹性是人体组织的重要物理特性，正常组织、异常组织、病变组织之间弹性存在显著差异。如从某一方向对组织施加一定应力，组织内会产生相应应变。相同应力作用下，组织应变越大，位移越大，组织硬度越小。高强度聚焦超声治疗后，焦域坏死表现为细胞核破裂或固缩、蛋白质变性、细胞结构破坏等，直接导致组织弹性的改变。
超声弹性成像的原理是对组织施加一个内部或外部的动态或者静态/准静态的压力，通过采集组织施压前后的射频回波信号，基于信号间的相关性，对其进行信号处理提取组织位移信息，可得到组织应变信息，从而间接反映组织弹性模量。
(3)基于磁共振弹性成像进行损伤检测：磁共振成像(Magnetic ResonanceImaging，MRI)由于其多方位立体成像、特有的温度监控、多种组织对比等特点，已成为聚焦超声治疗监控的最佳手段之一。
磁共振弹性成像的步骤主要是：1、需要一套磁共振兼容的激励源，在组织内产生频率50Hz至500Hz的剪切波；2、通过添加运动敏感梯度，将组织位移在磁共振相位图上进行编码，相位图与组织位移关系如下：
其中，φ为组织运动及施加梯度场在位置处产生的相位，θ为组织运动与运动敏感梯度之间的相位差，γ为磁旋比γ/2π＝42.57MHz/T，N为使用的运动敏感梯度的对数，T、分别为运动敏感梯度的时长与幅值，ξ0为最大位移，为波数；3、基于上述公式，对数据进行处理得到组织弹性的定量分布结果。
上述三种方法具有如下缺点：
基于组织T1、T2、质子密度进行损伤检测：受温度影响，需要等到组织内温度回复到热平衡温度，才可以快速准确定性组织是否已发生坏死，还是温度影响引起的图像信号改变反馈；
基于超声弹性成像进行损伤检测：超声图像本身分辨率明显低于磁共振图像，会直接影响判断的准确性。此外，高频超声的图像分辨率相对较高，但高频超声在组织内的传播距离有限，导致超声弹性成像更适用于浅表器官诊断；
基于磁共振弹性成像进行损伤检测：需要复杂的外部激励源，且激励源产生的剪切波存在与超声弹性成像相同的问题：空间分辨率会随着频率的增加而增加，但高频剪切波衰减速度快，传播距离有限，此外磁共振弹性成像的重建需要求解一个Helmholtz(赫姆霍兹)问题，极容易受到噪声的影响。
发明内容
本发明实施例提供一种高强度聚焦超声损伤判定方法，该方法不受温度的影响，可实现快速准确定性组织是否已发生坏死；无需复杂的外部激励源，不易受噪声影响；操作方法简单，不会对组织产生附加的热量及损伤，检测方法安全。该方法包括：
获得采用高强度聚焦超声对病变组织治疗前和治疗后的磁共振声辐射力位移检测序列；
在治疗前的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗前的病变组织内位移；
在治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗后的病变组织内位移；
当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。
本发明实施例还提供一种高强度聚焦超声损伤判定装置，采用本发明装置进行高强度聚焦超声损伤的判定，不受温度的影响，可实现快速准确定性组织是否已发生坏死；无需复杂的外部激励源，不易受噪声影响；操作简单，不会对组织产生附加的热量及损伤，检测方法安全。该装置包括：
序列获得模块，用于获得采用高强度聚焦超声对病变组织治疗前和治疗后的磁共振声辐射力位移检测序列；
第一位移确定模块，用于在治疗前的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗前的病变组织内位移；
第二位移确定模块，用于在治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗后的病变组织内位移；
损伤判定模块，用于当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。
在本发明实施例中，通过在治疗前和治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后的磁共振声辐射力位移检测序列的相位变化，确定治疗前和治疗后的病变组织内位移；当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。该方法不受温度的影响，可实现快速准确定性组织是否已发生坏死；无需复杂的外部激励源，不易受噪声影响；操作方法简单，不会对组织产生附加的热量及损伤，检测方法安全。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中高强度聚焦超声损伤判定方法流程图；
图2为本发明实施例中治疗前的组织内位移检测示意图；
图3为本发明实施例中治疗后(即刻)的组织位移检测示意图；
图4为本发明实施例中治疗中的温度监控曲线；
图5为本发明实施例中采用的声辐射力位移检测序列时序图；
图6为本发明实施例中无损伤产生时高强度聚焦超声加热前后组织位移曲线；
图7为本发明实施例中无损伤产生时温度变化曲线；
图8是本发明实施例中猪肉组织内产生损伤时高强度聚焦超声加热前后组织内位移曲线；
图9是本发明实施例中焦点处最大温度变化曲线；
图10是本发明实施例中牛肝组织内产生损伤时高强度聚焦超声加热前后组织内位移曲线；
图11是本发明实施例中焦点处最大温度变化曲线；
图12是本发明实施例中高强度聚焦超声损伤判定装置结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
现有的高强度聚焦超声治疗效果评判方法会受温度的影响，需要等到组织内温度恢复到热平衡温度才可以进行疗效的判断；需要外加激励源才可对组织弹性改变进行检测；极容易受到噪声的影响等等。本发明提出了一种高强度聚焦超声损伤判定方法及装置，可以克服现有技术中存在的缺点。
图1是本发明实施例中高强度聚焦超声损伤判定方法流程图，如图1所示，该方法包括如下步骤：
步骤101：获得采用高强度聚焦超声对病变组织治疗前和治疗后的磁共振声辐射力位移检测序列；
步骤102：在治疗前的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗前的病变组织内位移；
步骤103：在治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗后的病变组织内位移；
步骤104：当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。
具体实施时，在磁共振声辐射力位移检测序列中添加运动编码梯度，将组织在超声作用下引起的位移转换为磁共振图像相位变化，基于以下公式实现组织内位移的定量：
其中，Δx为病变组织内位移，Δφ为添加运动编码梯度前和后两次采集的磁共振图像的相位差，Ge为运动编码梯度的强度，τ为运动编码梯度的时长，γ为磁旋比。
图2和图3分别为治疗前和后的组织内位移检测示意图，可以看出位移有明显的变化。图4是治疗中的温度曲线，温度是持续升高的。
具体实施是，磁共振声辐射力位移检测序列可以采用分段读出的平面回波序列、一维线性扫描序列、二位自旋回波序列、单次激发平面回波序列等等。
位移编码梯度可以采用重复双极位移编码梯度、单极位移编码梯度、双极反向位移编码梯度等等。
具体实施是，在高强度聚焦超声治疗前后采集组织中的位移分布，根据需要可采集多组位移分布结果，保证检测的准确性。高强度聚焦超声治疗过程中可辅以梯度回波序列进行温度监控。
不同组织位移改变量可能存在一定的差异，所以根据前期实验选取合适的阈值，当治疗前后组织位移存在显著差异，且位移差大于该阈值时，可判定该点组织已坏死，从而实现对损伤区域的描绘。
为了验证本发明结果，所有实验在西门子3T磁共振成像系统(Siemens TIM Trio，Erlangen，German)上进行。本实验所使用的声辐射力位移检测序列为分段读出平面回波序列，序列中添加双极性运动编码梯度对组织位移进行编码，序列时序如图5所示，序列中所使用运动编码梯度强度为32mT/m，运动编码梯度时长10ms，其他成像参数为：TR＝200ms，TE＝36ms，带宽＝1149Hz/Pixel，层厚＝5mm，分辨率2.2*1.6mm2，EPI factor＝9；实验中所使用GRE序列TR/TE＝29/10ms，带宽＝390Hz/Pixel。加热前后分别采集20组组织内位移变化。
采用两种组织进行实验：猪肉组织和牛肝组织。
如图5所示，实验时共采集两组图像，第二组所用位移编码梯度(实线)与第一组极性相反(虚线)，以提高检测的灵敏度。
由图6和图7可以看出，尽管高强度聚焦超声治疗过程中组织内存在明显温度升高，但是加热前后的组织位移并无明显的变化，组织内并无热损伤产生。
由图8可以看出，当猪肉组织内产生损伤时，加热前后组织内位移存在显著差异，20组测量结果的平均最大位移之差为2.33μm，加热后组织内位移增大。由图9可以看出，焦点处最大温升为37.6℃，在对猪肉组织切开后，猪肉组织内可见明显热损伤。
由图10可以看出，当牛肝组织内产生损伤时，加热前后组织内位移存在显著差异，20组测量结果的平均最大位移之差为2.19um，加热后组织内位移减小。由图11可以看出，焦点处最大温升为40.5℃，在对牛肝组织切开后，牛肝组织内可见明显热损伤。
尽管在猪肉内加热后组织位移增大，而在牛肝内加热后组织位移减小(加热时间相同，功率不同)，但从实验结果可见当组织中产生热损伤时，加热前后的声辐射力位移检测结果存在显著差异。
基于同一发明构思，本发明实施例中还提供了一种高强度聚焦超声损伤判定装置，如下面的实施例所述。由于高强度聚焦超声损伤判定装置解决问题的原理与高强度聚焦超声损伤判定方法相似，因此高强度聚焦超声损伤判定装置的实施可以参见高强度聚焦超声损伤判定方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图12是本发明实施例的高强度聚焦超声损伤判定装置的一种结构框图，如图12所示，包括：
序列获得模块1201，用于获得采用高强度聚焦超声对病变组织治疗前和治疗后的磁共振声辐射力位移检测序列；
第一位移确定模块1202，用于在治疗前的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗前的病变组织内位移；
第二位移确定模块1203，用于在治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗后的病变组织内位移；
损伤判定模块1204，用于当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。
具体实施时，所述第一位移确定模块1202和第二位移确定模块1203具体用于：
按如下公式确定病变组织内位移：
其中，Δx为病变组织内位移，Δφ为添加运动编码梯度前和后两次采集的磁共振图像的相位差，Ge为运动编码梯度的强度，τ为运动编码梯度的时长，γ为磁旋比。
具体实施时，所述磁共振声辐射力位移检测序列为分段读出的平面回波序列、一维线性扫描序列、二位自旋回波序列或单次激发平面回波序列。
具体实施时，所述位移编码梯度为重复双极位移编码梯度、单极位移编码梯度或双极反向位移编码梯度。
综上，本发明具有以下优点：
1、该方法不需外加激励源，即可对组织弹性改变进行检测；
2、该方法操作简单，不需等组织内温度恢复到热平衡温度即可进行疗效判断；
3、该方法可兼容多种声辐射力位移检测序列；
4、该方法只需毫秒级高强度聚焦超声脉冲作用，不易引起组织内温升，检测方法安全。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种高强度聚焦超声损伤判定装置，其特征在于，包括：
序列获得模块，用于获得采用高强度聚焦超声对病变组织治疗前和治疗后的磁共振声辐射力位移检测序列；
第一位移确定模块，用于在治疗前的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗前的病变组织内位移；
第二位移确定模块，用于在治疗后的磁共振声辐射力位移检测序列中添加运动编码梯度，根据添加运动编码梯度前后获得的磁共振相位图的相位变化，确定治疗后的病变组织内位移；
损伤判定模块，用于当治疗前和治疗后的病变组织内位移的位移差大于预设阈值时，判定病变组织损伤。
2.如权利要求1所述的高强度聚焦超声损伤判定装置，其特征在于，所述第一位移确定模块和第二位移确定模块具体用于：
按如下公式确定病变组织内位移：
其中，Δx为病变组织内位移，Δφ为添加运动编码梯度前和后两次采集的磁共振图像的相位差，Ge为运动编码梯度的强度，τ为运动编码梯度的时长，γ为磁旋比。
3.如权利要求1所述的高强度聚焦超声损伤判定装置，其特征在于，所述磁共振声辐射力位移检测序列为分段读出的平面回波序列、一维线性扫描序列、二位自旋回波序列或单次激发平面回波序列。
4.如权利要求1所述的高强度聚焦超声损伤判定装置，其特征在于，所述位移编码梯度为重复双极位移编码梯度、单极位移编码梯度或双极反向位移编码梯度。
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2015/099843 WO2017113177A1 (zh) | 2015-12-30 | 2015-12-30 | 高强度聚焦超声损伤判定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107249690A CN107249690A (zh) | 2017-10-13 |
| CN107249690B true CN107249690B (zh) | 2019-05-31 |
Family
ID=59224245
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201580001252.9A Active CN107249690B (zh) | 2015-12-30 | 2015-12-30 | 高强度聚焦超声损伤判定方法及装置 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN107249690B (zh) |
| WO (1) | WO2017113177A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107362466B (zh) * | 2017-06-29 | 2019-12-10 | 哈尔滨医科大学 | 一种用于声动力治疗的超声聚焦方法 |
| CN113917562B (zh) * | 2021-09-27 | 2023-02-28 | 中国科学院武汉岩土力学研究所 | 深埋软弱夹层宏细观结构表征及三维空间构造方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102247163A (zh) * | 2010-04-12 | 2011-11-23 | 西门子公司 | 磁共振引导高强度聚焦超声聚焦的方法和装置 |
| CN104602761A (zh) * | 2012-07-09 | 2015-05-06 | 皇家飞利浦有限公司 | 声辐射力磁共振成像 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5679988B2 (ja) * | 2008-12-22 | 2015-03-04 | コーニンクレッカ フィリップス エヌ ヴェ | 印加された力に対する組織変位のリアルタイムの監視のためのアブレーション制御装置 |
| WO2011161559A1 (en) * | 2010-06-24 | 2011-12-29 | Koninklijke Philips Electronics N.V. | Real-time monitoring and control of hifu therapy in multiple dimensions |
| US20130296743A1 (en) * | 2012-05-02 | 2013-11-07 | Siemens Medical Solutions Usa, Inc. | Ultrasound for Therapy Control or Monitoring |
| CN104257357B (zh) * | 2014-10-14 | 2016-05-25 | 深圳先进技术研究院 | 磁共振温度成像方法和系统 |
-
2015
- 2015-12-30 CN CN201580001252.9A patent/CN107249690B/zh active Active
- 2015-12-30 WO PCT/CN2015/099843 patent/WO2017113177A1/zh active Application Filing
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102247163A (zh) * | 2010-04-12 | 2011-11-23 | 西门子公司 | 磁共振引导高强度聚焦超声聚焦的方法和装置 |
| CN104602761A (zh) * | 2012-07-09 | 2015-05-06 | 皇家飞利浦有限公司 | 声辐射力磁共振成像 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2017113177A1 (zh) | 2017-07-06 |
| CN107249690A (zh) | 2017-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kaye et al. | Rapid MR‐ARFI method for focal spot localization during focused ultrasound therapy | |
| US8725232B2 (en) | Therapeutic apparatus | |
| US8229544B2 (en) | Detecting temperature and protein denaturation during thermal therapy | |
| US9968333B2 (en) | Method and apparatus for analyzing elastography of tissue using ultrasound waves | |
| US10271890B2 (en) | High intensity focused ultrasound enhanced by cavitation | |
| Solovchuk et al. | Temperature elevation by HIFU in ex vivo porcine muscle: MRI measurement and simulation study | |
| EP2559375A1 (en) | Method and device for monitoring temperature of treatment site by using ultrasound, and system for treatment and diagnosis using ultrasound | |
| JP2011507585A5 (zh) | ||
| de Bever et al. | Simultaneous MR thermometry and acoustic radiation force imaging using interleaved acquisition | |
| Mougenot et al. | High intensity focused ultrasound with large aperture transducers: a MRI based focal point correction for tissue heterogeneity | |
| Zhang et al. | Enhanced lesion‐to‐bubble ratio on ultrasonic Nakagami imaging for monitoring of high‐intensity focused ultrasound | |
| Guertler et al. | Estimation of anisotropic material properties of soft tissue by MRI of ultrasound-induced shear waves | |
| Todd et al. | Model predictive filtering for improved temporal resolution in MRI temperature imaging | |
| Kujawska et al. | Determination of tissue thermal conductivity by measuring and modeling temperature rise induced in tissue by pulsed focused ultrasound | |
| CN110742645B (zh) | 多模成像系统、多模成像方法及存储介质 | |
| Allen et al. | MR‐based detection of individual histotripsy bubble clouds formed in tissues and phantoms | |
| Liu et al. | Concurrent visualization of acoustic radiation force displacement and shear wave propagation with 7T MRI | |
| CN107249690B (zh) | 高强度聚焦超声损伤判定方法及装置 | |
| Boulos et al. | Weighting the passive acoustic mapping technique with the phase coherence factor for passive ultrasound imaging of ultrasound-induced cavitation | |
| Kothapalli et al. | Acoustic field characterization of a clinical magnetic resonance‐guided high‐intensity focused ultrasound system inside the magnet bore | |
| JP2020527723A (ja) | 媒質の粘弾性の定量化方法及び装置 | |
| WO2019119536A1 (zh) | 一种同步检测组织位移和t2的方法及装置 | |
| Li et al. | Exploring potential mechanisms responsible for observed changes of ultrasonic backscattered energy with temperature variations | |
| US11406361B2 (en) | Mapping shear wave velocity and shear modulus in biological tissues | |
| WO2017123536A1 (en) | Elastography imaging with magnetic resonance imaging guided focused ultrasound |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20250118 Address after: 518000 101, 201, 301, Building 38, Xialang Industrial Zone, Heshuikou Community, Matian Street, Guangming District, Shenzhen, Guangdong Patentee after: UNIMICRO MEDICAL SYSTEMS (SHENZHEN) Co.,Ltd. Country or region after: China Address before: 1068 No. 518055 Guangdong city in Shenzhen Province, Nanshan District City Xili University School Avenue Patentee before: SHENZHEN INSTITUTES OF ADVANCED TECHNOLOGY Country or region before: China |