森林土木memoある自治体職員が、仕事で使うGIS、GPS、CADなどについてメモします。知床半島の立体図(CS立体図)を作成してみた2016年11月19日放送の「ブラタモリ」で、知床の地形について説明していたので、知床の立体図を作成してみた。知床の羅臼岳周辺の立体図また、2016年8月24日に大雨の影響で、羅臼において土砂崩れが発生しています。この箇所の地形を立体図でアップにしてみました。10mメッシュですので少し荒いですが、地形がよくわかります。土砂崩れした箇所は、マウスカーソルの十字がある場所です。corvature」コマンドで作成できる様々な地図を傾斜図に重ねてみました。Curvatureどうでしょうか?様々な立体図で見ると、土砂崩れした斜面の上流部には、沢地形があり、崖の手前で大きく南に曲がっています。水の集まりやすい地形ですね。このような地形は、大規模な土砂崩れが発生する可能性が高いので気をつけなければいけません。北海道の海岸沿いで、隆起した地形はこのような地形が多いかもしれませんね。立体図で見ると、このような地形がよく分かるのでいいですね。Curvature」で見ると、隠れた水道(みずみち)も地図で見ることができるようですね。このような立体図は、10〜30分程度で簡単に作成することができます。災害後の復旧計画にも、立体図を参考にすることができそうです。基盤地図情報の10mメッシュでも十分使えそうですね。|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立体図のうちの立体図の作成方法を紹介します。Curvature」という計算方法を使用していますが、これが本当に正しいのか、著者はわかっていません。他にいい方法があれば情報をください。Curvature」で作成した地図は、等高線上の地形のシワが強調されて表現されます。標高DEMを準備する標高DEMを基盤地図情報などからダウンロードします。株式会社エコリスの変換ツールを使って標高DEM(Tiffファイル)を作成します。作成したファイル(merge.tif)は、日本語を含まないフォルダ、ファイル名にして保存します。QGISのレイヤに標高DEMを追加する標高DEMをQGISに追加します。標高DEMのファイルをマップキャンバスにドラッグ&ドロップするとレイヤに追加することができます。プロセッシングツールで傾斜と曲率を計算プロセッシングツールを使って、標高DEMから傾斜図と曲率図を作成します。メニューに「プロセッシング」が表示されていない場合は、「プラグイン」→「プラグインの管理とインストール」で「Processing」を有効にしてください。メニュー「プロセッシング」→「ツールボックス」を選択して、プロセッシングツールボックスを表示します。Curvature」を選択します。表示されたダイアログの「Elevation」に標高DEMレイヤを選択します。Units」は、デフォルト設定で構いません。Curvature」の2つだけですので、それ以外の「アルゴリズムの実行後に出力ファイルを開く」のチェックを外します。Curvature」が作成されます。レイヤを名前をつけて保存する作成されたレイヤは、一時レイヤに作成されます。そのため、QGISを終了したら削除されてしまいます。レイヤを右クリックして、「名前をつけて保存する」を選択します。「出力モード」を「生データ」にします。「形式」を「GTiff」にします。「パス」の「参照」ボタンをクリックして、ファイルの保存先を設定します。座標参照形を設定する必要がある場合は、「CRS」で選択します。「OK」ボタンをクリックすると、ファイルが保存されて、レイヤに追加されます。Curvature)は削除します。傾斜レイヤのスタイルを設定する傾斜レイヤのプロパティを開いて、「スタイル」を選択します。「レンダータイプ」を「単バンド疑似カラー」にします。「色」を「Reds」にします。自動的にリストに数値が分類されます。「カラーレンダリング」の「混合モード」を「乗算」にします。「OK」ボタンをクリックするとスタイルの設定ができます。曲率レイヤのスタイルを設定する曲率レイヤのプロパティを開いて、「スタイル」を選択します。「レンダータイプ」を「単バンド疑似カラー」にします。「色」を「Blues」にします。自動的にリストに数値が分類されます。今回は、マイナス側を青にしたいので、「反転」にチェックを付けます。「カラーレンダリング」の「混合モード」を「乗算」にします。「OK」ボタンをクリックするとスタイルの設定ができます。必要であれば、地理院地図などを背景に追加する必要であれば、地理院地図や等高線などを地図に追加すると、地図がわかりやすくなるかもしれません。長野県林業総合センターの資料では、標高DEMも重ねていますが、標高DEMは重ねなくても立体的でわかりやすいと思います。図:地理院地図を重ねた立体図CS立体図は、赤色立体地図よりも地すべりなどの地形が見やすいように感じます。長野県では、1mメッシュDEMを利用していますので、かなりの微地形を確認することができるようですが、基盤地図情報の10mメッシュDEMでも十分地形判読を行うことができます。地質図や地すべり位置図などを重ねることで、更に地形判読を行いやすくなると思います。〜2016年11月24日追記〜長野県林業総合センターのCS立体図の中の人から、ご連絡をいただきました。CS立体図の作成方法は、オープンなので自由に紹介して使ってほしいとのことです。CS立体図の普及によって、林業がもっと発展すれば嬉しいとのことでした。素晴らしい考え方です。長野県では全国の10mメッシュDEMを使って、CS立体図を作成済みとのことです。そのうち、自由にダウンロードできるようになると思いますので、非常に楽しみです。ダウンロードできるようになれば、また紹介したいと思います。本ブログでは、QGISにおける赤色立体地図の作成方法を公開していましたが、この度赤色立体地図の生みの親である千葉達朗氏から、公開をやめてほしい旨の連絡をいただきました。その為、赤色立体地図の作成方法の公開についてはやめさせて頂きます。個人的に研究する分には、問題ないと思いますので、QGISの操作方法について確認したい方には、個別に対応させていただきたいと思います。質問のあるかたは、「qgisshitumonアットgmail.com(アットは@に置き換えてください)」にメールを頂ければ、返答します。ご迷惑をおかけして申し訳ありません。今後も「森林土木memo」を、よろしくお願いします。2.0業務で使う林業QGIS徹底使いこなしガイドブログ主が書きましたQGISのマニュアル本です。2017年6月発売!「林業」とありますが、林業だけでなくQGISの基本から応用まで幅広く網羅しています。御覧ください!Amazonでも販売中!下のリンクから、本で利用できるサンプルファイルをダウンロードできます。業務で使う林業QGIS徹底使いこなしガイド本の紹介ページ(全林協)公開しているテキスト当ブログ管理人が作成したQuantumGISのテキストなどをダウンロード出来ます。北海道Shapeファイル北海道全域の基盤地図情報をShapeに変換しました。)、「道路縁」があります。(04/19)ファンメッセージを送るこのブログの読者になる更新情報をチェックするブックマークする友達に教えるカテゴリ日記(14)その他(2)CAD(1)GIMP(1)林道(1)QGIS(50)GPS(4)リンク(1)DBF編集(2)GISその他(22)魚道(0)治山ダム(1)Android(10)カシミール3D(6)vb.net(2)Web地図(5)過去ログ2019年10月(2)2019年02月(1)2018年09月(3)2018年08月(1)2018年03月(1)2017年12月(1)2017年08月(1)2017年07月(3)2017年06月(2)2017年05月(2)2017年02月(1)2017年01月(1)2016年12月(1)2016年11月(4)2016年09月(1)2016年08月(3)2016年07月(3)2016年03月(1)2015年12月(1)2015年09月(2)Seesaaブログ。