[发明专利]一种基于边界裁剪的高速高精度矢栅叠置分析方法

权利要求书

1.一种基于边界裁剪的高速高精度矢栅叠置分析方法，用于提取地表面积的计算，其特征在于：包括以下步骤：

步骤S1：收集待测面积的栅格数据集A和矢量数据集B的数据；

步骤S2：遍历所述的矢量数据集B中的多边形b获取四至坐标极值；

步骤S3：遍历所述的多边形b获取完整覆盖多边形的栅格数据集A的最小栅格切片a；

具体步骤为：

步骤S31：读取所述的栅格数据集A的起始点坐标X0、Y0，以及像元尺寸Xs、Ys；

步骤S32：读取所述的多边形b的四至坐标极值XE、XW、YN、YS；

步骤S33：计算所述的栅格切片a的起始行号Rstart、起始列号Cstart、终止行号Rstop、终止列号Cstop，并读取栅格切片a的数据：

步骤S4：按行遍历所述的栅格切片a获取像元行的外包矩形R；

步骤S5：用所述的外包矩形R裁切所述的多边形b，得到结果点集多边形S；将遍历结果点集多边形S得到的结果存入结果集C，按照面积占比分配像元属性值并计入结果集C；

步骤S6：统计所述的结果集C中的属性值作为所述的栅格数据集A与所述的矢量数据集B进行叠置分析的结果。

2.根据权利要求1所述的一种基于边界裁剪的高速高精度矢栅叠置分析方法，其特征在于：所述的步骤S2中，具体步骤为：

步骤S21：遍历所述的矢量数据集B的多边形b；

步骤S22：初始化所述的多边形b的四至坐标极值分别为XE、XW、YN、YS；

步骤S23：遍历所述的多边形b的节点，设为(Xb，Yb)；

步骤S24：判断执行边界，若XbXE则XE＝Xb；若XbXW且XW＝Xb；若YbYN则YN＝Yb；若YbYS则YS＝Yb；

步骤S25：待所述的多边形b的节点遍历完成后，将所述的四至坐标极值XE、XW、YN、YS计入所述的多边形b的属性列表，从步骤S21循环执行，直到遍历完所述的矢量数据集B中的所有多边形。

3.根据权利要求1所述的一种基于边界裁剪的高速高精度矢栅叠置分析方法，其特征在于：所述的步骤S4中，具体步骤为：

步骤S41：遍历所述的栅格切片a的数据行，记为row；

步骤S42：读取所述的数据行row的行号Ra、起始列号Cstart、终止列号Cstop；

步骤S43：创建所述的数据行row的外包矩形R，分别计算R的顶点坐标(Xmin，Ymax)、(Xman，Ymax)、(Xmax，Ymin)、(Xmin，Ymin)：

Xmin＝Cstart·Xs，

Xmax＝(Cstop+1)·Xs，

Ymin＝(Ra+1)·Ys，

Ymax＝Ra·Ys。