一种弹性减影成像方法技术

本发明专利技术涉及一种弹性减影成像方法，属于超声弹性成像技术领域，该方法包括：对损伤引入前的组织施加一定的压缩，并采集压缩前后的超声射频数据；对损伤引入后的组织施加一定的压缩，并采集压缩前后的超声射频数据；根据所采集的超声射频数据计算损伤引入前的组织的位移场；所采集的超声射频数据计算损伤引入后的组织的位移场；根据上述两步骤得到的结果，计算弹性减影图像。本发明专利技术将弹性成像和减影成像这两种技术结合，能够实现将弹性图像中的差异突显，同时抑制相同的背景结构。同时，本发明专利技术给出了基于应变场的弹性减影成像实施例和基于规整化位移场的弹性减影成像实施例。本发明专利技术提出的方法还可以用在其他弹性成像方法上。

【技术实现步骤摘要】

本专利技术属于医学成像的超声弾性成像
，特别涉及。
技术介绍
弾性成像作为ー种成像模式，能够提供成像组织的硬度信息。在过去十几年中，弹性成像在一些微创消融治疗监控应用中得到了研究。在这些消融治疗中，通过ー些物理或者化学方法，产生局部组织坏死，这个过程通常伴随着组织硬度的变化。为达到治疗效果， 一方面要求完全摘除病灶组织，另一方面要尽可能少的损害周围健康组织，因此准确评估已经凝固坏死的区域尤其重要。弾性成像发展至今，其在HIFU治疗、射频消融、无水酒精注射治疗中的监控成像的可行性已经被验证。其它成像方式，比如磁共振成像(MRI)、增强型 CT也能够用于这些治疗的监控，但是它们比较昂贵，移动性不够，同时并非适用于所有病人 (MRI不适用帯心脏起搏器的病人，增强型CT不适用于对造影剤有严重过敏反应的病人)。 超声成像也被尝试用于上述消融治疗的监控中，但是由于组织的声特性同组织活性并不相关，因此效果并不理想。相比上述方法，弾性成像能直接提供硬度信息，而且具有实时、低成本、安全的特点，因此有望成为ー种消融治疗监控的理想方法。传统实现弹性成像方法步骤如下1)采集组织压缩前的超声射频数据；2)从外界对组织施加一定的压缩量，并获取組织压缩后的超声射频数据；3)根据步骤幻，对压缩前后得到的数据利用本专业人员所知的方法计算得到应变场，从而获得反应组织硬度分布的弹性图像。但是，目前的弾性成像有一定的局限性。如一些弹性成像中的伪影(机械伪影、声学伪影、以及信号处理方法导致的伪影)可能影响弹性成像的最终結果。而且，对于ー些病灶，如硬癌，本身硬度就大，治疗中损伤虽然会引起硬度的变化，但是由于背景硬度很大，对比度就会变小，从而不易检测治疗引起的硬度变化。在医学成像领域，减影技术已经被广为应用。通过从活板图像(live image)中减去蒙板图像(mask image)，使得两幅图像中的差异得到突显，同时相同的背景结构能够得到抑制。减影技术的应用可以极大的抑制伪影，提高成像结果的信号噪声比和对比度噪声比。现在，减影技术中最广为人知的应用是数字血管减影成像(digital subtraction angiography, DSA)，在DSA流程中，在感兴趣血管弾丸注射造影剤，并且获取注射前后的数字X光投影图像，之后将两图像相减获得没有背景信息的血管减影成像。虽然，弾性成像广为人们研究，减影技术在医学成像中广为应用，但是并没有将减影技术用于弾性成像的系统研究以及将弾性成像和减影技术这两种技术进行结合的相关报道
技术实现思路
本专利技术的目的是克服已有弹性成像技术的不足之处，提出，该方法结合弹性成像技术和减影成像技木，使成像质量得到提高，能够比传统的弹性成像在损伤监控中得到更好的效果。本专利技术提出的结合弹性成像技术和减影成像技术的弹性减影成像方法，包括以下步骤1)对损伤引入前的组织施加一定的压縮，并采集压缩前后的超声射频数据；2)对损伤引入后的组织施加一定的压縮，并采集压缩前后的超声射频数据；3)根据步骤1)所采集的超声射频数据计算损伤引入前的组织的位移场(计算方法可以采用互相关法，绝对值差和法，多普勒法，自相关法等本专业人员所知的位移估计方法)；4)根据步骤2~)所采集的超声射频数据计算损伤引入后的组织的位移场(计算方法可以采用互相关法，绝对值差和法，多普勒法，自相关法等本专业人员所知的位移估计方法)；5)根据步骤3)、4)得到的結果，计算弹性减影图像。本专利技术的原理弹性图像能够提供成像组织的硬度信息，减影技木通过从活板图像(live image) 中减去蒙板图像(mask image)，能够将活板图像和蒙板图像中的差异突显，同时抑制相同的背景结构。本专利技术正是基于弹性图像和减影成像技术的上述特点，通过将损伤引入后的弾性图像减去损伤引入前的弹性图像从而得到弾性减影图像，从而获得比普通弾性成像对于治疗引入的损伤的更好的检测。本专利技术的特点1)将弹性成像技术与减影技术相结合。2)能够将图像中的差异突显，同时抑制相同的背景结构，如背景硬物引入的应カ 集中导致的伪影可以得到有效抑制。3)能够获得较高的信号噪声比和对比度噪声比，成像质量得到提高，对组织中的损伤的检测能力得到提升。4)不容易受到硬组织背景的影响，而且能够提供更可靠的损伤边界信息。5)可以采用图像融合技术，将弹性减影成像结果和其他成像图进行融合显示，不仅可以突显两幅弾性图间的硬度差异，而且可以提供图像背景信息。附图说明图1为本专利技术的实施例1基于应变场的弹性减影成像流程；图2为本专利技术的实施例2基于规整化位移场的弹性减影成像流程；图3为实施例1与传统方法结果对比图；图4为实施例2结果效果图。具体实施例方式本专利技术提出的结合具体实施例和附图详细说明如下实施例1是基于应变场的弹性减影成像的实施例，该实施例包含下述几个步骤1)对损伤引入前的组织施加一定的压縮，并采集压缩前后的超声射频数据，具体步骤如下先根据如图1(a)所示，对损伤引入前的组织沿纵向施加总长度的压縮，井根据图1(b)所示采集压缩前后的超声射频数据。假设所采集超声射频数据包含N条扫描线， 每条扫描线包含L点数据，相邻两点数据对应的实际距离为Δζ。2)对损伤引入后的组织施加一定的压縮，并采集压缩前后的超声射频数据，具体步骤如下如图1(d)所示，对损伤引入后的组织纵向施加总长度的压縮，井根据图1(d) 所示采集压缩前后的超声射频数据。假设所采集超声射频数据包含N条扫描线，每条扫描线包含L点数据，相邻两点数据对应的实际距离为Δζ。3)根据步骤1)所采集的超声射频数据计算损伤引入前的组织的位移场，具体步骤如下3-1)将步骤1)中的每条压缩前后的扫描线划分为长度为2Μ+1、中心间隔为1个数据点的L-2M段数据。这L-2M段数据的中心分别为扫描线上的Μ+1、Μ+2、…、L-M点。则压缩前第k条扫描线上以ρ点为中心的数据段和压缩后第k条扫描线上以q点为中心的数据段的相关系数C；,,为MLrk(P + i)Rk(q + i)_5]= I MMX(rk(p + i))2X(Rk(q + i))2V i=-Mi=-M其中rk、Rk分別代表压缩前第k条扫描线数据和压缩后第k条扫描线数据；3-2)通过计算压缩前第k条扫描线上以ρ为中心的数据段和压缩后第k条扫描线上所有数据段的相关系数，并且搜索得到相关系数最大值cp，qmax如下Cpq = max{Cpq, q = M +1, M + 之…,L - Μ}qmax为Cp,%ax取最大值时候的q值，由此就可以得到中心处于ρ · Δ ζ深度，长度为QM+1) · Δζ的组织片段在压缩前后的位移如下d^z, k) |z = p#^,k = (qmax_p) · Δζ其中，ζ表示组织片段中心所处的深度；3-3)对第k条扫描线上的每一段数据重复上述过程，就可以获得第k条扫描线上的所有数据段对应的组织片段压缩前后的位移，这些位移組成ー个序列Cl1 (z，k)し=(M+1) · Δζ、(Μ+2) · Δζ、…、(L-M) · Δζ, k ‘3-4)对所有扫描线重复上述过程，就可以获得损伤引入前，压缩前后组织的位移场 Id1 {Z, k/ } Iz= (Μ+1) φ Δζ, (Μ+2) · Δζ,…、(L-M) · Δζ,本文档来自技高网...

【技术保护点】

【技术特征摘要】

【专利技术属性】
技术研发人员：邵金华，刘丹，孙锦，段后利，白净，
申请(专利权)人：北京索瑞特医学技术有限公司，清华大学，
类型：发明
国别省市：