白内障图像的颜色校正方法技术

本发明专利技术涉及内障图像的颜色校正方法，包括：A.设定像素光源估计层数w，每层设定不同聚类种子点的个数k，并将每一层的图像初始化；B.利用简单线性迭代聚类将每一层的图像分别划分为图像块；C.计算出每个图像块的光源估计，将其反投影到原始图像上，得到每一层的像素照明光源估计；D.将所有层的像素照明光源估计结果进行线性融合，得到最终的像素照明光源估计结果，根据像素照明光源估计对原始图像进行恢复，得到颜色校正后的结果图。克服了现有方法中受边缘影响导致白内障图像颜色校正中颜色偏差大和边缘颜色漂移的问题，棉线提高白内障图像颜色校正的准确率，降低了边缘影响。

Color Correction Method for Cataract Images

【技术实现步骤摘要】
白内障图像的颜色校正方法
本专利技术涉及医学图像处理的方法，具体的讲是白内障图像的颜色校正方法。
技术介绍
白内障是一种常见的致盲性眼部疾病，当晶状体中蛋白质积聚时会引发白内障，出现晶状体混浊，这会阻挡光线透过晶状体，导致视力下降。根据白内障出现的区域，可以分为将白内障分为三种类型：核性白内障(N型)，皮质性白内障(C型)和后囊型白内障(P型)。其中N型白内障是最常见的白内障类型，它在晶状体的焦点区域(核)深处形成，通常与细胞衰老有关。C型白内障是由于晶状体皮层(外层)变得不透明引起的，当晶状体周边的水含量改变时穿线裂缝时会发生这种情况。P型白内障发生在晶状体后部，患有糖尿病的人或服用高剂量类固醇药物的人患这种类型的白内障的风险更大。大多数老年人患有这三种类型的白内障中的一种。因此，如果眼科疾病能够被自动系统大规模筛选，对医生和病人来说都意义重大。皮质型白内障图像包括全部角膜、虹膜、睫状体、晶状体等多个眼睛结构，其中白内障病变区域仅出现在晶状体；图像的采集仅由眼科医生完成，故图像的质量受医生习惯、拍摄环境等影响较大，所以图像质量差别较大，一般可以通过颜色恒常性算法将图像校准为标准光照射下采集图像。目前，常用的颜色恒常性算法分为以下几类：(1)统计方法：统计类颜色恒常性方法不使用任何先验信息，仅利用图像的底层特征信息去估计光源，通常使用的图像底层特征信息包括边缘、纹理、对比度信息等等。灰度世界法、白斑法、灰度阴影法、灰度边缘法等均属于统计方法的范畴。尽管这类方法实现简单，但对于单幅输入图像，由于不知道场景中入射光的信息，以及场景中物体表面的反射特性，因此这类方法在没有合理的假设条件下问题是无法有效解决的。(2)监督学习方法：监督学习颜色恒常性算法需要学习一些先验信息来促进问题的解决，如光源的色度分布，或者场景中物体表面的颜色分布，然后用这些学习到的有用的先验信息去估计输入图像中的未知光源的颜色。常用方法有色域映射算法、支持向量回归法、语义分析法等，该类方法准确率很高，但基于监督学习的方法需要知道训练数据的光源信息，且每个学习模型只对特定类型图像有效。(3)基于生物医学的方法：计算机视觉的颜色恒常性模型就是以模拟人眼视觉实现为目的的，深入研究人类视觉系统的恒常性机制将进一步促进计算机视觉颜色恒常性问题的解决。不过目前生物医学对于人眼的颜色恒常性机制还不能非常清晰地从细胞级别进行解译，因此目前研究仅初步掌握人眼获取视觉信息的基本模型。目前应用于图像预处理阶段的颜色恒常性算法，主要为统计方法，但这一类方法均需假设所研究场景中只有单个光源且光照均匀分布，现实中的自然图像很难满足该条件，现实场景中一般都会有多个光源的存在，例如室外场景中阴影区域仅受到天空光线的照亮，而非阴影区域同时受到太阳光和天空光线的照亮；室内场景会同时受室内和室外光源的影响；物体表面间的相互反射也会引入一个“虚光源”。针对多光源的颜色恒常性方法尚处于起步阶段，现有方法有基于物理的方法、基于生物启发模型的方法、基于样本学习的方法；其中基于物理的方法是针对室外拍摄图像的方法；基于生物启发模型的方法，通过模拟人眼的拮抗和双拮抗特性来实现对光源的估计；基于样本学习的方法，需要利用特征从训练集中找到待处理图像中的所有表面最接近的表面，从而估计出照射未知表面的光源值。为了提高白内障图像颜色校正的准确率，考虑到常规颜色恒常性方法需假设场景光照均匀一致的特点，显然，多光源假设类基于样本学习的颜色恒常性方法更符合皮质型白内障图像特点，而该方法需要知道训练集图像光源信息，且算法复杂，容易受到表面边界影响，导致场景边缘颜色漂移；由于无法测量皮质型白内障图像的光源信息，同时为了消除边界效应，需要一种多层融合颜色恒常性的C型白内障图像颜色校正方法。
技术实现思路
本专利技术提供了一种白内障图像的颜色校正方法，以克服现有方法中受边缘影响导致白内障图像颜色校正中颜色偏差大和边缘颜色漂移的问题，并提高白内障图像颜色校正的准确率，降低边缘影响。本专利技术的内障图像的颜色校正方法，包括：A.设定像素光源估计层数w，每层设定不同聚类种子点的个数k，并将每一层的图像初始化；B.利用简单线性迭代聚类(simplelineariterativeclustering，SLIC)将每一层的图像分别划分为图像块；C.计算出每个图像块的光源估计，将其反投影到原始图像上，得到每一层的像素照明光源估计；D.将所有层的像素照明光源估计结果进行线性融合，得到最终的像素照明光源估计结果，根据像素照明光源估计对原始图像进行恢复，得到颜色校正后的结果图。本专利技术采用多光源颜色恒常性假设，将图像划分为多个图像块，假设每个图像块受单光源照射且光照均匀，更加符合皮质型白内障图像光照条件，解决了现有方法中由于采用单光源假设而导致光源估计与真实光源偏差大和采用多光源假设而需要一直光源信息数据集的问题，有效提高了光源估计准确率，降低了算法复杂度。而且引入了专门针对图像块边缘的多层融合，多层融合每一层采用不同图像块划分方法，每层块边缘不同，通过线性融合平滑块边缘，更好的描述了皮质型白内障图像光照信息，使算法的复杂度降低、使用条件降低，大幅度提高了算法的效率和算法的鲁棒性。具体的，步骤A包括：A1.获取白内障图像其中m和n分别表示白内障图像的长和宽，表示矩阵空间；A2.初始化所有层的白内障图像每层的初始图像均为原始图像，即I(t)＝I,t＝1,2…,w。具体的，步骤B包括：B1.将第t层的初始图像I(t)从RGB颜色空间转换到Lab颜色空间，得到图像为图像L(t)均匀分配k(t)个种子点，划分为k(t)个超像素，计算每个相邻种子点之间的步长S；B2.在每个种子点的s×s邻域内计算该邻域内所有像素点的梯度值，并将种子点移到该邻域内梯度最小的地方，其中s为所述领域的边长值，s可以根据实际情况按经验取值；B3.以每个种子点的2S×2S邻域内为搜索范围，对每个搜索到的像素点，分别计算该像素点和种子点之间的距离D；B4.由于每个像素点都会被多个种子点搜索到，所以每个像素点都会有一组与周围种子点的距离。每个像素点被所有种子点搜索完后，分别得到一组与周围种子点的距离，取各组里距离最小值对应的种子点作为其对应像素点的聚类中心；B5.计算每个类内所有像素的均值和坐标重心，将每个类作为一个超像素，将计算得到的坐标重心作为其对应的超像素的新种子点位置；B6.重复迭代步骤B2到B5得聚类结果和像素标签kLabel(t)∈Rm×n，以及图像分块结果。进一步的，步骤B1中计算所述步长S为：进一步的，步骤B3中计算所述距离D为：其中，dc表示像素间颜色的距离，ds表示像素间空间位置的距离，i和j分别表示第i个种子点和其搜索范围内的第j个像素，xi、xj、yi、yj分别表示像素点i和j的横坐标、纵坐标，li、lj、ai、aj、bi、bj像素点i和j的在Lab空间三通道的值。具体的，步骤C包括：C1.计算每个图像块的光源估计计算如下：其中，表示求梯度，Ω表示图像块内所有像素，I(t)为第t层的初始图像；C2.定义第t层像素照明光源估计的初始值初始值为：其中，xu表示第u个像素的横坐标，yu表示第u个像素的纵坐标；C3.更新第t层像素照明光源估计值其中dj(xu本文档来自技高网...

【技术保护点】
1.白内障图像的颜色校正方法，其特征包括：A.设定像素光源估计层数w，每层设定不同聚类种子点的个数k，并将每一层的图像初始化；B.利用简单线性迭代聚类将每一层的图像分别划分为图像块；C.计算出每个图像块的光源估计，将其反投影到原始图像上，得到每一层的像素照明光源估计；D.将所有层的像素照明光源估计结果进行线性融合，得到最终的像素照明光源估计结果，根据像素照明光源估计对原始图像进行恢复，得到颜色校正后的结果图。

【技术特征摘要】
1.白内障图像的颜色校正方法，其特征包括：A.设定像素光源估计层数w，每层设定不同聚类种子点的个数k，并将每一层的图像初始化；B.利用简单线性迭代聚类将每一层的图像分别划分为图像块；C.计算出每个图像块的光源估计，将其反投影到原始图像上，得到每一层的像素照明光源估计；D.将所有层的像素照明光源估计结果进行线性融合，得到最终的像素照明光源估计结果，根据像素照明光源估计对原始图像进行恢复，得到颜色校正后的结果图。2.如权利要求1所述的白内障图像的颜色校正方法，其特征为：步骤A包括：A1.获取白内障图像其中m和n分别表示白内障图像的长和宽，表示矩阵空间；A2.初始化所有层的白内障图像每层的初始图像均为原始图像，即I(t)＝I,t＝1,2…,w。3.如权利要求1所述的白内障图像的颜色校正方法，其特征为：步骤B包括：B1.将第t层的初始图像I(t)从RGB颜色空间转换到Lab颜色空间，得到图像为图像L(t)均匀分配k(t)个种子点，划分为k(t)个超像素，计算每个相邻种子点之间的步长S；B2.在每个种子点的s×s邻域内计算该邻域内所有像素点的梯度值，并将种子点移到该邻域内梯度最小的地方，其中s为所述领域的边长值；B3.以每个种子点的2S×2S邻域内为搜索范围，对每个搜索到的像素点，分别计算该像素点和种子点之间的距离D；B4.每个像素点被所有种子点搜索完后，分别得到一组与周围种子点的距离，取各组里距离最小值对应的种子点作为其对应像素点的聚类中心；B5.计算每个类内所有像素的均值和坐标重心，将每个类作为一个超像素，将计算得到的坐标重心作为其对应的超像素的新种子点位置；B6.重复迭代步骤B2到B5得聚类结果和像素标签kLabel(t)∈Rm×n，以及图像...

【专利技术属性】
技术研发人员：范文澜，彭真明，王慧，赵学功，王卓然，贺晨，袁国慧，曲超，彭凌冰，孔轩，蒲恬，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51