CN110969618B - 一种基于动态超声造影的肝肿瘤血管生成定量分析方法 - Google Patents
一种基于动态超声造影的肝肿瘤血管生成定量分析方法 Download PDFInfo
- Publication number
- CN110969618B CN110969618B CN201911308530.1A CN201911308530A CN110969618B CN 110969618 B CN110969618 B CN 110969618B CN 201911308530 A CN201911308530 A CN 201911308530A CN 110969618 B CN110969618 B CN 110969618B
- Authority
- CN
- China
- Prior art keywords
- perfusion
- tumor
- matching
- frame
- focus
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 18
- 230000005747 tumor angiogenesis Effects 0.000 title claims abstract description 16
- 238000002601 radiography Methods 0.000 title claims abstract description 10
- 206010019695 Hepatic neoplasm Diseases 0.000 title claims abstract description 7
- 208000014018 liver neoplasm Diseases 0.000 title claims abstract description 7
- 230000010412 perfusion Effects 0.000 claims abstract description 133
- 206010028980 Neoplasm Diseases 0.000 claims abstract description 46
- 238000010586 diagram Methods 0.000 claims abstract description 23
- 230000033001 locomotion Effects 0.000 claims abstract description 14
- 230000033115 angiogenesis Effects 0.000 claims abstract description 11
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000004458 analytical method Methods 0.000 claims description 16
- 238000002604 ultrasonography Methods 0.000 claims description 15
- 230000000241 respiratory effect Effects 0.000 claims description 14
- 230000004044 response Effects 0.000 claims description 14
- 230000002093 peripheral effect Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 10
- 230000007246 mechanism Effects 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- 230000010361 irregular oscillation Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 238000002607 contrast-enhanced ultrasound Methods 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 6
- 239000002872 contrast media Substances 0.000 claims description 5
- 230000007797 corrosion Effects 0.000 claims description 4
- 238000005260 corrosion Methods 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 3
- 238000012285 ultrasound imaging Methods 0.000 claims description 2
- 230000002440 hepatic effect Effects 0.000 claims 2
- 208000011580 syndromic disease Diseases 0.000 claims 1
- 238000011002 quantification Methods 0.000 abstract description 6
- 238000003745 diagnosis Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 4
- 239000000284 extract Substances 0.000 abstract description 2
- 238000004393 prognosis Methods 0.000 abstract description 2
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 238000005065 mining Methods 0.000 abstract 1
- 238000011158 quantitative evaluation Methods 0.000 abstract 1
- 230000003902 lesion Effects 0.000 description 7
- 230000017531 blood circulation Effects 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 230000036770 blood supply Effects 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 230000008081 blood perfusion Effects 0.000 description 2
- 238000003759 clinical diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012847 principal component analysis method Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 208000006678 Abdominal Neoplasms Diseases 0.000 description 1
- 206010006272 Breast mass Diseases 0.000 description 1
- 206010073069 Hepatic cancer Diseases 0.000 description 1
- 206010027476 Metastases Diseases 0.000 description 1
- 208000009453 Thyroid Nodule Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000009400 cancer invasion Effects 0.000 description 1
- 230000003778 catagen phase Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000001802 infusion Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000002601 intratumoral effect Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000036210 malignancy Effects 0.000 description 1
- 230000009401 metastasis Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000001991 pathophysiological effect Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000004614 tumor growth Effects 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
- G06T7/248—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments involving reference images or patches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30056—Liver; Hepatic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Multimedia (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种基于动态超声造影的肝肿瘤血管生成定量分析方法，通过局部灌注模式差异量化，构建反映肿瘤内及周边正常组织血管分布结构的网络图，从而表征肿瘤区域功能和形态的空间异质性，为临床影像医生提供新的肿瘤血管生成定量评估方法。本发明在抑制运动干扰的基础上提取局部灌注模式，并基于灌注相位约束量化局部灌注模式差异，然后通过复杂网络表征肿瘤局部血管生成差异性，从而挖掘与临床肿瘤诊断、预后评估相关联的异质性特征。
Description
技术领域
本发明涉及超声造影分析技术领域，尤其是一种基于动态超声造影的肝肿瘤血管生成定量分析方法。
背景技术
超声造影(Contrast-enhanced Ultrasound,CEUS)作为一种分子影像技术，能够准确地显示组织中微小动静脉和血流灌注情况，相比于传统灰度超声和彩色多普勒超声成像技术，能够有效地提高诊断的敏感性和特异性。同时，超声造影具有其安全、廉价、无创等特点，近10年来在我国临床肿瘤诊断中取得了很大的发展。肿瘤血管的生成是肿瘤生长和侵袭的基础，因此，肿瘤内及周边血管生成定量分析是评估肿瘤恶性程度、转移复发的重要方向。超声造影能够以较高的时间分辨率、空间分辨率实时动态地显示肿瘤滋养血管及血流灌注情况，更客观地辅助医生判断病灶大小及周围浸润状况。超声造影的研究已有20多年的历史，主要应用于腹部肿瘤(肝脏占位性病变等)，乳腺及甲状腺结节诊断。
目前，临床影像医生应用低机械指数连续扫描技术，观察病灶及正常实质内的血供。完整的灌注相位通常可分为灌入相位和消退相位，在肝脏占位性病变中，又具体分为动脉期、门脉期、实质期和延迟期。临床上，超声造影分析分为定性和定量两种方法。超声造影的定性分析是指观察在病灶和周边正常实质像素强度相对变化，例如，“快进快退”是恶性肝肿瘤的典型特征；增强程度，增强均匀性等。定量分析是指，通过时间-强度曲线(Time-intensity Curve,TIC)描绘感兴趣区内(region of interest,ROI)血供变化，能同时获得峰值强度、达峰时间、曲线下面积等客观参数。通过测定各种灌注参数的变化，可对活体肿瘤组织的血流灌注进行判断，为肿瘤治疗术后的疗效评估提供更加丰富的病理学及病理生理学信息。然而，超声造影评价肿瘤血流灌注还有诸多问题亟待解决，包括检查的可重复性差，呼吸运动干扰，区域灌注信息丢失，肿瘤血管生成均匀性、浸润性评估的主观性强等。
发明内容
为了解决现有技术中的问题，本发明提供一种基于灌注相位约束匹配的超声造影定量分析方法，实现基于病灶跟踪的呼吸干扰抑制、稳定灌注模式提取、像素级TIC差异定量表征，为肿瘤形状不规则度、肿瘤内部血管生成不均匀性及对正常组织侵袭程度评估提供新的定量分析指标。
为了实现上述目的，本发明采用以下技术方案，
一种基于动态超声造影的肝肿瘤血管生成定量分析方法，包括如下步骤：
(1)通过基于空间特征重标定的孪生卷积网络跟踪病灶区域，从而消除呼吸运动对超声造影定量分析带来的干扰；
(2)利用步骤(1)获取的无呼吸干扰的超声造影子序列，通过稳定成分分析算法提取局部灌注模式；
(3)通过基于灌注相位约束匹配的相似性度量量化肿瘤内部及周边正常组织的局部灌注模式差异；
(4)以局部组织为网络节点，基于步骤(3)计算得灌注相似性为边生成表征肿瘤内血管生成差异性的网络图。
所述步骤(1)中，利用与超声造影模态(Contrast-enhanced ultrasound,CEUS)同步成像的灰度模式超声(B-mode ultrasound,BUS)序列跟踪肿瘤区域；给定动态BUS序列
t表示某一帧，T表示输入序列总帧数，将首帧
中包含病灶的感兴趣区域作为模板Z，在后续动态序列中，以前t-1帧的病灶跟踪结果作为初始区域，在当前t帧中确定待搜索区域Ct；通过参数共享的空间特征重标定卷积网络fθ提取模板Z和待搜索区域Ct对应的深度卷积特征张量FZ＝fθ(Z)和
fθ表示特征学习网络，fθ(Z)为模板Z的特征学习网络，fθ(Ct)为待搜索区域Ct的特征学习网络，并基于特征张量计算匹配响应图
特征学习网络fθ引入了空间特征重标定机制，该机制通过二维卷积学习对应的权重矩阵W，用于增加网络对病灶内关键组织结构的关注度，从而提高模板匹配的准确性；重标定后的模板Z的卷积特征为F′z＝fθ(z)·Wz，待搜索区域Ct的卷积特征的
Wz,
分别对应于不同尺寸的权重矩阵；然后，将模板Z的卷积特征F′z作为卷积核，与待搜索区域Ct的卷积特征
作二维滑动卷积，得到匹配响应图
考虑到BUS序列连续运动特点，融合前t-1帧的病灶跟踪结果作为初始匹配估计图It＝[l1,l2,l3,…,lt-1]，lt-1为第t-1帧的病灶跟踪结果，并将匹配响应图St与其点乘得到最终匹配响应图
从而消除误匹配，
中最大值位置(xt,yt)，xt,yt为t帧中病灶所在的位置，即为第t帧病灶跟踪结果lt；最终，根据连续灰度超声序列中病灶位置匹配结果，提取与第一帧具有相同病灶位置的造影子序列
t为某一帧，K为子序列长度，从而消除呼吸运动带来的干扰。
所述步骤(2)中，采用稳定主成分分析法提取局部灌注模式，该算法将观测到的局部灌注强度变化分解为稳定的微泡浓度变化和微泡随机破裂引起的不规则震荡；将病灶及周边区域等间隔划分，得到所有的灌注模式U∈RN×K，R为实数，N为局部灌注区域个数，其中，稳定、平滑的微泡浓度变化对应于U中的低秩成分L，不规则震荡对应于U中的稀疏成分S，考虑到超声造影成像中固有噪声ε，采用不等式约束的成分分析算法：
s.t.||L+S-U||2≤ε
其中，||·||*，||·||1和||·||2分别表示矩阵的核范数，1范数和2范数，α为权重系数。
所述步骤(3)中，量化基于步骤(2)中提取的局部灌注模式间差异，给定任意局部灌注模式p＝[p1,p2,p3,…pt…pK]和q＝[q1,q2,q3,…qt′…qK]，其中，t为灌注模式p的某一帧，pt是灌注模式p在t帧的灌注强度，t′为灌注模式q的某一帧，qt′是灌注模式q在t′帧的灌注强度，K为灌注模式的长度，pK是灌注模式p在K帧的灌注强度，qK是灌注模式q在K帧的灌注强度，基于灌注相位约束的灌注差异度量PCD(p,q)为，
其中，M为匹配矩阵，Mtt′表示灌注模式p中t帧和灌注模式q中t′帧的灌注强度匹配权重，dtt′＝(pt-qt′)2表示不同帧灌注强度差异，Q为匹配范围约束，KL(M||Q)为矩阵间Kullback-Leibler散度，用于给序列间元素匹配增加相位约束，λ为灵活度约束权重参数；其中，为了给不同灌注相位赋予不同的匹配灵活度，Q分解为
其中，
和
分别表示造影剂灌入相位和灌出相位的匹配范围约束，π为圆周率，σwash-in和σwash-out分别表示灌入相位和灌出相位的高斯带宽，用于调节匹配灵活度大小，通常设置灌入相位更低的匹配灵活度，σwash-in<σwash-out。因此，该模块实现了像素级灌注模式差异度量。
所述步骤(4)包括以下步骤：
步骤(401)，以病灶及周边正常实质的局部灌注区域为节点E，计算任意节点Ei,Ej间灌注模式差异PCD(Ei,Ej)，得到距离矩阵D＝{Di,j|Di,j＝PCD(Ei,Ej)}，Di,j为节点Ei和节点Ej的距离；
步骤(402)，确定每个节点Ei的周边节点，并基于步骤(401)计算的距离矩阵D推算节点间测地距离，构造反映肿瘤局部区域灌注相似性的邻接矩阵A；
步骤(403)，统计病灶内所有节点的度中心性DC、聚类系数CC作为节点特征，并基于该特征对肿瘤内网络节点进行聚类分组，然后，统计每个聚类中心的度中心性DC及聚类系数CC，均值用于表征肿瘤内血管生成的异质性，并将聚类结果的彩色编码图覆盖于原造影帧，增强异质性分析的可视化效果；
步骤(404)，对病灶边界进行腐蚀膨胀，获取周边正常实质节点，统计正常实质节点的度中心性DC及聚类系数CC，均值用于表征肿瘤血管对周边的正常组织的浸润性，并将数值的彩色编码图覆盖于原造影帧肿瘤区域，增强浸润性分析的可视化效果。
与现有技术相比，本发明具有以下有益效果：
本发明在消除呼吸运动干扰基础上，提取像素级灌注模式，结合灌注相位量化像素级灌注模式差异并生成结点图，基于结点的方式表征肿瘤内灌注差异性，进一步，为医生量化肿瘤血管生成提供定量参数。
附图说明
图1为本发明的系统结构示意图。
图2为本发明的病灶跟踪示意图。
具体实施方式
下面结合实施例对本发明作更进一步的说明。
如图1和2所示，一种基于灌注相位约束匹配的超声造影定量分析方法，包括：基于病灶跟踪的呼吸干扰抑制模块，通过基于空间特征重标定的孪生卷积网络跟踪灰度模式窗口中病灶区域，根据超声造影的双模同步成像特点，我们可以根据病灶跟踪结果选择具有相同病灶位置的造影子序列，从而消除序列中呼吸运动带来的干扰；单像素灌注模式提取模块，提出基于稳定主成分分析算法提取单像素稳定灌注模式，消除由于微泡随机破裂或非线性震荡带来局部抖动干扰；基于灌注相位约束匹配的TIC差异性度量模块，提出基于灌注相位约束匹配的移土距离用于量化单像素灌注模式间差异，该度量在灌注序列差异比较中提供了灵活的局部灌注强度匹配方式，同时结合临床诊断经验，为不同的灌注相位(灌入相位和灌出相位)赋予不同的匹配自由度，从而为肿瘤及周边实质血流灌注提供像素级差异性度量。基于图的肿瘤血管生成量化模块，该模块以单像素为节点，像素级灌注模式差异为边，构造结点图表征肿瘤内局部灌注差异性或相似性大小，进一步，基于病灶内节点聚类系数向量量化肿瘤血管生成异质性，通过病灶边界腐蚀膨胀定量病灶边界与周边实质血供差异，从而量化肿瘤血管浸润性。
实施例1
一种基于动态超声造影的肝肿瘤血管生成定量分析方法，包括如下步骤：
(1)通过基于空间特征重标定的孪生卷积网络跟踪病灶区域，从而消除呼吸运动对超声造影定量分析带来的干扰；
(2)利用步骤(1)获取的无呼吸干扰的超声造影子序列，通过稳定成分分析算法提取局部灌注模式；
(3)通过基于灌注相位约束匹配的相似性度量量化肿瘤内部及周边正常组织的局部灌注模式差异；
(4)以局部组织为网络节点，基于步骤(3)计算得灌注相似性为边生成表征肿瘤内血管生成差异性的网络图。
所述步骤(1)中，提出基于空间特征重标定的孪生卷积网络对病灶区域进行跟踪，并基于跟踪结果选择具有相同病灶区域的造影子序列。一方面，由于超声造影成像是双模态同步成像，即病灶在传统灰度模式视图中和造影视图中的运动完全同步。另一方面，在造影视图中，由于基波抑制作用，窗口仅显示造影增强部分的灰阶图像，而在造影起始阶段病灶特征不明显。因此，该模块病灶跟踪是在与造影同步成像的灰度模式序列中完成的。连续超声序列中病灶跟踪可以简化为，在相邻图像帧中，以像素为单位在搜索区域内匹配病灶目标，得到与搜索区域相同尺寸的相似度响应图，相似图中得分最高的区域即为病灶位置。该模块中，我们采用结合空间特征重标定机制的孪生卷积网络实现病灶跟踪，通过端到端的层次化灰度超声图像特征提取和相似性度量学习，确定与目标帧有相同病灶位置的造影子序列。相较于一般目标跟踪场景，灰度模式超声图像具有空间分辨率低，对比度低和结构特征不明显等特点，病灶跟踪中特征提取难度大，误匹配可能性高。因此，区别于已有的基于孪生网络的目标跟踪算法，我们考虑到病灶位置在超声序列中连续性特点，将连续帧病灶检测信息融合作为当前帧病灶位置的初始估计，而不是独立地在每帧图像上匹配模板，最终位置估计应该在初始估计点的指定范围内，从而减少无匹配的可能性。同时，为了增加模型对目标区域关键特征的关注度的和模型的可解释性，我们在全卷积网络中引入空间特征重标定机制，调整目标区域内不同位置语义特征权重。具体地讲，利用与超声造影模态(Contrast-enhanced ultrasound,CEUS)同步成像的灰度模式超声(B-mode ultrasound,BUS)序列跟踪肿瘤区域；给定动态BUS序列
t表示某一帧，T表示输入序列总帧数，将首帧
中包含病灶的感兴趣区域作为模板Z，在后续动态序列中，以前t-1帧的病灶跟踪结果作为初始区域，在当前t帧中确定待搜索区域Ct；通过参数共享的空间特征重标定卷积网络fθ提取模板Z和待搜索区域Ct对应的深度卷积特征张量Fz＝fθ(Z)和
fθ表示特征学习网络，fθ(Z)为模板Z的特征学习网络，fθ(Ct)为待搜索区域Ct的特征学习网络，并基于特征张量计算匹配响应图
特征学习网络fθ引入了空间特征重标定机制，该机制通过二维卷积学习对应的权重矩阵W调整局部特征权重，增加网络对病灶内关键组织结构的关注度，从而提高模板匹配的准确性；重标定后的模板Z的卷积特征为F′z＝fθ(z)·Wz，待搜索区域Ct的卷积特征的
Wz,
分别对应于不同尺寸的权重矩阵；Wz,
分别对应于不同尺寸的权重矩阵；然后，将模板Z的卷积特征F′z作为卷积核，与待搜索区域Ct的卷积特征
作二维滑动卷积，得到匹配响应图
考虑到BUS序列连续运动特点，融合前t-1帧的病灶跟踪结果作为初始匹配估计图It＝[l1,l2,l3,…,lt-1]，lt-1为第t-1帧的病灶跟踪结果，并将匹配响应图St与其对应位置元素点乘得到最终匹配响应图
从而消除误匹配，
中最大值位置(xt,yt)，xt,yt为t帧中病灶所在的位置，即为第t帧病灶跟踪结果lt；最终，根据连续灰度超声序列中病灶位置匹配结果，提取与第一帧具有相同病灶位置的造影子序列
t为某一帧，K为子序列长度，从而消除呼吸运动带来的干扰。
所述步骤(2)中，尽管通过呼吸运动校正，像素级灌注模式仅包含固位置造影剂的反射信号，但信号强度仍然受到微泡的不规则震荡和破裂的影响，导致造影帧序列中同一像素位置的强度存在不规则抖动，给后续像素级灌注模式的差异定量带来误差。因此，采用稳定主成分分析法提取局部灌注模式，该方法将观测到的局部灌注强度变化分解为稳定的微泡浓度变化和微泡随机破裂引起的不规则震荡；将病灶及周边区域等间隔划分，得到所有的灌注模式U∈RN×K，R为实数，N为局部灌注区域个数，其中，稳定、平滑的微泡浓度变化对应于U中的低秩成分L，不规则震荡对应于U中的稀疏成分S，考虑到超声造影成像中固有噪声ε，采用不等式约束的成分分析算法：
s.t.||L+S-U||2≤ε
其中，||·||*，||·||1和||·||2分别表示矩阵的核范数，1范数和2范数，α为权重系数。
所述步骤(3)中，提出灌注相位约束的移土距离用于量化像素级灌注模式差异。单像素灌注模式差异比较本质上是一种时间序列距离度量，该度量能够刻画造影剂在灌入期，达峰期和消退期不同像素位置灌注强度差异。时间序列距离度量是基于不同时间点元素对齐匹配实现的，灌注模式距离即不同时间点灌注强度在对齐匹配后差异的和。不同时间点的强度匹配是必要的，以相邻像素的灌注达峰值比较为例，即使对于相邻像素，其灌注强度序列的达峰时间点也不会完全对应于同一时刻，点对点的像素差异比较显然不符合临床灌注模式差异比较需求，局部时间范围内灌注强度匹配具有重要的临床意义。同时，为了保证差异性刻画，一定时间范围内匹配灵活度是有限制的，结合临床医生病灶诊断经验，我们在不同的灌注相位设置不同的匹配灵活度，从而使得距离度量更符合超声造影临床诊断实践。具体地讲，量化基于步骤(2)中提取的局部灌注模式间差异，给定任意局部灌注模式p＝[p1,p2,p3,…pt…pk]和q＝[q1,q2,q3,…qt，…qk]，其中，t为灌注模式p的某一帧，pt是灌注模式p在t帧的灌注强度，t’为灌注模式q的的某一帧，qt是灌注模式q在t’帧的灌注强度，k为灌注模式的长度，pk是灌注模式p在k帧的灌注强度，qk是灌注模式q在k帧的灌注强度，基于灌注相位约束的灌注差异度量PCD(p,q)为，
其中，M为匹配矩阵，Mtt’表示灌注模式p中t帧和灌注模式q中t’帧的灌注强度匹配权重，dtt’＝(pt-qt’)2表示不同帧灌注强度差异，Q为匹配范围约束，KL(M||Q)为矩阵间Kullback-Leibler散度，用于给序列间元素匹配增加相位约束，λ为灵活度约束权重参数；其中，为了给不同灌注相位赋予不同的匹配灵活度，Q分解为
其中，
和
分别表示造影剂灌入相位和灌出相位的匹配范围约束，π为圆周率，σwash-in和σwash-out分别表示灌入相位和灌出相位的高斯带宽，用于调节匹配灵活度大小，通常设置灌入相位更低的匹配灵活度，σwash-in<σwash-out。因此，该模块实现了像素级灌注模式差异度量。
所述步骤(4)中，提出了基于复杂网络的肿瘤血管生成量化方法。超声造影序列中像素强度大小反映了肿瘤及周边实质实时的血流灌注量，临床上，为肿瘤血管生成提供了定性分析工具。基于以上模块获取到的像素级灌注模式差异，生成表征局部灌注差异的肿瘤网络，从而提取反映肿瘤内血管生成异质性和浸润性的定量特征。相比较传统基于纹理特征的异质性分析方法，基于复杂网络的表征方法对肿瘤内微环境刻画具有更好的解释性和可视化效果，不同时间点的复杂网络变化为肿瘤治疗、预后提供了新的评估方法。具体地讲，包括以下步骤：
步骤(401)，以病灶及周边正常实质的局部灌注区域为节点E，计算任意节点Ei,Ej间灌注模式差异PCD(Ei,Ej)，得到距离矩阵D＝{Di,j|Di,j＝PCD(Ei,Ej)}，Di,j为节点Ei和节点Ej的距离；
步骤(402)，确定每个节点Ei的周边节点，并基于步骤(401)计算的距离矩阵D推算节点间测地距离，构造反映肿瘤局部区域灌注相似性的邻接矩阵A；
步骤(403)，统计病灶内所有节点的度中心性DC、聚类系数CC作为节点特征，并基于该特征对肿瘤内网络节点进行聚类分组，然后，统计每个聚类中心的度中心性DC及聚类系数CC，均值用于表征肿瘤内血管生成的异质性，并将聚类结果的彩色编码图覆盖于原造影帧，增强异质性分析的可视化效果；
步骤(404)，对病灶边界进行腐蚀膨胀，获取周边正常实质节点，统计正常实质节点的度中心性DC及聚类系数CC，均值用于表征肿瘤血管对周边的正常组织的浸润性，并将数值的彩色编码图覆盖于原造影帧肿瘤区域，增强浸润性分析的可视化效果。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种基于动态超声造影的肝肿瘤血管生成定量分析方法，其特征在于，包括如下步骤：
(1)通过基于空间特征重标定的孪生卷积网络跟踪病灶区域，从而消除呼吸运动对超声造影定量分析带来的干扰；
(2)利用步骤(1)获取的无呼吸干扰的超声造影子序列，通过稳定成分分析算法提取局部灌注模式；
所述步骤(2)中，在步骤(1)获取的运动校正子序列基础上，采用稳定成分分析算法提取局部灌注模式，将观测到的局部灌注强度变化分解为稳定的微泡浓度变化和微泡随机破裂引起的不规则震荡；将病灶及周边区域等间隔划分，得到所有的灌注模式U∈RN×K，R为实数，N为局部灌注区域个数，其中，稳定、平滑的微泡浓度变化对应于U中的低秩成分L，不规则震荡对应于U中的稀疏成分S，考虑到超声造影成像中固有噪声ε，采用不等式约束的成分分析算法：
s.t.||L+S-U||2≤ε
其中，‖·‖*，‖·‖1和‖·‖2分别表示矩阵的核范数，1范数和2范数，α为权重系数；
(3)通过基于灌注相位约束匹配的相似性度量量化肿瘤内部及周边正常组织的局部灌注模式差异；
(4)以局部组织为网络节点，基于步骤(3)计算得灌注相似性为边生成表征肿瘤内血管生成差异性的网络图；
所述步骤(4)包括以下步骤：
步骤(401)，以病灶及周边正常实质的局部灌注区域为节点E，计算任意节点Ei,Ej间灌注模式差异PCD(Ei,Ej)，得到距离矩阵D＝{Di,j|Di,j＝PCD(Ei,Ej)}，Di,j为节点Ei和节点Ej的距离；
步骤(402)，确定每个节点Ei的周边节点，并基于步骤(401)计算的距离矩阵D推算节点间测地距离，构造反映肿瘤局部区域灌注相似性的邻接矩阵A；
步骤(403)，统计病灶内所有节点的度中心性DC、聚类系数CC作为节点特征，并基于该特征对肿瘤内网络节点进行聚类分组，然后，统计每个聚类中心的度中心性DC及聚类系数CC，均值用于表征肿瘤内血管生成的异质性，并将聚类结果的彩色编码图覆盖于原造影帧，增强异质性分析的可视化效果；
步骤(404)，对病灶边界进行腐蚀膨胀，获取周边正常实质节点，统计正常实质节点的度中心性DC及聚类系数CC，均值用于表征肿瘤血管对周边的正常组织的浸润性，并将数值的彩色编码图覆盖于原造影帧肿瘤区域，增强浸润性分析的可视化效果。
2.根据权利要求1所述的基于动态超声造影的肝肿瘤血管生成定量分析方法，其特征在于，所述步骤(1)中，利用与超声造影模态CEUS同步成像的灰度模式超声BUS序列跟踪肿瘤区域；给定动态BUS序列
t表示某一帧，T表示输入序列总帧数，将首帧
中包含病灶的感兴趣区域作为模板Z，在后续动态序列中，以前t-1帧的病灶跟踪结果作为初始区域，在当前t帧中确定待搜索区域Ct；通过参数共享的空间特征重标定卷积网络fθ提取模板Z和待搜索区域Ct对应的深度卷积特征张量FZ＝fθ(Z)和
fθ表示特征学习网络，fθ(Z)为模板Z的特征学习网络，fθ(Ct)为待搜索区域Ct的特征学习网络，并基于特征张量计算匹配响应图
特征学习网络fθ引入了空间特征重标定机制，空间特征重标定机制通过二维卷积学习对应的权重矩阵W调整局部特征权重；重标定后的模板Z的卷积特征为F′z＝fθ(z)·Wz，待搜索区域Ct的卷积特征的
Wz，
分别对应于不同尺寸的权重矩阵；然后，将模板Z的卷积特征F′z作为卷积核，与待搜索区域Ct的卷积特征
作二维滑动卷积，得到匹配响应图
融合前t-1帧的病灶跟踪结果作为初始匹配估计图It＝[l1，l2，l3，...，lt-1]，lt-1为第t-1帧的病灶跟踪结果，并将匹配响应图St与其对应位置元素点乘得到最终匹配响应图
中最大值位置(xt，yt)，xt，yt为t帧中病灶所在的位置，即为第t帧病灶跟踪结果lt；最终，根据连续灰度超声序列中病灶位置匹配结果，提取与第一帧具有相同病灶位置的造影子序列
1＜t＜K，t为某一帧，K为子序列长度。
3.根据权利要求1所述的基于动态超声造影的肝肿瘤血管生成定量分析方法，其特征在于，所述步骤(3)中，量化基于步骤(2)中提取的局部灌注模式间差异，给定任意局部灌注模式p＝[p1，p2，p3，…pt…pk]和q＝[q1，q2，q3，…qt，…qk]，其中，t为灌注模式p的某一帧，pt是灌注模式p在t帧的灌注强度，t’为灌注模式q的的某一帧，qt是灌注模式q在t’帧的灌注强度，k为灌注模式的长度，pk是灌注模式p在k帧的灌注强度，qk是灌注模式q在k帧的灌注强度，基于灌注相位约束的灌注差异度量PCD(p，q)为，
其中，M为匹配矩阵，Mtt’表示灌注模式p中t帧和灌注模式q中t’帧的灌注强度匹配权重，dtt’＝(pt-qt’)2表示不同帧灌注强度差异，Q为匹配范围约束，KL(M||Q)为矩阵间Kullback-Leibler散度，用于给序列间元素匹配增加相位约束，λ为灵活度约束权重参数；其中，为了给不同灌注相位赋予不同的匹配灵活度，Q分解为
其中，
和
分别表示造影剂灌入相位和灌出相位的匹配范围约束，π为圆周率，σwash-in和σwash-out分别表示灌入相位和灌出相位的高斯带宽，用于调节匹配灵活度大小，通常设置灌入相位更低的匹配灵活度，σwash-in＜σwash-out。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911308530.1A CN110969618B (zh) | 2019-12-18 | 2019-12-18 | 一种基于动态超声造影的肝肿瘤血管生成定量分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911308530.1A CN110969618B (zh) | 2019-12-18 | 2019-12-18 | 一种基于动态超声造影的肝肿瘤血管生成定量分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110969618A CN110969618A (zh) | 2020-04-07 |
| CN110969618B true CN110969618B (zh) | 2023-05-12 |
Family
ID=70034776
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911308530.1A Active CN110969618B (zh) | 2019-12-18 | 2019-12-18 | 一种基于动态超声造影的肝肿瘤血管生成定量分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110969618B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112790741B (zh) * | 2021-01-27 | 2021-12-28 | 复旦大学 | 脑血管形态特征定量分析方法 |
| CN113808090A (zh) * | 2021-09-08 | 2021-12-17 | 中国科学院深圳先进技术研究院 | 一种磁共振图像处理方法、终端设备及计算机存储介质 |
| CN116612193B (zh) * | 2023-07-20 | 2023-09-15 | 北京大学第三医院（北京大学第三临床医学院） | 一种基于深度学习的血管造影图像的处理及标定方法 |
| CN116993628B (zh) * | 2023-09-27 | 2023-12-08 | 四川大学华西医院 | 一种用于肿瘤射频消融引导的ct图像增强系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101756713A (zh) * | 2009-09-09 | 2010-06-30 | 西安交通大学 | 超声造影成像、灌注参量估计和灌注参量功能成像及其集成方法 |
| CN104688272A (zh) * | 2015-03-06 | 2015-06-10 | 西安交通大学 | 一种基于单像素tic源的超声造影灌注参量成像方法 |
| CN106056589A (zh) * | 2016-05-24 | 2016-10-26 | 西安交通大学 | 一种呼吸运动补偿的超声造影灌注参量成像方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140071125A1 (en) * | 2012-09-11 | 2014-03-13 | The Johns Hopkins University | Patient-Specific Segmentation, Analysis, and Modeling from 3-Dimensional Ultrasound Image Data |
| CN108697354B (zh) * | 2016-02-23 | 2022-12-06 | 梅约医学教育与研究基金会 | 超声血流成像 |
| CN109952063B (zh) * | 2016-11-14 | 2022-04-29 | 皇家飞利浦有限公司 | 用于表征造影剂流的肝脏灌注的系统和方法 |
-
2019
- 2019-12-18 CN CN201911308530.1A patent/CN110969618B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101756713A (zh) * | 2009-09-09 | 2010-06-30 | 西安交通大学 | 超声造影成像、灌注参量估计和灌注参量功能成像及其集成方法 |
| CN104688272A (zh) * | 2015-03-06 | 2015-06-10 | 西安交通大学 | 一种基于单像素tic源的超声造影灌注参量成像方法 |
| CN106056589A (zh) * | 2016-05-24 | 2016-10-26 | 西安交通大学 | 一种呼吸运动补偿的超声造影灌注参量成像方法 |
Non-Patent Citations (1)
| Title |
|---|
| 吴玉泉 等.《感兴趣区的描绘对肝细胞肝癌超声造影灌注参数的影响》.《广西医科大学学报》.2014,第208-212页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110969618A (zh) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110969618B (zh) | 一种基于动态超声造影的肝肿瘤血管生成定量分析方法 | |
| CN107230206B (zh) | 一种基于多模态数据的超体素序列肺部图像的3d肺结节分割方法 | |
| EP3132419B1 (en) | Systems and methods for image-based object modeling using multiple image acquisitions or reconstructions | |
| US7925064B2 (en) | Automatic multi-dimensional intravascular ultrasound image segmentation method | |
| CN102763135B (zh) | 用于自动分割和时间跟踪的方法 | |
| US6674879B1 (en) | Echocardiography workstation | |
| CN105719324B (zh) | 图像处理装置以及图像处理方法 | |
| De Craene et al. | 3D strain assessment in ultrasound (straus): A synthetic comparison of five tracking methodologies | |
| US10376236B2 (en) | Ultrasound diagnostic apparatus, image processing apparatus, and image processing method | |
| Haas et al. | Segmentation of 3D intravascular ultrasonic images based on a random field model | |
| US20060078182A1 (en) | Methods and apparatus for analyzing ultrasound images | |
| US20080095417A1 (en) | Method for registering images of a sequence of images, particularly ultrasound diagnostic images | |
| KR101697880B1 (ko) | 적응적 임계값 영상분할을 이용한 ct 영상 동맥경화반 용적 측정방법 및 장치 | |
| WO2017148332A1 (zh) | 一种灌注分析方法与设备 | |
| JP7041147B2 (ja) | 造影剤流の肝灌流を特徴付けるシステム及び方法 | |
| Heyde et al. | Anatomical image registration using volume conservation to assess cardiac deformation from 3D ultrasound recordings | |
| De Luca et al. | Estimation of large-scale organ motion in B-mode ultrasound image sequences: a survey | |
| AU2017358321A1 (en) | Method and system for modelling a human heart and atria | |
| EP1090372A2 (en) | Echocardiography workstation | |
| US6909914B2 (en) | Method for generating time independent images of moving objects | |
| TWI788629B (zh) | 圖像處理方法、裝置及系統、電子設備及電腦可讀儲存媒體 | |
| Lin et al. | Artificial intelligence in cardiovascular imaging: enhancing image analysis and risk stratification | |
| Babin et al. | Robust segmentation methods with an application to aortic pulse wave velocity calculation | |
| Ma et al. | Left ventricle segmentation from contrast enhanced fast rotating ultrasound images using three dimensional active shape models | |
| Noble et al. | Automated, nonrigid alignment of clinical myocardial contrast echocardiography image sequences: comparison with manual alignment |
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 |