パプアニューギニアの分解能が30mのデータ値(255)のない森林(値1)と森林破壊されたピクセル(値0)のラスタ(fcc.tif
)を持っています。これは、60459 x 40811ピクセルの大きなラスタです。gdal_proximityは大きなラスタファイルにゼロを生成します
私は、森林破壊されたピクセルまでの距離(m)を計算するためにgdal_proximity.py(GDALバージョン2.1.2)を使用しようとしています。私は、次のコマンドを使用しています:
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