バッファー解析

筑波大学農林工学系 西田顕郎

筑波大学の周辺に、どのくらい都市があるかを調べてみよう。

まずregionを初期設定に戻す。

d.region -d
d.erase

筑波大学付近を示すラスターマップを作る：

r.circle output=UT coordinate=140.05,36.09 min=0 max=200

そのまわりに、半径5kmから40kmの領域を設定：

r.buffer input=UT output=UT_buff distances=5,10,20,40 units=kilometers

それぞれの領域内の植生を調べてファイルに書き出す：

for i in 2 3 4 5; do
r.mapcalc "t=if(UT_buff==$i,LC_MEJ_MOD12cats,null())"
r.stats t -p -n > UT_${i}.txt
done

水田の割合：

cat UT_*.txt | grep ^253

都市の割合：

cat UT_*.txt | grep ^13