この投稿は、地震情報アプリ界隈 Advent Calendar 2021 の3日目の記事です。
昨日は、soshi1822さんのリアルタイムに防災情報を配信をするためにでした。

はじめに

GISデータは、防災情報を扱う上で欠かせないデータのひとつです。
そのなかでも、国土地理院や気象庁などが提供するGISデータは有難いことにとても親切で、港湾構造物など非常に面積の小さいポリゴンについてもしっかりとデータ収録が行われています。ただ、府県単位や地方単位で地図を作りたいといった場合、平方メートルオーダーのポリゴンまでは地図に載せなくても良かったりします。
そこで今回は、指定した面積未満のポリゴンを削除する方法を紹介します。

なお以下では、QGIS 3.12.1 を使用しています。

微小ポリゴンの削除

単位の設定

まず、ウィンドウ左上の プロジェクト メニューからプロパティを開きます。
そのなかで真ん中やや上にある 計測 の 面積計測の単位 を確認しておきます。今回は 平方メートル になっていますが、選択により 平方キロメートル などを指定できます。

単位の確認後は、プロパティウィンドウを閉じます。

シェープファイルの読み込み

今回は、気象庁GISデータのうち、地震情報／都道府県等のデータを使います。
データは、気象データ高度利用ポータルサイトにある予報区等GISデータの一覧からダウンロードできます。

ダウンロードしたshpファイルについて、ベクタレイヤの追加を行うと以下のように日本地図が出てきます。
（地図データの出典：気象庁ホームページ）

ここで、読み込んだレイヤ名（上の画像中の左側で選択されている部分）を右クリックしてプロパティを開きます。
気象庁のGISデータ特有の現象かもしれませんが、初期設定のままだと地物データが文字化けしているため、エンコーディングを UTF-8 に設定します。
また、その下にある座標参照系の設定について、 プロジェクトCRS: EPSG:4326 - WGS 84 を選択します。

設定後は、ウィンドウを閉じます。

マルチパートをシングルパートに変換

今回は、海ほたるPAのある人工島を例として、面積の小さいポリゴンを削除することにしましょう。

だめな例

海ほたるにクローズアップして、

黄色い四角のあるメニューを押してから地物を選択、

鉛筆のメニューから編集モードを切り替えて、

赤いゴミ箱メニューを押して地物を削除します。

ズームアウトすると…千葉県が消滅しています。悲しいですね。

これは、海ほたるPAが千葉県に属しているために、地物選択によって千葉県全体が選択されたのが原因です。
したがって、微小ポリゴンの削除をする際は、同じグループに属するポリゴン同士をそれぞれ別のものとして認識させてあげる必要があります。

よい例

まず、 ベクタ メニューの ジオメトリツール から マルチパートをシングルパートに変換… をクリックします。

右下の 実行 を押すと、変換が実行されます。

新しく 出力レイヤ という名前のレイヤが生成されました。このレイヤでは、各ポリゴンが別グループとして扱われています。したがって、海ほたるは海ほたるであり、他の千葉県の領域に影響を及ぼしません。
なお、既存の 地震情報／都道府県等 レイヤはチェックボックスを外して非表示にしておきましょう。また、レイヤ名は各自で分かりやすいように変えても構いませんが、ここではそのままにしておきます。