基于多通道分解的屏摄图像摩尔纹消除方法技术

本发明专利技术属于数字图像处理领域，为提出现有技术不存在消除屏摄摩尔纹的方法，实现对屏摄图像中摩尔纹图案的消除。本发明专利技术采用的技术方案是，基于多通道分解的屏摄图像摩尔纹消除方法，给出多相位图层分解的模型，将该模型用在Y通道，实现摩尔纹彩色条纹和高频部分的淡化；将得到的结果输入一个改进的拉普拉斯滤波器中，实现彩色条纹的进一步消除；最后对于上一步得到的结果在其RGB通道内分别应用多相位分解图层模型，得到最终结果。本发明专利技术主要应用于数字图像处理场合。

【技术实现步骤摘要】
基于多通道分解的屏摄图像摩尔纹消除方法
本专利技术属于数字图像处理领域，具体讲，涉及基于多通道分解的屏摄图像摩尔纹消除算法。
技术介绍
摩尔纹现象是一种常见的视觉现象，当拍摄一些比较密集的物体，比如：建筑装饰物、纺织物、屏幕等时，或者通过扫描仪扫描印刷品时，会出现大面积的彩色以及不规则的条纹，也就是摩尔纹现象。然而摩尔纹现象的出现严重影响了图像质量，所以消除图像中的摩尔纹是人们广泛关注的问题。在已存在的消除摩尔纹的方法中主要针对的是纺织物以及扫描图像上的摩尔纹。这些方法主要分为两大类。首先第一类就是“预处理”，顾名思义就是在拍摄图像之前对相机进行处理。在相机的镜头前放置光学低通滤波器，或者通过调整镜头角度找到最合适的角度此时拍摄的图像中拥有最少的摩尔纹。但是该方法是以降低图像锐度为代价来消除摩尔纹，且耗费时间。另一类方法就是“后处理”，也就是当图像被拍摄之后，对损坏的图像进行处理。通过基于滤波的方法或者基于字典学习和先验的方法实现摩尔纹的消除。然而这些方法主要应用于纺织物和扫描图像上的摩尔纹，对于屏摄摩尔纹的消除效果并不明显。相对图纺织物和扫描图像上的摩尔纹，屏摄摩尔纹结果更加复杂。其产生原因主要分为两种：一是屏幕的网格和相机传感器的网格的混叠；二是图像自身的二次采样。所以屏摄摩尔纹涉及的频率范围很广，既存在相对来说频率较低的大面积的彩色条纹，有存在频率较高的摩尔纹。消除屏摄摩尔纹的难度在于：1)彩色条纹类的摩尔纹处于较低频位置，且与图像背景信息有很高的相似性，二者很难区分；2)由于屏摄图像上的条纹分布并没有规律，很难找到其特征实现摩尔纹的消除；3)摩尔纹的高频部分与图像细节所处频率十分相近，故利用简单的低通滤波器虽然可以消除高频类摩尔纹，但是图像细节也会被平滑掉，最后得到的结果会过平滑。目前，对消除屏摄摩尔纹效果比较好的就是Photoshop，但是到目前为止仍旧没有提出消除屏摄摩尔纹较为高效的算法，保证既消除摩尔纹又保留图像细节。
技术实现思路
为克服现有技术的不足，本专利技术旨在提出现有技术不存在消除屏摄摩尔纹的方法，实现对屏摄图像中摩尔纹图案的消除。本专利技术采用的技术方案是，基于多通道分解的屏摄图像摩尔纹消除方法，给出多相位图层分解的模型，将该模型用在Y通道，实现摩尔纹彩色条纹和高频部分的淡化；将得到的结果输入一个改进的拉普拉斯滤波器中，实现彩色条纹的进一步消除；最后对于上一步得到的结果在其RGB通道内分别应用多相位分解图层模型，得到最终结果。具体步骤如下：1)多相位分解模型：首先将输入的图像通过二倍下采样的方式，分为四个多相位部分，也就是四个子采样图像；2)基于块先验的图层分解模型屏摄图像I为背景图层B和摩尔纹图层M以及方差为σ2的高斯噪声n的叠加，即得加性模型：I＝B+M+n(1)采用两个先验分别对摩尔纹图层和背景图层进行规则化，公式如下：其中，为保真项，Φ(B)和Ψ(M)定义了两个先验以及约束项对进行规则化；定义Φ(B)：Φ(B)对应的公式如下：其中，Pi为一个运算符，表示提取第i个块并将其变为向量的形式，块的大小自定义，K表示高斯部分的数量，πk表示混合权重，分别表示对应第k个部分的均值和协方差矩阵，N(·)表示64维的高斯分布，b则指示该参数对应的是背景图层的参数，因为自然图像的梯度较为稀疏，因此用于自然图像的边缘平滑，表示梯度运算符，·||1表示1范数,α是可自定义的参数；定义Ψ(M)：具体公式如下：Ψ(M)＝-∑ilogρm(PiM)(4)其中，ρm(·)是提前被只有摩尔纹的块训练出来，m则指示该参数对应的是摩尔纹图层的参数；因为基于块的高斯混合模型先验的存在，所以(2)式是非凸的，通过半二次分裂，得到如下优化方程：其中D,和分别是PiB和PiM的替换变量；3)局部拉普拉斯滤波器：首先建立一个高斯金字塔，并对每一个尺度进行点对点的操作从而产生部分频段的联合拉普拉斯金字塔，对于高斯金字塔的某一特定尺度l的每一个系数η，子区域Ω被提取出来之后，该区域内的每一个像素对应的摩尔纹都会被抑制，p对应Ω内的某一像素，其被修正之后得到新的像素值p，p由如下公式得到：其中是p关于η的正规化变化，参数σm是用来区分图像细节和图像边缘，参数σd进一步区分摩尔纹和图像细节，(△)γ是一个平滑方程，γ≥0是一个自定义的参数，当它的值大于1实现平滑，当它的值小于1时实现增强功能，当|p-η|<σm时，图像像素被认为是摩尔纹，实现平滑功能；当σm<|p-η|<σd时，图像像素被认为是图像细节，实现增强功能；修改后的子区域转化到一个拉普拉斯金字塔里，η所对应的每一个系数都被复制进同尺度同位置的联合金字塔里，当所有的系数都被处理完之后，产生的金字塔恢复出想要的结果，也就是将摩尔纹进一步淡化的图像。一个实例中的具体步骤如下：1)构造初始数据：11)采用拍摄屏幕的图像作为实验数据；12)拍摄大量白色背景的屏幕，并用其训练摩尔纹图层的先验，将图像分割成8*8的块，将这些块以向量形式放在一个矩阵里进行训练；同时收集大量自然图像，并用其训练自然图像背景图层，将图像分割成8*8的块，将这些块以向量形式放在一个矩阵里进行训练；2)首先对输入图像I的Y通道使用多相位分解模型LDPC，接下来介绍一下该模型，其主要分别两个步骤：21)多相位分解模型：首先将输入的图像通过二倍下采样的方式，分为四个多相位部分，也就是四个子采样图像。这些子采样图像不仅打乱了摩尔纹的结果特征，同时保留了自然图像的结构特征，接下来，对每个子采样图像应用图层分解模型实现摩尔纹的消除；22)基于块先验的图层分解模型屏摄图像I可以看成背景图层B和摩尔纹图层M以及方差为σ2的高斯噪声n的叠加，即得加性模型：I＝B+M+n(1)采用两个先验分别对摩尔纹图层和背景图层进行规则化，公式如下：其中，为保真项，Φ(B)和Ψ(M)定义了两个先验以及约束项对进行规则化。定义Φ(B)：Φ(B)对应的公式如下：其中，Pi一个运算符，表示提取第i个块并将其变为向量的形式，块的大小可以自定义，K表示高斯部分的数量，πk表示混合权重，分别表示对应第k个部分的均值和协方差矩阵，N(·)表示64维的高斯分布，b则指示该参数对应的是背景图层的参数，因为自然图像的梯度较为稀疏，因此用于自然图像的边缘平滑，表示梯度运算符，||·||1表示1范数,α设置为0.025；定义Ψ(M)：具体公式如下：Ψ(M)＝-∑ilogρm(PiM)(4)其中，ρm(·)是提前被只有摩尔纹的块训练出来，m则指示该参数对应的是摩尔纹图层的参数；因为基于块的高斯混合模型先验的存在，所以(2)式是非凸的，通过半二次分裂，得到如下优化方程：其中D,和分别是PiB和PiM的替换变量，β是一个正参数初始值设为200，每一次迭代后都将变为原来的2倍，σ＝5×10-3；3)局部拉普拉斯滤波器：起作用是平滑摩尔纹且增强图像细节，首先建立一个高斯金字塔，并对每一个尺度进行点对点的操作从而产生部分频段的联合拉普拉斯金字塔，对于高斯金字塔的某一特定尺度l的每一个系数η，子区域Ω被提取出来之后，该区域内的每一个像素对应的摩尔纹都会被抑制，p对应Ω内的某一像素，其被修正之后得到新的像素值p′，p′由如下公式本文档来自技高网...

【技术保护点】
一种基于多通道分解的屏摄图像摩尔纹消除方法，其特征是，给出多相位图层分解的模型，将该模型用在Y通道，实现摩尔纹彩色条纹和高频部分的淡化；将得到的结果输入一个改进的拉普拉斯滤波器中，实现彩色条纹的进一步消除；最后对于上一步得到的结果在其RGB通道内分别应用多相位分解图层模型，得到最终结果。

【技术特征摘要】
1.一种基于多通道分解的屏摄图像摩尔纹消除方法，其特征是，给出多相位图层分解的模型，将该模型用在Y通道，实现摩尔纹彩色条纹和高频部分的淡化；将得到的结果输入一个改进的拉普拉斯滤波器中，实现彩色条纹的进一步消除；最后对于上一步得到的结果在其RGB通道内分别应用多相位分解图层模型，得到最终结果。2.如权利要求1所述的基于多通道分解的屏摄图像摩尔纹消除方法，其特征是，具体步骤如下：1)多相位分解模型：首先将输入的图像通过二倍下采样的方式，分为四个多相位部分，也就是四个子采样图像；2)基于块先验的图层分解模型屏摄图像I为背景图层B和摩尔纹图层M以及方差为σ2的高斯噪声n的叠加，即得加性模型：I＝B+M+n(1)采用两个先验分别对摩尔纹图层和背景图层进行规则化，公式如下：其中，为保真项，Φ(B)和Ψ(M)定义了两个先验以及约束项对进行规则化；定义Φ(B)：Φ(B)对应的公式如下：其中，Pi为一个运算符，表示提取第i个块并将其变为向量的形式，块的大小自定义，K表示高斯部分的数量，πk表示混合权重，分别表示对应第k个部分的均值和协方差矩阵，N(·)表示64维的高斯分布，b则指示该参数对应的是背景图层的参数，因为自然图像的梯度较为稀疏，因此用于自然图像的边缘平滑，表示梯度运算符，||·||1表示1范数,α是可自定义的参数；定义Ψ(M)：具体公式如下：Ψ(M)＝-∑ilogρm(PiM)(4)其中，ρm(·)是提前被只有摩尔纹的块训练出来，m则指示该参数对应的是摩尔纹图层的参数；因为基于块的高斯混合模型先验的存在，所以(2)式是非凸的，通过半二次分裂，得到如下优化方程：其中D,和分别是PiB和PiM的替换变量；3)局部拉普拉斯滤波器：首先建立一个高斯金字塔，并对每一个尺度进行点对点的操作从而产生部分频段的联合拉普拉斯金字塔，对于高斯金字塔的某一特定尺度l的每一个系数η，子区域Ω被提取出来之后，该区域内的每一个像素对应的摩尔纹都会被抑制，p对应Ω内的某一像素，其被修正之后得到新的像素值p'，p'由如下公式得到：其中是p关于η的正规化变化，参数σm是用来区分图像细节和图像边缘，参数σd进一步区分摩尔纹和图像细节，(△)γ是一个平滑方程，γ≥0是一个自定义的参数，当它的值大于1实现平滑，当它的值小于1时实现增强功能，当|p-η|<σm时，图像像素被认为是摩尔纹，实现平滑功能；当σm<|p-η|<σd时，图像像素被认为是图像细节，实现增强功能；修改后的子区域转化到一个拉普拉斯金字塔里，η所对应的每一个系数都被复制进同尺度同位置的联合金字塔里，当所有的系数都被处理完之后，产生的金字塔恢复出想要的结果，也就是将摩尔纹进一步淡化的图像。3.如权利要求1所述的基于多通道分解的屏摄图像摩尔纹消除方法，其特征是，一个实例中的具体步骤如下：1)构造初始数据：11)采用拍摄屏幕的图像作为实验数据；12)拍摄大量白色背景的屏幕，并用其训练摩尔纹图层的先验，将图像分割成8*8的块，将这些块以向量形式放在一个矩阵里进行训练；同时收集大量自然图像，并用其训练自然图像背景图层，将图像分割成8*8的块，将这些块以向量形式放在一个矩阵里进行训练；2)首先对输入图像I的Y通道使用多相位分解模型LDPC，接下来介绍一下该模型，其主要分别两个步骤：21)多相位分解模型：首先将输入的图像通过二倍下采样的方式，分为四个多相位部分，也就是四个子采样图像。这些子采样图像不仅打乱了摩尔纹的结果特征，同时保留了自然图像的结构特征，接下来，对每个子采样图像应用图层分解模型实现摩尔纹的消除；22)基于块先验的图层分解模型屏摄图像I可以看成背景图层B和摩尔纹图层M以及方差为σ2的高斯噪声n的叠加，即得加性模型：I＝B+M+n(1)采用两个先验分别对摩尔纹图层和背景图层进行规则化，公式如下：...

【专利技术属性】
技术研发人员：杨敬钰，张雪，蔡常瑞，李坤，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津,12