[发明专利]一种地理栅格数据任选多边形区域的直方图实时生成方法

权利要求书

1.一种地理栅格数据任选多边形区域的直方图实时生成方法，其特征在于包括以下步骤：

1)进行预处理

首先将原始栅格数据均匀分割成一定尺寸的瓦片,然后分别统计原始栅格数据中各瓦片的 直方图并将其以文件方式存储到服务器；

2)判定各瓦片与用户选定多边形区域的关系，分别确定完全处于选定区域内部的瓦片和 处于选定区域边界的瓦片；

3)对于步骤2)确定的处于选定区域内部的瓦片，读入步骤1)中生成的文件，读取完全 处于选定区域内部的瓦片的直方图数据并将各瓦片的直方图数据汇总求和，得到直方图 inner_hist；

4)对于步骤2)确定的处于选定区域边界的瓦片，采用并行方式读取原始栅格数据并统计 直方图，最终将统计的直方图数据汇总求和得到直方图bound_hist；

5)将步骤3)中得到的直方图数据bound_hist和步骤4)中得到的直方图数据bound_hist 求和得到选定区域的直方图。

2.根据权利要求1所述的地理栅格数据任选多边形区域的直方图实时生成方法，其特征 在于：步骤2)的方法为：

首先计算能够将用户选定多边形区域完全包容进去的最小矩形即最小矩形外包框，将其表 示为(x1，y1，x2，y2)，其中(x1，y1)和(x2，y2)分别为最小矩形外包框的左下角和右 上角的顶点坐标；

基于客户端显示的缩放级别确定分辨率R，那么客户端显示的一个像素点就代表R×R个 栅格单元，从而可利用下取整函数floor得到纵向填充的范围跨度是从floor(y1/R)到 floor(y2/R)+1，横向填充的范围跨度是从floor(x1/R)到floor(x2/R)+1；

然后利用扫描线填充算法在分辨率R下对选定区域进行填充，假设一个瓦片有N×N个栅 格单元，那么为每个瓦片建立一个(N/R)×(N/R)的填充矩阵表示瓦片的所有点；在填充过程 中，对于在选定区域内的点，将填充矩阵中与该点对应的位置置为1，对于选定区域外的点， 填充矩阵中与该点对应的位置置为0；同时记录各填充矩阵中被置为1的点的个数；对于一个 填充矩阵，如果填充矩阵中的每一个点都被置为1，那么可以判断这个填充矩阵对应的瓦片完 全在选定区域内；那些未全部被置为1的填充矩阵对应的瓦片则处在选定区域边界。

3.根据权利要求2所述的地理栅格数据任选多边形区域的直方图实时生成方法，其特征 在于，步骤4)的方法为：

对于步骤2)确定的处于选定区域边界的瓦片，根据瓦片对应的填充矩阵确定瓦片处于选 定区域内的部分；将每个处于选定区域边界的瓦片作为一个子任务，程序的主进程P0将子任 务分配给不同的工作进程，将处在选定区域边界的瓦片对应的填充矩阵发送给不同的工作进程， 各工作进程分别使用开源栅格空间数据转换库GDAL在原始栅格数据中提取瓦片并根据瓦片对 应的填充矩阵来确定瓦片处于选定区域内的部分，然后统计这部分的直方图，最终将各工作进 程求得的结果进行求和得到处于选定区域边界处瓦片对应的直方图。