基于拉普阿斯金字塔重构的Retinex图像增强算法的制作方法

基于拉普阿斯金字塔重构的retinex图像增强算法
技术领域：
：1.本发明涉及图像增强
技术领域：
：，包括retinex算法，傅里叶变换，高斯高通滤波，锐化算法，小波变换。
背景技术：
：：2.随着社会的进步和科技的发展，数字图像处理技术已经逐步应用在其他学科领域并且在人们的日常生活中起着重要性的作用，但在恶劣天气影响的情况下所采集的图像会出现对比度低、细节丢失等因素，严重影响信息的传达，因此对图像增强的研究有着重要的意义。3.当前，图像增强可以分为空间域算法和频率域算法两大类；4.而当下传统的图像增强方法有：线性变换法、直方图均衡法、单尺度retinex算法、多尺度retinex算法。空间域算法是直接对图像进行运算以达到增强效果，如图像平滑、图像锐化等。频率域算法是将空间域图像变换为频率域图像，在频率域对图像进行处理以达到增强效果，如高通滤波、低通滤波等。这些算法虽然能解决一些图像增强问题，但往往只能适用于特定的情况，没用良好的适应性。5.图像融合是计算机视觉中的重要分支。通过图像融合，可以将同一场景的不同图像组合成单一融合图像。融合后的图像能提供更好的符合人们客观要求的信息，对人与机器的感知都更有帮助。例如，融合多光谱遥感图像[可以提高特征提取算法的性能。多曝光图像的融合可用于摄影。往往好的融合图像具有以下特点。首先，它可以保留不同图像的大部分有用信息。第二，它不产生工件。第三，该算法对误配和噪声等不完善条件具有较强的鲁棒性。[0006]现有的图像融合存在以下问题：[0007]图像失真，融合算法不能从两或多幅图像中提取有效信息，反而引入噪声。[0008]受低质量图像影响较大，当被融合图像中存在低质量图像时，融合效果会受低质量图像影响，使图像质量降低。[0009]块效应，由于视频及图像采用基于块的编码方式和量化造成相邻块之间存在明显差异的现象，在视频编码中人眼察觉到的小块边界处的不连续。技术实现要素：[0010]本发明要解决的技术问题是提供提出一种基于自适应retinex和小波融合的图像增强算法，该算法可以针对雾天，关照不均匀，低光照等环境恶劣图像进行边界压缩以及锐化操作，利用小波函数分解重构为一幅高信息熵与梯度的图像。[0011]为解决上述技术问题，本发明采用如下技术手段：[0012]基于自适应retinex和小波融合的图像增强算法，包括以下步骤：[0013](1)开始，输入原始图像s(i，j)，利用多尺度retinex估计照度图像用其得到反射图像，其表达式为：[0014][0015]其中si(x，y)为原始图像，ii(x，y)为光照图像，s′(x，y)为反射图像，n为执行单尺度retinex次数，ωk为每个尺度的权重，fk(x，y)为对原图像执行的模糊操作；[0016](2)利用步骤(1)中的估计照度图像求得表达物体本质的反射图像s′(x，y)，其表达式为：[0017]s′(x，y)＝logs(x，y)‑log[f(x，y)*s(x，y)]ꢀꢀꢀꢀ(2)[0018]其中，*表示卷积运算，f(x，y)为高斯核函数，即：[0019][0020]其中，e为常量且e≈2.71828，σ为高斯环绕尺度，ρ为归一化系数，需要满足：[0021]∫∫f(x，y)dxdy＝1ꢀꢀꢀ(4)[0022](3)对于边缘区域的像素点，往往会使图像过亮或者过暗，需要对图像像素边界进行压缩；根据对数域图像分布规律，亮暗点区域像素往往小于0点数量，故，以0对应的像素数量个数为参考值，对数值数量分布分别设置两个边界值分别为low与high，压缩对数域的数值分布；设logr(x，y)中的数值有x＝{x1，x2，x3……xn}，其对应的数量为y＝{y1，y2，y3......yn}，其中xi＝0，yi为0的数量，x1为最小值，xn为最大值，i为数值大小的总数；[0023](4)计算因子s1，s2；对于暗区域，当有x1≤xa≤xi(其中xa∈x)，且ya＝yi*s2时，选取low＝xa；对于亮区域，当xi≤xb≤xn(其中xb∈x)，且ya＝yi*s2时，选取high＝xb；以信息熵为s1，s2的选取标准；s1，s2在0.05到0.15选取能得到符合主观与客观评价标准的图像；[0024](5)将logr(x，y)中小于low的数值赋值为low；大于high的数值赋值为high；根据上述步骤进行处理，提高量化后的图像增强效果；[0025](6)量化后对数据进行线性拉伸使图片像素点转化为0‑255，得到自适应图像。线性拉伸公式如下：[0026]r(x，y)＝(value‑min)/(max‑min)*(255‑0)ꢀꢀꢀꢀꢀ(5)[0027](7)将(6)中得到的自适应图像进行中心化处理，其表达式为：[0028]g(x，y)＝f(x，y)*(‑1)x yꢀꢀꢀ(6)[0029]其中f(x，y)为反射图像，g(x，y)为中心化图像；[0030](8)将(7)中得到的中心图像进行傅里叶变换，其表达式为：[0031][0032]其中图像长m，高n，f(x，y)表示始于图像u的范围为[0，m‑1]，v的范围为[0，n‑1][0033](9)采用高斯高通滤波器对图像进行滤波，其公式为：[0034][0035]其中，h(u，v)表示输出图像，d0表示理想高通滤波截止频率，d(u，v)表示频率点(u，v)到频率中心的距离；[0036](10)对高通滤波后的二维图像用傅里叶逆变换得到频率域增强图像，其表达式为：[0037][0038](11)给出一种n*n的边缘锐化算子，其中n为奇数，且n大于等于5；公式表达为：[0039]k0＝(‑1)×20[0040]ki＝1×2i[0041]kmax＝2×2max(max＝(n‑1)/2)[0042]其中k0表示算子最外层，ki从外层到内层数的第j‑1层，kmax表示中心区域；n为5时，有k0＝‑1，k1＝2，kmax＝8。n＝5时，下表中算子权值之和为8，对掩模的结果除以8即为空间域锐化图像；[0043]k0k0k0k0k0k0k1k1k1k0k0k1kmaxk1k0k0k1k1k1k0k0k0k0k0k0[0044](12)将将步骤(6)中得到的自适应图像，步骤(10)中得到的频率域增强图像，步骤(11)中得到的空间域增强图像分别进行二维小波分解，其表达式为：[0045][0046][0047]其中，wφ(j0，m，n)为分解后低频系数，为分解后高频系数，表示为h、v、d，j0是一个任意开始尺度，通常令其为0，f(x，y)为离散函数，x和y表示离散变量；m，n表示为相对的相对偏移量，m与n代表图像由m×n个像素组成，φj0，m，n(x，y)为二维尺度函数，为i所对应的水平、垂直、对角方向的小波函数；[0048](13)对低频，高频分别进行加权融合，权值方程表达为：[0049]w1 w2 w3＝1[0050](12)[0051]其中wi为图片所占的权重，wi范围为0到1增长步长为0.1，因此可得w3＝1‑w2‑w1，即第三幅图片权值可由可由前两幅图片决定；所以可以使w1＞w2不断便利求出信息熵最高的权重比，再令w2＞w1不断遍历同样求出信息熵最高的权重比即遍历全部参数可求出最优结果；(14)图像融合公式为：[0052]j(x，y)＝w0×d0(i，j) w1×d1(i，j) w2×d2(i，j)[0053](13)[0054]其中d0(i，j)，d1(i，j)，d2(i，j)分别表示步骤(6)中的自适应图像，步骤(10)中的频率域增强图像，步骤(11)中的空间域增强图像中对应小波域的分量，j(x，y)表示融合后的分量；[0055](15)对步骤(13)中各个频率分量进行重构得到最终图像，二维离散小波重构公式如下：[0056][0057]其中图像长m，高n，f(x，y)表示始于图像u的范围为[0，m‑1]，v的范围为[0，n‑1][0058](16)重构后的图像即为最终图像，输出结果图像，结束。附图说明[0059]图1为本发明自适应retinex算法效果图；[0060]图2为本发明的算法流程图；[0061]图3为本发明实施例中5种算法对光照不均匀图像gril增强的结果图；其中，图(a)为原图，图(b)为经过bssr算法处理后的效果图，图(c)为经过bift算法处理后的效果图，图(d)为经过rvrg处理后的效果图，图(e)为经过bfr处理后的效果图，图(f)为本发明算法处理后的效果图；[0062]图4为本发明实施例中5种算法对光照不均匀图像boat增强的结果图；其中，图(a)为原图，图(b)为经过bssr算法处理后的效果图，图(c)为经过bift算法处理后的效果图，图(d)为经过rvrg处理后的效果图，图(e)为经过bfr处理后的效果图，图(f)为本发明算法处理后的效果图；[0063]图5为本发明实施例中5种算法对低照度图像da增强的结果图；其中，图(a)为原图，图(b)为经过bssr算法处理后的效果图，图(c)为经过he算法处理后的效果图，图(d)为经过ear处理后的效果图，图(e)为bfr算法处理后的效果图，图(f)为本发明算法处理后的效果图；[0064]图6为本发明实施例中5种算法对低照度图像ca增强的结果图；其中，图(a)为原图，图(b)为经过bssr算法处理后的效果图，图(c)为经过he算法处理后的效果图，图(d)为经过ear处理后的效果图，图(e)为bfr算法处理后的效果图，图(f)为本发明算法处理后的效果图；具体实施方式[0065]下面结合实施例，进一步说明本发明算法。[0066]结合图2，本发明基于自适应retinex和小波融合的图像增强算法，包括以下步骤：[0067]基于自适应retinex和小波融合的图像增强算法，包括以下步骤：[0068](1)开始，输入原始图像s(i，j)，利用多尺度retinex估计照度图像，其表达式为：[0069][0070]其中si(x，y)为原始图像，ii(x，y)为光照图像，s′(x，y)为反射图像，n为执行单尺度retinex次数，wk为每个尺度的权重，fk(x，y)为对原图像执行的模糊操作；[0071](2)利用步骤(1)中的估计照度图像求得表达物体本质的反射图像s′(x，y)，其表达式为：[0072]s′(x，y)＝logs(x，y)‑log[f(x，y)*s(x，y)]ꢀꢀꢀ(2)[0073]其中，*表示卷积运算，f(x，y)为高斯核函数，即：[0074][0075]其中，e为常量且e≈2.71828，σ为高斯环绕尺度，ρ为归一化系数，需要满足：[0076]∫∫f(x，y)dxdy＝1ꢀꢀꢀ(4)[0077](3)对于边缘区域的像素点，往往会使图像过亮或者过暗，需要对图像像素边界进行压缩；根据对数域图像分布规律，亮暗点区域像素往往小于0点数量，故，以0对应的像素数量个数为参考值，对数值数量分布分别设置两个边界值为low与high，压缩对数域的数值分布；设logr(x，y)中的数值有x＝{x1，x2，x3……xn}，其对应的数量为y＝{y1，y2，y3……yn}，其中xi＝0，yi为0的数量，x1为最小值，xn为最大值，i为数值大小的总数；[0078](4)计算因子s1，s2；对于暗区域，当有x1≤xa≤xi(其中xa∈x)，且ya＝yi*s2时，选取low＝xa；对于亮区域，当xi≤xb≤xn(其中xb∈x)，且ya＝yi*s2时，选取high＝xb；以信息熵为s1，s2的选取标准；s1，s2在0.05到0.15选取能得到符合主观与客观评价标准的图像；[0079](5)将logr(x，y)中小于low的数值赋值为low；大于high的数值赋值为high；根据上述步骤进行处理，提高量化后的图像增强效果；[0080](6)量化后对数据进行线性拉伸使图片像素点转化为0‑255，得到自适应图像。线性拉伸公式如下：[0081]r(x，y)＝(value‑min)/(max‑min)*(255‑0)[0082](5)[0083](7)将(6)中得到的自适应图像进行中心化处理，其表达式为：[0084]g(x，y)＝f(x，y)*(‑1)x yꢀꢀꢀꢀ(6)[0085]其中f(x，y)为反射图像，g(x，y)为中心化图像；[0086](8)将(7)中得到的中心图像进行傅里叶变换，其表达式为：[0087][0088]其中图像长m，高n，f(x，y)表示始于图像u的范围为[0，m‑1]，v的范围为[0，n‑1][0089](9)采用高斯高通滤波器对图像进行滤波，其公式为：[0090][0091]其中，h(u，v)表示输出图像，d0表示理想高通滤波截止频率，d(u，v)表示频率点(u，v)到频率中心的距离；[0092](10)对高通滤波后的二维图像用傅里叶逆变换得到频率域增强图像，其表达式为：[0093][0094](11)给出一种n*n的边缘锐化算子，其中n为奇数，且n大于等于5；公式表达为：[0095]k0＝(‑1)×20[0096]ki＝1×2i[0097]kmax＝2×2max(max＝(n‑1)/2)[0098]其中k0表示算子最外层，ki从外层到内层数的第j‑1层，kmax表示中心区域。n为5只是一个示例，有k0＝‑1，k1＝2，kmax＝8。n＝5时，下表中算子权值之和为8，对掩模的结果除以8即为空间域锐化图像；[0099]k0k0k0k0k0k0k1k1k1k0k0k1kmaxk1k0k0k1k1k1k0k0k0k0k0k0[0100](12)将步骤(6)中得到的自适应图像，步骤(10)中得到的频率域增强图像，步骤(11)中得到的空间域增强图像分别进行二维小波分解，其表达式为：[0101][0102][0103]其中，wφ(j0，m，n)为分解后低频系数，为分解后高频系数，表示为h、v、d，j0是一个任意开始尺度，通常令其为0，f(x，y)为离散函数，x和y表示离散变量；m，n表示为相对的相对偏移量，m与n代表图像由m×n个像素组成，φj0，m，n(x，y)为二维尺度函数，为i所对应的水平、垂直、对角方向的小波函数；[0104](13)对低频，高频分别进行加权融合，权值方程表达为：[0105]w1 w2 w3＝1[0106](12)[0107]其中wi为图片所占的权重，wi范围为0到1增长步长为0.1，因此可得w3＝1‑w2‑w1，即第三幅图片权值可由可由前两幅图片决定；所以可以使w1＞w2不断便利求出信息熵最高的权重比，再令w2＞w1不断遍历同样求出信息熵最高的权重比即遍历全部参数可求出最优结果；(14)图像融合公式为：[0108]j(x，y)＝w0×d0(i，j) w1×d1(i，j) w2×d2(i，j)[0109](13)[0110]其中d0(i，j)，d1(i，j)，d2(i，j)分别表示步骤(6)中的自适应图像，步骤(10)中的频率域增强图像，步骤(11)中的空间域增强图像中对应小波域的分量，j(x，y)表示融合后的分量；[0111](15)对步骤(13)中各个频率分量进行重构得到最终图像，二维离散小波重构公式如下：[0112][0113]其中图像长m，高n，f(x，y)表示始于图像u的范围为[0，m‑1]，v的范围为[0，n‑1][0114](16)重构后的图像即为最终图像，输出结果图像，结束。[0115]实验分析：[0116]本实验以光照不均图像、低照度图像分别为增强的目标图像，对本文提出的基于自适应retinex和小波融合的图像增强算法(以下简称为本发明算法)的增强结果分别进行论证，为了证明本发明算法的有效性，分别通过2组光照不均匀图像、2组低照度图像进行增强。为验证上一章提出算法模型的有效性。将本发明算法实验结果与基于双边滤波的retinex算法(bssr)、直方图均衡化方法(he)、常戬等人提出的基于图像融合技术的retinex图像增强算法(下面简称bift)、智宁等人提出的基于retinex变分模型(下面简称rvrg)、李宏宇等人提出的基于retinex的夜间增强算法(下面简称ear)、任营等人提出的改进双边滤波retine×的多聚焦图像融合算法(下面简称bfr)相比较。[0117]参照图1，对于边缘区域的像素点，往往会使图像过亮或者过暗，因此本发明算法为解决这一问题对图像像素边界进行压缩。根据对数域图像分布规律，可以发现亮暗点区域像素往往小于o点数量。可以以0对应的像素数量个数为参考值，对数值数量分布分别设置两个边界值。压缩对数域的边界值。从图(a)，图(b)看出图片整体偏暗，像素集中在0‑50区域内，视觉效果极差，经过本发明的自适应retinex算法增强后像素分布相对均匀，车窗细节清晰可见。[0118]参照图3、图4，选取2组光照不均匀图像girl、boat进行图像增强处理，对比实验为传统经典算法bssr和3组改进算法rvrg、bift、bfr。图像增强效果如图3，图4所示，以增强结果图像的均值、标准差、信息熵和平均梯度来对比和分析。图3和图4的参数对比如表1和表2所示。[0119]参照图4、图5，选取2组低照度图像ca、da进行图像增强处理，对比实验为2组传统经典算法bssr、he和2组改进算法ear、bfr。图像增强效果如图5，图6所示，以增强结果图像的均值、标准差、信息熵和平均梯度来对比和分析。图5和图6的参数对比如表3和表4所示。[0120]gril效果如图3，(a)为原图，图片整体昏暗区域较多，图中女孩和车子难以看见其细节，远方背景细节也较为模糊，图(b)为bssr算法增强效果示意图，可见昏暗部分得到初步处理，车窗和小女孩部分变得清晰，然而远处车辆与树木部分变得泛灰。图(c)和图(d)为bift和rvrg算法增强效果，和图(b)相比对比度不高，车窗反射细节亮度不足。图(e)为bfr算法，本文主要以bfr算法作为参考，以rvrg和bfr算法作为对比。bfr算法信息熵得到了显著提高，图片边缘信息保存良好，女孩与车窗细节清晰可见。本发明算法，对其采取明暗区域最优解及其融合算法的思路进行了改进，剔除明暗区域避免了图片过亮过暗，采用图片锐化的思路进一步提升细节，车窗反射细节更加明显，女孩细节更加突出，bfr算法头发区域呈现过亮效果，本发明算法则相对均匀。提高了信息熵的同时显著提高了图片梯度。对比rvrg算法，图片评价指标得到了全面提升。对比bfr算法，girl图片亮度降低了7.5％，方差降低了0.5％，信息熵提高了0.02，梯度提高了69％。虽然降低了方差和亮度但显著提高了梯度，因此可以说明本发明算法优于bfr算法。表1光照不均图像增强结果对比情况[0121]表1girl图像增强结果数据对比情况[0122]table1comparisoofgirlimageenhancementresultdata[0123][0124]表2boat图像增强结果数据对比情况[0125]table2comparisonofdataofboatimageenhancementresults[0126][0127]表3da图像增强结果数据对比情况[0128]table4datacomparisonofdaimageenhancementresults[0129][0130]表4ca图像增强结果数据对比情况[0131]table4datacomparisonofcaimageenhancementresults[0132][0133]关于boat图片对比rvrg算法，各个指标都得到了提升。对比bfr算法，boat图片亮度降低了7.5％，方差提升了23.0％，信息熵提高了0.05，梯度提高了54.5％。虽然亮度有所降低，但梯度和方差得到了显著提升，信息熵优于bfr算法，本发明算法得到更好的效果。boat的增强效果如图18所示。其中图(a)为原始图像，观察可见图片只有船帆比较清晰，图片上部和左下部以及中间人昏暗明显。由于该区域灰度值低，主观感觉较差，肉眼难以观察其细节。图(b)为bssr算法增强效果示意图，图像整体亮度明显提升，但细节仍然模糊，图片泛灰。图(c)为bift算法增强效果示意图，与图(b)相比图片细节明显增强，抑制了泛灰效果，图片对比度明显增强，但草地，波浪，以及左下角石头部分细节出现摸出；图(d)和图(e)为rvrg算法和bfr增强效果示意图，相对于前两种算法图像整体轮廓和细节增强效果明显，但图(e)视觉效果过亮，以及左下角石头颜色过白，这也能解释实验数据中，图(e)均值偏大的情况。本发明算法能有效避免图像过亮过暗的情况。虽然亮度有所降低，但梯度和方差得到了显著提升，信息熵也优于bfr，可以证明本发明算法得到更好的效果。[0134]对于da图片。图像整体昏暗，右上部区域存在雾气，经bssr算法，ear算法，bfr算法处理后虽然一定程度上缓解了图片昏暗的情况，但雾气效果明显，给图片一种朦胧感。从数值上看对然提高了信息熵，但梯度方面并未显著提高，方差也相对较低。he算法在低照度算法中主观感觉不错，但与本发明算法比较并未突出细节，在远处高楼与下部草丛中，不如本算法清晰。对比ear和bfr算法，图片评价指标得到了全面提升。对比两种算法。da图片亮度分别提升了24.8％和25.4％，方差分别提高了133％和131.6％，信息熵分别提高了14.6％和14.6％，梯度分别提高了228.7％和226.6％。[0135]对于ca图片。对比ear和bfr算法，图片评价指标得到了全面提升。对比两种算法。ca图片亮度分别提升了47.1％和46.5％，方差分别提高了41.9％和41.2％，信息熵分别提高了5.9％和5.8％，梯度分别提高了108.2％和99.6％。图像昏暗区域集中在左部楼房，下部街道，右中部楼房。经bssr算法处理后，图像得到了初步增强。he算法则使图片部分过亮。ear与bfr效果极为相似，评价数值也十分接近。本发明算法与ear与bfr相比，提升了亮度，对于低照度图片来说，应适当提高亮度，而对于高曝光图片来说则应该适当降低亮度。bfr算法中左侧窗户仍呈现昏暗状态，下侧街道较黑，右侧从下数第三个与第四个窗户细节不明显，本发明算法则可以清晰看到细节。[0136]综上所述，本发明对实际应用中常见的光照不均匀的图像、低照度图像均取得了较好的增强效果，可以适用于更多场合。[0137]以上所述仅为本发明较佳可行的实施例而已，并非因此局限本发明的权利范围，凡运用本发明说明书及附图内容所作的等效结构变化，均包含于本发明的权利范围之内。当前第1页12当前第1页12