[发明专利]一种基于土地覆盖信息的中分影像合成方法

权利要求书

1.一种基于土地覆盖信息的中分影像合成方法，其特征在于包括以下步骤：

步骤1、时间序列影像预处理；

步骤2、对目标区域划分格网，然后用格网裁切覆盖目标区域的时间序列影像得到分块数据；

步骤3、对分块数据进行NDVI最大值影像合成得到对应分块合成影像A；

步骤4、对分块数据进行蓝波段反射率最小值影像合成得到对应分块合成影像B；

步骤5、基于土地覆盖类型合成影像A和B，得到分块最终影像C，拼接分块影像，得到最终的合成影像D，具体如下：

基于目标区的土地覆盖先验数据将目标区分为植被区和非植被区，所述植被区包括农田、森林、草地、灌木，非植被区包括水体、不透水地表、裸地、雪、冰；对于每个像元所在位置，如果土地覆盖先验数据为植被区，选取影像A对应的像元作为合成后影像的像元值；反之，如果土地覆盖先验数据为非植被区，选取影像B对应的像元作为合成后影像的像元值；影像A和B合成得到分块最终影像C；拼接分块影像得到最终的合成影像。

2.根据权利要求1所述的一种基于土地覆盖信息的中分影像合成方法，其特征在于：所述步骤1的具体过程为：对固定时间段内时间序列影像做三步预处理工作，包括大气层顶反射率计算、正射校正和NDVI计算。

3.根据权利要求2所述的一种基于土地覆盖信息的中分影像合成方法，其特征在于：所述大气层顶反射率计算具体为：

使用辐射定标系数把影像各个波段的DN值转为辐亮度，具体公式如下：

L＝Gain*DN+Bias

其中L是卫星载荷波段入瞳处等效辐射亮度，Gain和Bias分别为定标系数增益、偏移量；

计算各波段的大气顶层反射率B，具体公式如下：

其中B为大气顶层反射率，i为第i波段，π为常量，L为辐亮度，D为日地距离，ESUN为大气层顶的平均太阳光谱辐照度，θ为太阳的天顶角。

4.根据权利要求2所述的一种基于土地覆盖信息的中分影像合成方法，其特征在于：所述正射校正具体为：对时间序列影像进行正射校正，消除原始影像的几何形变，并保证时间序列影像之间的几何偏差要小于一个像元。

5.根据权利要求2所述的一种基于土地覆盖信息的中分影像合成方法，其特征在于：所述NDVI计算具体为：

单景影像的NDVI计算需要红波段和近红外波段的反射率，具体公式如下：

NDVI＝(NIR–R)/(NIR+R)

其中R和NIR分别为红波段和近红外波段的大气层顶反射率。

6.根据权利要求1所述的一种基于土地覆盖信息的中分影像合成方法，其特征在于：所述步骤2的具体过程为：采用分块的方式进行影像合成；预先划定方形的矢量格网覆盖目标区域，每块大小根据计算能力确定。