Open Kanahiro opened 3 months ago
csmap-pyのなかでの挙動
process.py
with rasterio.open(input_dem_path) as dem:
...
with rasterio.open(
output_path,
"w",
driver="GTiff",
dtype=rasterio.uint8,
count=4,
width=out_width,
height=out_height,
crs=dem.crs,
transform=transform,
compress="LZW",
) as dst:
crs=dem.crs,
のように、CRSがそのまま通るように書かれている
CRS定義済み+末指定→もとのCRS
もとの通りの挙動(dem.crs
をそのままcrsとして採用)
CRS未定義+未指定→CRS未定義
もとの通りの挙動(dem.crs
をそのままcrsとして採用、つまりCRSなし)
CRS定義済み+指定→指定したCRS DEMが指定しているCRSを指定したものに上書きする(再投影はしない)。
CRS未定義+指定→指定したCRS DEMにCRSが設定されていない場合にCRSを設定する。
ラスターの領域が示されないと(コード中のtranform)、正しい位置を参照することが出来ないため、gdal_edit.pyに倣って
~ファイルを入力した時点で、CRSを取得して、ウィンドウに表示できるのが理想?~ これはQGISプラグインにした時の話だった
参考
gdalwarpに倣って、-te -te_srs、としようとしていましたが、teは「一部を切り抜くオプション」であり、誤解をうむ可能性があるため、やめます。
CRSを上書きするものなので、gdal_edit.pyに倣って:
-a_srs
としましょう。-a_ullr
としましょう(ullrってなんだろうと考えたら、upper left, lower rightですね)reference: https://gdal.org/programs/gdal_edit.html