日本一水部面積の割合が低い都道府県はどこ？

「日本一水部^＊面積の割合が高い都道府県はどこでしょう？」と聞かれたら、日本最大の湖沼面積を持つ琵琶湖（滋賀県）を思い浮かべる方が多いのではないでしょうか。では、日本一水部面積の割合が低い都道府県はどこでしょうか。

^＊河川・湖沼の総称

この答えを ArcGIS for Desktop とデータ製品「ArcGIS データコレクションスタンダードパック 2016 公共地図（以下「公共地図」）」に収録されているデータを用いて調べてみたいと思います。

ArcMap に「公共地図」に収録されている「WATER_SUIBU」フィーチャクラス（水部ポリゴン）と「BOUNDARY_PREF」フィーチャクラス（都道府県界ポリゴン）を ArcMap に追加します。

※「WATER_SUIBU」フィーチャクラスと「BOUNDARY_PREF」フィーチャクラスのソースは数値地図（国土基本情報）です。

ステップ1 : 陸域内の水部ポリゴンに都道府県情報を付与

■各都道府県フィーチャの面積を計算します。

1. 「BOUNDARY_PREF」フィーチャクラスの属性テーブルを開きます。[テーブルオプション] → [フィールドの追加] をクリックします。

2. 任意の名前で新しいフィールド（この例では「AREA_PREF」としています）を追加します。[タイプ] は、「double」を選択し、[OK] をクリックします。

3. 追加したフィールド名で右クリック → [ジオメトリ演算] をクリックします。

4. [プロパティ:] で「面積」が選択されていることを確認したら、出力する面積の単位を [単位:] から選択します。ここでは「平方キロメートル」を選択し、[OK] をクリックします。

※今回の例では便宜上、アルベルス正積円錐図法に投影して面積を計算しています。本例で実際に使用したアルベルス正積円錐図法の座標系ファイルは、こちらからダウンロードできます。

各都道府県フィーチャの面積が計算できました。

■[インターセクト] ツールを用いて、水部ポリゴンを出力します。

5. [ArcToolbox] ウィンドウから、[解析ツール] → [オーバーレイ] → [インターセクト] ツールをダブルクリックして起動します。

6. [入力フィーチャ] に「WATER_SUIBU」フィーチャクラスと「BOUNDARY_PREF」フィーチャクラスを設定し、[出力フィーチャクラス] に任意の名前（この例では「SUIBU_INTERSECT」としています）を設定して [OK] をクリックします。

各水部ポリゴンに都道府県情報を付与することができました。

※実行環境によっては1時間程度の処理時間を要します。

ステップ2 :都道府県ごとに水部ポリゴンの面積の合計値を算出

■各水部フィーチャの面積を計算します。

1. 「SUIBU_INTERSECT」フィーチャクラスの属性テーブルを開きます。[テーブルオプション] → [フィールドの追加] をクリックします。

2. 任意の名前で新しいフィールド（この例では「AREA_SUIBU」としています）を追加します。[タイプ] は、「double」を選択し、[OK] をクリックします。

3. 追加したフィールド名で右クリック → [ジオメトリ演算] をクリックします。

※今回の例では便宜上、アルベルス正積円錐図法に投影して面積を計算しています。

各水部フィーチャの面積が計算できました。

■[要約統計量] ツールを用いて、都道府県ごとに水部ポリゴンの面積の合計値を算出します。

5. [ArcToolbox] ウィンドウから、[解析ツール] → [統計情報] → [要約統計量] を起動します。

6. [入力フィーチャ] に「SUIBU_INTERSECT」フィーチャクラス、 [出力フィーチャクラス] に任意の名前（この例では「SUIBU_INTERSECT_FREQ」としています）を設定して、 [統計フィールド] に「AREA_PREF」フィールド（統計の種類：FIRST）と「AREA_SUIBU」フィールド（統計の種類：SUM）を設定します。また、[ケースフィールド] に「PREF」フィールドを入力して [OK] をクリックします。

都道府県ごとに水部ポリゴンの面積の合計値を算出できました。