基于超声造影的时空灌注异质性特征计算系统及其应用技术方案

本发明专利技术提供了一种基于超声造影的时空灌注异质性特征计算系统及其应用，系统涉及数据获取、关键帧的选择、超声造影图像的运动矫正、感兴趣区间的定义、像素级时间强度曲线提取以及时空灌注抑制性参数计算，基于感兴趣区间内超声造影信号强度数据在空间上的数据分布离散度，在时间强度曲线图中划定两条分位数曲线，求取分位数曲线之间围成的面积，通过该曲线间面积反映感兴趣区血管灌注的时空异质性特征。本发明专利技术借助曲线间面积定量化表征超声造影中肿瘤的造影剂灌注时空异质性特征，可用于评估肿瘤内血管时空异质性。相较于现有方案，保留了肿瘤区域内部血管分布的差异性特征，避免对肿瘤灌注信息平均化，有效评估恶性肿瘤血管的时空异质性。管的时空异质性。管的时空异质性。

【技术实现步骤摘要】
基于超声造影的时空灌注异质性特征计算系统及其应用


[0001]本专利技术涉及医学数据处理领域，具体而言，涉及一种基于超声造影的时空灌注异质性特征计算系统及其应用。

技术介绍

[0002]超声造影(Contrast Ultrasound)又称声学造影，是利用造影剂使后散射回声增强，明显提高超声诊断的分辨力、敏感性和特异性的技术。双幅模式的超声造影可以同时显示灰阶超声图像以及超声造影，帮助医师在疾病的超声造影检查中实时定位病灶的位置，提高对结果的判读与诊断能力。
[0003]由于超声造影剂为纯血池造影剂，造影剂在感兴趣区间的信号就代表着感兴趣区间内的血流信号。因此超声造影可以观察感兴趣区内的血流情况，包括血流速度、方向、到达先后时间、血供丰度等，从而了解肿瘤内部的血管分布和血流灌注模式情况。
[0004]造影剂的灌注异质性是指在病变的情况下，感兴趣区间内的造影剂信号可表现为多种形式，可呈均匀或不均匀、正态分布或偏态分布，尤其是在恶性肿瘤的生长过程中，肿瘤在不断增殖和分化过程中存在血管结构紊乱、扭曲、动静脉短路等差异性的血液灌注现象被称为血管灌注时空异质性，在超声造影下表现为造影剂灌注的时空异质性。灌注时空异质性可以更好地了解良恶性肿瘤内部的血流灌注情况，用来判断和评估肿瘤的增殖活性和恶性程度，有望用于指导放疗和化疗的治疗计划。
[0005]然而，临床上超声造影血流灌注的评估存在较强的主观性，对诊断医生经验依赖性较高，不同医生诊断结论间的一致性虽然尚可，但并非完全一致。并且，超声造影视频虽然含有大量的信息，但提供的帮助诊断的指标较少，目前主要围绕定性指标，如动脉期高增强、快速消退、不均匀灌注等，而定量指标如时间强度曲线相关参数，提供的是整个感兴趣区域的信息，相当于对感兴趣区域做了平均化，未曾表征区域内部造影增强强度的分布，不足以准确描述病灶的时空异质性。
[0006]鉴于此，本专利技术提出了一种基于超声造影的时空灌注异质性特征计算系统及其应用。

技术实现思路

[0007]基于现有技术存在的问题，本专利技术提供了一种基于超声造影的时空灌注异质性特征计算系统及其应用。具体方案如下：
[0008]一种基于超声造影的时空灌注异质性特征计算系统，包括如下：
[0009]数据获取单元，用于获取包含灰阶图像和造影图像的超声造影视频；
[0010]关键帧单元，用于将所述造影图像的某一帧，分别与其他帧进行相似度计算，筛选出在结节区域与其他帧相似度最高的帧，并将其作为关键帧；
[0011]运动矫正单元，用于裁剪相同大小、位置的灰阶图像和造影图像，基于关键帧对所述灰阶图像和所述造影图像进行运动矫正，得到第一图像集；
[0012]区间划分单元，用于在所述第一图像集中定义感兴趣区间，确保所述结节区域处于所述感兴趣区间，并以预设方形像素块为单位将所述感兴趣区间划分成多个像素单元；
[0013]曲线构建单元，用于计算各个像素单元的信号强度，并构建出各像素单元的信号强度随时间变化的、关于超声造影灌注曲线的曲线图；
[0014]特征计算单元，用于对各个像素单元随时间变化的信号强度进行离散度计算，以进行定量化时间
‑
离散度，并构建两条分位数曲线，通过两条分位数曲线之间的曲线间面积反映感兴趣区间内血管分布的时空异质性特征。
[0015]在一个具体实施例中，在曲线构建单元中，构建曲线图之前，对时间和各像素单元的信号强度进行归一化处理，使曲线间面积范围在0
‑
1之间。
[0016]在一个具体实施例中，取每一个时间点下各个像素单元的信号强度值，求取可表征数据离散度的统计学指标。
[0017]在一个具体实施例中，所述统计学指标包括四分位数间距、95％置信区间、变异系数、标准差或标准误差。
[0018]在一个具体实施例中，当以四分位数曲线为所述预设规则中的统计学指标时：
[0019]分别取所述超声造影灌注曲线上每个时间点的信号强度的25％和75％分位数做成两个曲线，得到一条表示四分之一分位的分位数曲线和一条表示四分之三分位的分位数曲线，通过求取两条分位数曲线与纵坐标之间的面积得到曲线间面积。
[0020]在一个具体实施例中，所述方形像素块的大小包括2*2、3*3或4*4。
[0021]在一个具体实施例中，所述运动矫正单元所采用的运动矫正算法包括滤波去运动、时间延迟去运动、空间滤波去运动、运动补偿去运动或基于深度学习的运动矫正。
[0022]在一个具体实施例中，在所述曲线构建单元中，对像素单元中各信号强度去零后取中位数即为该像素单元的信号强度。
[0023]一种计算机设备，所述计算机设备包括：
[0024]一个或多个处理器；
[0025]存储器，用于存储一个或多个程序；
[0026]当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如下处理：
[0027]获取包含灰阶图像和造影图像在内的超声造影视频；
[0028]将造影图像的某一帧，分别与其他帧进行相似度计算，筛选出在结节区域与其他帧相似度最高的帧，并将其作为关键帧；
[0029]裁剪相同大小、位置的灰阶图像和造影图像，基于关键帧对所述灰阶图像和所述造影图像进行运动矫正，得到第一图像集；
[0030]在所述第一图像集中定义感兴趣区间，确保所述结节区域处于所述感兴趣区间，以预设方形像素块为单位将感兴趣区间划分成多个像素单元；
[0031]计算各个像素单元的信号强度，并构建出各像素单元的信号强度随时间变化的、关于超声造影灌注曲线的曲线图；
[0032]对各个像素单元随时间变化的信号强度进行离散度计算，以进行定量化时间
‑
离散度，并构建两条分位数曲线，通过两条分位数曲线之间的曲线间面积反映感兴趣区间内血管分布的时空异质性特征。
[0033]一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如下处理：
[0034]获取包含灰阶图像和造影图像在内的超声造影视频；
[0035]将造影图像的某一帧，分别与其他帧进行相似度计算，筛选出在结节区域与其他帧相似度最高的帧，并将其作为关键帧；
[0036]裁剪相同大小、位置的灰阶图像和造影图像，基于关键帧对所述灰阶图像和所述造影图像进行运动矫正，得到第一图像集；
[0037]在所述第一图像集中定义感兴趣区间，确保所述结节区域处于所述感兴趣区间，以预设方形像素块为单位将感兴趣区间划分成多个像素单元；
[0038]计算各个像素单元的信号强度，并构建出各像素单元的信号强度随时间变化的、关于超声造影灌注曲线的曲线图；
[0039]对各个像素单元随时间变化的信号强度进行离散度计算，以进行定量化时间
‑
离散度，并构建两条分位数曲线，通过两条分位数曲线之间的曲线间面积反映感兴趣区间内血管分布的时空异质性特征。
[0040]有益效本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于超声造影的时空灌注异质性特征计算系统，其特征在于，包括如下：数据获取单元，用于获取包含灰阶图像和造影图像的超声造影视频；关键帧单元，用于将所述造影图像的某一帧，分别与其他帧进行相似度计算，筛选出在结节区域与其他帧相似度最高的帧，并将其作为关键帧；运动矫正单元，用于裁剪相同大小、位置的灰阶图像和造影图像，基于关键帧对所述灰阶图像和所述造影图像进行运动矫正，得到第一图像集；区间划分单元，用于在所述第一图像集中定义感兴趣区间，确保所述结节区域处于所述感兴趣区间，并以预设方形像素块为单位将所述感兴趣区间划分成多个像素单元；曲线构建单元，用于计算各个像素单元的信号强度，并构建出各像素单元的信号强度随时间变化的、关于超声造影灌注曲线的曲线图；特征计算单元，用于对各个像素单元随时间变化的信号强度进行离散度计算，以进行定量化时间
‑
离散度，并构建两条分位数曲线，通过两条分位数曲线之间的曲线间面积反映感兴趣区间内造影剂灌注的时空异质性特征。2.根据权利要求1所述的时空灌注异质性特征计算系统，其特征在于，在曲线构建单元中，构建曲线图之前，对时间和各像素单元的信号强度进行归一化处理，使曲线间面积范围在0
‑
1之间。3.根据权利要求1所述的时空灌注异质性特征计算系统，其特征在于，取每一个时间点下各个像素单元的信号强度值，求取可表征数据离散度的统计学指标。4.根据权利要求3所述的时空灌注异质性特征计算系统，其特征在于，所述统计学指标包括四分位数间距、95％置信区间、变异系数、标准差或标准误差。5.根据权利要求4所述的时空灌注异质性特征计算系统，其特征在于，当以四分位数曲线为所述预设规则中的统计学指标时：分别取所述超声造影灌注曲线上每个时间点的信号强度的25％和75％分位数做成两个曲线，得到一条表示四分之一分位的分位数曲线和一条表示四分之三分位的分位数曲线，通过求取两条分位数曲线与纵坐标之间的面积得到曲线间面积。6.根据权利要求1所述的时空灌注异质性特征计算系统，其特征在于，所述方形像素块的大小包括2*2、3*3或4*4。7.根据权利要求1所述的时空灌注异质性特征计算系统，其特征在于，所述运动矫正单元所采用的运动矫正...

【专利技术属性】
技术研发人员：黄品同，陈继繁，王亚华，林子梅，傅钰涵，
申请(专利权)人：浙江大学，
类型：发明
国别省市：