[发明专利]整合模式识别与上下文知识的光学遥感图像云判别方法

说明书

技术领域

本发明涉及一种光学遥感图像云判别方法，特别涉及一种整合统计模式识别与上下文知识的光学遥感图像云判别方法，属于图像识别技术领域。

背景技术

对于光学遥感图像，我们感兴趣的更多是地面目标，然而在很多情况下是存在云的遮挡的，这样云的部分是没有任何地面信息，却占用了处理系统中大量的存储空间处理能力和传输带宽。人工筛选和剔除这些数据通常需要花费大量的时间与精力，如果这些“无效区域”可以通过某种方法找出来，就可以将这部分数据去掉，进而节省数据存储空间和数据处理人员的工作量，特别对卫星上的系统，还可以节省下行带宽以传送更重要的数据。

因此，一种行之有效的自动云判技术就变得非常有必要。目前的云判别的方法的主要是利用模式识别的方法进行处理，通常包括对原图的抽样分块、特征提取、分类判决以及云区域剔除。这些方法通常只利用某一个分块本身的信息，一般用来处理分块中只有单一解译对象时效果较好。但当处理云和地物交界地方的分块时，由于其含有两个或两个以上解译对象，使得模式识别的判决效果变差，甚至，变为无效。例如，当一个分块中云和地物分布集中且含量各是一半（上半个分块是云且下半个分块是地物），这个分块本身就是既不是云也不是陆地，所以模式识别的结果就没有什么意义。

发明内容

有鉴于此，本发明提供了一种整合模式识别与上下文知识的光学遥感图像云判别方法，能够解决现有技术中只利用某一个分块本身的信息进行模式识别所带来的弊端。

为了解决上述技术问题，本发明是这样实现的：

一种整合模式识别与上下文知识的光学遥感图像云判别方法，包括：

第一步，对原图进行抽样及分块：读入一幅光学遥感图像并进行抽样和分块；所述分块采用有重叠的分块方式，即采用L×L的正方形在光学遥感图像中进行水平和竖直方向的步进，步进量为L/2，获得多个分块，称为Tile，每四个上下左右相邻相互重叠的Tile分块可以看成是九个不重叠的大小为L/2×L/2的小分块构成，这种小分块称为Block；

第二步，对各个Tile分块进行模式识别：利用模式识别依次对每个Tile分块进行判决，并且将判决值存入Tile索引矩阵，判决值包括属于云的置信度、属于地物的置信度和不确定，这里将判决值不确定赋值为基准B，置信度间隔为D，云的判决值应该是B+mD，对应地物的判决值是B-nD，其中m和n取正整数，云或地物的置信度越高m或n的值就越大；

第三步、确定不同类型对象特征的阈值范围；所述不同类型对象包括不同置信度的云、不确定和不同置信度的地物；

第四步，根据Tile索引矩阵确定Block的判决值并进行微调：

步骤S41)对于每一个Tile分块，分别计算组成该Tile分块的4个Block的灰度均值，基于第三步中得到的阈值范围采用门限判决法对Block进行云判决，判决值采用第二步中所用的形式，即不确定时为B，判为云时B+mD，判为地物时B-nD；

步骤S42)利用Tile分块及其内部4个Block的判决值计算Tile的微调因子Factor：

Factor=Sum-(4-n)×Predictn]]>

其中，Sum指Tile中判为云或地物Block判决值之和，n指Tile中判为不确定Block个数，Predict指Tile的判决值；

步骤S43)对于Tile中判决为不确定的Block，将步骤S42)中得到的微调因子与Tile的判决值相加作为该Block的临时判决值；

第五步，每一Block属于多个Tile分块，因此每一Block具有多个临时判决值，根据同一Block的临时判决值综合确定Block的判决值，Block判决值组成Block索引矩阵；利用上下文知识，对Block索引矩阵中的判决值进行修正；

第六步，利用修正后的Block判决值所标记的云位置，对光学遥感图像进行云剔除。

优选地，所述第三步为：利用Tile的模式识别结果确定不同类型对象特征的阈值范围；具体为：