2016年11月20日

QGISで「CS立体図」の「立体図」を作ってみた

 「CS立体図」という地図をご存知でしょうか?

 CS立体図は、傾斜図と曲率図と標高図を組み合わせた地図で、長野県林業総合センターで開発した、微地形を判読しやすくされた立体地図です。


・地形判読を行いやすくする立体図法

http://www.pref.nagano.lg.jp/ringyosogo/documents/iku_cs.pdf

・数値地形データを用いた「微地形図」の作成方法

http://www.pref.nagano.lg.jp/ringyosogo/seika/documents/bichikei.pdf


私も最近知ったのですが、林業の路網選定などには、非常に役に立つ地図のようです。


 長野県林業総合センターでは、ArcGISでの作成方法が紹介されています。この中では、「立体図」と「曲率図」を作成して重ね合わせる手法が紹介されていますが、「立体図」だけでもかなり立体的に見えるので、ここではQGISで作る、CS立体図のうちの立体図の作成方法を紹介します。



 Image 2016_11_20_205806.png

  図:ニセコアンヌプリ周辺の立体図


※今回作成する立体図の、曲率については、「Profile Curvature」という計算方法を使用していますが、これが本当に正しいのか、著者はわかっていません。もし他にいい方法があれば情報をください。

「Profile Curvature」で作成した地図は、等高線上の地形のシワが強調されて表現されます。




1. 標高DEMを準備する


  • 標高DEMを基盤地図情報などからダウンロードします。株式会社エコリスの変換ツールを使って標高DEM(Tiffファイル)を作成します。
  • 作成したファイル(merge.tif)は、日本語を含まないフォルダ、ファイル名にして保存します。



2. QGISのレイヤに標高DEMを追加する


  • 標高DEMをQGISに追加します。標高DEMのファイルをマップキャンバスにドラッグ&ドロップするとレイヤに追加することができます。

Image 2016_11_20_214521.png 

  図:レイヤに追加した標高DEM



3. プロセッシングツールで傾斜と曲率を計算


  • プロセッシングツールを使って、標高DEMから傾斜図と曲率図を作成します。メニューに「プロセッシング」が表示されていない場合は、「プラグイン」→「プラグインの管理とインストール」で「Processing」を有効にしてください。
  • メニュー「プロセッシング」→「ツールボックス」を選択して、プロセッシングツールボックスを表示します。
  • 「SAGA」→「Terrain Analysis - Morphometry」→「Slope aspect Curvature」を選択します。
  • 表示されたダイアログの「Elevation」に標高DEMレイヤを選択します。「Method」「Slope Units」「Aspect Units」は、デフォルト設定で構いません。
  • このコマンドは、たくさんのレイヤを作成する事ができますが、必要なのは「Slope」「Profile Curvature」の2つだけですので、それ以外の「アルゴリズムの実行後に出力ファイルを開く」のチェックを外します。
  • 「Run」ボタンをクリックすると、処理が開始され、新しいレイヤ「Slope」「Profile Curvature」が作成されます。

Image 2016_11_20_214422.png 

  図:「Slope aspect Curvature」のダイアログ



4. レイヤを名前をつけて保存する


  • 作成されたレイヤは、一時レイヤに作成されます。そのため、QGISを終了したら削除されてしまいます。
  • レイヤを右クリックして、「名前をつけて保存する」を選択します。
  • 「出力モード」を「生データ」にします。
  • 「形式」を「GTiff」にします。
  • 「パス」の「参照」ボタンをクリックして、ファイルの保存先を設定します。
  • 座標参照形を設定する必要がある場合は、「CRS」で選択します。
  • 「OK」ボタンをクリックすると、ファイルが保存されて、レイヤに追加されます。
  • 元のレイヤ(Slope、Profile Curvature)は削除します。

Image 2016_11_20_220440.png 

  図:「名前をつけて保存する」のダイアログ



5. 傾斜レイヤのスタイルを設定する


  • 傾斜レイヤのプロパティを開いて、「スタイル」を選択します。
  • 「レンダータイプ」を「単バンド疑似カラー」にします。
  • 「色」を「Reds」にします。自動的にリストに数値が分類されます。
  • 「カラーレンダリング」の「混合モード」を「乗算」にします。
  • 「OK」ボタンをクリックするとスタイルの設定ができます。

Image 2016_11_20_222358.png  Image 2016_11_20_222446.png

  図:傾斜レイヤのプロパティのスタイル設定



6. 曲率レイヤのスタイルを設定する


  • 曲率レイヤのプロパティを開いて、「スタイル」を選択します。
  • 「レンダータイプ」を「単バンド疑似カラー」にします。
  • 「色」を「Blues」にします。自動的にリストに数値が分類されます。
  • 今回は、マイナス側を青にしたいので、「反転」にチェックを付けます。
  • 「カラーレンダリング」の「混合モード」を「乗算」にします。
  • 「OK」ボタンをクリックするとスタイルの設定ができます。

Image 2016_11_20_222820.png Image 2016_11_20_222856.png 

  図:曲率レイヤのプロパティのスタイル設定




7. 必要であれば、地理院地図などを背景に追加する


 必要であれば、地理院地図や等高線などを地図に追加すると、地図がわかりやすくなるかもしれません。

 長野県林業総合センターの資料では、標高DEMも重ねていますが、標高DEMは重ねなくても立体的でわかりやすいと思います。


Image 2016_11_20_223435.png 

  図:地理院地図を重ねた立体図



 CS立体図は、赤色立体地図よりも地すべりなどの地形が見やすいように感じます。

 長野県では、1mメッシュDEMを利用していますので、かなりの微地形を確認することができるようですが、基盤地図情報の10mメッシュDEMでも十分地形判読を行うことができます。

 地質図や地すべり位置図などを重ねることで、更に地形判読を行いやすくなると思います。


〜2016年11月24日追記〜

 長野県林業総合センターのCS立体図の中の人から、ご連絡をいただきました。

CS立体図の作成方法は、オープンなので自由に紹介して使ってほしいとのことです。CS立体図の普及によって、林業がもっと発展すれば嬉しいとのことでした。大変素晴らしい考え方です。

 また、長野県では全国の10mメッシュDEMを使って、CS立体図を作成済みとのことです。そのうち、自由にダウンロードできるようになると思いますので、非常に楽しみです。


 ダウンロードできるようになれば、また紹介したいと思います。


 

ラベル:QGIS CS立体図
posted by kouichi at 22:47| Comment(0) | QGIS | このブログの読者になる | 更新情報をチェックする

2016年11月07日

赤色立体地図の作成方法の公開停止

 本ブログでは、QGISにおける赤色立体地図の作成方法を公開していましたが、この度赤色立体地図の生みの親である千葉達朗氏から、公開をやめてほしい旨の連絡をいただきました。

 その為、赤色立体地図の作成方法の公開についてはやめさせて頂きます。

 個人的に研究する分には、問題ないと思いますので、QGISの操作方法について確認したい方には、個別に対応させていただきたいと思います。

質問のあるかたは、「qgisshitumonアットgmail.com(アットは@に置き換えてください)」にメールを頂ければ、返答します。

ご迷惑をおかけして申し訳ありません。

今後も「森林土木memo」を、よろしくお願いします。

posted by kouichi at 21:18| Comment(0) | QGIS | このブログの読者になる | 更新情報をチェックする

2016年11月03日

林業QGIS本の出版が遅れてます


 昨年、現代林業で連載していた「QGISを林業で使う」を、単行本として出版するということをお知らせしていました。以前の記事では、10月ころに出版できるのではと書いていたかもしれませんが、出版が大幅に遅れて、現在「初校」という著者がチェックする段階になっています。それも11月段階で半分しかできておらず、全部の初校が出るのは、12月末か1月はじめになりそうだとのことです。


 IMG_20161103_151614.jpg

  写真:初校の原稿


 本として発売されるのは、2月以降になるかもしれません。

以前の記事を見て、楽しみにしていた方がいらっしゃいましたら、申し訳ありません。


 何日か前に「QGIS初心者質問グループ」回答者の方が紹介されていたのですが、「QGIS逆引きガイド」という本がKindle版で発売されているようです。

全2巻で、目次を見るとかなりいろいろな機能の紹介がされているようです。

 2巻買っても900円ということですが、私の本は安くても3000〜4000円にはなると思いますので、こんな安く最新版の紹介をされた本が発売されると、私の本がヤバイですね。


 私の本の発売はまだ未定ですが、発売日等が決まりましたらお知らせしますので、ぜひ手にとって見てください。

林業でQGISを使う実際の事例などもたくさん紹介しています。

よろしくお願いします。



ラベル:QGIS
posted by kouichi at 23:07| Comment(0) | QGIS | このブログの読者になる | 更新情報をチェックする