Unity上でNays2DHの計算結果を3D表示する【開発記録 vol.1】

Last updated at 2023-04-27Posted at 2023-04-06

はじめに

表題の通り、iRIC Nays2DHの計算結果をUnity上で3D表示できるようにしよう！
という個人的なプロジェクトの開発記録です。
（進捗報告的な記事なので詳細な手順までは説明しない旨、ご了承ください）

まずは3か月で、iRICのように時間毎の計算結果を切り替えられるようにすることが目標です。

なお、開発記録自体は開発の習慣化&自身へのプレッシャーとして
毎週アップロードする予定です。

Nays2DHの水位と地形をOBJ形式に変換し、Unity上で可視化するところまでは
できています。以下のような感じです。

まず、計算結果が無いと何も始まらないので作ります。
今回は北海道豊平川（kp16.0～17.0）を対象に一定流量を流した計算を実行しました。
（ここは3D関係なく、ただのiRICの操作なので割愛します）

なお、計算に用いる河床のデータはほとんど公開されていないのですが
2023年3月末に北海道開発局が保持するデータが公開されました。
（お陰で色々と開発がしやすくなりました。開発局の方々には感謝です。）

ただし、用途は研究開発を目的とした使用に限られており、
結果を公開する場合は北海道開発局河川管理課への確認が必要である点は注意してください。

ちなみに、今回の開発記録の公開に上記のデータを利用すると
確認が必要になって大変なので、ここでは自前のダミーデータを使用しています。

計算結果の値は以下のどちらかの方法で取得することができます。
（一般的にとっつき易いのはcsvのエクスポートかと）

上記いずれかの方法で、計算格子の座標値、水位、地形（標高）を取得し
Pythonを使ってOBJ形式に変換しています。

OBJ形式を作るモジュールみたいな都合のよいものは無いので
以下のサイトを参考にしながら変換するプログラムを作りました。

以上の流れで作成したプログラムがあるのですが、
ここで解説するとそれだけで2～3記事のボリュームになって（僕が）辛いので
ここではexe化したものを置いておくことにします。

cgns2obj_nays2dh.zip　※Windows11での使用を想定

解凍したフォルダ内のcgns2obj_nays2dh.exeを起動すると、
以下の画面が表示されるので何も考えずに地形と水深の両方で
同じCase1.cgn（iproを解凍すると入っています）を選択して、変換を押してください。

cgns2obj_nays2dh.exeと同じ階層にoutputフォルダが作成され、
中に地形（初期）のOBJファイルと水位のOBJファイルが格納されています。

上述の地形ファイルに貼り付けるテクスチャ画像として
国土地理院の航空写真を用意します。

ここはスマートなやり方が分かっていないが故に忘れやすいので
備忘録としてやや詳細に書きます。

手順の流れとしては

iRIC上で計算格子範囲をshpで出力

　
QGISで読み込み、緯度経度の範囲を確認
（xyzタイルをエクスポートしようとした時に表示される以下のウィンドウを利用する）

　
地理院地図で、指定範囲の画像をダウンロード
https://maps.gsi.go.jp/help/intro/kinolist/1-insatu.html

テキストエディタでland.mtlを作る

land.mtl

newmtl land
map_Kd img20230404163427947.png # ダウンロードしたpngファイル
d 1.0
illum 2

landフォルダを作成し、以下3つのファイルを格納する
- 地形のOBJファイル（land_step=1.obj）
- land.mtl
- ダウンロードしたpngファイル

Unityにてプロジェクトを立ち上げて、作成したlandフォルダや
水位のOBJファイルを読み込めば可視化できます。

ただし水面をそれっぽくしたかったので、今回はレンダリングパイプラインをURPとし、
Unity Asset StoreからSimple Water Shader URPをダウンロードして
水面のマテリアルにあてました。

vol.1は「地形と水面の3DモデルをUnity上で可視化できた」という報告でした。
今回の記事は2か月くらいの作業内容を凝縮してお届けしたのでボリュームが多いですが
次回からは1週間分の作業内容なので少ないボリュームでより詳細な記事になるかと思います。