07QGISで遺跡立地分析
1 like857 views
2019年9月24日〜27日開催奈良文化財研究所主催文化財担当者研修「遺跡GIS課程」の資料です。 DEMデータから各種地形指標を算出し、それらの地形指標をもとに遺跡立地可能性にかかるラスタ演算を行います。
![dependencies: [ { src: 'reveal.js-2.6.1/plugin/markdown/marked.js', condition: function() { return !!document.querySelector( '[data-markdown]' ); } }, { src:
'reveal.js-2.6.1/plugin/markdown/markdown.js', condition: function() { return !!document.querySelector( '[data-markdown]' ); } } ]](https://image.slidesharecdn.com/07predictive-190923082927/85/07QGIS-50-320.jpg)
07QGISで遺跡立地分析
- 1. 遺跡立地分析遺跡立地分析 2019文化財担当者研修 遺跡GIS課程 CREATED BY ISHII JUNPEI (CREATED BY ISHII JUNPEI ( ))@ISHIIJUNPEI@ISHIIJUNPEI
- 7. 傾斜方位ラスタ
- 8. 傾斜角度ラスタ
- 10. 日射量を算出日射量を算出
- 13. No. of day of the year
- 14. 実行すると40分くらいかかる
- 23. Proximity
- 26. 「ラスタ」→「ラスタ計算機」
- 27. ラスタ計算機
- 30. Data/Rule/Aspect_Reclass.txt に上記の計算式があ るのでコピー・ペーストで使用する。 #Aspectレイヤのバンド1 "Aspect@1" # 真（0より大きい）なら計算機は「1」を返し、偽なら「0」を返す "(Aspect@1">0) # 0より大きく45以下の値は「1」を、それ以外はすべて0 ("Aspect@1"<0)*("Aspect@1"<=45) # 「0以上45以下」という条件を満たすピクセルには「10」が代入される ("Aspect@1"verb|<|0)*("Aspect@1"<=45)*10 # 同様に45〜135（北）では20が代入され、135〜225（西）では30が代入され、 # 225〜315（南）では40が代入され、315〜（東）は10が代入される # 計算機が「1」を返す項は一つしかないので、全部の項を足し合わせると # 真となる項の数字だけが該当するピクセルに代入される
- 31. カテゴリ化された傾斜方位ラスタ
- 35. 標高評価式 # 標高5m未満=1、それ以外は0 ("DEM_utm54@1"<5) # 10を乗じているため、最初の括弧が「1」なら標高4mの地点は「1*4*10」で40が出 # 同様に、標高1mなら10が出力 ("DEM_utm54@1"<5)*("DEM_utm54@1"*10) # 標高が5m以上=1を、それ以外は0 ("DEM_utm54@1">=5) # 500を標高で除す。標高5mで最高点の100点が出力され、標高500mでは1が出力 ("DEM_utm54@1">=5)*(500/"DEM_utm54@1")
- 37. 傾斜の評価式 # 傾斜が0度の場合 100-0/45*100 = 100 # 傾斜が45度の場合 100-45/45*100 = 0 # 傾斜が30度の場合 100-30/45*100 = 33.33
- 40. 河川からの距離河川からの距離 河川からの距離は近いほうが遺跡立地に有利そうだ が あまりにも河川に近いところは不適切ではないか 河川から100mを最高得点とし、近くても遠くても 点数が下がるように配点
- 41. 河川からの距離評価式 ("WL_buffer@1"<50)*("WL_buffer@1"/50*100)+ ("WL_buffer@1">=50)*(50/"WL_buffer@1"*100)
- 45. 傾斜方位の評価式 # 傾斜方位が40（北）の場合に100点を付与 ("Aspect_reclass@1"=40)*100 # 傾斜方位が10（東）の場合に80点を付与 ("Aspect_reclass@1"=10)*80 # 傾斜方位が30（西）の場合に70点を付与 ("Aspect_reclass@1"=30)*70 # 傾斜方位が20（北）の場合に50点を付与 ("Aspect_reclass@1"=20)*50
- 46. すべての計算式すべての計算式 Data/Rule/Predictive_calc.txtの計算式をコピー・ ペースト # 標高 ("DEM_utm54@1"<5)*("DEM_utm54@1"*10)+ ("DEM_utm54@1">=5)*(500/"DEM_utm54@1")+ # 傾斜 ("Slope@1<=45)*(100-"Slope@1"/45*100)+ # 日射量 10*("Irradiation@1"-7067.65)/689.64+50+ # 河川からの距離 ("WL_buffer@1"<50)*("WL_buffer@1"/50*100)+ ("WL_buffer@1">=50)*(50/"WL_buffer@1"*100)+ # 傾斜方位 ("Aspect_reclass@1"=40)*100+("Aspect_reclass@1"=10)*80+ ("Aspect_reclass@1"=30)*70+("Aspect_reclass@1"=20)*50
- 48. 遺跡立地評価ラスタ
- 49. まとめまとめ 実際に遺跡予測地図を作成する場合には、それぞれ の地形指標と遺跡立地の関係について実証的な評価 が必要 統計モデルへの当てはめ 各地形指標の係数について「あてはまりの良さ」を 評価
- 50. dependencies: [ { src: 'reveal.js-2.6.1/plugin/markdown/marked.js', condition: function() { return !!document.querySelector( '[data-markdown]' ); } }, { src: 'reveal.js-2.6.1/plugin/markdown/markdown.js', condition: function() { return !!document.querySelector( '[data-markdown]' ); } } ]