2012年6月7日木曜日

花崗岩の分布図を作成したい

20万分の1日本シームレス地質図から花崗岩類を抽出して花崗岩類の分布図を作成する方法です.

特定の岩体を抽出する方法は,下記の通りです.
http://geologist-gis.blogspot.jp/2010/08/blog-post.html



1.ArcMapでシームレス地質図のポリゴンレイヤーとラインレイヤーを表示します。そして、ArcToolBox>>>解析ツール>>>抽出>>>選択を選びます。




2.入力フィーチャーで、地質図のポリゴンレイヤー(all_japan_d_a)を選択します。出力フィーチャークラスは、任意の出力ファイル名を入力 します。条件式(オプション)の右のボタンをクリックして、検索条件設定を指定します。ここでは、属性テーブルのうち、 ”Legend_J"をダブルク リックして、”=”ボタンを押し、個別値を取得ボタンを押してでてきたリストから、表示したい岩体名を選びました。今回は、”暁新世-前期始新世 (PG1)の花崗岩”を指定しています。

"Legend_J" = '暁新世-前期始新世(PG1)の花崗岩'

この条件式を組み合わせることで、複数の岩体を選択することが可能です。凡例番号がわかっていれば、”New No” = 52のように番号でも指定もできます。

3.OKを押して実行すると、その属性を持つ岩体だけのレイヤー(all_japan_Select2)ができます。この岩体だけのレイヤーとラインのレイヤーを表示すると下記のようになります。

4.岩体のレイヤーを右クリックして、属性テーブルを表示すると全部で297の岩体(ポリゴン)があることがわかります。


5.各レイヤーをダブルクリックすると、該当する岩体が拡大表示されます。下記では、北海道日高地域の暁新世-前期始新世(PG1)の花崗岩の分布が表示されています。





6.ここでは,シームレス地質図から花崗岩類すべてを抽出したいので,下記のように条件設定を行います.シームレス地質図の凡例表をみると,花崗岩類は,1260番〜1540番ですので,下記のように設定します.



7.抽出したgraniteレイヤのプロパティを開いて,シンボルタブ からインポートをクリックし,all_japan_d_aを選択すると下記のように色が付きます.



8.ラインについても海岸線のみを抽出すれば,下記のような花崗岩の分布図が完成します.
 マップのエクスポートから,PNGやKML, AI形式等で出力すればさまざまな用途に使えます.