CN113034426A - 超声图像病灶描述方法、装置、计算机设备和存储介质 - Google Patents
超声图像病灶描述方法、装置、计算机设备和存储介质 Download PDFInfo
- Publication number
- CN113034426A CN113034426A CN201911358149.6A CN201911358149A CN113034426A CN 113034426 A CN113034426 A CN 113034426A CN 201911358149 A CN201911358149 A CN 201911358149A CN 113034426 A CN113034426 A CN 113034426A
- Authority
- CN
- China
- Prior art keywords
- lesion
- contour
- region
- interest
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000003902 lesion Effects 0.000 claims abstract description 210
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 62
- 238000003709 image segmentation Methods 0.000 claims abstract description 30
- 230000006870 function Effects 0.000 claims description 66
- 238000002604 ultrasonography Methods 0.000 claims description 37
- 230000011218 segmentation Effects 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 14
- 230000002308 calcification Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- 230000012010 growth Effects 0.000 claims description 9
- 238000011478 gradient descent method Methods 0.000 claims description 6
- 238000012805 post-processing Methods 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 abstract description 4
- 230000002829 reductive effect Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 21
- 238000001514 detection method Methods 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 10
- 230000001788 irregular Effects 0.000 description 8
- 208000026310 Breast neoplasm Diseases 0.000 description 7
- 206010006187 Breast cancer Diseases 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 4
- 238000007493 shaping process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000035040 seed growth Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008571 general function Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 210000000952 spleen Anatomy 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001685 thyroid gland Anatomy 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- 239000013598 vector Substances 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/155—Segmentation; Edge detection involving morphological operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30068—Mammography; Breast
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Geometry (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Quality & Reliability (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Analysis (AREA)
Abstract
本申请涉及一种超声图像病灶描述方法、装置、计算机设备和存储介质。所述方法包括:识别超声图像中的病灶,以得到对应的感兴趣区域;通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓;根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。其方法通过对超声图像中的病灶进行识别,以定位对应的感兴趣区域,进而通过图像分割算法对感兴趣区域进行检测,以得到病灶轮廓,并根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征自动生成病灶描述信息,从而有效的减少了医生的工作量,且提高了诊断效率。
Description
技术领域
本申请涉及医学图像处理技术领域,特别是涉及一种超声图像病灶描述方法、装置、计算机设备和存储介质。
背景技术
乳腺癌是女性疾病中常见的恶性肿瘤,已成为严重威胁女性健康的病症之一。早发现、早诊断、早治疗是目前医学上对防治乳腺癌采取的基本原则。超声成像凭借其无创伤、无辐射、费用低廉等优势,已成为乳腺肿瘤临床诊断的主要手段之一。
然而,由于受成像设备的影响,超声图像常常具有较大的噪声、低对比度、灰度不均匀、不同程度的衰减以及浸润效应等,使得乳腺肿瘤在表面上与周围正常组织较为相似,即超声图像对人体器官形态的表达能力较弱,或者说器官在图像中的呈现形式较为模糊与抽象;此外,不同个体间乳腺肿瘤差异也较大。因此,乳腺超声图像中对于病灶区域的判断和阅读需要临床医生具有较高的专业水平和较丰富的经验,一般医生较难在超声图像中用肉眼准确、快速地将乳腺肿瘤区域与其周围的正常组织区分开来并做出病灶描述结果。尤其是当需要对大批量的患者进行病灶分析时,医生既要对病灶进行分析,同时还要手动填写病灶分析结果。当患者人数较多时,医生的工作量会急剧加大。
发明内容
基于此,有必要针对上述医生对超声图像的病灶进行分析时工作量较大的问题,提供一种超声图像病灶描述方法、装置、计算机设备和存储介质。
为了实现上述目的,一方面,本申请实施例提供了一种超声图像病灶描述方法,所述方法包括:
识别超声图像中的病灶,以得到对应的感兴趣区域;
通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓;
根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。
另一方面,本申请实施例还提供了一种超声图像病灶描述装置,包括:感兴趣区域识别模块,用于识别超声图像中的病灶,以得到对应的感兴趣区域;病灶分割模块,用于通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓;病灶描述信息生成模块,用于根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。
又一方面,本申请实施例还提供了一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现如上所述方法的步骤。
再一方面,本申请实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如上所述方法的步骤。
上述超声图像病灶描述方法、装置、计算机设备和存储介质,通过对超声图像中的病灶进行识别,以定位对应的感兴趣区域,进而通过图像分割算法对感兴趣区域进行检测,以得到病灶轮廓,并根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征自动生成病灶描述信息,从而有效的减少了医生的工作量,且提高了诊断效率。
附图说明
图1为一个实施例中超声图像病灶描述方法的应用环境图;
图2为一个实施例中超声图像病灶描述方法的流程示意图;
图3为一个实施例中超声图像示意图;
图4为对图3进行目标检测得到的感兴趣区域的示意图;
图5为对图4采用图像分割算法进行分割后得到的病灶轮廓示意图;
图6为在原图上显示病灶轮廓的示意图;
图7为一个实施例中通过图像分割算法得到病灶轮廓步骤的流程示意图;
图8A为另一个实施例中超声图像示意图;
图8B为对图8A进行目标检测得到的感兴趣区域的示意图;
图8C为对图8B进行阈值分割后得到的初始二值图像的示意图;
图9为对图8C进行形态学开运算后得到的新的二值图像的示意图;
图10为对图9进行分析后粗略估计出病灶区域的示意图;
图11为创建的与图8B尺寸相同的二值图像的示意图;
图12为基于距离正则化分割后得到的病灶轮廓示意图;
图13为另一个实施例中通过图像分割算法得到病灶轮廓步骤的流程示意图;
图14A为又一个实施例中超声图像示意图;
图14B为对图14A进行目标检测得到的感兴趣区域的示意图;
图15为初始零水平集函数对应的二值化图像的示意图;
图16为像素点x及其邻域像素y的示意图;
图17为目标水平集函数对应的目标零水平集二值图像的示意图;
图18为对图17进行处理后得到的新的目标零水平集二值图像的示意图;
图19为对图18进行分析后确定的病灶区域的示意图;
图20A为低回声灰度直方图的示意图;
图20B为无回声灰度直方图的示意图;
图21为映射的复数平面示意图;
图22A为在图6中拟合椭圆的示意图;
图22B为图22A中椭圆长轴与水平坐标轴的夹角示意图;
图23为一个实施例中描述病灶是否钙化的步骤的流程示意图;
图24A为再一个实施例中超声图像示意图;
图24B为对图24A进行目标检测得到的感兴趣区域的示意图;
图24C为采用LBF病灶分割算法对图24B分割后得到的有关病灶轮廓的二值图像示意图;
图24D是基于图24B和图24C得到的示意图;
图24E是对图24D做形态学闭运算后得到的示意图;
图24F是对图24E进行种子生长后得到的二值图像的示意图;
图25为一个实施例中超声图像病灶描述装置的结构框图;
图26为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
本申请提供的超声图像病灶描述方法,可以应用于如图1所示的应用环境中。其中,终端102与服务器104通过网络进行通信,在本实施例中,终端102可以是具有超声图像采集功能的设备,也可以是对采集的超声图像进行存储的设备,服务器104可以用独立的服务器或者是多个服务器组成的服务器集群来实现。具体的,终端102用于采集或存储超声图像,并将采集或存储的超声图像通过网络发送至服务器104,服务器104则对超声图像中的病灶进行识别,以定位对应的感兴趣区域,进而通过图像分割算法对感兴趣区域进行检测,以得到病灶轮廓,并根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征自动生成病灶描述信息,从而有效的减少了医生的工作量,且提高了诊断效率。
在一个实施例中,如图2所示,提供了一种超声图像病灶描述方法,以该方法应用于图1中的服务器为例进行说明,包括以下步骤:
步骤202,识别超声图像中的病灶,以得到对应的感兴趣区域。
其中,超声图像为临床上有关乳腺、甲状腺、肝脏、肾脏、脾脏等的超声图像。感兴趣区域(region of interest,简称ROI)是从超声图像中选择的一个需要处理的图像区域,这个区域是进行图像分析所关注的重点。具体的,可以通过目标检测对超声图像中的病灶进行自动识别、定位,从而得到超声图像中对应的感兴趣区域。其中,目标检测可以是基于图像识别算法的目标自动检测,也可以是基于人工框选的目标手动检测。如图3所示,为采集的超声图像示意图,图4则为对图3进行目标检测得到的感兴趣区域的示意图。
步骤204,通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓。
其中,图像分割算法包括基于距离正则化水平集演化(Distance RegularizedLevel Set Evolution,简称DRLSE)的病灶分割算法以及基于局部二值拟合(Local BinaryFitting,简称LBF)演化的病灶分割算法。具体的,通过上述图像分割算法对如图4所示的感兴趣区域进行检测,从而得到对应的分割结果,即如图5所示所示的二值图像,该二值图像的白色区域就是分割出的病灶区域,即图像分割的结果。图3的医学图像术语是蒙板(mask),其表示了被分割目标物的边界轮廓,因此分割结果就是一系列图像点的坐标构成的集合。将这些点的坐标在原图像(即如图3所示的原始超声图像)中呈现,即得到病灶轮廓,如图6所示,图6中的矩形框表示ROI(即感兴趣区域),不规则曲线框则表示病灶轮廓,不规则曲线框内部则为病灶区域。
步骤206,根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。
其中,病灶描述信息包括病灶的生长方向、回声类型、是否存在钙化现象以及病灶轮廓的形状、光整度、清晰度等。在本实施例中,基于上述步骤得到的病灶轮廓以及感兴趣区域,并根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征而自动生成病灶描述信息。
上述超声图像病灶描述方法,通过对超声图像中的病灶进行识别,以定位对应的感兴趣区域,进而通过图像分割算法对感兴趣区域进行检测,以得到病灶轮廓,并根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征自动生成病灶描述信息,从而有效的减少了医生的工作量,且提高了诊断效率。
在一个实施例中,如图7所示,以图像分割算法为基于距离正则化水平集演化的病灶分割算法为例,说明对感兴趣区域进行检测,得到病灶轮廓的具体实现方式,包括:
步骤702,对感兴趣区域进行阈值分割,以得到对应的初始二值图像,并根据初始二值图像估计病灶区域的质心。
具体的,对原超声图像(如图8A)进行目标检测,得到对应的感兴趣区域(如图8B),基于大津法(OTSU)对感兴趣区域进行阈值分割,从而得到对应的初始二值图像(如图8C所示),采用设定掩模尺寸对初始二值图像进行形态学开运算,得到运算后的新的二值图像(如图9所示)。根据邻域连通准则采用8邻域在新的二值图像中提取各邻域的连通分量,此时图9中有8个连通分量,即8个彼此独立的轮廓(每个白色区域为一个轮廓)。进而计算各连通分量所在区域的面积(即该连通分量区域所占的像素数)。将各连通分量所在区域的面积中面积最大的连通分量所在区域粗略估计为病灶区域(如图10所示)。并采用图像矩算法计算病灶区域的质心,其中,图像矩算法具体可以采用Hu矩算法,即通过Hu矩算法得到病灶区域的质心坐标O(X0,Y0)。
步骤704,根据感兴趣区域以及病灶区域的质心,构造初始零水平集函数。
具体的,根据ROI的尺寸创建一幅与ROI尺寸相同的初始图像,进而基于病灶区域的质心以及感兴趣区域的最小边长,在初始图像中确定水平集初始区域。即在初始图像中构造一个像素灰度值为-2(显示为黑色)的圆形区域,该圆形区域以外的像素灰度值为2(显示为白色),且圆形区域的圆心坐标为病灶区域的质心,半径为ROI较小一边的边长除以5,该圆形区域的圆形边界即水平集演化的初始轮廓,也即为水平集初始区域。
由于上述初始图像中水平集初始区域的像素灰度值为-2,而其他区域的像素灰度值为2,因此,得到初始图像的二值化图像,即图11所示的图像,该二值化图像即为初始零水平集图像,也是演化算法的初始对象。又由于数字图像可以通过二元函数表示,因此,对于图11所示的二值化图像可通过如下函数表示:
步骤706,利用距离正则化水平集演化算法基于感兴趣区域对初始零水平集函数进行迭代计算,得到演化终止时刻的目标水平集函数。
具体的,在DRLSE算法中,定义图像信息的能量泛函数为:
E(φ)=μRp(φ)+Eext(φ) (1.2),式中μ>0,为常数,Eext(φ)为外部能量泛函,使得零水平集朝目标边界演化,Rp(φ)为目标水平集函数的正则化项。
具体的,Eext(φ)=λLg(φ)+αAg(φ)=λ∫Ωgδε(φ)|▽φ|dxdy+α∫ΩgHε(-φ)dxdy (1.3);
式中,α是任意实数,即常数,λ为正实数,分别是式(1.3)中右侧长度项和面积项的权重,本实施例中可以取λ=4、α=3。▽φ是对目标水平集函数φ(x,y)求梯度,x,y为图像的横纵坐标;δε(·)和Hε(·)是一维正则化的Dirac函数和Heaviside函数。则有:
g是边缘停止函数,定义其中,Gσ表示标准差为σ的函数,即如图8B所示的感兴趣区域对应的函数,I表示如图11所示的二值图像,*表示卷积运算符。在本实施例中,则可以将DRLSE算法中的边缘停止函数改写为:
公式(1.2)中,定义Rp(φ)=∫Ωp(|▽φ|)dx (1.6);其中,p为势能,定义如下:
将公式(1.3)、(1.6)代入公式(1.2)中,即可得到下式:
E(φ)=μ·∫Ωp(|▽φ|)dx+λ·∫Ωgδε(φ)|▽φ|dxdy+α·∫ΩgHε(-φ)dxdy (1.8),该式是要求解的能量泛函,很显然,能量泛函E(φ)的自变量是φ,而φ是个函数,因此,对式(1.8)等号两边同时做微分,则得到如下式(1.9):
通过梯度下降法求解公式(1.9),以达到求解能量泛函的最小值的目的。将公式(1.9)中的偏微分方程近似转化为离散的有限差分形式,则得到水平集演化方程,即DRLSE模型:
其中,式(1.10)中的就是式(1.9)中的等式的右侧表达式,k为迭代次数iter,Δt为步长step。具体的,步长和迭代次数可根据病灶的回声类型而确定,在本实施例中,当确定回声类型为无回声时,则对应的迭代次数为120~260,时间步长为1.0;若回声类型为低回声时,则对应的迭代次数为650~950,时间步长为1.5。将设定的步长Δt和迭代次数k代入式(1.10)进行计算,从而得到演化终止时刻的目标水平集函数φ(x,y)。
步骤708,基于目标水平集函数得到对应的目标零水平集二值图像。
由于一幅数字图像可以理解为一个二元函数,则一个二元函数在笛卡尔三维坐标系中可以表示为一个连续的曲面。因此,在本实施例中,基于目标水平集函数可以得到对应的目标零水平集二值图像(如图12所示),其中,目标零水平集二值图像中前景与背景的边界像素则为病灶的轮廓数据。
在一个实施例中,如图13所示,以图像分割算法为基于局部二值拟合演化的病灶分割算法为例,说明对感兴趣区域进行检测,得到病灶轮廓的具体实现方式,包括:
步骤1302,根据感兴趣区域,构造初始零水平集函数。
具体的,对原超声图像(如图14A)进行目标检测,得到对应的感兴趣区域(如图14B),由于一幅数字图像可以理解为一个二元函数,在本实施例中,基于感兴趣区域的各顶点坐标,得到病灶区域的初始轮廓,如将感兴趣区域的各顶点坐标确定为参考点,根据设定的平移量对参考点进行平移,以得到相对于各参考点的新的顶点坐标,根据新的顶点坐标得到病灶区域的初始轮廓。
进而由感兴趣区域以及病灶区域的初始轮廓确定对应的初始零水平集函数,在本实施例中,初始零水平集函数用于表示病灶区域的初始轮廓,其也是演化算法的初始对象。其中,水平集方法的主要思想是将曲线作为零水平集嵌入到更高一维的曲面上,通过曲面的演化方程得到函数的演化方程。具体的,根据感兴趣区域(ROI)以及病灶区域的初始轮廓创建一幅与ROI尺寸相同的初始图像(如图15所示),其中,图15所示的黑色区域表示的是初始轮廓的内部,即为乳腺肿瘤病灶的初始区域,设定其像素值均为-2;白色区域表示的是初始轮廓的外部,设定其像素值均为2。
由于上述初始图像中水平集初始区域的像素灰度值为-2,而其他区域的像素灰度值为2,因此,该初始图像为二值化图像,该二值化图像即为初始零水平集图像,也是演化算法的初始对象。又由于数字图像可以通过二元函数表示,因此,对于图15所示的二值化图像可通过如下函数(即初始零水平集函数)表示:
式中r为图像中任意一个像素点的行坐标,c为列坐标,R0表示病灶的初始区域。
步骤1304,基于初始零水平集函数利用局部二值拟合演化算法定义能量泛函。
具体的,假设x是原图像中的任意一个像素点,y则是与像素点x相邻的任意一个像素点(称为x的邻域像素),如图16所示,其中,x和y均为二维向量,可表示为x(c,r),y(c,r)。则定义能量泛函为:
F(φ,f1,f2)=ELBF(φ,f1,f2)+μP(φ)+υL(φ) (2.2)
其中,等式右侧的第一项是能量泛函的主体项,第二项中的P是惩罚项,第三项中的L是水平集函数的零水平曲线的长度,μ,υ则是正常数。
步骤1306,通过梯度下降法求解能量泛函的最小值,以得到演化终止时刻的目标水平集函数。
在本实施例中,通过梯度下降法使得能量泛函最小化,得到水平集活动轮廓演化方程。
具体的,上式(2.2)中:
L(φ)=∫Ωδ(φ(x))|▽φ(x)|dx (2.5)
其中,式(2.3)中的H是Heaviside函数,本文中用正则化的Heaviside函数Hε,近似Heaviside函数,Kσ(x)是标准差为σ的高斯核函数(该函数基于对感兴趣区域进行高斯处理后获得)。I(y)表示如图15所示的二值化图像中任意像素点x的邻域像素y的像素灰度值,λ1、λ2是正的常数,为对应积分项的权重,在本实施例中,λ1=1(即恒为1),λ2则可以根据感兴趣区域的回声类型确定。式(2.4)中▽φ是对零水平集函数φ(c,r)求梯度。式(2.5)中的Dirac函数δ是Heaviside函数的一阶导数。正则化的Dirac函数表示为δε。则有:
将方程式(2.2)右侧的第一项和第三项正则化,则可以近似表达为:
Fε(φ,f1,f2)=Eε LBF(φ,f1,f2)+μP(φ)+υLε(φ) (2.9)
式(9)中f1(x)和f2(x)恒大于零,其中:
通过梯度下降法计算能量泛函的最小值,具体的,保持f1和f2固定不变,使用标准的梯度下降法将关于φ的能量泛函Fε(φ,f1,f2)最小化,从而得到水平集活动轮廓演化方程:
式(2.12)中,
e1(x)=∫ΩKσ(y-x)|I(x)-f1(y)|2dy (2.13)
e2(x)=∫ΩKσ(y-x)|I(x)-f2(y)|2dy (2.14)
基于水平集活动轮廓演化方程采用设定的步长和迭代次数迭代计算,以得到演化终止时刻的目标水平集函数。具体的,将式(2.12)中的偏微分方程近似转化为离散的有限差分形式:
式中的就是式(2.12)中的等式的右侧表达式,采用设定的步长Δt和迭代次数k对式(2.15)进行迭代计算,以得到演化终止时刻的目标水平集函数,即得到当能量泛函Fε(φ,f1,f2)取得最小值时对应的水平集函数φ的零水平集轮廓即为最终结果。其中,步长和迭代次数可根据病灶的回声类型而确定,在本实施例中,当确定回声类型为无回声时,则对应的迭代次数为80~260,时间步长为0.1~1.0,且对应式(2.3)中的λ2取值为2.0~3.3,对应式(2.1)中的V取值为0.003*255*255~0.008*255*255;若回声类型为低回声时,则对应的迭代次数为280~320,时间步长为0.1~1.0,且对应式(2.3)中的λ2取值为1.5~2.2,对应式(2.1)中的V取值为10~8。基于回声类型代入对应的参数进行计算,从而得到演化终止时刻的目标水平集函数,该函数对应如图17所示的目标零水平集二值图像。
步骤1308,基于目标水平集函数得到对应的目标零水平集二值图像,对目标零水平集二值图像进行后处理,以得到病灶区域的目标轮廓。
具体的,对如图17所示的目标零水平集二值图像进行反色处理,以得到反色处理后的多个待筛选的前景区域,其中,反色处理后的前景区域的像素灰度值为255(即白色),背景像素灰度值为0(即黑色)。基于上述反色处理后的多个待筛选的前景区域,对多个待筛选的前景区域进行孔洞填充,从而得到填充后的新的目标零水平集二值图像,如图18所示。
进而根据邻域连通准则采用8邻域在新的目标零水平集二值图像中提取各邻域的连通分量,此时图18中有4个连通分量,很显然,图18中的各白色区域与其它白色区域没有连通,即4个彼此独立的白色区域。进而计算各连通分量所在区域的面积(即该连通分量区域所占的像素数)。将各连通分量所在区域的面积中面积最大的连通分量所在区域确定为病灶区域,如图19所示,则病灶区域的边界像素为对应的目标轮廓,即图19中前景与背景的边界像素为病灶区域的目标轮廓。
在一个实施例中,病灶描述信息包括病灶的回声类型,通常,回声类型包括低回声和无回声,如图20A和图20B所示,分别表示低回声和无回声的灰度直方图。其横轴表示0到255共256个区间(即灰度可能的取值),纵轴为每个灰度值出现的频数。很显然,无回声的直方图中出现频数最多的灰度值的频数与其它灰度值出现的平均频数相差很大。因此可用直方图的这种分布特征来区分病灶的回声类型。具体的,可以采用如下公式来定量表示:ratio=maxFrequency/mean_num,其中,maxFrequency是直方图中出现频数最多的灰度值的频数,mean_num是其它灰度值出现的平均频数。
因此,在本实施例中,根据感兴趣区域得到对应的直方图,进而采用上述公式计算直方图中出现频数最多的灰度值的频数与其他灰度值的平均频数的比值,其中,其他灰度值是指直方图中除出现频数最多的灰度值以外的灰度值,基于比值的大小生成病灶的回声类型,比值越大,则越可能是无回声病灶。具体的,当比值大于7时,则可以确定病灶为无回声病灶,当比值小于7时,则可以确定病灶为低回声病灶。
在一个实施例中,病灶描述信息还包括病灶的形状,通常包括椭圆形、类圆形以及不规则形。则根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息,包括:
采用Fitzgibbon提出的最小二乘拟合椭圆算法将病灶轮廓点集拟合为一个椭圆(如图22A),比较病灶轮廓与该椭圆的相似性,进而基于两者的相似性生成病灶的形状描述,如生成病灶为椭圆形、类圆形或不规则形的描述。
具体的,分别计算病灶轮廓的圆形度、凸包度、紧致度,其计算公式如下:
紧致度公式:式中L是病灶轮廓的周长,S是病灶轮廓所包围的面积。紧致度越接近0,则形状越规则。综合考虑圆形度、凸包度、紧致度,则得到如下公式:式中,0<metric<1。在本实施例中,称metric为几何测度,其值越小,病灶轮廓越规则。具体的,可以设定几何测度阈值为0.4,即如果metric小于0.4,则认为病灶轮廓相对规则,否则认为病灶轮廓不规则。
进而计算病灶轮廓的傅里叶描绘子(Fourier Descriptors),即将图像平面坐标系中的轮廓点集映射到复数平面中。如图21所示的,设病灶轮廓有K个点,以逆时针在该轮廓上行进时,轮廓上的像素点可以看做是坐标对(x0,y0),(x1,y1),(x2,y2),...,(xK-1,yK-1)。这些点的坐标可以表示为x(k)=xk,y(k)=yk的形式,则病灶轮廓可以表示为坐标序列s(k)=[x(k),y(k)],k=0,1,...,K-1。每个坐标对可以当做一个复数来处理:
s(k)=x(k)+jy(k),k=0,1,2,...,K-1
则S(k)的离散傅里叶变换为:
其中,a(u)就是傅里叶描绘子,为了保持轮廓的平移、旋转、缩放的不变性,将其归一化为:
其中,anorm(u)为轮廓上每个点的傅里叶描绘子,其也可以由一个复数表示:
anorm(k)=x(k)+jy(k),k=0,1,2,...,K-1
先分别计算病灶轮廓与椭圆轮廓的傅里叶描绘子(这两个轮廓点的数量一般是不同的,所以需要分别在2个轮廓点集中做均匀采样,即每2个采样点的间隔相等,使得在2个轮廓中各采集K个点,然后采用欧氏距离计算两个轮廓的距离d:
式中,x、y即为复数anorm的实部虚部,距离d可以用来当做判断病灶轮廓和椭圆的相似度,若设定相似度阈值为0.4,则具体判断步骤为:如果d小于等于阈值,则病灶轮廓是“规则形”;如果d大于阈值,进一步用几何测度判断;如果几何测试metric<=0.4,则病灶轮廓是“规则形”,如果metric>0.4,则病灶轮廓是“不规则形”。
若判断为不规则形,则可以直接生成不规则形的形状描述,若判断为“规则形”,则进一步判断具体形状(如椭圆形、类圆形都属于规则形)。具体的,计算上述拟合椭圆的长轴与短轴的比值ratio=w/h,若ratio>=1.5,则确定该病灶的形状为“椭圆形”,若ratio<1.5,则确定该病灶的形状为“类圆形”。
在一个实施例中,病灶描述信息包括病灶的生长方向,通常,生长方向包括平行和非平行。则根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息,包括:采用Fitzgibbon提出的最小二乘拟合椭圆算法将病灶轮廓点集拟合为一个椭圆(如图22A),计算椭圆长轴与图像水平坐标轴的夹角θ,其中,图22B中直线1表示图像坐标系的水平坐标轴的方向,直线2表示椭圆的长轴的方向,直线1与直线2的夹角即θ,根据该夹角θ的大小生成病灶的生长方向。具体的,当00≤θ≤200时,则病灶的生长方向为“平行”,当θ>200时,则病灶的生长方向“非平行”。
在一个实施例中,病灶描述信息包括病灶轮廓的光整度,通常,光整度包括光整和不光整。则根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息,包括:由轮廓点的夹角均值和边界粗糙度来衡量病灶轮廓的光整度。
具体的,令病灶轮廓的轮廓点集为n个像素点,则轮廓点的夹角均值的计算方法:计算病灶轮廓上每个像素点及其相邻像素点的夹角之和,再除以轮廓上像素点的个数,可通过如下公式进行计算:
式中θi,j为像素点i与像素点j之间的夹角。
由于轮廓的质心到轮廓上某一点的距离称作该点的径向长度d(i),则在下式(2)中,令质心坐标为(xc,yc),轮廓上某一点的坐标是(xi,yi)。分母max d(i)表示径向长度最大的值,则归一化的径向长度是dnorm(i),则有:
则边界粗糙度的计算方法为:顺时针依次计算病灶轮廓上相邻像素点的归一化径向长度之差,然后对所有差值求和,再取平均数,具体可通过如下公式进行计算:
(1)和(3)的值越小,则轮廓越光整。利用轮廓点的夹角均值和边界粗糙度的加权值作为判断轮廓光整度的度量,可采用如下公式进行计算:
val=valθ+valrough (4)
基于公式(4)的计算结果生成病灶轮廓的光整度,具体的,可以设定一个阈值与式(4)的值val进行比较,如果val小于该阈值,则确定病灶轮廓为“光整”,反之则确定病灶轮廓为“不光整”,其中,阈值可以设定为3.6,则val小于3.6时确定病灶轮廓光整。
在一个实施例中,病灶描述信息还包括病灶轮廓的清晰度,具体包括清晰和不清晰。则根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息,包括:根据病灶轮廓的边缘点锐度值判断病灶轮廓的清晰度。
具体的,针对病灶轮廓的像素点依次选取8个邻域点求差值,根据8个点的距离设置权重并计算权重和,例如45、135度权重设置为最后对轮廓所有像素点的权重和计算平均值,即得到病灶轮廓的边缘点锐度值,如下公式(5)所示:
其中,m、n表示图像的宽、高,df表示像素点的灰度值变化,dx表示像素点之间的距离,df/dx是灰度变化的幅度。
进而基于病灶轮廓的边缘点锐度值的大小生成病灶轮廓的清晰度,由于边缘点锐度值越小,边缘越模糊。因此,可以设定一个阈值,将上述计算的边缘点锐度值与该阈值比较,若边缘点锐度值小于阈值,则确定病灶轮廓为“不清晰”,反之则确定病灶轮廓为“清晰”。其中,该阈值可以设置为130,则若边缘点锐度值小于130,则表示病灶轮廓不清晰,若边缘点锐度值大于130,则表示病灶轮廓清晰。
在一个实施例中,如图23所示,病灶描述信息还包括病灶是否钙化,则根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息,包括如下步骤:
步骤2302,根据病灶轮廓对感兴趣区域进行图像处理,得到处理后的感兴趣区域。
具体的,设原超声图像为图24A,则图24B是对图24A进行目标检测后得到的ROI图像,图24C是采用上述LBF病灶分割算法分割后得到的有关病灶轮廓的二值图像。将图24C的白色区域作为模板(即病灶轮廓内部区域),在图24B中选取与图24C中白色区域坐标相对应的像素,其余坐标的像素灰度值均设置为0(即黑色),如图24D所示。目的是减少ROI图像中病灶区域以外的干扰,并且减少计算量,提升运算速度。对图24D做形态学闭运算,得到如图24E(即处理后的感兴趣区域),此时,病灶区域的像素灰度变得均匀,从而有利于后续的处理。
步骤2304,遍历处理后的感兴趣区域中所有像素的灰度值,确定灰度值最大的像素的坐标及对应的灰度值。
具体的,对于图24E,首先遍历图中所有灰度值非零的像素,并统计选取其灰度值最大的像素,包括该像素的灰度值valmax及其像素坐标p(x0,y0)。
步骤2306,当灰度值大于设定阈值时,则基于灰度值最大的像素的坐标并采用区域生长算法对处理后的感兴趣区域进行图像分割,以得到对应的二值图像。
具体的,判断valmax是否大于阈值110。需要说明的,本实施例中的阈值110是经验值,由于一般钙化区域的像素灰度较高,本实施例以8位图像的灰度范围为0~255而选取。在本实施例中,如果valmax<110,则表示病灶“无钙化”,如果valmax≥110,则表示病灶可能存在钙化,此时需要进一步进行处理。即以上述灰度值最大的像素的坐标p(x0,y0)为种子点,采用种子生长(又称区域生长)算法进行图像分割,从而得到对应的二值图像,如图24F。
步骤2308,计算二值图像中疑似钙化区域的像素数,基于疑似钙化区域的像素数确定病灶是否钙化。
计算二值图像如图24F中白色区域(即疑似钙化区域)的面积S(像素个数),并基于疑似钙化区域的像素数确定病灶是否钙化。具体的,如果S>=9,则该区域是钙化区域,从而确定病灶“钙化”;如果S<9,则该区域不是钙化区域,从而确定病灶“无钙化”,对于此种情况,可能是由于图像中的高灰度的噪音像素点作为种子点而生成的区域。
应该理解的是,虽然图1-24F的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1-24F中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中,如图25所示,提供了一种超声图像病灶描述装置,包括:感兴趣区域识别模块2501、病灶分割模块2502和病灶描述信息生成模块2503,其中:
感兴趣区域识别模块2501,用于识别超声图像中的病灶,以得到对应的感兴趣区域;
病灶分割模块2502,用于通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓;
病灶描述信息生成模块2503,用于根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。
在一个实施例中,病灶描述信息包括病灶的生长方向;则病灶描述信息生成模块2503具体用于:采用最小二乘拟合椭圆算法将病灶轮廓的轮廓点集拟合为椭圆;计算椭圆的长轴与水平坐标轴的夹角,根据夹角的大小生成病灶的生长方向。
在一个实施例中,病灶描述信息包括病灶的回声类型;则病灶描述信息生成模块2503具体用于:根据感兴趣区域得到对应的直方图;计算直方图中出现频数最多的灰度值的频数与其他灰度值的平均频数的比值,其他灰度值是指直方图中除出现频数最多的灰度值以外的灰度值;基于比值的大小生成病灶的回声类型。
在一个实施例中,病灶描述信息包括病灶轮廓的光整度;则病灶描述信息生成模块2503具体用于:基于病灶轮廓的轮廓点集分别计算轮廓点的夹角均值和边界粗糙度;根据轮廓点的夹角均值和边界粗糙度的加权值大小生成病灶轮廓的光整度。
在一个实施例中,病灶描述信息包括病灶轮廓的清晰度;则病灶描述信息生成模块2503具体用于:基于病灶轮廓的轮廓点集分别计算各轮廓点相对于邻域点的权重和;根据各轮廓点相对于邻域点的权重和计算病灶轮廓的边缘点锐度值;基于病灶轮廓的边缘点锐度值的大小生成病灶轮廓的清晰度。
在一个实施例中,病灶描述信息包括病灶是否钙化;则病灶描述信息生成模块2503具体用于:根据病灶轮廓对感兴趣区域进行图像处理,得到处理后的感兴趣区域;遍历处理后的感兴趣区域中所有像素的灰度值,确定灰度值最大的像素的坐标及对应的灰度值;当灰度值大于设定阈值时,则基于灰度值最大的像素的坐标并采用区域生长算法对处理后的感兴趣区域进行图像分割,以得到对应的二值图像;计算二值图像中疑似钙化区域的像素数,基于疑似钙化区域的像素数确定病灶是否钙化。
在一个实施例中,病灶描述信息还包括病灶的形状;则病灶描述信息生成模块2503具体用于:采用最小二乘拟合椭圆算法将病灶轮廓的轮廓点集拟合为椭圆;基于病灶轮廓与椭圆的相似性生成病灶的形状描述。
在一个实施例中,图像分割算法包括基于距离正则化水平集演化的病灶分割算法;则病灶分割模块2502具体用于:对感兴趣区域进行阈值分割,以得到对应的初始二值图像,并根据初始二值图像估计病灶区域的质心;根据感兴趣区域以及病灶区域的质心,构造初始零水平集函数;利用距离正则化水平集演化算法基于感兴趣区域对初始零水平集函数进行迭代计算,得到演化终止时刻的目标水平集函数;基于目标水平集函数得到对应的目标零水平集二值图像,其中,目标零水平集二值图像中前景与背景的边界像素为病灶轮廓。
在一个实施例中,图像分割算法包括基于局部二值拟合演化的病灶分割算法;则病灶分割模块2502具体用于:根据感兴趣区域,构造初始零水平集函数,其中,初始零水平集函数表示病灶区域的初始轮廓;基于初始零水平集函数利用局部二值拟合演化算法定义能量泛函;通过梯度下降法求解能量泛函的最小值,以得到演化终止时刻的目标水平集函数;基于目标水平集函数得到对应的目标零水平集二值图像,对目标零水平集二值图像进行后处理,以得到病灶区域的目标轮廓。
关于超声图像病灶描述装置的具体限定可以参见上文中对于超声图像病灶描述方法的限定,在此不再赘述。上述超声图像病灶描述装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图26所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储超声图像数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种超声图像病灶描述方法。
本领域技术人员可以理解,图26中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
识别超声图像中的病灶,以得到对应的感兴趣区域;
通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓;
根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
识别超声图像中的病灶,以得到对应的感兴趣区域;
通过图像分割算法对感兴趣区域进行检测,得到病灶轮廓;
根据病灶轮廓的几何拓扑结构以及感兴趣区域的灰度特征生成病灶描述信息。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (12)
1.一种超声图像病灶描述方法,其特征在于,所述方法包括:
识别超声图像中的病灶,以得到对应的感兴趣区域;
通过图像分割算法对所述感兴趣区域进行检测,得到病灶轮廓;
根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息。
2.根据权利要求1所述的方法,其特征在于,所述病灶描述信息包括病灶的生长方向;则根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息,包括:
采用最小二乘拟合椭圆算法将所述病灶轮廓的轮廓点集拟合为椭圆;
计算所述椭圆的长轴与水平坐标轴的夹角,根据所述夹角的大小生成所述病灶的生长方向。
3.根据权利要求1所述的方法,其特征在于,所述病灶描述信息包括病灶的形状;则根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息,包括:
采用最小二乘拟合椭圆算法将所述病灶轮廓的轮廓点集拟合为椭圆;
基于所述病灶轮廓与所述椭圆的相似性生成所述病灶的形状描述。
4.根据权利要求1所述的方法,其特征在于,所述病灶描述信息包括病灶的回声类型;则根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息,包括:
根据所述感兴趣区域得到对应的直方图;
计算所述直方图中出现频数最多的灰度值的频数与其他灰度值的平均频数的比值,所述其他灰度值是指所述直方图中除出现频数最多的灰度值以外的灰度值;
基于所述比值的大小生成所述病灶的回声类型。
5.根据权利要求1所述的方法,其特征在于,所述病灶描述信息包括病灶轮廓的光整度;则根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息,包括:
基于所述病灶轮廓的轮廓点集分别计算轮廓点的夹角均值和边界粗糙度;
根据所述轮廓点的夹角均值和边界粗糙度的加权值大小生成所述病灶轮廓的光整度。
6.根据权利要求1所述的方法,其特征在于,所述病灶描述信息包括病灶轮廓的清晰度;则根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息,包括:
基于所述病灶轮廓的轮廓点集分别计算各轮廓点相对于邻域点的权重和;
根据所述各轮廓点相对于邻域点的权重和计算所述病灶轮廓的边缘点锐度值;
基于所述病灶轮廓的边缘点锐度值的大小生成所述病灶轮廓的清晰度。
7.根据权利要求1所述的方法,其特征在于,所述病灶描述信息包括病灶是否钙化;则根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息,包括:
根据所述病灶轮廓对所述感兴趣区域进行图像处理,得到处理后的感兴趣区域;
遍历所述处理后的感兴趣区域中所有像素的灰度值,确定灰度值最大的所述像素的坐标及对应的灰度值;
当所述灰度值大于设定阈值时,则基于灰度值最大的所述像素的坐标并采用区域生长算法对所述处理后的感兴趣区域进行图像分割,以得到对应的二值图像;
计算所述二值图像中疑似钙化区域的像素数,基于所述疑似钙化区域的像素数确定所述病灶是否钙化。
8.根据权利要求1至7任一项所述的方法,其特征在于,所述图像分割算法包括基于距离正则化水平集演化的病灶分割算法;则通过图像分割算法对所述感兴趣区域进行检测,得到病灶轮廓,包括:
对所述感兴趣区域进行阈值分割,以得到对应的初始二值图像,并根据初始二值图像估计病灶区域的质心;
根据所述感兴趣区域以及所述病灶区域的质心,构造初始零水平集函数;
利用距离正则化水平集演化算法基于所述感兴趣区域对所述初始零水平集函数进行迭代计算,得到演化终止时刻的目标水平集函数;
基于所述目标水平集函数得到对应的目标零水平集二值图像,所述目标零水平集二值图像中前景与背景的边界像素为所述病灶轮廓。
9.根据权利要求1至7任一项所述的方法,其特征在于,所述图像分割算法包括基于局部二值拟合演化的病灶分割算法;则通过图像分割算法对所述感兴趣区域进行检测,得到病灶轮廓,包括:
根据所述感兴趣区域,构造初始零水平集函数,所述初始零水平集函数表示病灶区域的初始轮廓;
基于所述初始零水平集函数利用局部二值拟合演化算法定义能量泛函;
通过梯度下降法求解所述能量泛函的最小值,以得到演化终止时刻的目标水平集函数;
基于所述目标水平集函数得到对应的目标零水平集二值图像,对所述目标零水平集二值图像进行后处理,以得到病灶区域的目标轮廓。
10.一种超声图像病灶描述装置,其特征在于,所述装置包括:
感兴趣区域识别模块,用于识别超声图像中的病灶,以得到对应的感兴趣区域;
病灶分割模块,用于通过图像分割算法对所述感兴趣区域进行检测,得到病灶轮廓;
病灶描述信息生成模块,用于根据所述病灶轮廓的几何拓扑结构以及所述感兴趣区域的灰度特征生成病灶描述信息。
11.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至9中任一项所述方法的步骤。
12.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至9中任一项所述方法的步骤。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911358149.6A CN113034426B (zh) | 2019-12-25 | 2019-12-25 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
PCT/CN2020/133026 WO2021129323A1 (zh) | 2019-12-25 | 2020-12-01 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911358149.6A CN113034426B (zh) | 2019-12-25 | 2019-12-25 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113034426A true CN113034426A (zh) | 2021-06-25 |
CN113034426B CN113034426B (zh) | 2024-03-08 |
Family
ID=76458528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911358149.6A Active CN113034426B (zh) | 2019-12-25 | 2019-12-25 | 超声图像病灶描述方法、装置、计算机设备和存储介质 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113034426B (zh) |
WO (1) | WO2021129323A1 (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113379739A (zh) * | 2021-07-23 | 2021-09-10 | 平安科技(深圳)有限公司 | 超声图像的识别方法、装置、设备及存储介质 |
CN113393462A (zh) * | 2021-08-18 | 2021-09-14 | 浙江杜比医疗科技有限公司 | 一种病灶区域分割方法、装置、电子设备及存储介质 |
CN113421293A (zh) * | 2021-06-30 | 2021-09-21 | 上海申瑞继保电气有限公司 | 变电站设备图像质心计算方法 |
CN113657553A (zh) * | 2021-09-01 | 2021-11-16 | 什维新智医疗科技(上海)有限公司 | 一种结节回声类型判断装置 |
CN113706473A (zh) * | 2021-08-04 | 2021-11-26 | 青岛海信医疗设备股份有限公司 | 超声图像中病灶区域的长短轴确定方法和超声设备 |
CN113781477A (zh) * | 2021-10-29 | 2021-12-10 | 平安科技(深圳)有限公司 | 基于人工智能的结石图像识别方法、装置及设备 |
CN114066781A (zh) * | 2022-01-18 | 2022-02-18 | 浙江鸿禾医疗科技有限责任公司 | 胶囊内窥镜肠道图像的识别定位方法、存储介质和设备 |
CN114757952A (zh) * | 2022-06-15 | 2022-07-15 | 深圳瀚维智能医疗科技有限公司 | 超声图像处理方法、装置、设备及存储介质 |
CN115063607A (zh) * | 2022-08-17 | 2022-09-16 | 腾讯科技(深圳)有限公司 | 切割轮廓的确定方法、装置、设备、存储介质及程序产品 |
CN115311268A (zh) * | 2022-10-10 | 2022-11-08 | 武汉楚精灵医疗科技有限公司 | 食管内窥镜图像的识别方法及装置 |
CN116309551A (zh) * | 2023-05-11 | 2023-06-23 | 浙江太美医疗科技股份有限公司 | 确定病灶取样区的方法、装置、电子设备及可读介质 |
CN116542942A (zh) * | 2023-05-10 | 2023-08-04 | 什维新智医疗科技(上海)有限公司 | 一种甲状腺无回声结节判定方法、系统、设备及介质 |
CN116862906A (zh) * | 2023-08-24 | 2023-10-10 | 武汉大学人民医院(湖北省人民医院) | 一种眼部检测装置及方法 |
CN117541580A (zh) * | 2024-01-08 | 2024-02-09 | 天津市肿瘤医院(天津医科大学肿瘤医院) | 一种基于深度神经网络的甲状腺癌图像比对模型建立方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102068281A (zh) * | 2011-01-20 | 2011-05-25 | 深圳大学 | 一种占位性病变超声图像的处理方法 |
CN103337074A (zh) * | 2013-06-18 | 2013-10-02 | 大连理工大学 | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 |
US20170270664A1 (en) * | 2016-03-21 | 2017-09-21 | The Board Of Trustees Of The Leland Stanford Junior University | Methods for characterizing features of interest in digital images and systems for practicing same |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4868843B2 (ja) * | 2005-01-26 | 2012-02-01 | 株式会社東芝 | 超音波診断装置及び超音波診断装置の制御プログラム |
CN105654490A (zh) * | 2015-12-31 | 2016-06-08 | 中国科学院深圳先进技术研究院 | 一种基于超声弹性图像的病变区域提取方法及装置 |
CN110223289A (zh) * | 2019-06-17 | 2019-09-10 | 上海联影医疗科技有限公司 | 一种图像处理方法和系统 |
-
2019
- 2019-12-25 CN CN201911358149.6A patent/CN113034426B/zh active Active
-
2020
- 2020-12-01 WO PCT/CN2020/133026 patent/WO2021129323A1/zh active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102068281A (zh) * | 2011-01-20 | 2011-05-25 | 深圳大学 | 一种占位性病变超声图像的处理方法 |
CN103337074A (zh) * | 2013-06-18 | 2013-10-02 | 大连理工大学 | 一种基于主动轮廓模型分割乳腺dce-mri病灶的方法 |
US20170270664A1 (en) * | 2016-03-21 | 2017-09-21 | The Board Of Trustees Of The Leland Stanford Junior University | Methods for characterizing features of interest in digital images and systems for practicing same |
Non-Patent Citations (2)
Title |
---|
杨谊 等: "基于改进水平集图像分割方法的乳腺超声病灶提取", 《计算机应用与软件》, vol. 31, no. 11, pages 218 - 221 * |
高东平 等: "乳腺肿瘤超声图像特征参数量化研究进展", 《北京生物医药工程》, vol. 30, no. 6, pages 656 - 659 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113421293B (zh) * | 2021-06-30 | 2023-12-29 | 上海申瑞继保电气有限公司 | 变电站设备图像质心计算方法 |
CN113421293A (zh) * | 2021-06-30 | 2021-09-21 | 上海申瑞继保电气有限公司 | 变电站设备图像质心计算方法 |
CN113379739B (zh) * | 2021-07-23 | 2022-03-25 | 平安科技(深圳)有限公司 | 超声图像的识别方法、装置、设备及存储介质 |
CN113379739A (zh) * | 2021-07-23 | 2021-09-10 | 平安科技(深圳)有限公司 | 超声图像的识别方法、装置、设备及存储介质 |
CN113706473A (zh) * | 2021-08-04 | 2021-11-26 | 青岛海信医疗设备股份有限公司 | 超声图像中病灶区域的长短轴确定方法和超声设备 |
CN113706473B (zh) * | 2021-08-04 | 2024-03-01 | 青岛海信医疗设备股份有限公司 | 超声图像中病灶区域的长短轴确定方法和超声设备 |
CN113393462A (zh) * | 2021-08-18 | 2021-09-14 | 浙江杜比医疗科技有限公司 | 一种病灶区域分割方法、装置、电子设备及存储介质 |
CN113657553A (zh) * | 2021-09-01 | 2021-11-16 | 什维新智医疗科技(上海)有限公司 | 一种结节回声类型判断装置 |
CN113657553B (zh) * | 2021-09-01 | 2023-12-26 | 什维新智医疗科技(上海)有限公司 | 一种结节回声类型判断装置 |
CN113781477A (zh) * | 2021-10-29 | 2021-12-10 | 平安科技(深圳)有限公司 | 基于人工智能的结石图像识别方法、装置及设备 |
WO2023071114A1 (zh) * | 2021-10-29 | 2023-05-04 | 平安科技(深圳)有限公司 | 基于人工智能的结石图像识别方法、装置及设备 |
CN114066781A (zh) * | 2022-01-18 | 2022-02-18 | 浙江鸿禾医疗科技有限责任公司 | 胶囊内窥镜肠道图像的识别定位方法、存储介质和设备 |
CN114066781B (zh) * | 2022-01-18 | 2022-05-10 | 浙江鸿禾医疗科技有限责任公司 | 胶囊内窥镜肠道图像的识别定位方法、存储介质和设备 |
CN114757952A (zh) * | 2022-06-15 | 2022-07-15 | 深圳瀚维智能医疗科技有限公司 | 超声图像处理方法、装置、设备及存储介质 |
CN115063607B (zh) * | 2022-08-17 | 2022-12-30 | 腾讯科技(深圳)有限公司 | 切割轮廓的确定方法、装置、设备及存储介质 |
CN115063607A (zh) * | 2022-08-17 | 2022-09-16 | 腾讯科技(深圳)有限公司 | 切割轮廓的确定方法、装置、设备、存储介质及程序产品 |
CN115311268A (zh) * | 2022-10-10 | 2022-11-08 | 武汉楚精灵医疗科技有限公司 | 食管内窥镜图像的识别方法及装置 |
CN116542942A (zh) * | 2023-05-10 | 2023-08-04 | 什维新智医疗科技(上海)有限公司 | 一种甲状腺无回声结节判定方法、系统、设备及介质 |
CN116309551B (zh) * | 2023-05-11 | 2023-08-15 | 浙江太美医疗科技股份有限公司 | 确定病灶取样区的方法、装置、电子设备及可读介质 |
CN116309551A (zh) * | 2023-05-11 | 2023-06-23 | 浙江太美医疗科技股份有限公司 | 确定病灶取样区的方法、装置、电子设备及可读介质 |
CN116862906A (zh) * | 2023-08-24 | 2023-10-10 | 武汉大学人民医院(湖北省人民医院) | 一种眼部检测装置及方法 |
CN116862906B (zh) * | 2023-08-24 | 2023-12-12 | 武汉大学人民医院(湖北省人民医院) | 一种眼部检测装置及方法 |
CN117541580A (zh) * | 2024-01-08 | 2024-02-09 | 天津市肿瘤医院(天津医科大学肿瘤医院) | 一种基于深度神经网络的甲状腺癌图像比对模型建立方法 |
CN117541580B (zh) * | 2024-01-08 | 2024-03-19 | 天津市肿瘤医院(天津医科大学肿瘤医院) | 一种基于深度神经网络的甲状腺癌图像比对模型建立方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113034426B (zh) | 2024-03-08 |
WO2021129323A1 (zh) | 2021-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113034426B (zh) | 超声图像病灶描述方法、装置、计算机设备和存储介质 | |
CN113112443B (zh) | 超声图像病灶的分割方法、装置和计算机设备 | |
Foruzan et al. | Improved segmentation of low-contrast lesions using sigmoid edge model | |
Villalobos-Castaldi et al. | A fast, efficient and automated method to extract vessels from fundus images | |
JP5864542B2 (ja) | 画像の異常を検出するための画像データの処理方法、システムおよびプログラム | |
US9230320B2 (en) | Computer aided diagnostic system incorporating shape analysis for diagnosing malignant lung nodules | |
WO2006024974A1 (en) | Feature weighted medical object contouring using distance coordinates | |
Luo et al. | A Novel Segmentation Approach Combining Region‐and Edge‐Based Information for Ultrasound Images | |
US8577104B2 (en) | Liver lesion segmentation | |
Liu et al. | An effective approach of lesion segmentation within the breast ultrasound image based on the cellular automata principle | |
Elizabeth et al. | A novel segmentation approach for improving diagnostic accuracy of CAD systems for detecting lung cancer from chest computed tomography images | |
WO2021129324A1 (zh) | 超声图像病灶的分割方法、装置和计算机设备 | |
CN109753997B (zh) | 一种ct图像中的肝脏肿瘤自动精确鲁棒分割方法 | |
CN112308846B (zh) | 血管分割方法、装置及电子设备 | |
CN116071355A (zh) | 一种用于外周血管影像的辅助分割系统及方法 | |
WO2011158135A1 (en) | Tissue classification | |
Iqbal et al. | Potential lung nodules identification for characterization by variable multistep threshold and shape indices from CT images | |
Faisal et al. | Computer assisted diagnostic system in tumor radiography | |
Masood et al. | Automated segmentation of skin lesions: Modified Fuzzy C mean thresholding based level set method | |
Fotin et al. | A multiscale Laplacian of Gaussian (LoG) filtering approach to pulmonary nodule detection from whole-lung CT scans | |
CN114187320B (zh) | 脊柱ct图像的分割方法、椎体成像的识别方法及装置 | |
Rawas et al. | HCET-G 2: dermoscopic skin lesion segmentation via hybrid cross entropy thresholding using Gaussian and gamma distributions | |
Vanmore et al. | Survey on automatic liver segmentation techniques from abdominal CT images | |
US20080107321A1 (en) | Spiculation detection method and apparatus for CAD | |
Kumar et al. | Semiautomatic method for segmenting pedicles in vertebral radiographs |
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 | ||
CB02 | Change of applicant information |
Country or region after: China Address after: 215000 5th floor, building a, 4th floor, building C, 27 Xinfa Road, Suzhou Industrial Park, Suzhou City, Jiangsu Province Applicant after: Feiyinuo Technology Co.,Ltd. Address before: 215000 5th floor, building a, 4th floor, building C, 27 Xinfa Road, Suzhou Industrial Park, Suzhou City, Jiangsu Province Applicant before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd. Country or region before: China |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |