[发明专利]一种光晕去除的多曝光图像融合方法

摘要

本发明公开了一种光晕去除的多曝光图像融合方法，其计算每幅曝光图像的亮度图像的局部梯度稀疏复杂度映射图和曝光质量映射图；然后融合得到每幅曝光图像的亮度图像的初步融合权值图像，进而得到每幅曝光图像的亮度图像的初步融合权值图像的归一化图像的优化图像的归一化图像；再获得所有曝光图像的最终融合图像；优点是通过针对亮度图像的梯度幅值图像的稀疏表示构建局部梯度稀疏复杂度映射图及曝光质量映射图进行图像融合，有效抑制了大幅值的无效梯度在图像融合过程中产生的影响，消除了融合图像中的光晕现象，且能更有效地提取多曝光图像序列中不同曝光条件下的图像纹理细节，因此融合图像中含有更丰富的纹理细节和具有更高的清晰度。

主权项

1.一种光晕去除的多曝光图像融合方法，其特征在于包括以下步骤：步骤一：选取不同场景的清晰图像共Nclear幅，选取的清晰图像包括室内场景的清晰图像和室外场景的清晰图像；然后获取每幅清晰图像的亮度图像；接着获取每幅清晰图像的亮度图像的梯度幅值图像；其中，清晰图像为RGB格式，Nclear≥2，室内场景的清晰图像的幅数与室外场景的清晰图像的幅数相同；步骤二：通过滑动窗口技术将每幅清晰图像的亮度图像的梯度幅值图像划分成Nblock个部分重叠的尺寸大小为Block×Block的图像块；然后将每幅清晰图像的亮度图像的梯度幅值图像中的每个图像块中的所有像素点的像素值按序排列构成该图像块对应的幅值列向量；接着将每幅清晰图像的亮度图像的梯度幅值图像中的所有图像块对应的幅值列向量构成该幅清晰图像的亮度图像的梯度幅值图像对应的幅值矩阵；再将所有清晰图像的亮度图像的梯度幅值图像对应的幅值矩阵构成总幅值矩阵，记为Sclear；其中，滑动窗口的窗口大小为Block×Block、滑动步长为1个像素点，Nblock＝(W‑Block+1)×(H‑Block+1)，W表示清晰图像的宽度，H表示清晰图像的高度，Block＝2n+1，n为正整数，1≤n≤5，幅值列向量的维数为(Block×Block)×1，幅值矩阵的维数为(Block×Block)×Nblock，Sclear的维数为(Block×Block)×(Nblock×Nclear)；步骤三：采用K‑SVD算法对Sclear进行处理，得到最优的过饱和字典，记为Dbest；其中，Dbest的维数为(Block×Block)×K，K表示Dbest中的原子的总个数；步骤四：选取待融合的同一场景而曝光程度不同的N_ex幅曝光图像，所有曝光图像中各自含有场景中的在不同曝光条件下所能呈现出来的纹理细节，将第n_ex幅曝光图像记为然后扩展每幅曝光图像的尺寸大小，将扩展后的图像作为待处理图像，每幅待处理图像的宽度为W'+2n且高度为H'+2n；接着获取每幅曝光图像的亮度图像和每幅待处理图像的亮度图像，将的亮度图像记为再获取每幅待处理图像的亮度图像的梯度幅值图像；其中，曝光图像为RGB格式，N_ex≥3，n_ex为正整数，1≤n_ex≤N_ex，W'表示曝光图像的宽度，H'表示曝光图像的高度；步骤五：通过滑动窗口技术将每幅待处理图像的亮度图像的梯度幅值图像划分成N'_block个部分重叠的尺寸大小为B_lock×B_lock的图像块；然后将每幅待处理图像的亮度图像的梯度幅值图像中的每个图像块中的所有像素点的像素值按序排列构成该图像块对应的幅值列向量；接着将每幅待处理图像的亮度图像的梯度幅值图像中的所有图像块对应的幅值列向量构成该幅待处理图像的亮度图像的梯度幅值图像对应的幅值矩阵，将第n_ex幅待处理图像的亮度图像的梯度幅值图像对应的幅值矩阵记为其中，滑动窗口的窗口大小为B_lock×B_lock、滑动步长为1个像素点，N'_block＝W'×H'，幅值列向量的维数为(B_lock×B_lock)×1，的维数为(B_lock×B_lock)×N'_block；步骤六：计算每幅待处理图像的亮度图像的梯度幅值图像对应的幅值矩阵的稀疏系数矩阵，将的稀疏系数矩阵记为然后计算每幅曝光图像的亮度图像的局部梯度稀疏复杂度映射图，将的局部梯度稀疏复杂度映射图记为将中坐标位置为(x_ex,y_ex)的像素点的像素值记为设中坐标位置为(x_ex,y_ex)的像素点为中的第i'个像素点，且第n_ex幅待处理图像的亮度图像的梯度幅值图像中的第i'个图像块的中心像素点与中的第i'个像素点对应，则有其中，的维数为K×N'_block，1≤x_ex≤W'，1≤y_ex≤H'，i'为正整数，1≤i'≤N'_block，符号“|| ||₁”为l₁范数求取符号，a_i'表示的第i'列的所有元素构成的向量，亦表示中的第i'个幅值列向量的稀疏系数向量；步骤七：对每幅曝光图像的亮度图像进行归一化处理，得到每幅曝光图像的亮度图像的归一化图像，将的归一化图像记为然后计算每幅曝光图像的亮度图像的曝光质量映射图，将的曝光质量映射图记为将中坐标位置为(x_ex,y_ex)的像素点的像素值记为其中，表示曝光质量好，表示曝光质量差，α表示曝光质量阈值，α∈[0,1]，表示中坐标位置为(x_ex,y_ex)的像素点的像素值；步骤八：计算每幅曝光图像的亮度图像的初步融合权值图像，将的初步融合权值图像记为将中坐标位置为(x_ex,y_ex)的像素点的像素值记为然后对每幅曝光图像的亮度图像的初步融合权值图像进行归一化处理，得到每幅曝光图像的亮度图像的初步融合权值图像的归一化图像，将的归一化图像记为将中坐标位置为(x_ex,y_ex)的像素点的像素值记为接着通过递归滤波技术对每幅曝光图像的亮度图像的初步融合权值图像的归一化图像进行优化处理，得到每幅曝光图像的亮度图像的初步融合权值图像的归一化图像的优化图像，将的优化图像记为将中坐标位置为(x_ex,y_ex)的像素点的像素值记为再对每幅曝光图像的亮度图像的初步融合权值图像的归一化图像的优化图像进行归一化处理，得到每幅曝光图像的亮度图像的初步融合权值图像的归一化图像的优化图像的归一化图像，将的归一化图像记为将中坐标位置为(x_ex,y_ex)的像素点的像素值记为其中，τ是一个很小的正数，F_c(x_ex,y_ex)为引入的中间变量，e表示自然基数，sigma_s表示递归滤波的标准差，sigma_r表示递归滤波的标准偏差范围，表示中坐标位置为(x_ex,y_ex)的像素点的像素值，若1≤y_ex+1≤H'则表示中坐标位置为(x_ex,y_ex+1)的像素点的像素值，若y_ex+1＞H'则令表示中坐标位置为(x_ex,H')的像素点的像素值，若1≤y_ex+1≤H'则表示中坐标位置为(x_ex,y_ex+1)的像素点的像素值，若y_ex+1＞H'则令表示中坐标位置为(x_ex,H')的像素点的像素值，F_b(x_ex,y_ex)为引入的中间变量，若1≤y_ex‑1≤H'则表示中坐标位置为(x_ex,y_ex‑1)的像素点的像素值，若y_ex‑1＜1则令表示中坐标位置为(x_ex,1)的像素点的像素值，F_a(x_ex,y_ex)为引入的中间变量，若1≤x_ex+1≤W'则表示中坐标位置为(x_ex+1,y_ex)的像素点的像素值，若x_ex+1＞W'则令表示中坐标位置为(W',y_ex)的像素点的像素值，若1≤x_ex+1≤W'则表示中坐标位置为(x_ex+1,y_ex)的像素点的像素值，若x_ex+1＞W'则令表示中坐标位置为(W',y_ex)的像素点的像素值，若1≤x_ex‑1≤W'则表示中坐标位置为(x_ex‑1,y_ex)的像素点的像素值，若x_ex‑1＜1则令表示中坐标位置为(1,y_ex)的像素点的像素值；步骤九：计算所有曝光图像的最终融合图像，记为I_F，将I_F的第m个通道中坐标位置为(x_ex,y_ex)的像素点的像素值记为I_F(x_ex,y_ex,m)，其中，I_F为RGB格式，m＝1,2,3，m＝1时表示R通道，m＝2时表示G通道，m＝3时表示B通道，表示的第m个通道中坐标位置为(x_ex,y_ex)的像素点的像素值。