专利名称：一种基于灰度分割的自适应矩匹配条带噪声去除方法
技术领域：
本发明属于遥感影像处理领域，涉及一种基于灰度分割的自适应矩匹配条带噪声去除方法。
背景技术：
条带噪声是影响光学卫星影像成像质量的一个重要因素，抑制或去除条带噪声是卫星地面预处理进行辐射处理的基本环节之一。光学卫星的成像系统由于内外部因素的干扰，如CCD器件(电荷耦合器件)的性能随时间发生变化、大气干扰等，在进行影像均一化相对辐射校正之后，影像中仍然会残留有条带噪声，这类噪声的存在，极大的降低了影像的清晰度，为影像的后续判读处理增加了难度，因此必须予以剔除。通过对影像中残留的条带噪声进行分析发现该类条带噪声存在有以下特点
I.噪声出现的空间位置是随机的。2.噪声与周围地物之间高度非线性相关。目前常用的去条带噪声的方法归纳起来可以分为两类一类是针对影像空间域特征提出的去噪方法；另一类是将空间域和频率域相结合，采用适当的滤波算子去除条带噪声的方法。其中空间域去噪的典型算法有Wegener提出的基于改进直方图匹配的多载荷成像的条带噪声去除法，该方法为了得到期望直方图，需要统计大量的有代表性的卫星影像，工作量巨大，不利于单景影像实时处理，而且影像中条带噪声去除的效果依赖于样本影像的选取，因此该算法的处理效果不稳定。于是，Gadallahd等人提出了基于矩匹配的多载荷影像的条带噪声去除法，该方法要求影像灰度均一，且条带噪声影像与其相邻影像之间存在较高的相关性，这类算法对于处理局部的地物单一的卫星影像具有很好的效果，但是处理整幅地物丰富的卫星影像，效果不佳。在此基础上，刘正军等人提出了基于均值补偿的矩匹配方法，该方法利用高光谱影像的多光谱信息来纠正矩匹配算法产生的灰度畸变，该方法很好的处理了标准矩匹配算法纠正影像时造成的探元间的灰度畸变，但是并没有纠正标准矩匹配算法造成的探元内的灰度畸变，且该方法对于单波段的全色影像并不适用。空间域和频率域相结合的去噪算法目前主要是利用小波变换的时频特点，通过对影像进行小波变换，研究噪声的小波系数变换规律，从而提炼出条带噪声的成分并对其进行剔出，但小波变换计算量大，并且会对光谱信息造成较大损失，对于整幅卫星影像的处理，计算速度和光谱保留度并不理想，同时由于影像中条带噪声空间位置是随机的，因此完全定位该类条带噪声在影像中的位置，可实现性尚存在质疑。
发明内容
本发明所要解决的问题是，针对全色影像中存在的随机条带噪声，提供一种有效的影像噪声去除方法。本发明的技术方案为一种基于灰度分割的自适应矩匹配条带噪声去除方法，包括以下步骤步骤1，获取原始含噪影像的亮度值曲线；
步骤2，通过采用不同大小的处理窗口对步骤I所得亮度值曲线进行中值滤波去除基线，获取处理单元中心；步骤3，根据两个相邻处理单元中心的间距进行调整，确定处理单元；步骤4，对单个处理单元采用标准矩匹配的处理方式进行影像处理。而且，设原始含噪影像的大小为MXN，(i, j)表示原始含噪影像中第i行第j列的像素，其中i的取值范围为0，I…M-l，j的取值范围为0，I…N-1，步骤I中获取原始含噪影像的亮度值曲线是根据如下公式计算得到，
i=M-Lj=N-IS(k) = {艺 //；(/,./) I Ij. {iJ) = k}
I=O^j=O 其中，k表示对应的灰度级，IE(i, j)表示影像中位于原始含噪影像(i，j)处像素的灰度值。而且，步骤2包括如下子步骤，步骤2. 1，设中值滤波的处理窗口大小记为WF，初始化设置Wf的值，设迭代次数Count=O ；步骤2. 2，执行以下计算公式获取去基线的亮度值曲线S(k)f ；S (k) f = S (k) -med (S (k_t) t G ff)其中，W={t I-(Wf-I)/2 ^ t ^ (Wf-I)/2}，med 是指中值滤波处理;步骤2. 3，对去基线的亮度值曲线3&\采用均值y (k)f和方差O &\进行归一化处理，处理公式如下，
…、S(k)f-^(k)fS(k)n =-—-L
O-(Zc)f其中，S(k)^|^是均一化后的结果；步骤2. 4，根据预设的阈值T选择候选处理中心，选择公式如下，
fl I S(k) \>TI\k), =J、
其中，m表示影像量化的位数。而且，步骤3中对处理单元中心进行调整，实现方式为当两个相邻处理单元中心的间距小于预设阈值h时，则将两个相邻处理单元中心合并为一个处理单元中心，合并后的处理单元中心的值为这两个相邻处理单元中心的均值；当两个相邻处理单元中心的间距大于或等于预设阈值h时，两个相邻处理单元中心之间分裂产生新的处理单元中心，分裂后新的处理单元中心的值为原始两相邻处理单元中心的平均值。而且，预设阈值h采用原始含噪影像的随机噪声灰度跃变最大值的4飞倍。本发明只通过全色影像本身的信息来抑制标准矩匹配算法在去噪过程中产生的灰度畸变。该方法能够有效地复原影像、改善影像质量，运行效率也很高。
图I为本发明实施例的流程示意图。
具体实施例方式根据影像处理的基本理论可以知道，影像直方图的不同峰值代表了影像中包含的不同地物类型，标准矩匹配方法虽然在处理地物丰富的影像时会造成灰度畸变，即在影像直方图具有多个峰值或者直方图峰值曲线跨度较大的情况下，采用标准矩匹配方法会造成灰度畸变现象，但是对于灰度范围跨度较小的单个峰值的曲线分布，标准矩匹配方法在对影像进行噪声去除处理时，并不会造成影像中产生灰度畸变现象。因此将单一探元内部的成像数据根据影像直方图分布的差异(即不同的灰度等级)进行区分，然后分别求取不同灰度区域成像数据的方差和均值，最后对不同的灰度区域的成像数据，分别进行标准矩匹配处理，这样就可以解决标准矩匹配由于均值改变导致的探元间的灰度畸变以及由于方差改变导致的探元内的灰度畸变。以下结合附图和实施例详细说明本发明技术方案。基于灰度分割的自适应矩匹配条带噪声去除方法，其原理公式和处理步骤如下所示
M
UM ^t
+产(I)
、^ _w=0m=0 ^mi
/ j mi
m=0式⑴中，x，y分别为影像第i个扫描列各像素校正前、后的灰度值；O ，U 为灰度等级为m时，无噪声参考列的方差和均值；O mi，ii mi为灰度等级为m时第i个扫描列的方差和均值。m的取值为0，I…M，灰度等级共有M+1个。对于每个灰度等级的成像数据，并非所有的探元的数据都需要参与相应的校正计算，只有当卫星载荷对应的探元在扫描成像时包含有对应灰度等级的成像数据，才参与相应的校正计算。本发明技术方案可由本领域技术人员采用计算机软件技术实现自动运行流程，关键在于分割处理，即处理单元中心的自适应获取。实施例规定灰度等级划分的原则如下I)处理单元的左边取本单元中心和左边相邻中心的均值为界限，处理单元的右边取本单元中心和右边相邻中心的均值为界限。2)灰度级跨度范围应小于某一阈值，以免处理单元的灰度跨度较大造成标准矩匹配处理时产生灰度畸变现象；3)依据影像直方图，处理单元中心对应的灰度值所包含的像素个数应该在单元内处最大值，且要大于一个域值以保证处理单元中大部分是地物反射成分而不是噪声成分，从而达到去噪的目的而并非噪声增强处理；4)处理单元的划分应该自动收敛，使得处理无需人工干预。获取处理单元的关键在于准确定位直方图的峰值点，为了提高直方图峰值点获取的准确性，本发明实施例采用多层滤波的方式，并设置限定条件，最后对于灰度跨度范围依据噪声的特点进行限定，从而在去除条带噪声的同时抑制标准矩匹配算法带来的灰度畸变 现象。通过对国产卫星全色影像中残留的条带噪声进行分析发现，随机条带噪声和其四周影像的能量相关性较高，因此随机条带噪声也可以根据成像数据灰度等级区分方式，被分配到不同的灰度等级之中，根据这个原理，本发明实施例的处理步骤如下所示，如图I :步骤1，获取原始含噪影像的亮度值曲线。实施例中，假设一幅影像的大小为MXN，(i, j)表示影像中第i行第j列的像素，其中i的取值范围为0，I…M-l，j的取值范围为0，I…N-1。那么亮度值分别的曲线根据如
下公式计算得到
S(k) = { Yj lE(iJ)\IF{iJ) = k\{2j
z=o,y=o式⑵中，k表示对应的灰度级，IE(i, j)表示影像中位于图像(i，j)处像素的灰度值。亮度值曲线S(k)通常由两部分组成变化较平缓的基线和一些峰值，这些变换较平缓的基线说明了影像包含地物的丰富程度；而一些峰值则说明了影像中一些反映纹理特征的地物类型，由于随机条带噪声和地物类型存在高非线性相关性，因此随机条带噪声的成分随机分布于这些峰值周围，它们共同组成影像中地表地物反射率。步骤2，通过采用不同大小的处理窗口对步骤I所得亮度值曲线进行中值滤波去除基线，获取处理单元中心。依据步骤I所得亮度值曲线S (k)的特点可知，获取处理单元中心首先要进行亮度值曲线中的基线成分去除处理。本实施例采用中值滤波的方法来求取基线，事实证明它能有效的抑制峰值的影响.为了去掉亮度值曲线中的基线成分，使用一个大小为Wf的窗口作为中值滤波的处理窗口，输出的去基线的亮度值曲线S(k)f为S (k) f = S (k) -med (S (k_t) t G ff)(3)式(3)中集合W = UI-(Wf-I)/2 ( t ( (Wf-I)/2}，med 是指中值滤波处理。将去基线的亮度值曲线3&八采用均值y (k)f和方差O仏八进行归一化处理，处理公式如式(4)所示
权利要求
1.一种基于灰度分割的自适应矩匹配条带噪声去除方法，其特征在于，包括以下步骤 步骤1，获取原始含噪影像的亮度值曲线； 步骤2，通过采用不同大小的处理窗口对步骤I所得亮度值曲线进行中值滤波去除基线，获取处理单元中心； 步骤3，根据两个相邻处理单元中心的间距进行调整，确定处理单元； 步骤4，对单个处理单元采用标准矩匹配的处理方式进行影像处理。
2.如权利要求I所述基于灰度分割的自适应矩匹配条带噪声去除方法，其特征在于设原始含噪影像的大小为MXN，(i, j)表示原始含噪影像中第i行第j列的像素，其中i的取值范围为O，Ρ··Μ-1，j的取值范围为O，Ρ··Ν-1，步骤I中获取原始含噪影像的亮度值曲线是根据如下公式计算得到， /WVr-I) = { X Z/: (Z，力 11/从./) = k * /=0,/=0 其中，k表示对应的灰度级，iE(i, j)表示影像中位于原始含噪影像(i，j)处像素的灰度值。
3.如权利要求2所述基于灰度分割的自适应矩匹配条带噪声去除方法，其特征在于步骤2包括如下子步骤， 步骤2. 1，设中值滤波的处理窗口大小记为Wf，初始化设置Wf的值，设迭代次数Count=O ； 步骤2. 2，执行以下计算公式获取去基线的亮度值曲线S(k)f ； S (k) f=S (k) -med (S (k_t) t e ff) 其中，W= {t|-(WF-l)/2 ( t ( (Wf-I)/2}，med 是指中值滤波处理； 步骤2. 3，对去基线的亮度值曲线3&\采用均值μ (k)f和方差σ (k)f进行归一化处理，处理公式如下， n ^U)/ 其中，s(k)n就是均一化后的结果； 步骤2. 4，根据预设的阈值T选择候选处理中心，选择公式如下，肌j 1聊>7’ 当去基线的亮度值曲线s(k)f上某点的候选值P(k)。为I时，该点是候选处理中心，否则不是； 步骤2. 5,令Count=Count+l,判断Count是否达到预设次数,否则进入步骤2. 6,是则减小Wf的值，返回步骤2. 2迭代执行中值滤波； 步骤2. 6，对迭代所得的所有候选处理中心进行过滤， 首先计算判断函数值，计算公式如下， I I\k) > Count 12将某候选处理中心的判断函数值P(k)带入到步骤I所得亮度值曲线S(k)中，若对应的像素个数大于预设阈值I，则确定为处理单元中心，否则去除。
4.如权利要求3所述基于灰度分割的自适应矩匹配条带噪声去除方法，其特征在于预设阈值I采用以下公式计算，
5.如权利要求I或2或3或4所述基于灰度分割的自适应矩匹配条带噪声去除方法，其特征在于步骤3中对处理单元中心进行调整，实现方式为， 当两个相邻处理单元中心的间距小于预设阈值h时，则将两个相邻处理单元中心合并为一个处理单元中心，合并后的处理单元中心的值为这两个相邻处理单元中心的均值； 当两个相邻处理单元中心的间距大于或等于预设阈值h时，两个相邻处理单元中心之间增加一个新的处理单元中心，新的处理单元中心的值为原始两相邻处理单元中心的平均值。
6.如权利要求5所述基于灰度分割的自适应矩匹配条带噪声去除方法，其特征在于预设阈值h采用原始含噪影像的随机噪声灰度跃变最大值的4 5倍。
全文摘要
一种基于灰度分割的自适应矩匹配条带噪声去除方法，包括获取原始含噪影像的亮度值曲线；通过采用不同大小的处理窗口对亮度值曲线进行中值滤波去除基线，获取处理单元中心；根据两个相邻处理单元中心的间距进行调整，确定处理单元；对单个处理单元采用标准矩匹配的处理方式进行影像处理。该方法针对数据量大的遥感全色影像，基于基本影像处理理论，提供了一种有效的影像噪声去除方法。该方法能够有效地复原影像、改善影像质量，算法效率也很高。
文档编号G06T5/40GK102819827SQ20121023711
公开日2012年12月12日 申请日期2012年7月10日 优先权日2012年7月10日
发明者王密, 张炳先, 潘俊, 李德仁 申请人:武汉大学