gdal_proximity在大型光栅文件上生成零点

Question

我有一个森林（值为1）的栅格（fcc.tif）和去掉像素（值为0）的数据值为255的巴布亚新几内亚分辨率为30m的像素（255）。这是一个60459 x 40811像素的大光栅。

我尝试使用gdal_proximity.py（GDAL版本2.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

在此区域中，砍伐像素都存在，但它们似乎没有被考虑用于计算距离。

光栅fcc.tif可here

Answer 1

这是由于一个整数溢出对大数据集距离的计算过程中发生的。已按照https://trac.osgeo.org/gdal/ticket/7102解决问题。感谢解决此问题的GDAL的开发人员和维护人员。