[发明专利]一种地表覆盖表面积栅格化计算方法及装置

说明书

本发明实施例提供一种地表覆盖表面积栅格化计算方法及装置，所述方法包括：将地表覆盖矢量数据栅格化至目标格网大小的栅格数据；将数字高程模型重采样至目标格网大小的栅格数据；根据重采样的数字高程模型数据进行坡度计算；计算表面积栅格化数据；将目标区域划分为多个子区域，将各个子区域对应的地表覆盖栅格化数据及表面积栅格化数据同步处理，获得各个子区域对应的各地类表面积；将所有子区域对应的各地类表面积汇总，获得目标区域对应的各地表覆盖地类的表面积。所述方法极大地降低了表面积的计算量，实现了地表覆盖表面积的轻量、并行计算，能够解决大数据时代的海量时序地表覆盖数据的表面积高效、动态计算技术难题。

技术领域

本发明涉及地理空间信息领域，具体而言，涉及一种地表覆盖表面积栅格化计算方法及装置。

背景技术

随着智慧城市大数据时代的到来，对地观测、卫星遥感、生态评估、国土监管逐渐向宏观、动态、精细化方向发展，诸多行业领域对于精确动态的地表覆盖表面积快速统计与动态对比具有较大潜在需求。

传统的地表覆盖表面积统计方法的计算方法复杂度高，采用传统的串行计算逻辑，难以实现大规模并行计算，总体计算效率低下，不满足大数据时代地表覆盖表面积大尺度、高性能统计的现实需求。

发明内容

有鉴于此，本发明实施例的目的在于提供一种地表覆盖表面积栅格化计算方法及装置，以解决上述问题。

为了实现上述目的，本发明实施例采用的技术方案如下：

第一方面，本发明实施例提供了一种地表覆盖表面积栅格化计算方法，该方法包括：将地表覆盖矢量数据栅格化至目标格网大小的栅格数据，获得地表覆盖栅格化数据；将数字高程模型重采样至所述目标格网大小的栅格数据，获得重采样的数字高程模型数据；根据所述重采样的数字高程模型数据进行坡度计算，获得坡度栅格数据；根据所述坡度栅格数据，计算表面积栅格化数据；将目标区域划分为多个子区域，将所述地表覆盖栅格化数据及所述表面积栅格化数据按照所述子区域进行划分，将各个子区域对应的地表覆盖栅格化数据及所述表面积栅格化数据同步处理，获得各个子区域对应的各地类表面积；将所有所述子区域对应的各地类表面积汇总，获得所述目标区域对应的各地表覆盖地类的表面积。

作为一种实施方式，所述将目标区域划分为多个子区域，包括：将所述目标区域按照经纬网划分为多个相同大小的子区域。。

作为一种实施方式，所述将各个子区域对应的地表覆盖栅格化数据及所述表面积栅格化数据同步处理，获得各个子区域对应的各地类表面积的步骤，包括：将各个所述子区域对应的地表覆盖栅格化数据及所述表面积栅格化数据进行叠加分析，并逐栅格汇总各个所述子区域对应各地类的表面积，获得各个所述子区域对应的各地类的表面积。

作为一种实施方式，根据所述坡度栅格数据，计算表面积栅格化数据，包括：将所述坡度栅格数据对应的格网面积除以所述坡度栅格数据对应坡度的余弦值获得的数值，作为所述表面积栅格化数据。

作为一种实施方式，所述目标格网大小为根据所述目标区域的大小进行设置。

第二方面，本发明实施例提供了一种地表覆盖表面积栅格化计算装置，所述装置包括：栅格化模块，用于将地表覆盖矢量数据栅格化至目标格网大小的栅格数据，获得地表覆盖栅格化数据；将数字高程模型重采样至所述目标格网大小的栅格数据，获得重采样的数字高程模型数据；计算模块，用于根据所述重采样的数字高程模型数据进行坡度计算，获得坡度栅格数据；根据所述坡度栅格数据，计算表面积栅格化数据；处理模块，用于将目标区域划分为多个子区域，将所述地表覆盖栅格化数据及所述表面积栅格化数据按照所述子区域进行划分，将各个子区域对应的地表覆盖栅格化数据及所述表面积栅格化数据同步处理，获得各个子区域对应的各地类表面积；汇总模块，用于将所有所述子区域对应的各地类表面积汇总，获得所述目标区域对应的各地表覆盖地类的表面积。

作为一种实施方式，所述处理模块，还用于将所述目标区域按照经纬网划分为多个相同大小的子区域。