グリッドの作成について

492 views
Skip to first unread message

S F

unread,
Feb 3, 2022, 1:26:06 AM2/3/22
to QGIS初心者質問グループ
いつもお世話になっております.
グリッドの作成関係について,2点質問させていただきます.

以前は,業務でArcGISをメインで使用していましたが,最近QGISを使用し始めたものです.

1点目はQGISとArcGISで作成したグリッドの形状の違いについてです.
QGIS上で地図のシェープファイルをもとにグリッドを作成ツールから4 km×4 kmのグリッドを作成したところ,以前ArcGISで作ったものと四角形の形状が異なってしまいました.QGISで作ったものの方がやや横長に表示されます.範囲の四隅を緯度経度などで指定しても同じ症状になるのですが,原因は何でしょうか.形状をそろえることはできますでしょうか.
どちらも,WGS1984で表示されていた地図シェープファイルをkm単位で作業するために,JGD2011の第12系に変換して格子線作成の作業を行いました.

2点目は,1つのグリッド内のラインの距離の算出についてです.
ある地域をグリッドで分け,1つずつのグリッド内に通っているラインの総距離を求めたい(等面積ごとの努力量を知りたい)と思い,以下の作業を行いました.
①ラインのシェープファイルを格子線(グリッド)のポリゴンで交差(intersect)
で分割し,ラインをグリッドで分割
②①でインターセクトされたラインのテーブルを開いてフィールド計算機により長さ($length)の列を追加
③②で作成した長さが追加されたラインのシェープファイルとグリッドのポリゴンを空間結合(集計つき)の交差するで結合し,グリッドのシェープファイルを作成(計算する集計関数は合計としました)

③で作成されたグリッドのポリゴンの1行ごとの長さは,グリッド内のラインの距離(長さ)として算出されていると考えていたのですが,いくつかのグリッドでグリッド内のラインを計測ツール(物差しのアイコン)で手動で計測したところ,2倍以上距離が増えていることもありました.
②のインターセクトの段階では,手動計測の結果と行ごとに計算された長さに違いは見られませんでしたので,③の段階が誤っていると考えられるのですが,原因がわかりません.
すべてのシェープファイルはJGD2011の第12系で作業しています.

PCはWindows10,QGISは3.16.16を使用しています.
2点について,ご教示いただけますと幸いです.どうぞよろしくお願いいたします.

キタ

unread,
Feb 3, 2022, 8:30:35 AM2/3/22
to QGIS初心者質問グループ

こんにちは
1つ目質問ですが、プロジェクトの座標参照系は何になっているでしょうか?
QGISの右下に「EPSG:xxxx」(xは数字)と表示されていると思いますが、4桁の数字はなんですか?
プロジェクトの座標参照系が地理座標系になっていると、地図の表示が横長になります。
プロジェクトの座標参照系もJGD2011の12系(EPSG:6680)にすると、正方形に近く表示されると思います。
プロジェクトの座標参照系をJGD2011の12系に変更しても横長になっている場合には、グリッドを作るときのレイヤの座標参照系が間違っている可能性があります。

キタ
2022年2月3日木曜日 15:26:06 UTC+9 S F:

キタ

unread,
Feb 3, 2022, 8:55:50 AM2/3/22
to QGIS初心者質問グループ
2つ目の質問ですが、ラインレイヤ、グリッドレイヤの両方の座標参照系がJGD2011の12系であるならば、長さが2倍などにあるというのであれば線が同じ場所で重なっているのではないでしょうか?
地物情報の表示で地物をクリックすると、複数の情報が表示されませんか?

地物がぴったり同じ位置で重なっていると仮定したら、次の方法で、グリッド内の線の長さの合計をできます。
1,グリッドを長方形(ポリゴン)で座標参照系JGD2011の12系(EPSG:6680)で作成します。
2,プロセッシングツールの「線をセグメントに分解」で線を屈曲点で分割します。
3,プロセシングツール「重複ジオメトリの削除」で重なっている線を1つにします。
4,プロセッシングツール「線長の合計」でポリゴン内の線の長さを合計します。

ラインが重なっていない場合は、1と4を実行することでポリゴン内の線の長さの合計を算出できます。

キタ

2022年2月3日木曜日 15:26:06 UTC+9 S F:
いつもお世話になっております.

S F

unread,
Feb 3, 2022, 9:42:03 PM2/3/22
to QGIS初心者質問グループ
キタ様

早速のご回答ありがとうございました.

1つ目質問ですが、プロジェクトの座標参照系は何になっているでしょうか?
QGISの右下に「EPSG:xxxx」(xは数字)と表示されていると思いますが、4桁の数字はなんですか?
プロジェクトの座標参照系が地理座標系になっていると、地図の表示が横長になります。
プロジェクトの座標参照系もJGD2011の12系(EPSG:6680)にすると、正方形に近く表示されると思います。
プロジェクトの座標参照系をJGD2011の12系に変更しても横長になっている場合には、グリッドを作るときのレイヤの座標参照系が間違っている可能性があります。

>プロジェクトの座標系がEPSG:4326のまま,作業をしていたようです.プロジェクトの座標座標系を6680に変更したところ問題なく表示されました.今度は,ArcGISの方で作った方が歪んでしまったため,こちらを作成したときに間違っていたのかもしれません.こちらを確認してみたいと思います.ありがとうございます.

2つ目の質問ですが、ラインレイヤ、グリッドレイヤの両方の座標参照系がJGD2011の12系であるならば、長さが2倍などにあるというのであれば線が同じ場所で重なっているのではないでしょうか?
地物情報の表示で地物をクリックすると、複数の情報が表示されませんか?

地物がぴったり同じ位置で重なっていると仮定したら、次の方法で、グリッド内の線の長さの合計をできます。
1,グリッドを長方形(ポリゴン)で座標参照系JGD2011の12系(EPSG:6680)で作成します。
2,プロセッシングツールの「線をセグメントに分解」で線を屈曲点で分割します。
3,プロセシングツール「重複ジオメトリの削除」で重なっている線を1つにします。
4,プロセッシングツール「線長の合計」でポリゴン内の線の長さを合計します。

ラインが重なっていない場合は、1と4を実行することでポリゴン内の線の長さの合計を算出できます。
>ラインが重なっていることはなかったので,1・4の処理をしたところ,無事に似た値が出力されました.ありがとうございました.

2022年2月3日木曜日 22:55:50 UTC+9 キタ:
Reply all
Reply to author
Forward
0 new messages