CRSを設定する

[geogra-geogra]

CRSを設定する

[Kanahiro]

で、ファイルの他にパラメータとしてCRSを渡して付与する機能はcsmap-pyには現状ないです

社内OSSなのでいじってもらって一向にかまわないです
https://github.com/MIERUNE/csmap-py

[geogra-geogra]

ありがとうございます

確かに、再投影と定義が曖昧になっていました。

その仕様でいきたいと思います

明日以降、まずはcsmap-pyの方でissue立ててブランチ切ります

[Kanahiro]

csmap-py自体は、入力データと同じCRSで出力するようにはなっています
入力データにCRS未設定の場合の挙動は…未定義 😢 エラーにはならないと思われ、多分単に位置情報を持たないデータとなる

[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指定するUIが必要になる)
• 再投影する(これまで同様？指定されてない場合どうする？)
• その他

変更するコードの方向性

process.pyのcrs=dem.crs,
の部分を書き換えることになる。
main.pyの部分で、他のパラメータと同様にQGISプラグイン上で動作するように行う

また、選択したCRSを出力する

[geogra-geogra]

@Kanahiro

-crs を加えて指定することで、元のファイルとは異なるCRSに変更できるようにしたいです。

未指定の場合は、もとのCRSが採用されるようにすれば良いと思います。
CRS未定義のもので、特に何も選択しなかった場合については、無理やりCRSをつけたりはしなくて良いかと思います。

CRSの定義/未定義はもとのファイルにCRSが設定されているかどうか
未指定/指定は-crsやQGIS上でCRSを指定するかどうかになります

CRS定義済み+末指定→もとのCRS
CRS定義済み+指定→指定したCRS
CRS未定義+未指定→CRS未定義
CRS未定義+指定→指定したCRS

以上で考えたのですがいかがでしょうか？

[Kanahiro]

未指定の場合は、もとのCRSが採用されるようにすれば良いと思います。

賛成です。従来と挙動が変わらないということなので。

ゆえに、以下の挙動を考えることになります。

CRS定義済み+指定→指定したCRS
CRS未定義+指定→指定したCRS

[Kanahiro]

CRS未定義+指定→指定したCRS

まずこちら。DEMにCRSが設定されていない場合、オプションがあったら設定する、のは正しいように思います。
ただしその場合、EPSGコードだけでなく、そのラスターがどの領域を示しているか（コード中ではtranformのこと）という情報がなければ、結局は正しい位置を参照することが出来ません。なので以下のような追加オプションが必要になるでしょう。

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

GDALを参考にすれば、-te -te_srsというのがあるので：

csmap input.tif output.tif -te_crs EPSG:6668 --te 140.1 40.2 140.3 41.2

となるでしょうか。

[Kanahiro]

CRS定義済み+指定→指定したCRS

一方こちらは、定義済みのCRSがある＝位置情報を持つということです。このデータの投影法を変更するということは、データの再投影（reprojection）を行うことを意味します。

再投影はコストの高い処理です。現状のcsmap-pyは、どれだけ大きいDEMであってもチャンクごとに処理しているので、メモリ効率などは悪くないはずです。
一方、再投影は完成したCS立体図に対して実行しなければならず、現状の変換処理が終わった「あと」に実行されることになります。

以上から、IMO：

• 入力のDEMを予め、望みの投影法にしておいてもらうか
• 出力したCS立体図を、gdalwarpなどで再投影してもらう

ということで良いのではないかと思います。（再投影はこのライブラリ内では行わない）

[Kanahiro]

以上から、私の意見は：

• -teオプション, -te_srsオプションを追加する
• これらが設定されている場合は、元のCRSがなんであれ、常に位置情報を「上書き」する（再投影はしない）

という仕様です。いかがでしょうか。 @geogra-geogra