[发明专利]栅格数据快速区域统计方法和系统

权利要求书

1.一种栅格数据快速区域统计方法，其特征在于，包括：

将目标矢量范围从第一坐标系转换到第二坐标系下，获取所述目标矢量范围在所述第二坐标系的竖直方向上的顶端坐标和底端坐标；所述目标矢量范围为待统计区域的边界围成的闭合曲线，所述第一坐标系为所述待统计区域的矢量数据所在的坐标系，所述第二坐标系为目标栅格数据所在的坐标系，所述目标栅格数据为包括所述待统计区域的栅格数据；

计算多条预设扫描线与所述目标矢量范围的交点的地理坐标；所述多条预设扫描线为所述第二坐标系中，分布在所述顶端坐标与所述底端坐标之间的、与水平方向平行的、以所述栅格数据在竖直方向的分辨率为间隔的多条直线；

基于所述地理坐标，将每条扫描线与所述目标矢量范围的交点两两分组，得到多个交点对；

将每个交点对的地理坐标映射到所述目标栅格数据的行列坐标，得到多个像素点对；

在所述目标栅格数据中，基于每个像素点对之间的像素点对所述待统计区域进行栅格数据区域统计，得到统计结果。

2.根据权利要求1所述的方法，其特征在于，在将目标矢量范围从第一坐标系转换到第二坐标系下之前，所述方法还包括：

基于待统计区域的矢量数据，获取目标矢量范围。

3.根据权利要求1所述的方法，其特征在于，基于所述地理坐标，将每条扫描线与所述目标矢量范围的交点两两分组，得到多个交点对，包括：

将每条扫描线与所述目标矢量范围的交点，按照地理坐标的水平坐标排序，再按照排序两两分组，得到多个交点对。

4.根据权利要求1所述的方法，其特征在于，将每个交点对的地理坐标映射到所述目标栅格数据的行列坐标，包括：

通过如下算式进行映射：col＝floor(x-x0)/xreso，row＝floor(y-y0)/yreso；其中，x和y分别为交点的地理坐标的X轴坐标和Y轴坐标；x0为所述目标栅格数据左端X轴的地理坐标值；y0为所述目标栅格数据的顶端Y轴地理坐标值；xreso和yreso分别为X和Y方向空间分辨率；floor为向下取整函数；col为所述目标栅格数据的列坐标，row为所述目标栅格数据的行坐标。

5.根据权利要求1所述的方法，其特征在于，在所述目标栅格数据中，基于每个像素点对之间的像素点对所述待统计区域进行栅格数据区域统计，得到统计结果，包括：

获取每个像素点对之间的所有像素点的像素值；

统计所述像素值中的有效像素值的最大值、最小值、累加值和平方累加值；

统计所述有效像素值所对应的像素点的数量值，得到有效像元数；

基于所述累加值、所述平方累加值和所述有效像元数，得到所述待统计区域的平均值和所述待统计区域的方差；

将所述最大值、所述最小值、所述待统计区域的平均值、所述待统计区域的方差和所述有效像元数，作为统计结果。

6.一种栅格数据快速区域统计系统，其特征在于，包括：转换模块，计算模块，分组模块，映射模块和统计模块，其中，

所述转换模块，用于将目标矢量范围从第一坐标系转换到第二坐标系下，获取所述目标矢量范围在所述第二坐标系的竖直方向上的顶端坐标和底端坐标；所述目标矢量范围为待统计区域的边界围成的闭合曲线，所述第一坐标系为所述待统计区域的矢量数据所在的坐标系，所述第二坐标系为目标栅格数据所在的坐标系，所述目标栅格数据为包括所述待统计区域的栅格数据；

所述计算模块，用于计算多条预设扫描线与所述目标矢量范围的交点的地理坐标；所述多条预设扫描线为所述第二坐标系中，分布在所述顶端坐标与所述底端坐标之间的、与水平方向平行的、以所述栅格数据在竖直方向的分辨率为间隔的多条直线；

所述分组模块，用于基于所述地理坐标，将每条扫描线与所述目标矢量范围的交点两两分组，得到多个交点对；

所述映射模块，用于将每个交点对的地理坐标映射到所述目标栅格数据的行列坐标，得到多个像素点对；

所述统计模块，用于在所述目标栅格数据中，基于每个像素点对之间的像素点对所述待统计区域进行栅格数据区域统计，得到统计结果。