CN115660990A - 基于亮度分类的内窥镜图像镜面反射检测与修复方法 - Google Patents

基于亮度分类的内窥镜图像镜面反射检测与修复方法 Download PDF

Info

Publication number
CN115660990A
CN115660990A CN202211345244.4A CN202211345244A CN115660990A CN 115660990 A CN115660990 A CN 115660990A CN 202211345244 A CN202211345244 A CN 202211345244A CN 115660990 A CN115660990 A CN 115660990A
Authority
CN
China
Prior art keywords
image
brightness
highlight
endoscope
classification
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
Application number
CN202211345244.4A
Other languages
English (en)
Inventor
徐超
聂超
李正平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui University
Original Assignee
Anhui University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Anhui University filed Critical Anhui University
Priority to CN202211345244.4A priority Critical patent/CN115660990A/zh
Publication of CN115660990A publication Critical patent/CN115660990A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Endoscopes (AREA)

Abstract

本发明公开了一种基于亮度分类的内窥镜图像镜面反射检测与修复方法,其应用于图像处理技术领域,将原始内窥镜图像根据平均亮度进行图像分类;对分类结果进行镜面反射检测,得到全局高光二值图像和局部高光二值图像,并进行融合,同时基于COV的判别方法去除饱和非亮点区域;基于改进的Criminisi算法来修复图像亮点区域。本发明不仅能准确检测出各种亮度的内窥镜图像中的高光区域,还能避免检测到中等亮度图像中大的饱和非亮点区域。对高光区域进行修复,不仅能够保留器官表面的纹理、结构和细节,修复时间也大大缩短。特别是对于低亮度图像,本方案不仅能实现亮度增强,还能准确检测出增强图像的镜面反射区域并修复,提升了图像质量。

Description

基于亮度分类的内窥镜图像镜面反射检测与修复方法
技术领域
本发明涉及图像处理技术领域,更具体的说是涉及基于亮度分类的内窥镜镜面图像反射检测与修复方法。
背景技术
目前对内窥镜图像镜面反射检测的方法主要分为基于颜色空间的方法、基于深度学习的方法、基于双色反射模型的方法。由于从不同类型的MIS程序或医疗扫描仪获得的图像数据的特征存在很大的差异,现有镜面反射检测方法的主要缺陷是普适性差、计算复杂度高,误检测率高。大多数检测方法局限于特定的图像,在低亮度图像、图像整体亮度较高、器官表面纹理复杂的情况下无法正确检测到镜面反射区域,这些方法可能更适用于检测某一类中等亮度内镜图像中高强度的镜面高光,但不能适应不同类型内窥镜图像序列的颜色变化。
1.基于颜色空间的方法;
基于颜色空间的方法与机器学习技术相比计算量更少,因此在实践中具有优势。它利用不同模式(如HSV、YUV、RGB)下的颜色分布特征来检测内窥镜图像中的高光,它们的成功实施通常涉及到经验参数的选择,参数不能在不同场景下进行自适应调整,并且对于低亮度图像和高亮度图像,检测性能会变差,检测效果不佳。
2.基于双色反射模型的方法;
双色反射模型广泛应用于自然图像中的镜面反射检测,具体来说,该方法利用强度比从图像中提取镜面反射和漫反射分量,漫反射分量表示真实的颜色信息,镜面反射分量表示高光。基于双色反射模型的方法之间的差异在于如何执行这种分离。该方法的性能受到许多因素的影响,如对象的材质、运动和结构,因此它适用于静态的自然图像。但是内窥镜图像内亮点为粘液的镜面反射,亮度过高,漫反射成分很低甚至没有,使得高亮区域检测及修复效果很差。
3.基于深度学习的方法;
相比于基于颜色空间的方法,使用深度学习对内窥镜图像中的镜面反射区域进行检测可以获得更高的精度。虽然这类方法不依赖于许多假设,但它们的计算复杂度很高,不适合实际应用。此外,这些方法大多需要大量完全标记的数据。
现有的研究以及相应的解决方案都只集中在一个具体的问题上,即要么只检测与修复图像中的亮点,要么只增强低亮度图像,现有镜面反射检测方法的主要缺陷是普适性差、计算复杂度高、误检测率高。大多数检测方法局限于特定的图像,在低照度图像、图像整体亮度较高、器官表面纹理复杂的情况下无法正确检测到镜面反射区域,这些方法可能更适用于检测某一类中等亮度内镜图像中高强度的镜面高光,但不能适应不同类型内窥镜图像序列的颜色变化。
发明内容
有鉴于此,本发明提供了基于亮度分类的内窥镜镜面图像反射检测与修复方法,以解决背景技术中提及的技术问题。
为了实现上述目的,本发明采用如下技术方案:
一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,包括如下具体步骤:
将原始内窥镜图像根据平均亮度进行图像分类;
对分类结果进行镜面反射检测,得到全局高光二值图像和局部高光二值图像,并进行融合,同时基于COV的判别方法去除饱和非亮点区域;
基于改进的Criminisi算法来修复图像亮点区域。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,图像分类具体步骤如下:基于图像平均亮度设定阈值来进行图像分类,分类规则如下:
Figure BDA0003918158830000031
Figure BDA0003918158830000032
其中,I为输入的原始内窥镜图像,la为图像的平均亮度,La为预期正常图像的全局平均亮度,Ih是I经过分类识别为高亮度图像,Im是I经过分类识别为中等亮度图像,Il是I经过分类识别为低亮度图像。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,得到局部高光二值图像的具体步骤如下:
对低亮度图像增强图像的亮度分量;
将待检测高光图像转化为灰度图,使用边缘检测算子Laplacian算子计算出灰度图的梯度图像,再剔除低亮度边缘;
细化检测得到局部高光二值图像。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,对低亮度图像增强图像的亮度分量具体步骤如下:
将低亮度图像从RGB颜色空间转到HSV颜色空间,提取V通道分量,对其进行自适应伽马校正,得到低频分量LFγ
用奇异值分解提取原始V通道的低频分量LF和伽马校正值通道Vγ的低频分量LFγ对应的奇异值矩阵Δ和Δγ
计算均衡化奇异值矩阵;
对均衡化奇异值矩阵执行逆奇异值分解,以生成均衡的低频分量LFeq
使用软阈值,对原始V通道高频分量进行去噪,利用奇异值均衡的低频分量LFeq和去噪后的高频分量HF执行逆小波变换,得到原始内窥镜图像增强后的V通道分量;再将增强后的HSV颜色空间图像转到RGB空间,得到亮度增强图像。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,得到全局高光二值图像的具体步骤如下:
采用高帽滤波和低帽滤波增强输入的待检测高光图像的对比度;
提取对比度增强图像的G、B通道分量,再分别计算G、B通道分量的标准差;
定义基于图像亮度的自适应阈值函数检测出全局高光二值图像。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,基于COV的判别方法去除饱和非亮点区域具体步骤如下:
利用高光区域轮廓长度是否大于c,判断是否存在潜在的饱和非高光区域,若存在,计算出潜在区域的cov值cv1、其他高光区域的cov值cv2,
C={C1,C2,C3……Ci……Cn} (2-8);
Figure BDA0003918158830000041
Figure BDA0003918158830000042
在式2-8中,C为亮点轮廓集合,ci为亮点轮廓Ci的长度,γ为调节参数,n为亮点个数,σ为标准差,μ为均值;
若cv1≥cv2,则为大亮点区域,无需去除,若cv1<cv2,则为非高光区域,则在检测到的高光结果中去除该大尺寸区域。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,改进的Criminisi算法来修复图像亮点区域具体步骤如下:
确定优先级计算范围为该单亮点边缘;
基于改进的优先级函数计算亮点边缘优先级,并确定优先级最大的像素点;
基于自适应搜索范围策略搜索最佳匹配块;
将最佳匹配块中的像素对应复制到目标块中的未知像素点,更新置信度和亮点边缘直至该亮点修复完成,最后依次修复每个亮点区域直至修复完成。
优选的,在上述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法中,改进的优先级函数:
Figure BDA0003918158830000051
曲率项:
Figure BDA0003918158830000052
上式中,w表示控制置信项、数据项曲线平滑度的正则化因子,0<w<1,
Figure BDA0003918158830000053
β是权重因子,
Figure BDA0003918158830000054
经由上述的技术方案可知,与现有技术相比,本发明公开提供了一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,本发明提出的检测算法不仅能准确检测出高光区域,还能避免检测到中等亮度图像中大的饱和非亮点区域。通过改进的Criminisi函数对高光区域进行修复,不仅能够保留器官表面的纹理、结构和细节,修复时间也大大缩短。特别是对于低亮度图像,本发明不仅能实现亮度增强,还能准确检测出增强图像的镜面反射区域并修复,提升了图像质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明提供的方法流程图;
图2为本发明提供的图像分类及低亮度图像亮度增强流程图;
图3为本发明提供的高光图像的R、G、B三通道直方图;
图4为本发明提供的去除错误检测示意图;
图5为本发明提供的亮点区域及其邻域示意图;
图6为本发明提供的镜面反射检测结果对比图;
图7为本发明提供的内窥镜图像镜面反射检测与修复结果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,提出的检测算法不仅能准确检测出高光区域,还能避免检测到中等亮度图像中大的饱和非亮点区域。通过改进的Criminisi函数对高光区域进行修复,不仅能够保留器官表面的纹理、结构和细节,修复时间也大大缩短。特别是对于低亮度图像,本发明不仅能实现亮度增强,还能准确检测出增强图像的镜面反射区域并修复,提升了图像质量。
本发明的实施例公开了一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,如图1所示。
步骤1图像分类及低亮度图像增强。首先,基于图像平均亮度设定阈值来进行图像分类,将原始内窥镜图像分为高亮度图像h1,中等亮度图像m1,低亮度图像l1三类。其次,对于低亮度图像l1,提取图像V通道分量,对V分量进行自适应伽马校正(AGC),对原始和伽玛校正后的V通道执行离散小波分解,对低频分量进行奇异值均衡,得到低亮度图像增强图像l2。
步骤2镜面反射检测。首先采用高帽滤波和低帽滤波增强图像对比度,提取增强图像的g,b通道,接着使用提出的自适应阈值函数检测出全局高光二值图像,再使用Laplacian算子对原始图像的灰度图进行滤波得到梯度图像,剔除低亮度边缘得到局部高光二值图像,结合全局高光二值图像、局部高光二值图像,对检测到的掩膜采用形态学膨胀操作来扩大与镜面反射对应的区域,得到边界扩充后的高亮点区域掩膜图像,再基于亮点轮廓长度判断是否存在潜在饱和非亮点区域,若存在,基于COV的判别方法以去除饱和非亮点区域。
步骤3镜面反射修复。基于改进的Criminisi算法来修复图像亮点区域。首先确定待修复亮点边缘,提出改进的优先级计算函数,计算单个亮点上的像素点优先级,基于自适应搜索范围策略搜索匹配块,将匹配准则修改为平均SSD距离,再将最佳匹配块中的像素点对应复制到目标块中的未知像素点,更新置信度和亮点边缘直至该亮点修复完成,最后依次修复每个亮点区域直至修复完成。
在步骤1中,图像分类及低亮度图像亮度增强算法:
由于从不同类型的MIS程序或医疗扫描仪获得的图像数据存在很大的亮度差异,会对后续的镜面反射检测算法产生较大影响,因此基于图像平均亮度设定阈值来进行图像分类,分类规则如下:
Figure BDA0003918158830000071
Figure BDA0003918158830000072
其中,I为输入的原始内窥镜图像,la为图像的平均亮度,La为预期正常图像的全局平均亮度,Ih是I经过分类识别为高亮度图像,Im是I经过分类识别为中等亮度图像,Il是I经过分类识别为低亮度图像,经过大量实验得到,当T1值选用112,t1值选用0.3时,分类效果最佳。
对于低亮度图像,需要充分增强图像的亮度分量。首先将低亮度图像从RGB颜色空间转到HSV颜色空间,提取V通道分量,对其进行自适应伽马校正,为了避免图像颜色失真,保留H、S分量。伽马修正函数定义为:
Vγ=cVγ (1-3)
其中Vγ为伽马校正值通道(V),c和γ为根据图像信息塑造变换曲线的自适应参数。
利用μ+σ与μ-σ的差值D来判断图像是低对比度还是中等对比度:
D=(μ+σ)-(μ-σ) (1-4)
Figure BDA0003918158830000081
时图像为低对比度,当
Figure BDA0003918158830000082
时图像为中等对比度,参数τ设为3。
对于低对比度的图像,计算γ的公式为(1-5):
γ=-log2(σ) (1-5)
c的值为:
Figure BDA0003918158830000083
k的值为:
k=Vγ+(1-Vγ)×μγ (1-7)
对于对比度中等的图像,用式(1-8)计算γ:
Figure BDA0003918158830000084
对于低亮度图像,c依赖于k,因此将较低的强度级扩展到更高的强度级,对于σ较低的图像,变换曲线更陡。
奇异值分解(SVD)被认为可以解决照明问题,奇异值矩阵表示输入图像的强度信息,对奇异值矩阵的任何修改都会影响图像的强度。用奇异值分解(SVD)提取原始V通道的低频分量LF和伽马校正值通道Vγ的低频分量LFγ对应的奇异值矩阵Δ和Δγ
LFSVD=UΔVT (1-9)
Figure BDA0003918158830000091
其中,U和V对应低频分量LF,Uγ和Vγ对应低频分量LFγ
增强因子∈为:
Figure BDA0003918158830000092
计算均衡化奇异值矩阵Δeq:
Δeq=∈Δ (1-12)
对均衡奇异值矩阵执行逆SVD,以生成均衡的低频分量LFeq:
LFeq=UΔeqVT (1-13)
使用软阈值,对原始V通道高频分量进行去噪,利用奇异值均衡的低频分量LF和去噪后的高频分量HF执行逆小波变换,得到原始内窥镜图像增强后的V通道分量。再将增强后的HSV颜色空间图像转到RGB空间,得到亮度增强图像lv。
在步骤2中,镜面反射检测算法:
通常情况下,在内窥镜图像非高光区域,由于血红蛋白的存在,大多数内窥镜图像是红色的,并且红色通道R的数值高于绿色G和蓝色B通道的数值。然而,在高光区域下,三个通道的数值几乎相同且很高,特别是对于较大的高光区域。图3为含有镜面反射区域的内窥镜图像的R、G、B三通道的直方图,由图3可见,G、B通道颜色在漫反射和镜面反射之间具有较大区分度。
考虑到上述信息,对于全局高光,本实施例提出的自适应阈值函数只考虑G、B通道的颜色变化,首先采用高帽滤波和低帽滤波增强输入的待检测高光图像的对比度,接着提取对比度增强图像的G、B通道分量,再分别计算G、B通道分量的标准差σgb:
Figure BDA0003918158830000093
Figure BDA0003918158830000101
Figure BDA0003918158830000102
Figure BDA0003918158830000103
其中,n是图像像素数量,xi是当前像素值,μgb分别是G、B通道分量的均值。
自适应阈值函数Th定义为:
Figure BDA0003918158830000104
Figure BDA0003918158830000105
式(2-6)中,max(g,b)为g,b通道强度最大值,τ为根据图像亮度自适应调整的参数。当图像整体亮度偏高时,固定参数的自适应阈值会导致阈值较小,将更多的亮度偏高的非高光区域误检测为高光;当图像整体亮度偏低时,固定参数的自适应阈值会导致阈值较大,将遗漏更多的正确的高亮像素,因此,对于高亮度图像,自适应参数应设置小一点以保证阈值大一点,对于低亮度图像,自适应参数应设置大一点以保证阈值小一点。经过大量统计实验,确定高亮度图像的τ为0.3,中等亮度图像的τ为0.8,亮度增强后的低亮度图像的τ为1.1最合适。
而对于局部高光二值图像,考虑了梯度信息的大小,首先将待检测高光图像转化为灰度图,再使用边缘检测算子Laplacian算子计算出灰度图的梯度图像,再剔除低亮度边缘,阈值为Th1,进一步细化检测得到局部高光二值图像。
Th1=Th-10 (2-7)
融合全局高光二值图像、局部高光二值图像,由于光晕效应导致亮点周围可能有黑圈或者彩圈,采用形态学膨胀操作扩大检测到的亮点区域,本实施例使用的是结构元素是大小为3的圆盘形状。对于中等亮度图像,图像中可能会存在明亮的饱和非高光区域,如图4(a)左下角区域,现有的算法将其视之为高光,本实施例基于亮点轮廓长度和变异系数(COV)去除这种错误分割。首先利用高光区域轮廓长度是否大于c这一信息判断是否存在潜在的饱和非高光区域,若存在,计算出该潜在区域(图4(b)蓝框)的cov值cv1、其他高光区域(图4(b)紫框)的cov值cv2,
C={C1,C2,C3……Ci……Cn} (2-8)
Figure BDA0003918158830000111
Figure BDA0003918158830000112
在式2-8中,C为亮点轮廓集合,ci为亮点轮廓Ci的长度,γ为调节参数,n为亮点个数,本实施例所有实验γ值均取4;σ和μ为区域内的标准差和均值。
变异系数(COV)反映了区域内的强度均匀性,一般用来定量评估无地面真相图像的修复结果,由于高光区域的cov值要远大于无高光区域,若cv1≥cv2,则为大亮点区域,无需去除,若cv1<cv2,则为非高光区域,则在检测到的高光结果中去除该大尺寸区域。
在步骤3中,改进的Criminisi图像修复算法
如图5中:Ω表示待修复的亮点区域,
Figure BDA0003918158830000113
表示亮点区域的边界线,Ψp是以位于亮点区域边界线上的以p点为中心的待修复样本块,np作为p点的法线矢量,
Figure BDA0003918158830000114
表示p点的等照度线方向,Φ表示信息完好区域。
Criminisi算法根据优先级函数来确定图像修复顺序,以保持结构完整连贯,但可能会出现两种现象,会对填充顺序产生负面影响,从而影响修复结果。第一个是置信项C(p)可能会呈指数下降,这种现象导致了C(p)和D(p)之间的不平衡,从而使得现有的修复算法不适合处理某些含有强结构区域的内窥镜图像。有两种方法可以避免这种现象。第一种是用于控制C(p)下降的正则化项,第二种是将优先项的乘法定义改为加权求和的方式。第二种现象是,当置信项可能很大时,数据项D(p)可能为零,为了避免产生错误的填充顺序,现有研究方法提出了一个曲率项与数据项相结合的方法。
综上所述,为了使优先级函数具有鲁棒性,本实施例提出了一个新的优先级函数计算公式:
Figure BDA0003918158830000121
曲率项:
Figure BDA0003918158830000122
上式中,w表示控制置信项、数据项曲线平滑度的正则化因子,0<w<1,
Figure BDA0003918158830000123
β是权重因子,
Figure BDA0003918158830000124
在本实施例中,w=0.7,
Figure BDA0003918158830000125
β=0.7,
Figure BDA0003918158830000126
小于β的目的在于优先修复结构信息较为明显的区域。
由于内窥镜图像中镜面反射区域的大小,位置都是随机的,当我们首次确定待修复单亮点边缘时,只需要计算该单个亮点边缘像素点的优先级大小,找到优先级最大的像素点,以其为中心确定一个9×9的待修复目标块,同时采用自适应搜索范围策略,根据亮点轮廓长度lc采取不同的扩张系数n来确定搜索范围:
扩张系数n:
Figure BDA0003918158830000127
其中a由图片分辨率决定,n1,n2,按照图像纹理结构复杂度不同而进行不同的调整,lc为高亮点轮廓长度。
搜索范围的上下左右边界us,ds,ls,rs
us=yup-nh (3-3)
ds=ydown+nh (3-4)
ls=xleft-nw (3-5)
rs=xright+nw (3-6)
其中,yup为亮点上极点的纵坐标,ydown为下极点的纵坐标,xleft为左极点的横坐标,xright为右极点的横坐标。
在一定程度上,平均SSD距离能更好地反映样本快的平均相似度。平均SSD距离定义为:
Figure BDA0003918158830000131
其中
Figure BDA0003918158830000132
为每个颜色通道的强度平均值。
在搜索范围内搜索出与目标块之间平均SSD距离最小的样本块即为最佳匹配块,将匹配块中的像素对应复制到目标块中的未知像素点,更新置信度和亮点边缘直至该亮点修复完成,最后依次修复每个亮点区域直至修复完成。
图6显示了每一种检测方法所识别出的镜面反射区域(绿色部分),第一行从左至右分别为亮图像1(h1),亮图像2(h2),中等图像1(m1),低亮度图像1(l1),第二行到第五行分别为其他算法的检测结果,第六行为提出算法的检测结果。相比于其他算法,结果的可视化表示清楚地表明提出的算法镜面反射区域识别结果更加精确。
图7显示了本实施例具体每一步的效果图。结果表明,本发明所提出的检测算法能够实现对不同亮度的内窥镜图像的分类,有效地增强低亮度图像的亮度,准确检测出不同亮度图像的镜面反射区域,有效地修复内窥镜图像镜面反射区域,在很大程度上保留了原图像纹理、细节和结构。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (8)

1.一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,包括如下具体步骤:
将原始内窥镜图像根据平均亮度进行图像分类;
对分类结果进行镜面反射检测,得到全局高光二值图像和局部高光二值图像,并进行融合,同时基于COV的判别方法去除饱和非亮点区域;
基于改进的Criminisi算法来修复图像亮点区域。
2.根据权利要求1所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,图像分类具体步骤如下:基于图像平均亮度设定阈值来进行图像分类,分类规则如下:
Figure FDA0003918158820000011
Figure FDA0003918158820000012
其中,I为输入的原始内窥镜图像,la为图像的平均亮度,La为预期正常图像的全局平均亮度,Ih是I经过分类识别为高亮度图像,Im是I经过分类识别为中等亮度图像,Il是I经过分类识别为低亮度图像。
3.根据权利要求2所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,得到局部高光二值图像的具体步骤如下:
对低亮度图像增强图像的亮度分量;
将待检测高光图像转化为灰度图,使用边缘检测算子Laplacian算子计算出灰度图的梯度图像,再剔除低亮度边缘;
细化检测得到局部高光二值图像。
4.根据权利要求3所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,对低亮度图像增强图像的亮度具体步骤如下:
将低亮度图像从RGB颜色空间转到HSV颜色空间,提取V通道分量,对其进行自适应伽马校正,得到低频分量LFγ
用奇异值分解提取原始V通道的低频分量LF和伽马校正值通道Vγ的低频分量LFγ对应的奇异值矩阵Δ和Δγ
计算均衡化奇异值矩阵;
对均衡化奇异值矩阵执行逆奇异值分解,以生成均衡的低频分量LFeq
使用软阈值,对原始V通道高频分量进行去噪,利用奇异值均衡的低频分量LFeq和去噪后的高频分量HF执行逆小波变换,得到原始内窥镜图像增强后的V通道分量;再将增强后的HSV颜色空间图像转到RGB空间,得到亮度增强图像。
5.根据权利要求3所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,得到全局高光二值图像的具体步骤如下:
采用高帽滤波和低帽滤波增强输入的待检测高光图像的对比度;
提取对比度增强图像的G、B通道分量,再分别计算G、B通道分量的标准差;
定义基于图像亮度的自适应阈值函数检测出全局高光二值图像。
6.根据权利要求1所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,基于COV的判别方法去除饱和非亮点区域具体步骤如下:
利用高光区域轮廓长度是否大于c,判断是否存在潜在的饱和非高光区域,若存在,计算出潜在区域的cov值cv1、其他高光区域的cov值cv2,
C={C1,C2,C3……Ci……Cn} (2-8);
Figure FDA0003918158820000021
Figure FDA0003918158820000022
在式2-8中,C为亮点轮廓集合,ci为亮点轮廓Ci的长度,γ为调节参数,n为亮点个数,σ为标准差,μ为均值;
若cv1≥cv2,则为大亮点区域,无需去除,若cv1<cv2,则为非高光区域,则在检测到的高光结果中去除该大尺寸区域。
7.根据权利要求1所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,改进的Criminisi算法来修复图像亮点区域具体步骤如下:
确定优先级计算范围为该单亮点边缘;
基于改进的优先级函数计算亮点边缘优先级,并确定优先级最大的像素点;
基于自适应搜索范围策略搜索最佳匹配块;
将最佳匹配块中的像素对应复制到目标块中的未知像素点,更新置信度和亮点边缘直至该亮点修复完成,最后依次修复每个亮点区域直至修复完成。
8.根据权利要求7所述的一种基于亮度分类的内窥镜镜面图像反射检测与修复方法,其特征在于,改进的优先级函数:
Figure FDA0003918158820000031
曲率项:
Figure FDA0003918158820000032
上式中,w表示控制置信项、数据项曲线平滑度的正则化因子,0<w<1,
Figure FDA0003918158820000033
β是权重因子,
Figure FDA0003918158820000034
CN202211345244.4A 2022-10-31 2022-10-31 基于亮度分类的内窥镜图像镜面反射检测与修复方法 Pending CN115660990A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211345244.4A CN115660990A (zh) 2022-10-31 2022-10-31 基于亮度分类的内窥镜图像镜面反射检测与修复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211345244.4A CN115660990A (zh) 2022-10-31 2022-10-31 基于亮度分类的内窥镜图像镜面反射检测与修复方法

Publications (1)

Publication Number Publication Date
CN115660990A true CN115660990A (zh) 2023-01-31

Family

ID=84992837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211345244.4A Pending CN115660990A (zh) 2022-10-31 2022-10-31 基于亮度分类的内窥镜图像镜面反射检测与修复方法

Country Status (1)

Country Link
CN (1) CN115660990A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116703798A (zh) * 2023-08-08 2023-09-05 西南科技大学 基于自适应干扰抑制的食管多模态内镜图像增强融合方法
CN117314839A (zh) * 2023-09-08 2023-12-29 深圳市汇天益电子有限公司 基于工业内窥镜的无损检测方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010007599A1 (en) * 1999-12-28 2001-07-12 Ryosuke Iguchi Image processing method and image processing apparatus
KR101875004B1 (ko) * 2017-01-04 2018-07-05 금오공과대학교 산학협력단 캡슐 내시경 영상을 이용한 자동 출혈 감지 방법 및 컴퓨터 프로그램
CN113989147A (zh) * 2021-10-26 2022-01-28 安徽大学 基于直方图均衡和改进反锐化掩模的内窥镜图像增强方法
CN114529475A (zh) * 2022-02-23 2022-05-24 安徽大学 基于二维伽马校正和色调映射的图像增强方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010007599A1 (en) * 1999-12-28 2001-07-12 Ryosuke Iguchi Image processing method and image processing apparatus
KR101875004B1 (ko) * 2017-01-04 2018-07-05 금오공과대학교 산학협력단 캡슐 내시경 영상을 이용한 자동 출혈 감지 방법 및 컴퓨터 프로그램
CN113989147A (zh) * 2021-10-26 2022-01-28 安徽大学 基于直方图均衡和改进反锐化掩模的内窥镜图像增强方法
CN114529475A (zh) * 2022-02-23 2022-05-24 安徽大学 基于二维伽马校正和色调映射的图像增强方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
聂超: "内窥镜图像镜面反射区域检测与修复方法研究", 万方数据库, 8 October 2022 (2022-10-08) *
雯霞: "基于样本块的数字图像修复算法研究", 中国优秀硕士学位论文全文数据库 (信息科技辑), 30 April 2019 (2019-04-30), pages 138 - 1097 *
马坤: "基于DWT-SVD和拐点伽马的木材彩色图像增强与识别", 中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑), 31 January 2020 (2020-01-31), pages 024 - 786 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116703798A (zh) * 2023-08-08 2023-09-05 西南科技大学 基于自适应干扰抑制的食管多模态内镜图像增强融合方法
CN116703798B (zh) * 2023-08-08 2023-10-13 西南科技大学 基于自适应干扰抑制的食管多模态内镜图像增强融合方法
CN117314839A (zh) * 2023-09-08 2023-12-29 深圳市汇天益电子有限公司 基于工业内窥镜的无损检测方法及系统

Similar Documents

Publication Publication Date Title
CN115660990A (zh) 基于亮度分类的内窥镜图像镜面反射检测与修复方法
CN109035253A (zh) 一种语义分割信息指导的深度学习自动图像抠图方法
CN111199524B (zh) 一种针对可调光圈光学系统的图像紫边校正方法
CN111161222B (zh) 一种基于视觉显著性的印刷辊筒缺陷检测方法
CN109978848B (zh) 基于多光源颜色恒常模型检测眼底图像中硬性渗出的方法
CN105046653B (zh) 一种视频雨滴去除方法及系统
CN110782407B (zh) 一种基于天空区域概率分割的单幅图像去雾方法
CN110717865B (zh) 图片检测方法及装置
CN116542982B (zh) 一种基于机器视觉的出发判断器缺陷检测方法及装置
CN115511907B (zh) 一种用于led屏幕的划痕检测方法
CN103729828B (zh) 视频去雨方法
CN110147816B (zh) 一种彩色深度图像的获取方法、设备、计算机存储介质
CN117274113B (zh) 基于图像增强的碎硅片清洗效果视觉检测方法
Trongtirakul et al. Single backlit image enhancement
CN115809966A (zh) 一种低照度图像的增强方法和系统
Yu et al. Efficient highlight removal of metal surfaces
Kuzovkin et al. Descriptor-based image colorization and regularization
JP6375138B2 (ja) パープルフリンジ除去処理方法及びその処理を遂行するパープルフリンジ除去処理装置
CN110930358B (zh) 一种基于自适应算法的太阳能面板图像处理方法
CN112270683B (zh) 一种ihc数字预览图识别和组织前景分割方法及系统
CN109801238B (zh) 一种基于简化大气散射模型的弱光照图像增强方法
CN110473190B (zh) 一种基于尺度的自适应织物疵点检测方法
CN109934190B (zh) 基于变形高斯核函数的自适应强光人脸图像纹理恢复方法
CN112949389A (zh) 一种基于改进目标检测网络的雾霾图像目标检测方法
Thayammal et al. Performance analysis of image denoising using deep convolutional neural network

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