5mのグリッド作成

[Aki]

QGIS初心者です。
QGIS(3.16または3.22.5)で、5mのグリッド（メッシュ）を作成したいのですが、

・ベクタ→調査ツール→グリッドを作成、

もしくは

・プロセシングツールボックス→ベクタ作成→グリッドを作成

の間隔の設定で、単位を(m)、x方向は0.0000625、y方向は0.00004166…を設定したくても、0.0000625が0.000063に丸められてしまいます。

小数点以下6桁の制限があると思い、単位をkmにして、x:0.062500、y:0.041667で実行すると、処理中に画面が黒くなり、復活したと思ったらQGISが終了しています。

試しに、単位(m)、x:0.000063、y:0.000042でやっても、同じでした。

2台のPCで同じ現象です。windows10、メモリ16MBは2台共通、CPUは、1台がcore-i7(8565U)、もう1台がcore-i7(10875H)です。

原因や解決法について、アドバイスを頂けると助かります。

よろしくお願い致します。

[福岡]

Aki 様

こんにちは

プロセシングツールの「ベクタ作成 － グリッドを作成」で動作を確認してみました。
まず単位を「m」に指定できるということは、CRSに投影系を設定しているようだと思いました。
そうすると「5m」のメッシュを作りたいということでしたら、x方向、y方向の数値は素直に「5」と入力すると良いと思います。

それと余計な心配かもしれませんが、数値からみると「0.15秒」メッシュを作ることがご希望の作業かと思いました。
以前、このことに関連した投稿がありましたのでリンクを張っておきます。
東雲さんのアドバイスが参考になると思います。

QGISによる洪水想定浸水図の作成について
https://groups.google.com/g/qgisshitumon01/c/FEREdYvGu-I

2022年4月28日木曜日 21:35:44 UTC+9 Aki:

[Aki]

福岡様

ありがとうございます。

アドバイス頂いたとおり、5(m)指定でグリッド(EPSG:6677)の作成ができました。

ただ、先に作成済み（知人に作成してもらったもの）の50mのグリッド(EPSG:6668)と、当然ながらずれてしまいます。

このように、座標系のこともまだ理解が浅いため、ご案内頂いた点も含めて勉強し、上手く行かない時には改めて質問させて頂きたいと存じます。

その際には、またよろしくお願いします。

2022年4月29日金曜日 6:57:42 UTC+9 福岡:

[Aki]

福岡様

先日はありがとうございました。
また、中途半端な質問で失礼しました。

ご心配頂いたとおり、緯度経度でメッシュを作成する必要があります。
目的は、5mメッシュに基盤地図の5mメッシュ標高データを持たせて、浸想図の作成です。
基盤地図は直角座標ということなので、DEMポイントを含まないメッシュには、ArcGISで近傍点の地盤高を持たせるようにします。（別途、ArcGISをお持ちの方にお願いします。）
5mのメッシュ作成も知人にお願いできるのですが、できるだけ自力でを目指して、ここで質問させて頂いている次第です。

ということで改めてのご相談なのですが、今現在、
・グリッドタイプ：長方形
・水平方向の間隔：0.0000625？
・垂直方向の間隔：0.00004166…？
　　※↑小数点以下6桁までなので、0.000063、0.000042になってしまう。

・グリッド範囲：同じCRSのレイヤ
・水平方向・垂直方向の重なり：0.000000のまま
・出力グリッドのCRS：EPSG:6668
という状況です。

先日ご案内頂いた「QGISによる洪水想定浸水図の作成について　https://groups.google.com/g/qgisshitumon01/c/FEREdYvGu-I」の質問者の方が「できた」と仰られていましたし、私の設定がおかしいのでしょうが、どうすればいいのかわからず、ご教授頂きければ幸いです。

よろしくお願い致します。

2022年4月29日金曜日 10:53:42 UTC+9 Aki:

[福岡]

Aki 様

こんにちは

私なりの理解では、次のような理由から「グリッドを作成」でメッシュを作らない方が良いと思っています。
確かな機関から提供されている３次メッシュを利用して、5mメッシュを作成してみてはいかがでしょうか？

まずメッシュは、全国標準地域メッシュ３次メッシュ（約 1km）を200等分したものを使うようです。
基準地域メッシュ（３次メッシュ）についてはP22～23、1/200メッシュについてはP24で解説されています。

浸水想定区域図データ電子化ガイドライン（第 3 版）
https://www.mlit.go.jp/river/shishin_guideline/bousai/saigai/tisiki/syozaiti/pdf/e-guideline_1910.pdf

1/200メッシュは、緯度0.15秒、経度0.225秒の区画のようです。
これを小数点表示の度になるように計算してみます。

緯度（1辺）＝ 0.15 / 3600 ＝ 0.000041666・・・（循環小数）
経度（1辺）＝ 0.225 / 3600 ＝ 0.0000625

このように小数点で正確に数値を入力するのには難があります。
さらにQGISの数値丸めの仕様もあるので、東雲さんは適当な３次メッシュを分割するという手法を提案したのだと思います。

３次メッシュはG空間情報センターから提供されていますので、それをお使いになると良いと思います。

全国標準地域メッシュ３次メッシュ（約1㎞四方）
https://www.geospatial.jp/ckan/dataset/biodic-mesh/resource/38bd3651-120e-480f-99cf-7bb89cad7a05?inner_span=True

QGISに読み込んでから関連する場所を選択し、レイヤパネルのレイヤ上で右クリックして「エクスポート － 新規ファイルに選択地物を保存」で書き出してください。
それを「split polygon」というプラグインを利用して、縦横200等分することによって、ご要望の成果が得られると思います。

2022年4月30日土曜日 18:53:43 UTC+9 Aki:

[Aki]

福岡様

詳細な解説を頂き、誠にありがとうございます。

DLした3次メッシュを対象範囲に切り取って保存し、プラグインからsplit polygonをインストール
→[ベクタ]→[split polygon]→[split polygon]

Input Layer：3次メッシュから切り出し保存したレイヤ
Split only selected featuresのチェックをオン（チェックなしだとQGISシャットダウン）
Angleは、とりあえず0°0′0″（次の段階で90°0′0″）
The way of splitting：Equal width：200
でOKすると、レイヤパネルに「splitted」レイヤが出現するのですが、メイン画面には変化が起きない状況です。

どこが間違っているでしょうか。選択や入力値を何通りか変えてみましたが、変化がありません。
何度もすみません。
よろしくお願い致します。

2022年4月30日土曜日 21:09:40 UTC+9 福岡:

[福岡]

Aki 様

こんにちは

思い込みで述べたところもあり、ご迷惑をおかけしました。
実際に手を動かして動作を確認しましたので、お知らせします。

１ 検証用のデータ作成
ダウンロードしたデータから１区画だけを別ファイルに保存して、検証用のデータにしました。

２ 入力するデータが区分数ではなく、角度で入力すべきことが判明
私のデスクトップPCでは「Equal width」を「200」にして実行すると、QGISが倒れてしまいます。
別環境ではどうかとノートPCで試したところ、同様に倒れてしまいます。
もしかしてということで、角度の場合と同桁程度の小さい数を入力してみたところ動作します。
ここで角度で入力すべきだったのか、と思いました。

３ 入力する数値を計算してみる
レイヤのプロパティの「情報」タブの「プロバイダからの情報」という項目の中に「領域」の数値が表示されています。
私の検証用データでは、次のようなものです。

140.1250000000000000,39.7083330000000032 : 140.1374999999999886,39.7166670000000011

緯度、経度で差を取って、200で除して得た数値は次のとおりです。

緯度：0.00004167
経度：0.0000625

４ Split polygonで試してみる
上記数値でSplit polygonは動きましたが、１つ問題が発生しました。
200で分割したポリゴンを別方向に200分割しようとすると、QGISが倒れてしまいます。

５ 代替手順
緯度、経度でそれぞれ分割します。
どちらかのポリゴンをメニューの「ベクタ － ジオメトリツール － ポリゴンを線に変換」でラインにします。
プロセッシングツールの「ベクタオーバーレイ － 線で分割」で、もう片方のポリゴンに対して上記ラインを使って分割します。

以上の手順で、私の環境では分割することができました。
ついでですが、緯度の数値を「0.000041667」でも試してみました。
結果からいうと、分割数が201になったので「0.00004167」で良い様です。

2022年4月30日土曜日 22:31:51 UTC+9 Aki:

[Aki]

福岡様

ありがとうございます。

こちらでも作成できました。

試しで1区画のみですので、次に対象範囲全体のものを作成しようと思います。

この度は、私のしつこい質問に丁寧に答えて下さり、深く感謝しております。

また質問させて頂くことがあると思います。

今後とも、よろしくお願い致します。

ありがとうございました。

2022年5月1日日曜日 8:37:49 UTC+9 福岡: