メッシュ・レイヤにポイント(位置)と凸包で作成したバッファを結合したい

[arakawa.sh...@gmail.com]

TO：福岡さま、　ご助力いただいております皆様

大変お世話になっております。　新川(あらかわ)です。

前回から引き続き、と言いますか度々ご質問をいたします。

Japanメッシュにてグリッド(shape)ファイルを作成し、それぞれのポイントとバッファが

どのグリッドに属しているかを処理(メッシュ・レイヤに結合)したいのですが、何か良い

方法はございますでしょうか。

※実際は、ポイント数やグループ(バッファ)も増えてしまい、重なってしまう状況です。

お忙しい中、申し訳ございません。　よろしくお願いいたします。

[福岡]

新川　様

こんばんは

メッシュもポイントもデータ数が多く、以前の方法を応用するのは手間が掛かりすぎるもですね？
それではメッシュがポリゴンになっていることが前提ですが、「Point sampling tool」というプラグインを使うと良いと思います。
これを使えば、ポイントにメッシュの属性を取り込むことができます。

日々是戯言・Geo出張版－QGISを使ってポリゴンの属性をポイントに読み込む－
http://d.hatena.ne.jp/wata909/touch/20120109/p2

Lutetiums－QGISのPoint sampling tool プラグインを使うときの注意－
http://lutetiums.koushijima.com/gis/qgis20151214

ちなみにメッシュをラインで作っていたときのために、地域メッシュを簡単に作るための「JapanMesh」というプラグインを紹介しておきます。

dayjournal blog－QGIS #033 – 地域メッシュ作成「JapanMesh」－
http://day-journal.com/blog/qgis-033-%E5%9C%B0%E5%9F%9F%E3%83%A1%E3%83%83%E3%82%B7%E3%83%A5%E4%BD%9C%E6%88%90%E3%80%8Cjapanmesh%E3%80%8D/

[きた]

新川さん

一対多でポイントとポリゴンなどの関係を確認したい場合、多の地物に一の属性データを結合することになります。

新川さんの例だと、ポイントの属性データにメッシュの属性データを結合して、メッシュポリゴンの属性に凸包ポリゴンの属性データを結合することになります。

福岡さんからは「Point sampling toolプラグイン」を紹介されましたが、

ベクタメニュー→「データマネジメントツール」→「場所で属性を結合する」でも行うことができます。

重なっている地物の属性データを、対象レイヤの属性データに結合することができます。

ポイントレイヤとメッシュレイヤの場合、「対象ベクタレイヤ」にポイントレイヤを選択し、「結合するベクタレイヤ」にメッシュレイヤを選択します。

2つ以上のレイヤが重なっている場合（一対多でなく多対多の場合）には、はじめに見つけた地物の属性データを結合します。

一対多の関係のレイヤの場合、リレーションを作ることもできます。

メニュー「プロジェクト」→「プロジェクトプロパティ」→「リレーション」で、

「リレーションの追加」をクリックして、「参照レイヤ（子）」にポイントレイヤ、「参照レイヤ（親）」にメッシュレイヤを選択し、「参照フィールド」をメッシュ番号など親レイヤでユニークな値を選択します。

親レイヤの属性データを開いて、フォームビューにすると、そのデータに関係する子レイヤが表示できます。

子レイヤはテーブルビューにすると見やすいと思います。

試してみてください。

きた

[福岡]

きた　様

投稿ありがとうございます。

「場所で属性を結合する」は、全然念頭にありませんでした。
ベクタ同士だと、こちらの方が使い勝手が良さそうな感じです。

勉強になりました。m(__)m

[arakawa.sh...@gmail.com]

TO：福岡さま

おはようございます。　新川(あらかわ)です。

返信ありがとうございます。

Point Samoling Toolを試してみます。　Japanメッシュのリンクもありがとうございます。

※実は、Japanメッシュで作成した属性テーブルを見て「何だこりゃ」ってなっております。

また、お助けをお願いするとおもいますが、よろしくお願いいたします。

[arakawa.sh...@gmail.com]

TO：きた　さま

初めまして、新川(あらかわ)です。

ご連絡ありがとうございます。

リレーションなるものがあるとは、処理の方法が広がるので

すごいですね。　早速試してみます。

以後、よろしくお願いします。

[arakawa.sh...@gmail.com]

TO：きた　さま

毎度お世話になっております。　新川(あらかわ)です。

一対多の関係レイヤを行う予定ですので、リレーションを試してみましたが、上手くできませんでした。

　教えていただきました、メッシュ・レイヤを親としてポイント・レイヤを子としております。　処理が完了しましてメッシュ(親)レイヤを

開きますとポイント情報が反映されない状況です。　※画面添付いたします。

上手く行けば、画面の黄色く選択しているメッシュ部分のポイントデータが子のテーブルに反映されるかと思っております。

情報が少なくてすいませんが、ご教授お願いいたします。

宜しくお願い致します。

[きた]

新川様

親レイヤと子レイヤには共通するキーフィールドが必要です。
多分、まだ属性データを結合する前のデータでリレーションしてないでしょうか?
福岡さんの紹介されていたプラグインか、私の紹介した場所で属性を結合するを実行してから、リレーションを設定してください。

[arakawa.sh...@gmail.com]

TO：きた　さま

尾根輪になっております。　新川(あらかわ)です。

すいません「属性を結合」を行って、再度ためします。

ありがとうございます。

[arakawa.sh...@gmail.com]

TO：きた　さま

毎度お世話になっております。　新川(あらかわ)です。

「属性の結合」にてポイント・レイヤにメッシュ・レイヤのIDを関連付けましたら

上手くリレーションできました。　ありがとうございます。

※画面を添付いたします。

メッシュを選択すると、そこに含まれるポイントがきちんと表示されました。

※さらに凸包で作成したバッファ・レイヤも併せてみます。

ご対応に感謝いたします。　今後もよろしくお願いいたします。

[arakawa.sh...@gmail.com]

TO：きた　さま、　福岡さま

毎度お世話になっております。　新川(あらかわ)です。

「属性の結合」について教えてほしいのですが、ポイントにメッシュのIDを反映させてリレーショナルを行う要領で

凸法で作成したバッファもメシュに結合(ID付け)を行おうとした時、一部のメュシュに重なるバッファが存在してしまい

最初に出てくるバッファIDのみが結合されてしまいます。

複数のグループ(バッファ)が重なる場合の「属性の結合」は、出きるのでしょうか。　もしくは、別の方法でしょうか？

※比較の為、メュシュの２か所を試した画面を添付いたします。

黄色く選択したメッシュの部分には、２つのグループが存在しているのですが、片側を結合してしまします。

結合のための「属性カラム」を複数作成するとか、いろいろ考えたのですが、ハマってしまいました。

ちなみに、お気づきかとおもいますが・・・メッシュが斜めになっているのはCRS設定で戦っておりました。

福岡さまよりジャパンメッシュのリンクなど勉強になりました。奥が深いなあ。。。

何度もお手数おかけいたしますが、よろしくお願いいたします。

[福岡]

新川　様

こんばんは

CRSですが、属性テーブルをみると結構細かいメッシュになっているようなので、平面直角座標系が良い様な気がします。
（距離や面積などを確認しているようなので、投影系の座標が良いと思いました）
19の地域に分かれていますので、該当する地域を確認してください。

わかりやすい平面直角座標系－19の平面直角座標系－
http://www.gsi.go.jp/sokuchikijun/jpc.html

使うCRSが決まったら、全てのレイヤとプロジェクトのCRSをそのCRSに揃えて、念のためオンザフライをオフにします。
（分析などを行うときは、CRSを揃えてオンザフライをオフにするのを習慣にしておくと、トラブルを未然に防げます）

次にポイントやメッシュが複数のグループに含まれているという問題です。
これは最終的にどのようにデータを整理したいか、ということによって解決策が違うと思います。
例えばポイントが属するメッシュとグループが分かるように整理したいということであれば、ポイントの属性テーブルで整理すると良いと思います。
まず、メッシュについてはきたさんが説明したような手法で１つのフィールドにメッシュ番号を使って整理します。
一方、グループに関してはグループ数だけフィールドを用意し、フィールド名をグループ名にしてタイプを整数値にします。
以前説明した方法（調査ツールとフィールド計算機を使う方法）の応用で、フィールド計算機では「"名称"」の替わりに「1」として実行します。
これでグループに含まれているところだけ「1」とフラグを立てることができます。

とりあえず思いついたものを投稿いたしますので、最終目的が違っていたらごめんなさい。
参考になるような内容であれば良いなと思っております。

[arakawa.sh...@gmail.com]

TO：福岡さま、　きた　さま

ご返信ありがとうございます。　新川(あらかわ)です。

　以前教えていただきました「調査ツール」と「計算機ツール」でメッシュレイヤにフラグ立てる

処理につきまして関連付けが出来るのですが、広範囲でグループ数が1,000個以上とか、種類が確定

していない場合では厳しいかと思います。

重なるバッファ(ポリゴン)の処理を「解析ツール」やプラグイン等で行い、属性テーブルに書き出せる

様な物がございますでしょうか？

面倒な質問ですいません、宜しくお願いいたします。

[福岡]

新川　様

なかなかの難題ですねー。
調べて試したりしているのですが、なかなかうまくいきません。

Adding attributes to datapoints from overlapping polygons in qgis
https://stackoverflow.com/questions/32394277/adding-attributes-to-datapoints-from-overlapping-polygons-in-qgis

最初の回答にある「"Clip points with polygons" tool」は、おそらくSAGAのプロセッシングツールだと思うのですが、
ポリゴンの縁辺部で重なっているポイントが抜けたりします。
２番目の回答にあるプロセッシングツールは、日本語表記でどれなのか特定できず、試せていません。

中途半端な情報ですが、他に検討なさっている方のためにも情報は多い方が良いと思いますので投稿いたします。

[arakawa.sh...@gmail.com]

TO：福岡さま

大変、お世話になっております。　新川(あらかわ)です。

お手数をおかけして、すいません。　ありがとうございます。

色々、調べたり試したりしてみます。　進展がありましたら投稿したいと

思いますので、今後ともよろしくお願いいたします。

[福岡]

新川　様

Excelで関係を整理することまではできましたのでお知らせします。

テスト用として用意したデータです。
後で判別しやすいように、大きくポリゴンを重ねています。
使うツールはプロセッシングツール － QGISジオアルゴリズム － ベクタオーバーレイツール － 交差　です。

レイヤの交差の右端にある「回転矢印」アイコンはクリックして有効にしてください。
これを有効にすることによって、ポリゴンの個数分の繰り返し処理が行われます。

レイヤパネルに一時レイヤが追加されますが、不要なので削除してください。

出力先として指定したフォルダ内に、ポリゴン別に含まれているポイントを抽出したShapeファイルが生成されました。

生成されたShapeファイルを統合します。
統合は、ベクタ － データマネジメントツール － 複数のシェープファイルを１つに統合する で行います。

出力されたShapeファイルの属性テーブルです。

Excelに貼り付けて、PointのIDで整理しました。ポリゴンのIDで整理することも可能です。

[きた]

新川様

きたです。

メッシュポリゴンと、凸包ポリゴンの数にもよりますが、数百程度であれば私なら少し手作業を入れます。

１．凸包ポリゴンの属性データには、IDなどの固有の番号が記録されているとします。

２．前に紹介した「場所で属性を結合する」で、対象ベクタレイヤにグリッドレイヤ、結合するベクタレイヤに凸包レイヤを指定します。

３．「属性の集計方法」で「交差する全ての・・・」にチェックを付けて、「最小」「最大」にチェックを付けます。

４．実行すると、新しいメッシュレイヤが出来ます。

５．属性テーブルを開き、「count」というフィールドを探して、降順に並び替えます。「count」フィールドに２以上が入力されているメッシュが、凸包ポリゴンが２つ以上重なっているメッシュです。

６．「count」に２が入力されている場合、最大と最小がその凸包ポリゴンのIDになります。

７．「count」が３以上のメッシュは、そのメッシュにズームして、地図で確認し、３番目のデータを入力できるフィールドを作成して、自分でデータを入力します。

数百くらいなら、この方が早いかなと思います。

参考まで。

きた

[arakawa.sh...@gmail.com]

TO：福岡さま、　

　　きた　さま、

　　サポートされてる一同　さま

ご連絡遅くなりましてすいません、大変お世話になっております。　新川(あらかわ)です。

　福岡さま、確認の為データまで仕込んでいただきまして、本当にありがとうございます。

きた様、色々とアドバイスありがとうございます。　確認しました内容をお知らせ致します。

また、私が使用しているWindows版　QGIS2.18.5(バージョン1個前)の為なのかメニューや設定画面が

違うようなので、手順が間違っているかもと不安に思いキャプチャを添付いたします。

こんな感じかなあーと思いつつ、テストの凸包グループが２個しかないので、グループを少し増やして試してみます。

今後もお世話になると思いますが、、、お世話になりますので宜しくお願いします。

以上、お知らせ致します。

[きた]

新川様

2.18では、ベクタメニューの一部が、プロセッシングツールと統合されたので、画面が違ってしまいましたね。

私の紹介した方法の2.18版が、新川さんの行った方法で正しいです。

とりあえず解決したようで良かったです。

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

きた