地理院地図のMBTilesの作り方メモ

注意！この記事はかなり古い情報です。現在のQGISとは操作方法が違いますので注意してください

タイル地図を1つのファイルで管理でき、QGISやGeopaparazziで利用できる、

「MBTiles」の作り方をメモ。

幾つかの手順があり、忘れそうなので。

尚この手順は、Windowsのみで有効です。ご了承ください。

１．QGISをインストールする

---

GISをインストールしておきます。
QGIS本体ももちろん使いますが、ここで重要なのは、一緒にインストールされる、
「OSGeo4w Sell」

２．地理院地図のタイルをダウンロードする

---

地理院地図のタイル画像をダウンロードします。

QGISで「タイルレイヤプラグイン」と「QMetaTilesプラグイン」を使ってもダウンロードできますが、

ズームレベルが場所によって合わないのか、モザイク状になってしまう箇所があるので、

ここでは、「GeopapaofflineTile」というソフトを使います。

何を隠そう、私が作ったソフトです。こちらからダウンロードできます。

QGISでタイルレイヤプラグインをつかって、XYZFreamを表示します。

必要なズームレベルのXY座標をメモして、

GeopapaOfflineTileのリストに入力し、タイルをダウンロードします。

GeopapaOffilenTileの詳しい使い方は、同梱の説明書を参照してください。

タイル地図は、Jpegでダウンロードすることをオススメします。

ファイルサイズが小さくなるので。

３．mbutilを使えるようにする

---

タイル地図をMBTilesに変換するためのソフトが、「mbutil」です。

コマンドラインで動くソフトですので、少し使い方が難しいです。

まずは、Githubからダウンロードします。

こちらのGithubページを開いて、「Clone or Download」ボタンをクリックし、

zipファイルをダウンロードします。

適当なフォルダに保存し解凍しますが、できるだけ簡単なフォルダのほうが後で楽です。

ここでは、「D:\mbutil」に解凍したとことにします。

次に、mbutilをpythonにインストールします。

この手順は、初めの１度だけでOKです。

QGISをインストールした段階で、Pythonは使えるようになっているはずです。

1. スタートメニューの「QGIS Essen」の中にある、「OSGeo4w Sell」を開きます。
2. 「d：」と入力します。
3. 「cd　d:\mbutil」と入力し、対象フォルダを「d:\mbutil」に変更します。
4. 「python setup.py install」と入力します。
5. メッセージが表示されて、mbutilがインストールされました。
6. 正常にインストールされている場合、「python mb-util -h」と入力すると、
   ヘルプが表示されます。

＝＝＝＝2017・8・15追記＝＝＝＝

※Setup.pyが更新されて、QGIS2.1xのOSGeo4w Sellではエラーが出る場合があります。

↑エラー画面

この場合、こちらの「Setup_old.zip」をダウンロードして、

解凍して「mbutil」フォルダにおいてください。

そして、上記４コマンドを「python setup_old.py install」としてください。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

４．mbutilでタイル地図をMBTilesに変換する

---

GeopapaOfflineTileでダウンロードした地理院地図のタイル画像を、

日本語の含まないフォルダに保存します。

ここでは、「d:\tile」に保存したとします。

「d:\tile」の直下にズームレベルごとのフォルダが有る状態です。

次の手順で、mbtilesに変換します。

dドライブの直下に、「chiriin.mbtiles」というファイルを作成します。

1. スタートメニューの「QGIS Essen」の中にある、「OSGeo4w Sell」を開きます。
2. 「d：」と入力します。
3. 「cd　d:\mbutil」と入力し、対象フォルダを「d:\mbutil」に変更します。
4. 「python mb-util --image_format=jpg d:\tile d:\chiriin.mbtiles」と入力します。
5. 変換が始まります。（コマンドがたくさん流れていく）
6. 変換が終了すると、MBTilesファイルが作成されています。

※今回のタイル画像がjpegなので、オプションに「--image_format=jpg」と入力しています。

作成されたMBTilesファイルをQGISにドラッグして、内容を確認します。

基盤地図情報から傾斜区分地図を作る（その２）

注意！この記事はかなり古い情報です。現在のQGISとは操作方法が違いますので注意してください

前の記事：基盤地図情報から傾斜区分地図を作る（その１）

＜QGISに標高データを表示する＞

QGISを入手していない場合は、右の「GIS関係リンク集」の「QGIS　WIKI」からダウンロードして下さい。

QGISは、Ver1.7.2で説明します。

QGISを起動し、「ラスタレイヤの追加」で前の手順で作成した「merge.tif」を選択します。

エクスプローラから「レイヤ」に直接ドロップしてもOKです。

　　↓

　　↓

開くと、グレーの四角が画面に表示されます。

　　↓

＜プロジェクトの測地系を設定する＞

表示された標高データは、変換時に指定した座標系になっています。

しかし、プロジェクト（地図全体）の測地系は設定されていないので、

標準では「WGS８４」になっています。

標高データの測地系と、プロジェクトの測地系を合わせます。

レイヤ名を右クリックし、「レイヤのCRSをプロジェクトに設定する」を選択します。

　　↓

これで、プロジェクトの測地系が標高データとおなじになりました。

一応「設定」→「プロジェクトのプロパティ」で確認しておきます。

　　↓

他の測地系のデータを同時に表示する場合には、

「オンザフライCRS変換を有効にする」にチェックを付けておきます。

データの測地系が正しく設定されていれば、変換されて同じ地図上に表示されます。

　　↓

＜傾斜区分図を作成する＞

「ラスタベース地形解析」プラグインを使用して、傾斜区分図を作成します。

「プラグイン」→「ラスタベース地形解析」を選択して下さい。

　　↓

このプラグインが表示されない場合は、

「プラグイン」→「プラグインの管理」で「ラスタ空間解析プラグイン」にチェックを付けてください。

　　↓

ダイアログが表示されたら、

解析手法を「傾斜」にして、入力レイヤに変換するレイヤ名を選択します。

出力レイヤに保存するファイル名を入力するのですが、

フォルダ名に日本語があると変換できません。

また、ファイル名に日本語を使うと文字化けします。

なので、出力するフォルダ名、ファイル名は、すべて半角英数文字で設定します。

「結果をプロジェクトに追加する」にチェックを付けておくと、

ファイル作成後、地図に追加されます。

　　↓

新しいレイヤが追加されたら、「merge」レイヤを非表示にします。

　　↓

新しく追加したレイヤ名をダブルクリックし、プロパティを表示します。

　　↓

「スタイル」タブの「カラーマップ」を「カラーマップ」にします。

　　↓

「カラーマップ」タブを選択し、「エントリを追加」をクリックします。

ここで、角度ごとに色を設定していきます。

　　↓

数字部分をクリックして、０から１０度刻みで数字を入力していきます。

　　↓

０は白、10は青、20は緑、30は黄色、40はピンク、50は赤、それ以上は黒

に設定します。

（角度や色は自由に変更しても構いません）

　　↓

「OK」ボタンをクリックすると、角度ごとの色分けがされています。

　　↓

起伏がわかりやすいように、陰影起伏図と重ねてみます。

前の手順で作成した「merge_sh_img.tif」を「ラスタレイヤの追加」で追加します。

　　↓

レイヤの順番を入れ替えます。

　　↓

傾斜区分図のレイヤをダブルクリックし、プロパティを表示します。

「透過性」タブの「全体の透過率」を５０％くらいにします。

　　↓

これで、陰影起伏図と重なって、地形がわかりやすいくなりました。

　　↓

この他、等高線や林小班などを重ねると、もっとわかりやすい地図になります。

今回作った傾斜区分図は、

０〜１０度までは「青」、

１０〜２０度までは「緑」

２０〜３０度までは「黄色」

３０〜４０度までは「ピンク」

４０〜５０度までは「赤」

５０〜９０度までは「黒」

に色分けしました。

林業専用道や森林作業道を計画する場合、２５度以下くらいの横断傾斜の箇所を選定するので、

今回の地図だと「青」、「緑」、「黄色」の部分を通る計画をすれば、安い経費で計画できるということです。

現地の踏査を行う前に、この地図で線形の候補を設定しておけば、

現地踏査も少しは楽になるかもれません。

QGISでは、レイヤをGPXファイルに変換できます。

GPSにQGISで設定したラインをいれて、現地に行けば更に踏査が捗るでしょう。

基盤地図情報から傾斜区分地図を作る（その１）

注意！この記事はかなり古い情報です。現在とは操作方法が違いますので注意してください

林業ニュースさんに紹介されたので、すこしヤル気を出して、

数回に分けて、「基盤地図情報から傾斜区分地図を作る」を説明していきたいと思います。

＜国土地理院から標高データDEMを入手する＞

国土地理院の「基盤地図情報ダウンロードサービス」から標高データを入手します。

以前は自由にダウンロードできたのですが、現在はユーザー登録を行う必要があります。

ユーザー登録を行い、ログインすると次の画面になります。

「JPGIS(GML)形式を選択します。

　　↓

地図の種類と、場所を選択します。

地図の種類は「10ｍメッシュ標高」を選択します。

場所は、必要なところを選択して下さい。複数選択することも出来ます。

選択したら、「次へ」をクリックします。

　　↓

そのまま、選択した都道府県全体をダウロードすることも出来ますが、

必要な部分のみ地図で選択することも出来ます。

北海道など広い都道府県は、一括のダウンロードはできないようです。

　　↓

「地図から選択」ボタンをクリックすると、

地図からダウンロードする場所を選択できます。

はじめに地図の上にあるコンボボックスで都道府県を選び、

ある程度地図を拡大すると、メッシュが表示されます。

「メッシュ選択」をクリックすると、

地図からダウンロードしたい場所を選択できます。

　　↓

「ダウンロード」をクリックすると、次の画面になります。

右の「ダウンロードボタンをクリックして、

データをダウンロードしてください。

　　↓

ダウンロードしたら、Zipファイルを解凍し、

１つのフォルダにファイルをまとめます。

＜標高DEMデータ変換ツールを入手し、データの変換を行う＞

基盤地図情報標高DEM変換ツールを、株式会社エコリスのページからダウンロードします。

便利なツールを開発して頂き、感謝です。

「ダウンロードはこちら」をクリックし、「demtool.zip」をダウンロードしたら、

解凍し、「変換結合.vbs」を実行します。

投影法、座標の系番号を入力します。

　　↓

　　↓

段彩図と陰影起伏図を作成しますか？と聞いてくるので、「はい」を選択します。

スケールは25000のままでいいです。

　　↓

　　↓

海域の標高は0なので「はい」を選択します。

　　↓

黒い画面が表示され、変換が始まります。

　　↓

ファイルの数にもよりますが、３０分から１時間くらいかかるときもありますので、

気長に待ってください。

終了するとメッセージが表示されます。

作成されたデータで傾斜区分図に必要なものは、

「merge.tif」（標高図）と「merge_sh_img.tif」（陰影起伏図）だけです。

それ以外のファイルは削除しても構いません。

（その２）につづく・・・・

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝