我有一个森林(值1)和无数据值(255)的森林像素(值0)的栅格(fcc.tif),其分辨率为巴布亚新几内亚的30m。这是一个60459 x 40811像素的大光栅。
我尝试使用gdal_proximity.py (GDALversion2.1.2)来计算到森林破坏像素的距离(以m为单位)。我使用以下命令:
gdal_proximity.py fcc.tif dist_defor.tif -values 0 -ot UInt32 -distunits GEO -use_input_nodata YES -nodata 4294967295计算对于大部分区域都工作得很好,但在具有某种磁盘形状(下图中用红色标记)的大区域上产生错误的零值。
Distance to deforestation in Papua New Guinea
在这一区域,森林被毁的像素是存在的,但它们似乎没有在计算距离时考虑在内。
here中提供了栅格fcc.tif
发布于 2017-10-19 15:03:25
这是由于在大型数据集上进行距离计算时发生的整数溢出。感谢GDAL的开发者和维护者,他们解决了这个问题。
https://stackoverflow.com/questions/46814362
复制相似问题