平均傾斜値の取得について

[熊沢]

お世話になります。

DEMファイルからポリゴンの平均傾斜の値を取得したいと思い、

地域統計を使用して取得したのですが、その値がちょっとおかしい？ようなのです。

手順としては、基盤地図情報サイトから10m標高データをDLし、

エコリスさんのソフトでDEMを作成しました（2＝平面直角座標系、系＝11）。

それを元にラスタ-地形解析-傾斜で傾斜区分図を作成。

地域統計でその区分図を選択。

傾斜値を追加したいポリゴンレイヤーにちゃんとフィールドは追加されます。

しかし、その平均傾斜（mean）の値が全て１～４の間になっています。

どう考えてももっと傾斜がきつい場所もあるのですが…？？

測地系はJGD2011の11系（EPSG6679）で、

DEMもポリゴンファイルもプロジェクトも揃えています。

QGISは2.14.18です。

お知恵を拝借できれば幸いです。よろしくお願いします。

[福岡]

熊沢 様

こんにちは

すぐに原因が思いつかないで申し訳ありません。m(__)m

ただ、傾斜のラスタが正しく出力されているかは確認した方が良いと思います。

プラグインの「Value Tool」を使うとラスタの値が確認できますので、試してみてはどうでしょうか？

[きた]

熊沢様

傾斜図が正しくできていない場合、

もしかしたら、JGD2011にしていることが問題かもしれません。

標高DEMを名前をつけて保存でJGD2000の11系で保存してから、

傾斜図を作るとうまくいくかもしれません。

きた

[熊沢]

>福岡様、きた様

早速のご返信、どうもありがとうございました。

DEMを最初から作り直し、JGD2000、2453のまま傾斜区分図を作成し、

プロジェクトCRSも2453（OTF)にして、Value Toolで確認してみました。

画像は地図上の469上にカーソルを置いた図です。標高が微妙に違いますが、誤差範囲内？

傾斜はやっぱり異常に小さい数字になります。

[福岡]

熊沢　様

思いつきで申し訳ありませんが、原因として思いつくことを２つほど書き並べてみます。

予想が外れていても、ご容赦くださいますようお願いいたします。m(__)m

１　DEMに欠損がある。

　欠損が確実に無いのは、10mメッシュ（地形図の等高線）だけだったと思います。

　5mメッシュでは計測範囲が限られていたり、河川・湖沼などはレーザーで計測できないため欠損している場合があります。

　もし5mメッシュを使いたくて一部にデータ欠損がある場合には、欠損値に10mメッシュの値を上書きするような処理が必要かと思います。

　手順としては以前の投稿が参考になると思いますので、そちらをご覧ください。

　データ無しの水面の標高を基盤地図情報数値標高モデルのGeotiffより推察して埋めたい

　https://groups.google.com/forum/#!topic/qgisshitumon01/mTnnGJGk0-U

２　CRSが揃っていない

　地理院地図を表示するためオンザフライを有効にしているようですが、オンザフライをオフにして確認してください。

　解析ではCRSが揃っていないと不具合が起こる可能性があります。オフにすると、解析の際のデータ位置を確認できます。

[熊沢]

https://groups.google.com/forum/#!topic/qgisshitumon01/FFm_hpyv7_0

上記にて福岡様が載せてくださった傾斜区分図の作成手順を参考に、

今一度、初心に戻ってやってみました。

すると、「Zファクタ」の数値がなぜか10になっていたので、1.0に戻して作成、

そして地域統計を実行してみると、まともな数字が出てきました。

傾斜を強調しようとして、適当な数字を入れたままになっていたようです（汗汗。

ちなみにJGD2011、OTFでも問題なく出力されました。

標高データは10mです。

念のためDLし直そうとしたら基盤地図情報サイトにアクセスできませんでした。

こんな感じのオチで何とも微妙ですが、丁寧に対応してくださった福岡様、きた様、

本当にどうもありがとうございました。