CRSを設定するオプションを追加する

[Kanahiro]

gdalwarpに倣って、-te -te_srs、としようとしていましたが、teは「一部を切り抜くオプション」であり、誤解をうむ可能性があるため、やめます。

CRSを上書きするものなので、gdal_edit.pyに倣って:

• CRSの指定は-a_srsとしましょう。
• 領域を、-a_ullrとしましょう（ullrってなんだろうと考えたら、upper left, lower rightですね）

reference: https://gdal.org/programs/gdal_edit.html

[geogra-geogra]

課題整理

現状把握

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を指定するパラメータ -a_srs
• 領域を指定するパラメータ -a_ullr
  とする

追加検討

ファイルを入力した時点で、CRSを取得して、ウィンドウに表示できるのが理想？
これはQGISプラグインにした時の話だった

参考

• https://gdal.org/programs/gdalwarp.html#cmdoption-gdalwarp-te
• https://gdal.org/programs/gdal_edit.html