基盤地図情報で発生する境界線?を消す方法

[yuto otsuka]

お世話になります｡

現在､基盤地図情報の建物データを用いて分析を行っているのですが､

建物データに関わらず基盤地図情報のデータは一定の緯度､経度で境界があり､

そこで別々のポリゴン(地物)として分割されてしまいます｡

その例が添付している画像なのですが､もとは一つの建物なのにも関わらず一定のラインで分割されています｡

分析の都合上､｢一つの建物で一つの地物｣として扱いたいのですが､いい方法はないでしょうか｡

単純にディゾルブしただけではすべてのポリゴン合わせて一つの地物となってしまうため､隣接する建物同士だけディゾルブするような機能があればよいと思っているのですが…

代わりにosmからダウンロードしようとしたのですが､osmではカバーしきれていない建物が意外と多く､地理院地図のデータを使いたかったため質問させていただきました｡

ご回答いただければ幸いです｡

[ありた]

こんにちは

一般にディゾルブする際の連結にはキーとなる属性を指定しますが

基盤地図情報は提供単位をまたぐ地物でも同一キーをもっていません。

単に、隣接したもののみにしたければ、マージ、ディゾルブを行い、

「マルチパートをシングルパートに変換」でシングルポリゴンに分割しては

いかがでしょうか。

[yuto otsuka]

返信ありがとうございます｡概ねうまくいきました｡

しかし､ところどころディゾルブできていないポリゴンがあったのですが､これは仕方ないものなのでしょうか｡

体感9割はうまくディゾルブ→マルチパートへの変換がうまくいっているのですが､のこりの1割､ディゾルブしきれていないポリゴンがあります｡

もしうまくいく方法などあれば､お教えいただけると幸いです｡

2022年10月9日日曜日 1:47:20 UTC+9 ありた:

[遠藤真]

私も以前同様の件で、悩んでいました。

その際に私なりに解決した手順を記します。

下で述べるように完璧ではありませんが、参考まで。

使用した地図は、数値地図2500を使用していますが、同じような処理ができるのでは無いかと思います。

■目的

関係する数値地図の家屋外形（BldA）を選択し、管轄全域の2次メッシュで分断された家屋を接合し、全ての家屋を抽出する。

■なぜ、回りくどい方法になったか。

接している家屋があった場合に結合してしまうため、極力その影響を排除するためメッシュ境界にある家屋のみ結合対象にする。
※このバッファで抽出された家屋が複数接していた場合、結合してしまう可能性がある。

■手順

1．領域のDKG-SHP-BldAを全てマージする。
2．領域の2次メッシュを線に変え、10cmバッファを生成する。（「japan-mesh-tool」を用いてメッシュを生成）
3．1)を2)で交差する家屋を場所による抽出をする。
4．1)を2)で離れている家屋を場所による抽出をする。
5．家屋の内部保証点を生成する。
6．3)を融合する。（基準となる属性：type）
7．マルチパートをシングルパートにする。
8．4）と7）をマージする。

■2次メッシュの生成はjapan-mesh-toolで作成しました。
python main.py 2 -e e0,n0 e1,n1 -d ./
(e0,n0:領域左下、e1,n1:領域右上)

2022年10月9日日曜日 22:11:22 UTC+9 take...@gmail.com: