一种融合梯度和灰度信息的遥感图像增强方法技术

本发明专利技术公开了一种融合梯度和灰度信息的遥感图像增强方法，主要步骤有：1、统计输入遥感图像的梯度灰度联合直方图；2、根据标准差和各灰度级的频数自适应地对梯度灰度联合直方图进行变换得到优化的直方图；3、基于直方图均衡化方法建立灰度映射关系，对输入遥感图像灰度重构生成全局增强遥感图像；4、计算全局增强图和原图的梯度差异，对增强图中梯度下降的区域进补偿。本发明专利技术能够有效克服现有直方图变换的图像结构细节信息的缺失，且对梯度下降区域进行补偿能够避免局部细节的丢失。使用该方法增强后的遥感图像视觉效果好、细节丰富，且整个过程无需人为干预。

【技术实现步骤摘要】
一种融合梯度和灰度信息的遥感图像增强方法
本专利技术涉及图像增强领域，特别是关于遥感图像整体对比度和局部细节的增强方法。
技术介绍
整体对比度和局部细节信息是衡量遥感图像质量的重要信息，也是影响后续遥感应用的重要因素。高质量的遥感图像是地物分类、目标识别、变换检测、地质勘探、环境监测中发挥着不可替代的作用。然而，由于天气影响、设备老化等不确定因素导致遥感图像整体对比度下降、细节信息模糊。因此，对低对比度图像进行增强处理，是遥感应用的前提。图像增强是为了提升图像的对比度和细节信息，在尽可能保持信息量的同时获得更好视觉效果和更高清晰度的图像，目前提出了多种类型的图像增强方法：1)基于频率域的图像增强，如基于DCT的图像增强、基于DWT的图像增强、基于SVD的图像增强等；这类方法在细节增强方面有较好的效果，但往往会使图像平坦区域产生伪细节，造成图像的细节失真，如导致平坦的地面出现波纹状的纹理。2)基于空间域的图像增强，如直方图均衡化、线性拉伸、线性截断拉伸等；这类方法常需要进行灰度合并，会导致增强效果生硬、大区域细节丢失，如导致红外波段图像的纹理信息丢失。3)其他类型的增强算法，如Retinex增强、空间域和频率域结合的图像增强等。这些算法也大都针对自然图像增强，在遥感图像增强中常常失效，无法满足实际应用。在实际遥感图像显示中，我们需要获得整体对比度较好且局部细节保持较好的图像，而现有方法难以满足该要求。
技术实现思路
本专利技术提出了一种融合梯度和灰度信息的遥感图像增强方法，充分考虑梯度信息在遥感图像增强中的重要性，在全局增强中能够克服直方图均衡化、线性拉伸等增强方法的在全局增强中效果不自然的不足，且通过局部细节补偿能够极大程度保持局部细节信息。本专利技术的技术方案如下：该融合梯度和灰度信息的遥感图像增强方法，主要包括以下步骤：步骤a、对输入原始遥感图像统计得出梯度灰度联合直方图；所述梯度灰度联合直方图包含原始遥感图像的灰度信息和梯度信息(即能够描述图像除灰度分布以外的地物结构信息，从而更加适合遥感图像的增强)；步骤b、对梯度灰度联合直方图进行优化：计算灰度梯度联合直方图的标准差，结合梯度灰度联合直方图各灰度级对应的归一化频数，构建梯度灰度联合直方图优化参数，对每个灰度级的直方图频数进行修正，得到优化的梯度灰度联合直方图；步骤c、对优化的灰度梯度联合直方图进行均衡化处理计算，建立原始遥感图像到整体增强遥感图像的灰度映射关系，得到整体增强遥感图像；步骤d、计算所述整体增强遥感图像与原始遥感图像的梯度差异，对整体增强遥感图像中梯度下降的区域进行梯度细节补偿，输出最终增强遥感图像。其中，进行梯度细节补偿的方法是：用整体增强遥感图像的梯度逐像素减去原始遥感图像的梯度，标记小于0的像素；将相连的标记像素按8邻域进行编号得到不同编号的标记区域；并用对应标记区域的原始遥感图像的线性拉伸值对整体增强遥感图像中的像素值进行改变，得到最终增强遥感图像。以上各步骤较佳的具体实现过程分别如下：步骤a具体按照以下公式计算输入的原始遥感图像f(i,j)的梯度灰度联合直方图GIH；G(k)＝sum(G(i,j)),iff(i,j)＝k其中：Gx(i,j)＝f(i+1,j-1)+2×f(i+1,j)+f(i+1,j+1)-f(i-1,j-1)-2×f(i-1,j)-f(i-1,j+1)Gy(i,j)＝f(i-1,j+1)+2×f(i,j+1)+f(i+1,j+1)-f(i-1,j-1)-2×f(i,j-1)-f(i+1,j-1)其中(i,j)为像元水平和垂直方向的位置；k＝0,1,2,…,K-1为图像的灰度级；K＝2B-1；B为输入原始遥感图像的位数；G(k)表示图像在灰度值为k的像素梯度值之和，GIH(k)表G(k)的归一化值，即梯度灰度联合直方图各灰度级的归一化频数。步骤b中所述梯度灰度联合直方图优化参数为各灰度级调节参数；步骤b具体是：首先，根据梯度灰度联合直方图的标准差和各灰度级的频数，计算梯度灰度联合直方图在各灰度级的调节参数T(k)；stdh＝std(GIH)式中，k＝0,1,2,…,K-1为图像的灰度级；K＝2B-1；B为输入原始遥感图像的位数；GIH为输入图像f(i,j)的梯度灰度联合直方图；std(·)为标准差函数；然后，利用调节参数T(k)对梯度灰度联合直方图进行调节，得到优化的梯度灰度联合直方图GIHrefine(k)，并进行归一化；GIHrefine(k)＝GIH(k)T(k)步骤c具体按照以下公式计算累积分布梯度灰度联合直方图Frefine(k)，并计算灰度映射函数将原始遥感图像灰度k变换为yrefine(k)，实现图像整体增强得到整体增强遥感图像fc(i,j)；其中yu，yd表示增强后输出图像的最小值和最大值。步骤d具体可以是：首先，设Pc(i,j)值为0，计算整体增强遥感图像fc(i,j)和原始遥感图像f(i,j)梯度的差异，找出整体增强后梯度下降的像素位置；Pc(i,j)＝1,if(Gc(i,j)-G(i,j)＜0)然后，对梯度下降的相连像素进行标记，得到梯度下降区域，并按照8邻域对梯度下降区进行编号，标记得到N个不同区域P1＝1，P2＝2，...PN＝N；最后，对梯度下降的每个区域的像素进行重新赋值；其中fliner(i,j)表示在线性拉伸增强的图像；tn为整体增强遥感图像fc(i,j)在区域Pn内外边界均值的差值减去原始遥感图像f(i,j)在区域Pn内外边界均值的差值；如此按照标记值从小到大，对整幅图中所有梯度下降区域进行更新得到fc_f(i,j)。本专利技术有如下效果：1、本专利技术中提出的融合梯度和灰度信息(梯度灰度联合直方图变换和梯度细节补偿)原理的图像自适应增强方法，能够有效地提升图像整体的对比度，同时也能够保持局部细节信息，对遥感图像适用性高，增强效果明显，适用于各种极端情况下(薄雾、低照度、过曝光、传感器退化等)遥感图像的对比度增强。2、本专利技术方法涉及的所有参数，都根据图像自身特性进行计算，避免了人为设置参数导致增强结果的不稳定，自动化程度高。3、本专利技术方法兼具图像整体增强和局部增强的优点，增强后的遥感图像视觉效果好、细节丰富，适用性较现有方法好。附图说明图1是本专利技术的遥感图像增强方法流程图。图2是按照8连通区域对梯度下降区进行编号的示意图。图3是一幅火星遥感图像。图4是火星遥感图像的直方图。图5是直接使用直方图均衡化的结果。图6是使用梯度灰度联合直方图均衡化处理的结果。图7是使用本专利技术中优化的梯度灰度联合直方图均衡化处理的结果。图8是优化的梯度灰度联合直方图均衡化处理后梯度细节减弱区域。图9是使用本专利技术方法处理的结果。图10是使用本专利技术方法处理的结果的直方图。具体实施方式以下结合附图对本专利技术的具体实施流程做进一步的说明。受光照、探测器、目标辐射特性等成像条件的影响，获取的遥感图像往往对比度较低、细节特征不明显，难以从中获取有效的信息。图3为一幅全色火星遥感图像，如图4所示，其灰度分布集中在[40,140]，该范围的像素大约有90％；在灰度区间[0,39]和[141,255]大约2/3范围中分布较少，像素个数大约10％左右。这就造成大量的灰度级浪费，导致图像对比度不强，细节纹理难以辨认。本专利技术的融合本文档来自技高网...

【技术保护点】
1.一种融合梯度和灰度信息的遥感图像增强方法，其特征在于，包括以下步骤：步骤a、对输入原始遥感图像统计得出梯度灰度联合直方图；所述梯度灰度联合直方图包含原始遥感图像的灰度信息和梯度信息；步骤b、对梯度灰度联合直方图进行优化：计算灰度梯度联合直方图的标准差，结合梯度灰度联合直方图各灰度级对应的归一化频数，构建梯度灰度联合直方图优化参数，对每个灰度级的直方图频数进行修正，得到优化的梯度灰度联合直方图；步骤c、对优化的灰度梯度联合直方图进行均衡化处理计算，建立原始遥感图像到整体增强遥感图像的灰度映射关系，得到整体增强遥感图像；步骤d、计算所述整体增强遥感图像与原始遥感图像的梯度差异，对整体增强遥感图像中梯度下降的区域进行梯度细节补偿，输出最终增强遥感图像。

【技术特征摘要】
1.一种融合梯度和灰度信息的遥感图像增强方法，其特征在于，包括以下步骤：步骤a、对输入原始遥感图像统计得出梯度灰度联合直方图；所述梯度灰度联合直方图包含原始遥感图像的灰度信息和梯度信息；步骤b、对梯度灰度联合直方图进行优化：计算灰度梯度联合直方图的标准差，结合梯度灰度联合直方图各灰度级对应的归一化频数，构建梯度灰度联合直方图优化参数，对每个灰度级的直方图频数进行修正，得到优化的梯度灰度联合直方图；步骤c、对优化的灰度梯度联合直方图进行均衡化处理计算，建立原始遥感图像到整体增强遥感图像的灰度映射关系，得到整体增强遥感图像；步骤d、计算所述整体增强遥感图像与原始遥感图像的梯度差异，对整体增强遥感图像中梯度下降的区域进行梯度细节补偿，输出最终增强遥感图像。2.根据权利要求1所述的融合梯度和灰度信息的遥感图像增强方法，其特征在于，步骤a具体按照以下公式计算输入的原始遥感图像f(i,j)的梯度灰度联合直方图GIH；G(k)＝sum(G(i,j)),iff(i,j)＝k其中：Gx(i,j)＝f(i+1,j-1)+2×f(i+1,j)+f(i+1,j+1)-f(i-1,j-1)-2×f(i-1,j)-f(i-1,j+1)Gy(i,j)＝f(i-1,j+1)+2×f(i,j+1)+f(i+1,j+1)-f(i-1,j-1)-2×f(i,j-1)-f(i+1,j-1)其中(i,j)为像元水平和垂直方向的位置；k＝0,1,2,…,K-1为图像的灰度级；K＝2B-1；B为输入原始遥感图像的位数；G(k)表示图像在灰度值为k的像素梯度值之和，GIH(k)表G(k)的归一化值，即梯度灰度联合直方图各灰度级的归一化频数。3.根据权利要求1所述的融合梯度和灰度信息的遥感图像增强方法，其特征在于，步骤b中所述梯度灰度联合直方图优化参数为各灰度级调节参数；步骤b具体是：首先，根据梯度灰度联合直方图的标准差和各灰度级的频数，计算梯度灰度联合直方图在各灰度级的调节参数T(k)；stdh＝std(GIH)式中，k＝0,...

【专利技术属性】
技术研发人员：陈铁桥，柳稼航，朱锋，刘佳，张航，王一豪，陈军宇，
申请(专利权)人：中国科学院西安光学精密机械研究所，
类型：发明
国别省市：陕西,61