基于矢量几何的地理空间栅格数据集综述
rasterstats的Python项目详细描述
rasterstats是一个python模块,用于根据矢量几何图形汇总地理空间光栅数据集。 它包括用于zonal statistics和插值的point查询的函数。命令行界面允许 与其他geojson工具的简单互操作性。
文档
有关安装和使用的详细信息,请访问http://pythonhosted.org/rasterstats上的文档。
它是做什么的?
给定一个矢量层和一个光栅带,计算每个矢量几何图形的汇总统计信息。 例如,使用多边形矢量图层和数字高程模型(DEM)光栅,计算 每个多边形的平均高程。
命令行快速启动
命令行与zonalstats和point_query的接口 是读写geojson特性的子命令
$ fio cat polygon.shp | rio zonalstats -r elevation.tif $ fio cat points.shp | rio pointquery -r elevation.tif
请参阅CLI Docs。更多细节。
python快速入门
对于区域统计
>>>fromrasterstatsimportzonal_stats>>>stats=zonal_stats("tests/data/polygons.shp","tests/data/elevation.tif")>>>stats[1].keys()['count','min','max','mean']>>>[f['mean']forfinstats][756.6057470703125,114.660084635416666]
对于点查询
>>>fromrasterstatsimportpoint_query>>>point="POINT(245309 1000064)">>>point_query(point,"tests/data/elevation.tif")[723.9872347624]
问题
找到虫子了吗?通过github问题报告
- 下载再现错误所需的最小可能光栅和矢量数据集的链接
- 再现错误的python代码或命令
- 有关您的环境的信息:python、gdal和numpy的版本以及系统内存