一种针对折射性色差的快速检测校正算法制造技术

本发明专利技术公开了一种针对折射性色差的快速检测校正算法，首先求出所得图像绿色通道在垂直方向和水平方向的梯度；然后将该梯度图像以阈值T进行过滤，得到较高梯度值的图像，即对应的包含折射性色差的区域；利用二值化分区方法将不同强光区域区分开来；依次提取每个区域，将该区域的边界进行适当扩充，从而得到折射性色差区域的边界；再利用现有的色差校正方法进行校正。本方法主要针对除了强光区域周围色差之外的普通折射性色差，只需简单利用绿色通道的梯度值就可以检测出主要色差区域边界，在图像处理领域具有非常重要的意义。

【技术实现步骤摘要】

本专利技术主要涉及到数字图像处理领域，特指针对折射性色差的快速检测校正算法。
技术介绍
目前，单反相机、摄像机等成像设备在人们的日常生活中发挥越来越大的作用。但是受到光学镜片自身特性的影响，实际光学系统与理想光学系统有很大的差异，即一个物点发出的光线经实际光学系统后，不再会聚于一点，而是一个弥散斑，弥散斑的大小与光学系统的像差有关。光学系统的像差具体可分为单色像差与色差。其中，色差是彩色成像设备中常见的一种畸变。色差产生的主要原因是光学系统的光学系统的色散效应，光学系统对穿过透镜的不同波长的色光具有不同的折射率，从而该成像系统不能将不同色光很好地聚焦在一起。根据在图像中的不同位置，色差又可以简单分为两类：(1)强光区域周围的色差。主要是由于强光包含大量白光，白光向周围扩散时不同波长的单色光出现误差，从而导致强光区域周围一般会出现较大的色差。(2)一般折射性色差。除了强光区域，如果某个边界两边颜色差别比较大，则由于色散效应，该边界处会出现色差。这种色差可以称为折射性色差，如图1所示。如果同时处理这两种不同位置处的色差，会影响算法的整体稳健性和速度。强光区域周围的色差可以通过计算彩色图像的纯度V而快速检测并校正，但是利用纯度V却无法检测出其他边界处的折射性色差。因此，如何有效去除折射性色差也是图像处理领域亟需解决的一个问题。色差具体可以分为两种类型：垂轴色差和轴向色差，如图2所示。目前色差校正的算法有单独针对轴向色差的校正算法，但是垂轴色差不能很好地消除。有的算法虽然针对两个类型的色差，但是对强光区域周围的色差，检测色差边界的方法比较复杂，而且检测不够准确。直接通过图像计算消除色差的方法较少，更多的方法是基于镜头和相机参数，这类方法需要的参数较多，计算过程复杂。
技术实现思路
针对除了强光区域周围色差之外的普通折射性色差，本专利技术提出一种针对折射性色差的快速检测校正算法，其是一种基于绿色通道梯度值的折射性色差快速检测校正算法。由绿色通道的梯度检测出包含折射性色差的主要边界区域，然后通过适当扩展得到折射性的色差区域的边界，并进行校正，可以有效地消除图像较大对比度边界处的折射性色差。本专利技术的技术方案是：一种针对折射性色差的快速检测校正算法，包括以下步骤：步骤一：利用相机拍摄含有边界的图像，首先求出其绿色通道G在垂直方向的梯度，采用垂直Sobel算子，如式(1)所示： S G _ v ( i , j ) = G ( i - 1 , j - 1 ) + 2 G ( i - 1 , j ) + G ( i - 1 , j + 1 ) - G ( i + 1 , j - 1 ) - 2 G ( i + 1 , j ) - G ( i + 1 , j + 1 ) - - - ( 1 ) ]]>其中，表示SG_v(i,j)表示坐标(i,j)处垂直方向的Sobel梯度，i表示图像中像素点的横坐标，j表示图像中像素点的纵坐标，G(i,j)表示坐标(i,j)处的像素值。然后求出其绿色通道G在水平方向的梯度，采用水平Sobel算子，如式(2)所示： S G _ h ( i , j ) = G ( i - 1 , j - 1 ) + 2 G ( i , j - 1 本文档来自技高网...

【技术保护点】
一种针对折射性色差的快速检测校正算法，其特征在于，包括以下步骤：步骤一：利用相机拍摄含有边界的图像，首先求出其绿色通道G在垂直方向的梯度，采用垂直Sobel算子，如式(1)所示：SG_v(i,j)=G(i-1,j-1)+2G(i-1,j)+G(i-1,j+1)-G(i+1,j-1)-2G(i+1,j)-G(i+1,j+1)---(1)]]>其中，表示SG_v(i,j)表示坐标(i,j)处垂直方向的Sobel梯度，i表示图像中像素点的横坐标，j表示图像中像素点的纵坐标，G(i,j)表示图像中坐标(i,j)处的像素值；然后求出其绿色通道G在水平方向的梯度，采用水平Sobel算子，如式(2)所示：SG_h(i,j)=G(i-1,j-1)+2G(i,j-1)+G(i+1,j-1)-G(i-1,j+1)-2G(i,j+1)-G(i+1,j+1)---(2)]]>其中，表示SG_h(i,j)表示坐标(i,j)处水平方向的Sobel梯度，i表示图像中像素点的横坐标，j表示图像中像素点的纵坐标，G(i,j)表示图像中坐标(i,j)处的像素值；步骤二：针对步骤一中所求得的图像绿色通道在垂直方向和水平方向的梯度，设置阈值T，然后利用该阈值T分别对垂直方向和水平方向梯度进行过滤，将垂直方向和水平方向梯度图像中小于该阈值T的像素都令为0，其余像素值不变，则过滤得到的区域包含存在折射性色差的区域；其操作方法如下式所示：Gg(i,j)new=Gg(i,j),ifGg(i,j)≥T0,ifGg(i,j)<T---(3)]]>其中，Gg(i,j)表示梯度图像中坐标为(i,j)的像素点的像素值，Gg(i,j)new表示过滤后对应梯度图像中坐标为(i,j)的像素点的像素值；步骤三：对步骤二中得到的包含折射性色差区域的图像利用matlab中的bwlabel函数进行二值化分区，bwlabel函数能将图像中不同的连通区域区分开来，且不同区域的像素值不同，该函数能将不同色差区域区分并进行标记；步骤四：依次提取步骤三中得到的每个色差区域，检测每个色差区域间的边界并进行扩展，从而得到折射性色差的基本区域范围；对所检测色差区域的边界进行扩展是指，以某个色差区域边界最先被检测到的一点为中心，该点坐标为k，分别向与该边界垂直的两个方向扩展；垂直方向的边界扩展公式如式(4)所示：up=kv-dvd1dp=kv+dvd2---(4)]]>其中，kv表示步骤四中所检测垂直梯度图像中色差区域边界最先被检测到的点的纵坐标，dvd1和dvd2分别表示从步骤四所检测出的边界向上和向下所扩展的宽度，up和dp分别表示所检测色差区域扩充之后的上边界和下边界的纵坐标；水平方向的边界扩展公式如式(5)所示：lp=kh-dhd1rp=kh+dhd2---(4)]]>其中，kh表示步骤四中所检测水平梯度图像中色差区域边界最先被检测到的点的横坐标，dhd1和dhd2分别表示从步骤四所检测出的边界向左和向右所扩展的宽度，lp和rp分别表示所检测色差区域扩充之后的左边界和右边界的横坐标；步骤五：由步骤四中所确定的扩展后的色差区域边界，按照图像边缘颜色差异理论进行校正，得到色差消除后的图像。...

【技术特征摘要】
1.一种针对折射性色差的快速检测校正算法，其特征在于，包括以下步骤：
步骤一：利用相机拍摄含有边界的图像，首先求出其绿色通道G在垂直方向的梯度，采
用垂直Sobel算子，如式(1)所示：
S G _ v ( i , j ) = G ( i - 1 , j - 1 ) + 2 G ( i - 1 , j ) + G ( i - 1 , j + 1 ) - G ( i + 1 , j - 1 ) - 2 G ( i + 1 , j ) - G ( i + 1 , j + 1 ) - - - ( 1 ) ]]>其中，表示SG_v(i,j)表示坐标(i,j)处垂直方向的Sobel梯度，i表示图像中像素点的横坐
标，j表示图像中像素点的纵坐标，G(i,j)表示图像中坐标(i,j)处的像素值；
然后求出其绿色通道G在水平方向的梯度，采用水平Sobel算子，如式(2)所示：
S G _ h ( i , j ) = G ( i - 1 , j - 1 ) + 2 G ( i , j - 1 ) + G ( i + 1 , j - 1 ) - G ( i - 1 , j + 1 ) - 2 G ( i , j + 1 ) - G ( i + 1 , j + 1 ) - - - ( 2 ) ]]>其中，表示SG_h(i,j)表示坐标(i,j)处水平方向的Sobel梯度，i表示图像中像素点的横坐
标，j表示图像中像素点的纵坐标，G(i,j)表示图像中坐标(i,j)处的像素值；
步骤二：针对步骤一中所求得的图像绿色通道在垂直方向和水平方向的梯度，设置阈值T，
然后利用该阈值T分别对垂直方向和水平方向梯度进行过滤，将垂直方向和水平方向梯度图
像中小于该阈值T的像素都令为0，其余像素值不变，则过滤得到的区域包含存在折射性色
差的区域；其操作方法如下式所示：
G g ( i , j ) n e w = G g ( i , j ) , i f G g ( i , j ) ≥ T 0 , i f G g ( i , j ) < T - - - ( 3 ) ]]>其中，Gg(i,j)表示梯度图像中坐标为(i,j)的像素点的像素值，Gg(i,j)new表示过滤后对应
梯度图像中坐标为(i,j)的像素点的像素值；
步骤三：对步骤二中得到的包含折射性色差区域的图像利用matlab中的bwlabel函数进行
二值化分区，bwlabel函数能将图像中不同的连通区域区分开来，且不同区域的像素值不同，
该函数能将不同色差区域区分并进行标记；
步骤四：依次提取步骤三中得到的每个色差区域，检测每个色差区域间的边界并进行扩
展，从而得到折射性色差的基本区域范围；对所检测色差区域的边界进行扩展是指，以某个
色差区域边界最先被检测到的一点为中心，该点坐标为k，分别向与该边界垂直的两个方向
扩展；
垂直方向的边界扩展公式如式(4)所示：
u p = k v - dv d 1 d p = k v + dv d 2 - - - ( 4 ) ]]>其中，kv表示步骤四中所检测垂直梯度图像中色差区域边界最先被检测到的点的纵坐标，
dvd1和dvd2分别表示从步骤四所检测出的边界向上和向下所扩展的宽度，up和dp分别表示所
检测色差区域扩充之后的上边界和下边界的纵坐标；
水平方向的边界扩展公式如式(5)所示：
l p = k h - dh d 1 ...

【专利技术属性】
技术研发人员：刘煜，李卫丽，王炜，张茂军，徐玮，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南;43