一种红外图像双向直方图均衡处理方法及系统技术方案

本发明专利技术提供的一种红外图像双向直方图均衡处理方法及系统，通过处理所输入的红外图像形成直方图，进而计算高、低边界，通过预设平台值及所述高低边界计算图像灰度重心位置，再将由所述重心位置划分的暗部及亮部分别进行均衡处理，如此，计算方法适合适应性更强，在处理大背景小目标细节的图像时具有明显优势，有效降低数据运算量及设计成本，同时，从本发明专利技术技术方案求得的图像重心分析，对红外图像亮部目标细节体现将更加明显，较之经更加符合人眼观察习惯，并且相对现有平台直方图均衡相比较，本发明专利技术技术方案处理后的亮部目标细节不至于太过夸张，有利于红外图像暗部其他细节的体现。

【技术实现步骤摘要】
一种红外图像双向直方图均衡处理方法及系统
：本专利技术涉及红外图像处理，具体来说，涉及红外图像直方图均衡处理技术。
技术介绍
：红外热像仪获取的大动态范围内的红外图像普遍存在对比度低和细节表现不足的问题，普通的线性拉伸或分段线性拉伸方法对于图像背景和图像细节一视同仁，在提高了图像对比度的同时也提高了背景图案对显示灰度级的占有，而真正目标细节只占有了少量的显示灰度级，造成目标细节分辨力不足。于是在红外热像仪图像处理时，直方图均衡技术特别是平台直方图技术得到了广泛应用。直方图均衡技术一大特点是原始图像的灰度重心必然映射到显示灰度中心，这在处理大背景图像情况时直方图均衡的缺点非常明显，使得目标细节非常聚集，甚至不如线性拉伸方法。于是对平台直方图均衡的研究和应用越来越多，平台直方图的目的就是要抑制大背景所占有的灰度比重，提高目标细节的比重，然而平台直方图对图像重心把握有很大偏差，造成处理后的图像目标细节部分过于夸张，图像其余部分灰度级过于聚集。自适应平台直方图均衡是解决此问题的一种方法，其实施方案是，对于不同图像输入灰度，不采用统一的平台值，而是根据图像直方图分布的普遍共性，多采用接近于“V”型平台值，提高分布在两端的平台值，降低分布在中心的平台值。然而，虽然自适应平台直方图均衡技术引进加权函数明显大大提高了运算量，但也不太适合红外热像仪这种在大动态范围内无固定图像灰度基准的情况，针对红外热像仪获取的不同场景，加权函数要不断地做伸缩处理，算法较为复杂。
技术实现思路
：本专利技术的目的之一在于提供红外图像处理效果较好、数据运算量大大下降、节约了设计成本的技术方案。本专利技术的目的之一在于解决平台直方图均衡在处理红外图像时图像灰度中心难以确定，导致图像灰度重心偏移造成图像视觉感光不协调以及目标细节不明显的问题。为此，本专利技术提供一种红外图像双向直方图均衡处理方法，包括：通过处理所输入的红外图像形成直方图，根据所述红外图像比特位数b来获取对应的像素灰度级值区间记为n、及与所述区间内各像素灰度级值对应的各直方图统计值记为HIST_STAT（n）；按所述像素灰度级值由小到大累加所对应的各直方图统计值直至大于阈值A时，将当前直方图统计值所对应像素灰度级值作为低边界；按所述像素灰度值由大到小累加所对应的各直方图统计值直至大于所述阈值A时，将当前直方图统计值所对应像素灰度级值作为高边界；对所述高、低边界之间的像素灰度级值对应的直方图统计值进行平台第一累加计算，以获取平台第一累加计算结果记为HIST_T，所述平台第一累加计算为将大于预设平台值T的所述直方图统计值替换为T进行累加；对所述高低边界之间的各像素灰度级值对应直方图统计值作平台第一累加计算，直至大于所述第一累加计算结果HIST_T的1/2时，以当前直方图统计值对应的像素灰度级值为图像重心记为HIST_M；分别对小于所述像重心HIST_M的各像素灰度级值、大于所述图像重心HIST_M的各像素灰度级值作均衡处理。所述分别对小于所述像重心HIST_M的各像素灰度级值、大于所述图像重心HIST_M的各像素灰度级值作均衡处理包括：（a）对小于图像重心HIST_M的各像素灰度级值作平台第二累加计算，所述平台第二累加计算公式为：得到平台第二累加结果，记为HIST_ADD(m)，m∈[0,HIST_M)，并且获取小于图像重心HIST_M区域的平台统计总和HIST_SUM_LOW=HIST_ADD（0）；（b）对大于图像重心HIST_M的各像素灰度级值作平台第三累加计算，所述平台第三累加计算公式为：得到平台第三累加结果，记为HIST_ADD(m)，m∈(HIST_M,2b-1]，并且获取大于图像重心各的平台统计总和HIST_SUM_HIGH=HIST_ADD(2b-1)；（c）根据公式：VOUT=128-128*HIST_ADD(m)/HIST_SUM_LOW,m∈[0,HIST_M)的计算结果，对小于图像重心HIST_M的各像素灰度级值对应像素点数量作均衡处理；（d）根据公式：VOUT=128+128*HIST_ADD(m)/HIST_SUM_HIGH,m∈(HIST_M,2b-1]的计算结果，对大于图像重心HIST_M的各像素灰度级值对应像素点数量作均衡处理。所述步骤（c）和（d）前还包括将图像重心HIST_M映射到所述直方图显示的步骤。所述直方图统计值为所述各像素灰度级值的像素点数量。进一步的，所述阈值A为所输入的红外图像的像素点总数的1/4。本专利技术还提供一种红外图像双向直方图均衡处理系统，包括：输入直方图统计单元，用于通过处理所输入的红外图像形成直方图，根据所述红外图像比特位数b来获取对应的像素灰度级值区间记为n、及与所述区间内各像素灰度级值对应的各直方图统计值记为HIST_STAT（n）；高低边界获取单元，用于按所述像素灰度级值由小到大累加所对应的各所述直方图统计值直至大于阈值A时，将当前直方图统计值所对应像素灰度级值作为低边界，以及用于按所述像素灰度值由大到小累加所对应的各所述直方图统计值直至大于所述阈值A时，将当前直方图统计值所对应像素灰度级值作为高边界；图像重心获取单元，用于对所述高低边界之间的各像素灰度级值对应的各所述直方图统计值进行平台累加，以获取总平台统计累加值HIST_T，其中，所述平台累加指的是将大于预设平台值T的所述直方图统计值替换为T进行累加，对所述高低边界之间的各像素灰度级值对应直方图统计值作所述平台累加，直至大于所述第一累加计算结果HIST_T的1/2时，以当前直方图统计值对应的像素灰度级值为图像重心HIST_M；双向均衡处理单元，用于分别对小于所述像重心HIST_M的各像素灰度级值、大于所述图像重心HIST_M的各像素灰度级值作均衡处理。所述双向均衡处理单元包括双向统计累加单元、双向均衡输出单元，所述双向均衡处理单元用于所述分别对小于所述像重心HIST_M的各像素灰度级值、大于所述图像重心HIST_M的各像素灰度级值作均衡处理指的是：所述双向统计累加单元，用于：对小于图像重心HIST_M的各像素灰度级值作平台第二累加计算，所述平台第二累加计算公式为：得到平台第二累加结果，记为HIST_ADD(m)，m∈[0,HIST_M)，并且获取小于图像重心HIST_M区域的平台统计总和HIST_SUM_LOW=HIST_ADD（0）；对大于图像重心HIST_M的各像素灰度级值作平台第三累加计算，所述平台第三累加计算公式为：得到平台第三累加结果，记为HIST_ADD(m)，m∈(HIST_M,2b-1]，并且获取大于图像重心各的平台统计总和HIST_SUM_HIGH=HIST_ADD(2b-1)；所述双向均衡输出单元，用于根据公式：VOUT=128-128*HIST_ADD(m)/HIST_SUM_LOW,m∈[0,HIST_M)的计算结果，对小于图像重心HIST_M的各像素灰度级值对应像素点数量作均衡处理；根据公式：VOUT=128+128*HIST_ADD(m)/HIST_SUM_HIGH,m∈(HIST_M,2b-1]的计算结果，对大于图像重心HIST_M的各像素灰度级值对应像素点数量作均衡处理。所述双向均衡处理单元还用于将图像重心HIST_M映射到所述直方图显本文档来自技高网...

【技术保护点】
一种红外图像双向直方图均衡处理方法，其特征在于，包括：通过处理所输入的红外图像形成直方图，根据所述红外图像比特位数b来获取对应的像素灰度级值区间记为n、及与所述区间内各像素灰度级值对应的各直方图统计值记为HIST_STAT（n）；按所述像素灰度级值由小到大累加所对应的各直方图统计值直至大于阈值A时，将当前直方图统计值所对应像素灰度级值作为低边界；按所述像素灰度级值由大到小累加所对应的各直方图统计值直至大于所述阈值A时，将当前直方图统计值所对应像素灰度级值作为高边界；对所述高、低边界之间的像素灰度级值对应的直方图统计值进行平台第一累加计算，以获取平台第一累加计算结果记为HIST_T，所述平台第一累加计算为将大于预设平台值T的所述直方图统计值替换为T进行累加；对所述高低边界之间的各像素灰度级值对应直方图统计值作平台第一累加计算，直至大于所述第一累加计算结果HIST_T的1/2时，以当前直方图统计值对应的像素灰度级值为图像重心记为HIST_M；分别对小于所述图像重心HIST_M的各像素灰度级值、大于所述图像重心HIST_M的各像素灰度级值作均衡处理。

【技术特征摘要】
1.一种红外图像双向直方图均衡处理方法，其特征在于，包括：通过处理所输入的红外图像形成直方图，根据所述红外图像比特位数b来获取对应的像素灰度级值区间记为n、及与所述区间内各像素灰度级值对应的各直方图统计值记为HIST_STAT(n)；按所述像素灰度级值由小到大累加所对应的各直方图统计值直至大于阈值A时，将当前直方图统计值所对应像素灰度级值作为低边界；按所述像素灰度级值由大到小累加所对应的各直方图统计值直至大于所述阈值A时，将当前直方图统计值所对应像素灰度级值作为高边界；对所述高低边界之间的像素灰度级值对应的直方图统计值进行平台第一累加计算，以获取平台第一累加计算结果记为HIST_T，所述平台第一累加计算为将大于预设平台值T的所述直方图统计值替换为T进行累加；对所述高低边界之间的各像素灰度级值对应直方图统计值作平台第一累加计算，直至大于所述第一累加计算结果HIST_T的1/2时，以当前直方图统计值对应的像素灰度级值为图像重心记为HIST_M；分别对小于所述图像重心HIST_M的各像素灰度级值、大于所述图像重心HIST_M的各像素灰度级值作均衡处理：(a)对小于图像重心HIST_M的各像素灰度级值作平台第二累加计算，所述平台第二累加计算公式为：得到平台第二累加结果，记为HIST_ADD(m)，m∈[0,HIST_M)，并且获取小于图像重心区域的平台统计总和值HIST_SUM_LOW＝HIST_ADD(0)；(b)对大于图像重心HIST_M的各像素灰度级值作平台第三累加计算，所述平台第三累加计算公式为：得到平台第三累加结果，记为HIST_ADD(m)，m∈(HIST_M,2b-1]，并且获取大于图像重心区域的平台统计总和值HIST_SUM_HIGH＝HIST_ADD(2b-1)；(c)根据公式：VOUT＝128-128*HIST_ADD(m)/HIST_SUM_LOW,m∈[0,HIST_M)的计算结果，对小于图像重心HIST_M的各像素灰度级值对应像素点数量作均衡处理；(d)根据公式：VOUT＝128+128*HIST_ADD(m)/HIST_SUM_HIGH,m∈(HIST_M,2b-1]的计算结果，对大于图像重心HIST_M的各像素灰度级值对应像素点数量作均衡处理。2.如权利要求1所述的红外图像双向直方图均衡处理方法，其特征在于，在所述步骤(c)和(d)前还包括将图像重心HIST_M映射到所述直方图显示的步骤。3.如权利要求1至2中任意一项所述的红外图像双向直方图均衡处理方法，其特征在于，所述直方图统计值为所述各像素灰度级值的像素点数量。4.如权利要求3所述的红外图像双向直方图均衡处理方法，其特征在于，所述阈值A为所输入的红外图像的像素点总数的1/4。5.一种红外图像双向直方图均衡处理系统，其特征在于，包括：输入直方图统计单元，用于通过处理所输入的红外图像形成直方图，根据所述红外图像比特位数b来获取...

【专利技术属性】
技术研发人员：梅平，吕超英，龚亮，
申请(专利权)人：无锡华润矽科微电子有限公司，
类型：发明
国别省市：江苏;32