[发明专利]一种基于Landsat8影像二维特征空间的芦笋种植面积提取方法

权利要求书

1.一种基于Landsat 8影像二维特征空间的芦笋种植面积提取方法，其特征在于，具体步骤包括：

(1)获取数据

获取芦笋种植区域的相关数据，芦笋种植区域的相关数据包括：芦笋种植区域经纬度、芦笋种植区域地表覆盖类型、芦笋种植区域地形、芦笋种植区域耕地参数；选取Landsat 8卫星数据作为芦笋种植区域遥感影像；

(2)对步骤(1)获取的芦笋种植区域遥感影像依次进行影像剪裁、几何校正、大气辐射校正处理；

(3)对通过步骤(2)获取的芦笋种植区域遥感影像求取归一化差值植被指数NDVI，进行阈值分割，掩膜去除除芦笋外的其它植被像元，除芦笋外的其它植被像元包括：水体像元、小麦像元及除芦笋、水体、小麦外的其它植被像元；

归一化差值植被指数NDVI是指植被生长状态及植被覆盖度的最佳指示因子；

植被指数NDVI的计算公式如式(Ⅰ)所示：

NDVI＝(NIR-R)/(NIR+R) (Ⅰ)

式(Ⅰ)中，NIR表示Landsat 8卫星第5波段即近红外波段的反射率；R表示Landsat 8卫星第4波段即可见光红波段的反射率；

求取的归一化差值植被指数NDVI中，NDVI<0的像元为水体像元，NDVI>0.3的像元为小麦像元及除芦笋、水体、小麦外的其它植被像元；

(4)通过步骤(3)获取的芦笋种植区域遥感影像包括建筑物、道路和芦笋，获取建筑物的像元、道路的像元和芦笋的像元在近红外波段和可见光红波段的二维特征空间，即NIR-Red特征空间的分布散点图；在获取的分布散点图中，道路的像元和建筑物的像元呈离散分布，芦笋的像元呈土壤线分布，对芦笋的像元拟合出土壤线y＝ax+b，a＝0.618,其中，x＝R，y＝NIR；

(5)在遥感图像处理软件envi中，采用波段运算的方法，将近红外波段和可见光红波段的值进行线性运算，所述波段运算的计算公式为：r＝NIR-aR；其中，r表示波段运算后的结果，NIR表示近红外波段，R表示可见光红波段；

(6)经过步骤(5)所述波段运算，确定合理阈值，对步骤(3)得到的掩膜结果图像进行阈值分割，掩膜去除建筑物及道路，提取芦笋种植区域的面积。

2.根据权利要求1所述的一种基于Landsat 8影像二维特征空间的芦笋种植面积提取方法，其特征在于，所述影像剪裁，具体是指：在envi软件中，选择主菜单中Open Image File，打开芦笋种植区域遥感影像，点击basic tool下的resize data，选择Spatial Subset；进入芦笋种植区域遥感影像裁剪界面，选择根据已有的ROI文件，进行影像剪裁，选择路径及文件名，完成影像剪裁。