专利名称：Modis遥感影像中海面条带噪声和条状云的过滤方法
技术领域：
本发明涉及MODIS遥感影像处理领域中一种图像预处理的方法，尤其是海面条带噪声和条状云的过滤方法。
背景技术：
MODIS (Moderate Resolution Imaging Spectroradiometer,中分辨率成像光谱辐射计)是搭载在terra和aqua卫星上的一种光学遥感仪器，是卫星上唯一将实时观测数据通过X波段向全世界直接广播，并可以免费接收数据并无偿使用的星载仪器，全球许多国家和地区都在接收和使用MODIS数据，由美国国家航空航天局提供免费数据，我国已经建立了多个MODIS遥感影像地面接收站。MODIS具有成像范围大，噪声低的优点，非常适合海洋内波的大范围研究，因为太阳光照条件对MODIS有很大的限制，云的干扰对内波的观测和识别也有很大的干扰。 海洋内波的最大振幅是发生在海洋内部的，与表面波相比，内波具有更大的能量和破坏力，从而会对水下航行的潜艇和海上石油平台的安全造成严重的威胁。因此，预报海洋内波，降低海洋内波带来的损失，早发现早预防，对于航海、海上科研和海上生产，特别是深海石油平台的安全是具有现实意义的。使用遥感卫星对海洋内波进行观测是当前海洋内波研究的主要方法之一。MODIS各级数据产品均是以HDF文件格式存储、处理和发布的。HDF文件格式是一种高效存储和分发科学数据的新型的数据格式，是美国国家高级计算应用中心(NationalCenter for Supercomputing Application)研制的，主要是为了满足各领域的研究需求。HDF文件能够由一些通用软件，如ENVI、MATLAB等直接读取。ENVI将所有数据直接读取出来后，不容易实现在软件基础上的再次开发，而MATLAB可以直接读取HDF文件中的原始数据，然后在读取的数据的基础上进行编程、开发，实现内波检测和信息提取等工作，且程序实现相对容易，程序运行效率较高。HDF文件中的原始数据灰度值偏低，视觉效果较差，很难直接进行图像分析，从而对于海洋内波的检测很困难，所以必须首先对获取的数据进行增强处理，使其更适合于海洋内波的检测算法。在MODIS遥感影像中海洋内波表现为一种明暗相间的条纹，由此出发，可以采用纹理分析的方法对海洋内波的特征进行提取，实现内波的自动检测。在MODIS遥感影像中，有两类因素会对内波条纹的纹理分析带来很大的影响一是MODIS遥感影像中呈条带状分布的噪声，这种条带噪声是由卫星传感器的光、电器件在反复扫描地物的成像过程中，由扫描探测元正反扫描响应差异、传感器的机械运动和温度变化等造成的；另一类就是海平面上空的云，特别是较细的条状的云。这两类因素也具有很强的纹理特征，从而会对内波的识别造成很大的影响。所以在提取海洋内波的纹理特征之前，为了减少这些噪声的干扰，首先就要对这些噪声进行滤除，尽量减少它们造成的干扰，目前都是采用在同一地区、不同时段的多幅遥感图像进行常规的噪声去云处理，处理效率低且处理的效果较差。
发明内容
本发明提供一种MODIS遥感影像中海面条带噪声和条状云的过滤方法，能够在短时间内快速有效过滤掉影像MODIS海面内波检测的条带状噪音和条状云，同时有效保留MODIS影像中海洋内波的基本信息，为海洋内波的自动检测提供了很好的基础。本发明采用的技术方案是包括以下步骤(I)由HDF文件中读取同一地区、同一时段的一幅MODIS遥感影像的原始图像；(2)对原始图像进行直方图均衡化处理；(3)将直方图均衡化处理后的图像分割为300X 300像素的各个子图；(4)求出各个子图的灰度平均值; (5)选取阈值= fl-10，小于 阈值的像素为条状云或条带噪声的像素，大于阈值A的像素则为海平面或海洋内波的像素。进一步地，本发明在小于阈值的被检测的像素点周围选取IOX 10的区域随机进行探测，若探测到的像素点的灰度值大于阈值A = +10，则由探测点的灰度值代替被检
测像素点灰度值；探测次数的上限1000次，若探测次数达到上限，则由值《3 = flf + 5代替被检测像素点灰度值。本发明的优点和效果是与现有技术中采用多幅同一地区、不同时段的遥感图像进行去云处理不同，本发明采用同一地区、同一时段的一幅MODIS遥感影像进行影响内波检测的条状云和条带噪音的过滤。由于光照条件的不同，同一幅图像上的不同的区域，区分海平面背景和云的阈值也不相同。基于此，采用了将整幅大图分割为多幅300 X 300的子图，在子图内光照条件基本相同，采用固定阈值的方法区分海平面背景和云，对MODIS遥感影像进行特殊的处理，能够在短时间内对某一地区某一时段的MODIS影像中海平面的条状云和条带状噪音进行有效的滤除，方法简单，运算效率高，对海平面的条状云和条带噪音的滤除效果好。下面结合附图和具体实施方式
对本发明作进一步的详细说明。
图I是对原始图像进行直方图均衡化后的MODIS影像 图2是图I中局部海洋内波与条带噪声 图3是图I中局部海面上的条状云 图4为本发明过滤方法的流程 图5为采用本发明方法的处理结果图。
具体实施例方式参见图4，本发明的具体实施步骤如下
一、由HDF文件中读取MODIS遥感影像的原始图像
MODIS各级数据均是以HDF文件格式存储、处理和发布的。HDF文件格式是一种高效存储和分发科学数据的新型的数据格式，是美国国家高级计算应用中心(National Centerfor Supercomputing Application)研制的,主要是为了满足各领域的研究需求。HDF文件能够由一些通用软件如ENVI、MATLAB等直接读取。由HDF文件中读取MODIS遥感影像的原始图像时，选择250m分辨率的原始数据。ENVI将所有数据直接读取出来后，不容易实现在软件基础上的再次开发，而MATLAB可以直接读取HDF文件中的原始数据，然后在读取的数据的基础上进行编程、开发，实现内波检测和信息提取等工作，且程序实现相对容易，程序运行效率较高。二、对原始图像进行直方图均衡化
HDF文件中的原始数据灰度值偏低，视觉效果较差，很难直接进行图像分析，从而对于内波的检测很困难，所以必须首先对获取的数据进行增强处理，使其更适合于内波的检测。本发明采用直方图均衡化的方法，可以有效地改善图像的质量。a、直方图的基本原理
灰度级为
范围的数字图像/(AJ)的直方图A(rfc)=fl.fc，这是一个离散函数，其中4表示第k级灰度，^表示图像中灰度级为h的像素个数。归一化的直方图是P{rk)^nk!n , * = 0,1,1 ,其中n表示图像中像素的总数，L为图像中可能的灰度级总
数，Pfrt)表示灰度级为％的像素出现的概率估计值，E =I。直方图增强就是一种仅仅依靠给定图像的直方图信息，通过变换函数来增强图像的方法。b、直方图均衡化
对于数字图像f(x,y)，灰度级I出现的概率为P(Tk) = nJn ’ ^ = I…，i-1,其
中，n为图像中像素的总数；％是图像中灰度级为rk的像素个数；L为图像中可能的灰度级总数。变换函数％ = 7(1), A = 0, I, ,￡-1,其中％是变换后的灰度级。变换函数W)须满足两个条件(I)当r被归一化到[W]时，W)在区间[CU]上单调递增；(2)当re
时，S=T(r) €
。条件(I)保证了输出图像仍然是从黑到白的顺序；条件(2)保证输出灰度级与输入灰度级有同样的范围。直方图的均衡化的变换函数可表示为
h = Tirk)^ Pirj) ^ ^U，使得直方图均衡化后的图像灰度级能够跨
JnmQ^一 * I, * ■ * , Lj — I
越更大的范围，直方图均衡化后图像灰度级为64。直方图均衡化后，MODIS遥感影像可以由肉眼清楚地进行识别。如果含有内波，则其中的内波影像比较清楚，这就具备了对MODIS遥感影像中是否含有内波进行自动检测的基础。对由MATLAB直接读取HDF文件获得的MODIS影像进行直方图均衡化后，结果如图I所示，图像灰度级为64。获取的图I中局部海洋内波与条带噪声如图2所示，获取的图I中局部海面上的条状云如图3所示。三、将经直方图均衡化后的图像分割为300 X 300像素的各个子图
由于光照的不同，在MODIS影像的不同区域，将云、条带噪声和海平面背景有效区分的阈值也不同，因此本发明采用固定阈值，固定阈值不能将整个的待检查区域中的条状云与海面背景以及内波有效地区分开来。本发明首先将待检查区域的图像分割成多个较小的子图像；然后在子图像内，光照情况基本相同，采用固定的阈值将海平面背景与条状云及条带状噪声进行分割。而针对不同的子图像，采用自适应的方法调整阈值；最后将条带状噪声和条状云的像素灰度值由海平面背景的灰度值代替，这样将条带状噪声和条状云部分与海平面背景融为一体，以消除它们对内波检测所造成的干扰。将大图分割成大小不同的子图的多次实验中发现在250m分辨率的MODIS遥感影像中，如果子图分割过大，当然内波信息很完整，但算法执行的时间将会很长，而且在较大的子图中，由于光照等条件的不同，对子图进行条状云与条带噪声过滤的阈值就非常难以确定，不利于条带噪声和条状云的滤除；如果子图分割过小，则子图中含有的内波信息不完整，而如果将大图分割为UU 500像素的子图，则对于阈值的选择较为有利，而且该子图的大小对于进一步的提取内波的纹理特征也非常有利。此时阈值可以由子图像素灰度的平 均值来确定。四、依次对每个子图进行如下操作
采用图像处理技术中常规的方法可求出各个子图的灰度平均值a。
对2000多幅300 300像素的子图的多次实验，找到将MODIS内波影像中背景与云区分开的阈值。选取阈值4=0-10 ,小于阈值巧的像素为云或者条带状噪声的像素，大于阈值巧的像素则为海平面或者海洋内波的像素。在小于阈值^的被检测的像素点周围，选取IOX 10的区域随机地进行探测，如果
探测到的像素点的灰度值大于阈值A + W，则由探测点的灰度值代替被检测像素点灰度值。为避免在大块云的附近(大于10X10)随机探测的时间过长，给定探测次数的上限，如果达到上限，则直接由背景值％ = + 5取代该被检测点像素的值。本发明提出的MODIS遥感影像中海面条带噪声和条状云的过滤方法，是针对同一地区、同一时段获取的一幅MODIS遥感影像进行条状云和条带噪音去除处理的方法。采用将原始数据进行直方图均衡化后，将整幅大图的海平面的区域分割为多幅300X300的子图，然后在每一幅子图内，基于图的灰度的平均值给出分割海面背景和条状云和条带噪音的阈值。在此基础上，将影响MODIS内波检测的条状云和条带状噪音进行有效的滤除，同时保留内波的条纹信息。图5是采用本发明方法的处理结果图，其中图5a、5b、5c为原始图像直方图均衡化后获得的三幅300 00的子图，图5A、5B、5C分别是图5a、5b、5c经过滤处理之后的结果图。
权利要求
1.一种MODIS遥感影像中海面条带噪声和条状云的过滤方法，其特征在于包括以下步骤 (1)由HDF文件中读取同一地区、同一时段的一幅MODIS遥感影像的原始图像； (2)对原始图像进行直方图均衡化处理； (3)将直方图均衡化处理后的图像分割为300' 300像素的各个子图； (4)求出各个子图的灰度平均值a； (5)选取阈值A= -10 ,小于阈值^的像素为条状云或条带噪声的像素，大于阈值 的像素则为海平面或海洋内波的像素。
2.根据权利要求I所述的过滤方法，其特征在于步骤(5)中，在小于阈值％的被检测的像素点周围选取10X10的区域随机进行探测，若探测到的像素点的灰度值大于阈值% = fl+10，则由探测点的灰度值代替被检测像素点灰度值；探测次数的上限1000次，若探测次数达到上限,则由值％ = a+ 5代替被检测像素点灰度值。
3.根据权利要求I所述的过滤方法，其特征在于步骤(I)中，原始图像的分辨率是250mo
4.根据权利要求I所述的过滤方法，其特征在于步骤(2)中，直方图均衡化后的图像灰度级为64。
5.根据权利要求I所述的过滤方法，其特征在于步骤(3)中，子图的个数大于2000。
全文摘要
本发明公开一种MODIS遥感影像中海面条带噪声和条状云的过滤方法，由HDF文件中读取同一地区、同一时段的一幅MODIS遥感影像的原始图像；对原始图像进行直方图均衡化处理；将直方图均衡化处理后的图像分割为300300像素的各个子图；求出各个子图的灰度平均值；选取阈值，小于阈值的像素为条状云或条带噪声的像素，大于阈值的像素则为海平面或海洋内波的像素；能够在短时间内对某一地区某一时段的MODIS影像中海平面的条状云和条带状噪音进行有效的滤除，同时有效保留MODIS影像中海洋内波的基本信息，方法简单，处理效率高。
文档编号G06T5/20GK102750683SQ201210200149
公开日2012年10月24日 申请日期2012年6月18日 优先权日2012年6月18日
发明者姜忠义, 康慧燕, 江兴方, 陈芳芳 申请人:常州大学