gis - シェープファイルと重ならない位置で geotiff ファイルをクリッピングする

Question

シェープファイルと重複する geotiff ファイルがあります。tif ファイルの重複部分を切り取るには、次のようにします。

gdalwarp -co compress=deflate -dstnodata 255 -cutline shapefile.shp  original.tif overlap.tif

しかし、交差しない部分を切り取るにはどうすればよいですか? つまり、"overlap.tif"wrtの補数を作成したいのです"original.tif"。

Answer 1

gdal_rasterize を使用して、シェープファイルがファイルと重なっている値を焼き付けることができます。既存のファイルで機能するため、必ずコピーを使用してください。

gdal_rasterize -burn 255 shapefile.shp copy_of_original.tif

これは 255 の値を-a_nodata 255書き込みます。私のバージョンの GDAL では設定が機能しません。実際の nodata 値にする必要がある場合は、後でgdal_translatewithを使用するとうまくいきます-a_nodata 255。

Gdal_rasterize に-iは、シェープファイルを反転する便利なフラグもあります。