GRASS GISその1

GRASSは"Geographic Resources Analysis Support System"の略。
最初にお世話になったオープンソースGISソフト。

公式サイトはここ。

もともとはラスタ解析用のソフトだが，最近はベクタ解析も充実しつつあるような感じ。
GUIも用意されているが，各機能は全てコマンドプログラムとなっており，やりたいことはコマンドを組み合わせたシェルスクリプトとして書いて自動化するっていうのが本来の使い方。

Unix流のシェルスクリプトや地形解析やGIS関連の知識をある程度持った人じゃないと使いにくい(使いこなせない)と思う。

最初に仕事で使ったのは，マルチバンドリモートセンシングデータの有効活用ネタを探すこと。

一番使ったのは，平成13年にできた土砂災害対策法関連で某県より発注された傾斜区分図作成。
これは，図化機やレーザープロファイラで取得した地形標高値より生成した10mメッシュ地形標高データから，2度以下と30度以上の傾斜のあるエリアを抽出してポリゴン化する作業。

ちなみに土砂災害防止法とは急傾斜地の崩壊(崖崩れ)，土石流の発生する可能性のおそれのある区域を指定し，都市計画・宅地建物取引業・建築基準・住宅金融公庫各法と連携して，人的物的被害を最小限に抑えて行こうというもの。
一般の人には今後ハザードマップという形で目にする機会があると思う。

傾斜区分図はエリア指定作業を行うための基本的な資料となるもの。

大体以下の様な感じで，1:2,500レベルの図郭100毎程度の量を自動化してました。
そのまま使っても多分動かないですが。

elevation_all_hosei.dem=min(elevation_all.dem[-1,-1],\
                            elevation_all.dem[0,-1], \
                            elevation_all.dem[1,-1], \
                            elevation_all.dem[-1,0], \
                            elevation_all.dem[0,0],  \
                            elevation_all.dem[1,0],  \
                            elevation_all.dem[-1,1], \
                            elevation_all.dem[0,1],  \
                            elevation_all.dem[1,1])
#傾斜計算
slope_ns_all_hosei=atan((elevation_all_hosei.dem[-1,0]-\
                         elevation_all_hosei.dem[1,0])/20)
slope_ew_all_hosei=atan((elevation_all_hosei.dem[0,-1]-\
                         elevation_all_hosei.dem[0,1])/20)
#絶対値変換
slope_ns_abs_all_hosei=abs(slope_ns_all_hosei)
slope_ew_abs_all_hosei=abs(slope_ew_all_hosei)
#最大勾配
slope_max_all_hosei=max(slope_ns_abs_all_hosei,\
                        slope_ew_abs_all_hosei)
#2度以下エリアの属性値を1にセット
deg2_all_hosei=(slope_max_all_hosei <= 2.0)
#30度以上エリアの属性値を1にセット
deg30_all_hosei=(slope_max_all_hosei >= 30.0)