一种基于肝癌免疫组化图像的血管角度自动测量方法技术

本发明专利技术公开了一种基于肝癌免疫组化图像的血管角度自动测量方法，包括：1）采集肝癌免疫组化的血管图像；2）对图像进行染色分离、阈值分割、平滑、填充处理，使血管清晰、完整；3）对目标血管进行特征点提取并将所提取的特征点分为两个等量的集合；4）根据所分的两个集合，拟合两根直线，并对这两根直线进行判断、调整，最后计算出血管角度。本发明专利技术的自动测量肝癌免疫组化图像的血管角度的方法，对免疫组化图像血管角度测量准确率高，对判断恶性实体肿瘤的生长、转移有很重要的意义。

【技术实现步骤摘要】
一种基于肝癌免疫组化图像的血管角度自动测量方法
本专利技术属于病理学
，具体涉及一种基于肝癌免疫组化图像(IHC)的血管角度自动测量方法。
技术介绍
现代医学认为肝癌的发生是一个多因素、多步骤的生物学过程，肝癌细胞的增殖、侵袭需要依赖血管。而血管的生成与肝癌的生长、浸润和复发均有明显的相关性。临床上，当癌肿增大明显时，受侵袭的血管段不规则、压迹、移位，直接影响肝脏毛细血管走向，从而提出了根据血管角度的大小来判别肝癌分化程度的方法。现阶段免疫组化技术的应用提高了临床病理诊断的准确性，已经成为判断肿瘤良恶性及生长程度的重要评价标准。然而目前免疫组化自动分析方法不足以解决血管的全部形态问题，只能获得血管长度、面积，长宽比等常规统计学信息。目前尚无研究人员提出血管角度的自动测量方法，仅靠人工测量角度工作量大、易受个人主观因素的影响。免疫组化是生物成像的核心技术，生物医学图像处理技术的提高促进了免疫组化图像自动分析工具的发展，这对辅助病理学家进行与疾病相关的生物学活性检测具有重要的意义。面对IHC存在的染色剂分布不均、成像有噪声等复杂情况，如何精确地提取血管结构以用于辅助诊断等有着关键性作用，它的处理结果将影响后续定量检测的精度。由于IHC图像使用多种染色剂进行着色，直接采用传统的方法不能达到对象准确分割的目的。比如：有学者采用R、G、B分量上的色差来作为分割依据，但由于R、G、B这三个颜色空间存在很强的相关性，直接在RGB颜色模型下很难实现免疫组化阳性区、阴性区和背景区的颜色分离。一些学者将RGB颜色模型转换为HSI模型或者CMYK模型进行多着色的图像分割，但是这种方法依赖于阈值的设置，容易导致大量信息的丢失，影响检测的精确度。
技术实现思路
专利技术目的：针对现有技术中存在的不足，本专利技术的目的是提供一种基于肝癌免疫组化图像的血管角度自动测量方法，对免疫组化图像血管角度测量准确率高，对判断恶性实体肿瘤的生长、转移有很重要的意义。技术方案：为了实现上上述专利技术目的，本专利技术采用的技术方案为：一种基于肝癌免疫组化图像的血管角度自动测量方法，包括如下步骤：1)采集肝癌免疫组化的血管图像；2)对图像进行染色分离、阈值分割、平滑、填充处理，使血管清晰、完整；3)对目标血管进行特征点提取并将所提取的特征点分为两个等量的集合；4)根据所分的两个集合，拟合两根直线，并对这两根直线进行判断、调整，最后计算出血管角度。步骤2)中，染色分离利用的是颜色去卷积算法(具体参考论文“FastBoundaryArtifactReductionAlgorithmforImageDeconvolution”)，颜色去卷积算法为了获取每种着色的独立信息，需对RGB信息进行正交变换，而且为了获得每种着色正确的吸收因子，需对变换进行标准化。通过临床上最常用的两种染色剂DAB与H&E的实验证明，颜色去卷积的方法可以成功地分割出不同的染色区域；步骤2)中，阈值分割法利用的是改进的Otsu算法，本专利技术在分析传统的Otsu算法原理的基础上，用后验概率信息代替传统的Otsu算法均值信息对于免疫组化图像来说具有更好的自适应性,能够有效地减少图像的灰度线性变化和平移变化对传统的Otsu算法的影响，在复杂背景下仍能获得最佳的分割阈值。本专利技术改进的Otsu算法步骤具体包括所获得的染色分离图像为灰度图像，其灰度直方图产生的灰度信息概率为u(a)(0＜a＜L-1)，同时满足染色分离图像均值当阈值t选择到最佳值时，图像则会被分类成两个部分：血管和基质。血管和基质的后验概率分别为pb(t)和pj(t)。血管和基质的灰度均值可表示为zb(t)和zj(t)。则归一化的Otsu阈值分割准则函数为：血管和基质的后验概率pb(t)和pj(t)对应的Rényi信息熵Rn为其中n为信息熵参数。将获得的信息熵Rn(p,t)作为Otsu阈值准则函数的约束条件，则获得改进的Otsu算法的正则化最佳阈值判别准则函数为：步骤2)中，平滑利用的是中值滤波方法，中值滤波是图像平滑的一种方法，它是一种非线性平滑滤波技术，通常采用一个含奇数个点的滑动窗口，用窗口中的灰度值的中值来代替中心点的灰度值；步骤2)中，填充利用的是四连通域种子填充算法，该算法从多边形区域的一个内点开始，由内向外用给定的颜色画点直到边界为止，如果边界是以一种颜色指定的，则种子填充算法可逐个像素地处理直到遇到边界颜色为止；步骤3)中，对目标血管进行特征点提取采用的是SURF算法(具体参考论文“ResearchonMedicalImageRegistrationAlgorithmBasedonSURFImprovement”)，SURF特征点提取，通过计算Hessian矩阵行列式的局部极值点,来确定特征点的位置。SURF算法利用盒子滤波器代替高斯二阶微分,并引入了积分图像,极大的提高了运算速度和精度。将盒子滤波器与图像进行卷积运算,其结果分别为Dxx、Dxy和Dyy。于是,Hessian矩阵的行列式可以简化为:Det(H)＝DxyDyy-(wDxy)2(式中w是权重系数)，按照式中近似Hessian矩阵行列式的计算方法，对图像中每一点求取响应并记录下来当响应大于阈值时，则标记为特征点；步骤3)中，将所提取的特征点分为两个等量的集合包括如下步骤：(1)对获取的N个特征点进行编号包括如下步骤：①根据特征点的坐标位置关系，采用线扫描的方式从上、下、左、右四个方向进行扫描，分别找出四个方向上经扫描后的第一个点，标记为Tup、Tdown、Tleft、Tright，分别计算出该点与最邻近的两个点的角度，记为∠Tup1、∠Tup2、∠Tdown1、∠Tdown2、∠Tleft1、∠Tleft2、∠Tright1、∠Tright2，再计算出同方向上两个角度之间的差值的绝对值：如∠Tup1-2＝|∠Tup1-∠Tup2|，取其中角度差值最小的点为血管起始点，记为P1；②根据欧几里得距离关系Oρ＝sqrt((x1-x2)2+(y1-y2)2)，分别计算P1与其他N-1个特征点的距离，取到P1距离最小的点记为P2，再分别计算P2到其他N-2个特征点的距离，当最小时，记i为P3，以此类推，直到标记出所有的特征点。(2)N为奇数时，取编号为(N+1)/2的特征点为标记点，N为偶数时，取编号为N/2与N/2+1之间的距离作为标记点，把血管中N个特征点分成等量的两个集合S1与S2，(其中N为奇数时，S1＝{P1，...,P(N-1)/2},S2＝{P(N+1)/2，...,PN},N为偶数时，S1＝{P1，P2，...,PN/2}，S2＝{PN/2+1，...，PN})；步骤4)中，根据所分的两个集合，拟合两根直线包括如下步骤：(1)根据两个集合中的特征点利用最小二乘法拟合两根直线L1、L2(具体参考论文“MultipleAdaptiveLeastSquareCurveFittingAlgorithmandApplications”)。步骤4)中，对所拟合成的两根直线进行判断、调整包括如下步骤：(1)由公式：(公式描述：公式中Cov(X,Y)为X,Y的协方差，D(X)、D(Y)分别为X、Y的方差。)可知，当|ρXY|＝1时，称X,Y完全相关，此时，X,Y本文档来自技高网...

【技术保护点】
一种基于肝癌免疫组化图像的血管角度自动测量方法，其特征在于，包括如下步骤：1)采集肝癌免疫组化的血管图像；2)对图像进行染色分离、阈值分割、平滑、填充处理，使血管清晰、完整；3)对目标血管进行特征点提取并将所提取的特征点分为两个等量的集合；4)根据所分的两个集合，拟合两根直线，并对这两根直线进行判断、调整，最后计算出血管角度。

【技术特征摘要】
1.一种基于肝癌免疫组化图像的血管角度自动测量方法，其特征在于，包括如下步骤：1)采集肝癌免疫组化的血管图像；2)对图像进行染色分离、阈值分割、平滑、填充处理，使血管清晰、完整；3)对目标血管进行特征点提取并将所提取的特征点分为两个等量的集合；4)根据所分的两个集合，拟合两根直线，并对这两根直线进行判断、调整，最后计算出血管角度。2.根据权利要求1所述的基于肝癌免疫组化图像的血管角度自动测量方法，其特征在于，步骤2)中，利用颜色去卷积算法进行染色分离，对RGB信息进行正交变换获取每种着色的独立信息，对变换进行标准化获得每种着色正确的吸收因子。3.根据权利要求1所述的基于肝癌免疫组化图像的血管角度自动测量方法，其特征在于，步骤2)中，阈值分割时，利用图像染色分离后所得到的DAB染色和H染色对象的后验概率为新的正则化约束，实现Otsu最佳阈值判别准则函数的正则化，正则化最佳阈值判别准则函数为：4.根据权利要求1所述的基于肝癌免疫组化图像的血管角度自动测量方法，其特征在于，步骤3)中，采用SURF算法对目标血管进行特征点提取，SURF特征点提取，通过计算Hessian矩阵行列式的局部极值点，来确定特征点的位置；将盒子滤波器与图像进行卷积运算，结果分别为Dxx、Dxy和Dyy；Hessian矩阵的行列式简化为：Det(H)＝DxyDyy-(wDxy)2，式中，w是权重系数；按照式中近似Hessian矩阵行列式的计算方法，对图像中每一点求取响应并记录下来，根据判断阈值判断哪些响应点作为特征点。5.根据权利要求1所述的基于肝癌免疫组化图像的血管角度自动测量方法，其特征在于，步骤3)中，将所提取的特征点分为两个等量的集合包括如下步骤：(1)对获取的N个特征点进行编号，包括如下步骤：①根据特征点的坐标位置关系，采用线扫描的方式从上、下、左、右四个方向进行扫描，分别找出四个方向上经扫描后的第一个点，标记为Tup、Tdown、Tleft、Tright，分别计算出该点与最邻近的两个点的角度，记为∠Tup1、∠Tup2、∠Tdown1、∠Tdown2、∠Tleft1、∠Tleft2、∠Tright1、∠Tright2，再计算出同方向上两个角度之间的差值的绝对值，取其中角度差值最小的点为血管起始点，记为P1；②根据欧几里得距离关系Oρ＝sqrt((x1-x2)2+(...

【专利技术属性】
技术研发人员：张堃，陆逸骋，陈莉，吴建国，张培建，张洪彬，
申请(专利权)人：南通大学，
类型：发明
国别省市：江苏,32