[发明专利]一种高分辨率建设用地图斑识别方法

权利要求书

1.一种基于PanTex和直线特征的高分辨率建设用地图斑识别方法，其特征在于：包括以下步骤：

步骤一、对高分辨率遥感影像和对应的土地利用图斑进行影像与矢量配准；

步骤二、对步骤一处理后的土地利用图斑矢量文件与高分辨率遥感影像进行叠加掩模操作，获取每个图斑多边形对应的影像区域，得到独立图斑影像；

步骤三、计算步骤一处理后高分辨率遥感影像的PanTex特征图像，统计每个图斑内PanTex指数的总和；

步骤四、对经过所述步骤二处理后的独立图斑影像，提取图斑影像内的直线；

步骤五、对经过所述步骤四处理后的直线计算图斑的直线特征；

步骤六、基于步骤三和步骤五得到的PanTex指数特征和直线特征，利用SVM两类分类器对所述步骤二得到的图斑进行分类，提取建设用地图斑。

2.如权利要求1所述的基于PanTex和直线特征的高分辨率建设用地图斑识别方法，其特征在于：步骤一中的高分辨率遥感影像与矢量文件配准的方法是利用ArcGIS软件的空间校准工具，通过仿射模型校正矢量文件，将矢量文件与基期影像配准对齐。

3.如权利要求1所述的基于PanTex和直线特征的高分辨率建设用地图斑识别方法，其特征在于：步骤二中叠加掩模操作是将步骤一处理后土地利用图斑矢量文件叠加至高分辨率遥感影像，通过矢量文件中每个多边形的边界，将边界内的影像区域裁剪出来，存为独立的图斑影像。

4.如权利要求1所述的基于PanTex和直线特征的高分辨率建设用地图斑识别方法，其特征在于：步骤三中PanTex特征图像的计算步骤如下：

1、选取如下10个方向

(shiftX,shiftY)＝{(1,-2),(1,-1),(2,-1),(1,0),(2,0),(0,1),(1,1),(2,1),(0,2),(1,2)}

其中，shiftX为像素间的X方向偏移量，shiftY为像素间的Y方向偏移量；

2、计算灰度共生矩阵

设f(x,y)为一幅二维数字图象，其大小为M×N，灰度级别为Z，则满足一定空间关系的灰度共生矩阵公式为：

P(i,j)＝#{(x₁,y₁),(x₂,y₂)∈M×N|f(x₁,y₁)＝i,f(x₂,y₂)＝j}

其中，P为Z×Z的矩阵，#{x}表示集合x中的元素个数，i、j、x₁、y₁、x₂和y₂的取值范围均为[0,Z-1]；

利用步骤1的10个方向计算各个方向的灰度共生矩阵P(i,j|shiftX,shiftY)；

3、计算PanTex特征图像

3.1、利用步骤2中的10个方向的灰度共生矩阵计算10个对比度测度；

对比度测度计算公式如下：

3.2、分别对步骤3.1中的10个对比度测度进行归一化处理；

3.3、将步骤3.2归一化处理后的10幅纹理测度图像进行逐像素取最小值，得到PanTex特征图像。

5.如权利要求1所述的基于PanTex和直线特征的高分辨率建设用地图斑识别方法，其特征在于：步骤三中统计每个图斑内PanTex指数的总和是将PanTex特征图像与土地利用图斑叠加，统计每个多边形内部PanTex特征值的总和。