CN109903293A - 一种利用血流灌注参数的超声映射图进行组织分割的方法 - Google Patents
一种利用血流灌注参数的超声映射图进行组织分割的方法 Download PDFInfo
- Publication number
- CN109903293A CN109903293A CN201910072643.XA CN201910072643A CN109903293A CN 109903293 A CN109903293 A CN 109903293A CN 201910072643 A CN201910072643 A CN 201910072643A CN 109903293 A CN109903293 A CN 109903293A
- Authority
- CN
- China
- Prior art keywords
- parameter
- matrix
- blood perfusion
- perfusion parameter
- degree
- 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.)
- Pending
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明实施例公开了一种利用血流灌注参数的超声映射图进行组织分割的方法,通过算法对超声造影视频进行计算,获取视频中各个像素点的达峰时间,峰值强度,上升速率,上升时间等血流灌注参数,同时依据参数在图像中的位置及参数大小构建参数矩阵。基于参数大小分布,利用模糊聚类算法将参数划分不同类的区间,每个区间对应一种映射颜色。由于参数矩阵中元素的位置与图像像素点位置是一一对应的,结合像素点对应位置的参数的大小及大小所处分类区间的对应颜色,对图像中所有像素点编码对应的颜色,从而实现了对原图像不同区域映射不同颜色,使得生物组织结构差异性大的地方颜色差异大。从而获得超声造影血流灌注参数彩色映射图,实现了病灶分割。
Description
技术领域
本发明涉及人工智能技术中图像分割领域,尤其涉及一种利用血流灌注参数的超声映射图进行组织分割的方法。
背景技术
目前国内外广泛应用的医学图像分割方法有很多种,许多学者应用数学、信号处理、人工智能等领域知识拓展医学图像分割的理论方法。常用的分割方法主要可以分为几类:基于阈值的分割、基于边缘的分割、基于运动的分割、基于活动轮廓模型的分割方法、基于小波变换的方法以及深度学习算法的分割等。这些分割方法均是从视频与图像本身具有的特性的角度去完成分割。
然而通过研究发现,组织结构相似位置的灌注参数值大小相近,而组织结构差异性大的地方,灌注参数值大小相差很大。反映到该参数矩阵中,就是矩阵中的数据具有分块性,某一块数据大小相近,可以聚成一类。但是在聚类过程中,某些区间块边缘数据,既可以属于这一聚类,亦可属于另一聚类,存在一个模糊的区间,难以实现对组织的精确分割。
发明内容
本发明实施例所要解决的技术问题在于,提供一种利用血流灌注参数的超声映射图进行组织分割的方法。可将模糊C均值聚类算法用于血流灌注参数的离散数据矩阵进行模糊聚类分区,从而获得了明晰的分割边界,很好地实现对不同组织的精确分割。
为了解决上述技术问题,本发明实施例提供了一种利用血流灌注参数的超声映射图进行组织分割的方法,包括以下步骤:
S1:设定聚类个数k(1<k<n),模糊指数m(m>1),最大迭代次数为T,收敛的误差为△;
S2:用随机数初始化隶属度矩阵UC;
S3:将UC代入到以下公式中,计算聚类中心矩阵X:
S4:将步骤S3计算得到的聚类中心矩阵X代入以下公式,重新计算隶属度矩阵UC+1:
S5:判断|UC+1-UC|≤△,满足条件,则输出U和X;否则,利用隶属度矩阵UC+1继续执行步骤S3~S5,在满足小于最大迭代次数条件下,直到满足步骤S5判断条件,输出U和X为止;
S6:对于最后一次输出的U矩阵中,每列的最大值即为该样本所属类别;
S7:计算每类所辖所有样本里面的的最大值dMAX和最小值dMIN,那么,对于每个类别的区间dMIN≤d≤dMAX,即为彩色映射的区间。
进一步地,还包括将血流灌注参数的离散数据矩阵转换为数据集D:
D={d1,d2,d3,……,dj,……dn}(1≤j≤n),n=p×q
将所述数据集D中的数据可以分成k类:
{S1,S2,S3,……,Si,……,Sk}(1≤i≤k)
定义隶属度矩阵为U,U中的任一元素uij是数据集D的第j个样本dj对于第i类Si的隶属度函数:
定义模糊C均值聚类目标函数及隶属度函数uij的约束条件分别如下:
其中,uij为隶属度值,m为模糊指数,||dj-xi||2指的就是欧氏距离下样本dj跟中心点xi之间的距离。
更进一步地,所述模糊指数m取值为2。
更进一步地,所述△设定为0.05。
实施本发明实施例,具有如下有益效果:从组织结构生物特性出发,利用不同组织在造影过程中血流灌注参数(起始时间,达峰时间,峰值强度,上升速率等)的差异特性来实现分割,将模糊C均值聚类算法用于血流灌注参数的离散数据矩阵进行模糊聚类分区,从而获得了明晰的分割边界,很好地实现对不同组织的精确分割。
附图说明
图1是造影过程中血流灌注参数时间-强度曲线;
图2是用超声造影颜色编码映射区域的外沿边界像素进行分割示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
超声造影(CEUS)是纯血池显像,通过向周围静脉内注射直径与红细胞相似的气体微泡造影剂,利用气体较强的散射性及与人体组织不同的声学特性,增大血流或病变与相邻组织间的声阻抗差异,从而获得反差较大的声像图,可动态直观地显示肿瘤微血管,已广泛应用于各种肿瘤癌症的诊断和疗效评估中。
对于人体某组织超声造影过程中造影强度的变化,反映到视频中,就是图像像素灰度的变化。采用时间-强度曲线(TimeIntensityCurve,TIC)定量分析法,如图1所示,可获得起始时间(造影剂微泡刚刚进入肿瘤内的时间)、峰值强度(TIC上显示的最大的增强强度)、达峰时间(从造影开始起至造影剂到达峰值强度的一段时间)、上升速率(造影过程中的造影强度升高的平均速率)等定量的血流灌注参数指标。由于异常组织与正常组织生物组织结构的差异性(如:恶性肿瘤相对于良性及正常区域来说,其血流灌注参数具有:达峰时间短,峰值强度高,上升速率大,上升时间短等特点),定量分析良恶性肿瘤TIC曲线的血流灌注参数指标,发现良恶性肿瘤血流灌注参数指标差异均有统计学意义。
超声造影(CEUS)是纯血池显像,通过向周围静脉内注射直径与红细胞相似的气体微泡造影剂,利用气体较强的散射性及与人体组织不同的声学特性,增大血流或病变与相邻组织间的声阻抗差异,从而获得反差较大的声像图,可动态直观地显示肿瘤微血管,已广泛应用于各种肿瘤癌症的诊断和疗效评估中。
对于人体某组织的超声造影过程,采用时间-强度曲线(Time Intensity Curve,TIC)定量分析法,可获得起始时间(造影剂微泡刚刚进入肿瘤内的时间)、峰值强度(TIC上显示的最大的增强强度)、达峰时间(从造影开始起至造影剂到达峰值强度的一段时间)、上升速率(从造影开始起至造影剂到达峰值强度的一段时间)等定量的血流灌注参数指标。由于异常组织与正常组织生物组织结构的差异性(如:恶性肿瘤相对于良性及正常区域来说,其血流灌注参数具有:达峰时间短,峰值强度高,流入速率大,上升时间短等特点),定量分析良恶性肿瘤TIC曲线的血流灌注参数指标,发现良恶性肿瘤血流灌注参数指标差异均有统计学意义。
具体实施步骤如下。
对视频分辨率大小为p×q的某块组织的超声造影视频进行计算,可以获得任一血流灌注参数的离散数据矩阵。
矩阵(1)中所示的v11-vqp为计算机分析计算得到的某一血流灌注参数,即造影过程中,该位置造影强度变化的的平均速率,反映到图像中为图像像素点灰度变化的平均速率。
将该血流灌注参数的离散数据矩阵转换为数据集D:
[v11,v12,v13……v1q,v21,v22,v23……,v2q,v31,v32,v33……,v3q,……,vp1,vp2,vp3……,vpq] (2)
该数据集D也可表征如下:
D={d1,d2,d3,……,dj,……dn}(1≤j≤n),n=p×q (3)
为了方便行文说明,将行向量(2)用行向量(3)表示,其中,d1=v11,d2=v22,……,dn=vpq,里面的元素是一一对应相等的。
假设该数据集中的数据可以分成k类:
{S1,S2,S3,……,Si,……,Sk}(1≤i≤k), (4)
其中,Si(1≤i≤k)表示为任意一个类别。k为不小于1的整数。
设X为所有类的聚类中心的集合
X={x1,x2,……,xi,……,xk) (5)
其中,x1,x2,……,xi,……,xk分别是S1,S2,S3,……,Si,……,Sk各个类的数据中心,简称聚类中心。
定义隶属度矩阵为U,如(6)所示。U中的任一元素uij是数据集D的第j个样本dj对于第i类Si的隶属度函数(一个样本属于某个类别的相似程度称为样本的隶属度,隶属度越大,说明样本更倾向于属于该类别)。
定义模糊C均值聚类目标函数及隶属度函数uij的约束条件分别如下:
其中,uij指的是隶属度函数(数据集D中的样本dj对于类别Si的隶属程度),m为模糊指数(一般取值>1,本文专利中设定为2),||dj-xi||2指的就是欧氏距离下样本dj跟中心点xi之间的距离。
具体分割步骤如下:
首先用随机数初始化隶属度矩阵U,用(9)式计算各个聚类的中心:
其中,xi表示的是任意类别Si的聚类中心,dj为数据集D中的样本,uij指的是隶属度函数(数据集D中的样本dj对于类别Si的隶属程度),m为模糊指数(一般取值>1,本文专利中设定为2)。
定义隶属度函数如(10)式,将(9)式的计算得到的聚类中心X代入(10)式,得到更新隶属度函数U:
反复迭代执行这两步操作,当更新的隶属度矩阵UC+1与上次迭代隶属度矩阵UC满足下式:
|UC+1-UC|≤△,本专利中△设定为0.05
则所有的样本隶属度值收敛趋于稳定,算法停止迭代。
此时,各样本对于各类的隶属度都具有确定的值,求取每个样本对于各个类隶属度中的最大值,该最大隶属度所在的类即为该样本数据所属的类,这样就可以实现所有样本的归类,此时每一个类别Si中都包含了该类别的样本数据;在此基础上,求取每个聚类Si所包含样本的最大值和最小值,作为该类别的聚类区间,由于每个样本在组织位置上都有相应的坐标,相近坐标的样本若处于同一聚类区间,则该处样本将聚类成一个区域,该区域会具有明确的边界,将每一个聚类区域标记一种映射的颜色,以该区域颜色的外沿边界像素作为分割的边界,完成对特定区域的分割。如图2所示。
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (4)
1.一种利用血流灌注参数的超声映射图进行组织分割的方法,其特征在于,包括以下步骤:
S1:设定聚类个数k(1<k<n),模糊指数m(m>1),最大迭代次数为T,收敛的误差为△;
S2:用随机数初始化隶属度矩阵UC;
S3:将UC代入到以下公式中,计算聚类中心矩阵X:
S4:将步骤S3计算得到的聚类中心矩阵X代入以下公式,重新计算隶属度矩阵UC+1:
S5:判断|UC+1-UC|≤△,满足条件,则输出U和X;否则,利用隶属度矩阵UC+1继续执行步骤S3 ̄S5,在满足小于最大迭代次数条件下,直到满足步骤S5判断条件,输出U和X为止;
S6:对于最后一次输出的U矩阵中,每列的最大值即为该样本所属类别;
S7:计算每类所辖所有样本里面的的最大值dMAX和最小值dMIN,那么,对于每个类别的区间dMIN≤d≤dMAX,即为彩色映射的区间。
2.根据权利要求1所述的利用血流灌注参数的超声映射图进行组织分割的方法,其特征在于,还包括将血流灌注参数的离散数据矩阵转换为数据集D:
D={d1,d2,d3,……,dj,……dn}(1≤j≤n),n=p×q
将所述数据集D中的数据可以分成k类:
{S1,S2,S3,……,Si,……,Sk}(1≤i≤k)
定义隶属度矩阵为U,U中的任一元素uij是数据集D的第j个样本dj对于第i类Si的隶属度函数:
定义模糊C均值聚类目标函数及隶属度函数uij的约束条件分别如下:
其中,uij为隶属度值,m为模糊指数,||dj-xi||2指的就是欧氏距离下样本dj跟中心点xi之间的距离。
3.根据权利要求2所述的利用血流灌注参数的超声映射图进行组织分割的方法,所述模糊指数m取值为2。
4.根据权利要求1-3任一项所述的利用血流灌注参数的超声映射图进行组织分割的方法,所述△设定为0.05。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910072643.XA CN109903293A (zh) | 2019-01-24 | 2019-01-24 | 一种利用血流灌注参数的超声映射图进行组织分割的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910072643.XA CN109903293A (zh) | 2019-01-24 | 2019-01-24 | 一种利用血流灌注参数的超声映射图进行组织分割的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109903293A true CN109903293A (zh) | 2019-06-18 |
Family
ID=66944135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910072643.XA Pending CN109903293A (zh) | 2019-01-24 | 2019-01-24 | 一种利用血流灌注参数的超声映射图进行组织分割的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109903293A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110458836A (zh) * | 2019-08-16 | 2019-11-15 | 深圳开立生物医疗科技股份有限公司 | 一种超声造影成像方法、装置和设备及可读存储介质 |
CN112807018A (zh) * | 2021-01-08 | 2021-05-18 | 武汉中旗生物医疗电子有限公司 | 一种超声造影剂到达时间成像方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105809670A (zh) * | 2016-02-29 | 2016-07-27 | 上海联影医疗科技有限公司 | 灌注分析方法 |
CN105976379A (zh) * | 2016-05-11 | 2016-09-28 | 南京邮电大学 | 一种基于布谷鸟优化的模糊聚类彩色图像分割方法 |
CN107203785A (zh) * | 2017-06-02 | 2017-09-26 | 常州工学院 | 多路径高斯核模糊c均值聚类算法 |
CN107301645A (zh) * | 2017-06-20 | 2017-10-27 | 上海联影医疗科技有限公司 | 一种数据处理方法和装置 |
-
2019
- 2019-01-24 CN CN201910072643.XA patent/CN109903293A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105809670A (zh) * | 2016-02-29 | 2016-07-27 | 上海联影医疗科技有限公司 | 灌注分析方法 |
CN105976379A (zh) * | 2016-05-11 | 2016-09-28 | 南京邮电大学 | 一种基于布谷鸟优化的模糊聚类彩色图像分割方法 |
CN107203785A (zh) * | 2017-06-02 | 2017-09-26 | 常州工学院 | 多路径高斯核模糊c均值聚类算法 |
CN107301645A (zh) * | 2017-06-20 | 2017-10-27 | 上海联影医疗科技有限公司 | 一种数据处理方法和装置 |
Non-Patent Citations (1)
Title |
---|
张军毅: "弥散敏感度b值在脑梗塞多谱MR图像分割中的应用研究", 《万方数据库》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110458836A (zh) * | 2019-08-16 | 2019-11-15 | 深圳开立生物医疗科技股份有限公司 | 一种超声造影成像方法、装置和设备及可读存储介质 |
CN112807018A (zh) * | 2021-01-08 | 2021-05-18 | 武汉中旗生物医疗电子有限公司 | 一种超声造影剂到达时间成像方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111539930B (zh) | 基于深度学习的动态超声乳腺结节实时分割与识别的方法 | |
CN107578416B (zh) | 一种由粗到精级联深度网络的心脏左心室全自动分割方法 | |
Shahangian et al. | Automatic brain hemorrhage segmentation and classification algorithm based on weighted grayscale histogram feature in a hierarchical classification structure | |
CN109215032A (zh) | 图像分割的方法及系统 | |
CN102800089B (zh) | 基于颈部超声图像的主颈动脉血管提取和厚度测量方法 | |
CN109087703B (zh) | 基于深度卷积神经网络的腹腔ct图像腹膜转移标记方法 | |
CN108257135A (zh) | 基于深度学习方法解读医学图像特征的辅助诊断系统 | |
CN107203989A (zh) | 基于全卷积神经网络的端对端胸部ct图像分割方法 | |
CN107977952A (zh) | 医学图像分割方法及装置 | |
CN112785598B (zh) | 基于注意力增强改进u型网络的超声乳腺肿瘤自动分割方法 | |
CN113570627B (zh) | 深度学习分割网络的训练方法及医学图像分割方法 | |
CN110473207A (zh) | 一种交互式分割肺叶的方法 | |
CN106780453A (zh) | 一种基于深度信任网络实现对脑肿瘤分割的方法 | |
CN109903293A (zh) | 一种利用血流灌注参数的超声映射图进行组织分割的方法 | |
Senthilkumaran et al. | Brain image segmentation | |
CN110473206B (zh) | 一种基于超体素与测度学习的弥散张量图像分割方法 | |
CN105869169B (zh) | 一种肿瘤组织微排列图像的自动分割方法 | |
CN114445429A (zh) | 基于多标签与多解码器的全心脏ct分割方法及装置 | |
CN111986216B (zh) | 一种基于神经网络改进的rsg肝脏ct图像交互式分割算法 | |
CN108876789A (zh) | 一种基于超像素和邻域块特征结合的连续图割方法 | |
Shen et al. | CT medical image segmentation algorithm based on deep learning technology | |
Patel et al. | Medical image enhancement using histogram processing and feature extraction for cancer classification | |
Cameron et al. | Patient-Specific Anatomic Models-Geometric Surface Generation from Three-Dimensional Medical Images Using a Specified Polygonal Budget | |
Peng et al. | Lung contour detection in chest X-ray images using mask region-based convolutional neural network and adaptive closed polyline searching method | |
Frannita et al. | Thyroid nodule classification based on characteristic of margin using geometric and statistical features |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190618 |
|
RJ01 | Rejection of invention patent application after publication |