おわり
2014年08月30日
QGIS入門書
最近出たこの本↓，いいと思います。
おわり
おわり
2014年08月10日
QGISで地理院オルソを表示
QGISで国土地理院のオルソ写真画像を表示できれば便利なのになあと思っていたら，
赤木 実さんがプラグインを開発しておられました。
プラグインの説明サイトはここです↓。
http://space.geocities.jp/bischofia_vb/qgis-plugins/TileLayerPlugin/
説明の通りにやってみたらオルソ画像が表示できたので，その手順をメモしておきます。
まずQGISを開き，プラグイン/プラグインの管理とインストールをクリック↓
アルファベット順に並んでいるプラグインから，TileLayer Pluginを選んでインストールをクリック↓
これでプラグイン自体のインストールは終了です。
前述の赤木さんのプラグイン説明サイトから，以下の部分を，テキストエディタ（例えばメモ帳）にコピペする↓。
スペースを，キーボードの[Delete]キーで削除し，代わりに[Tab]キーでタブを挿入する↓。
編集後のテキストファイルを，文字コード＝UTF-8，拡張子＝tsvで，適当なフォルダに保存する。
私の場合は，D:\TileLayer_Setting\GSI_Tile.tsv↓
日本語のフォルダ名やファイル名は使わないことをお勧めします。
QGISで，プラグイン/タイルレイヤプラグイン/設定をクリック↓
開く窓の外部レイヤ情報ディレクトリ右側のアイコンをクリック↓
先ほどのフォルダ（ディレクトリ）を選び，[フォルダーの選択]をクリック↓
[OK]をクリック↓
これで設定は完了ですが，いきなり地理院地図やオルソを選んでも，どこを表示すればいいのか，QGISが判断できないので，まず最初に，別の地図を開くことにします。
ここでは，以前このブログで取り上げた（http://gpsrsgis.seesaa.net/article/386679119.html），立山町の森林計画図No.1を開くことにします。
地図画像を選ぶ↓
狭い範囲の地図なら何でもいいです。
地図画像オリジナルの，Tokyo測地系第VII系（EPSG＝30167）を選んでOK↓
地図画像が開いたところ↓
プロジェクト/プロジェクトのプロパティをクリック↓
オンザフライCRS変換を有効にし，地理院タイル画像のCRSである，EPSG＝3857を選んでOK↓
QGIS画面右下のEPSGが3857になった↓
プラグイン/タイルレイヤプラグイン/タイルレイヤを追加するをクリック↓
とりあえず，標準地図2012を反転表示して，追加をクリック↓
画面左下の青いプログレスバー？が左右に動く（サーバにアクセス中）↓
地理院の地図画像が表示される↓
再び，タイルレイヤを追加するから，今度は電子国土オルソを選んで追加をクリック↓
地図画面右上に何やら警告が出る↓
オルソ画像のレベルは15〜17なのに，現在のズームレベルが14なので，表示できないという意味だと思われる。
＋アイコンをクリックして倍率を上げる（ズームレベルの数字を大きくする）とオルソが見える↓
もう少し倍率を上げたところ↓
はい。めでたしめでたし。
欲を言えば，オリジナルのオルソ画像は地上分解能（セルサイズ）が約20cmなのに，地理院サーバで提供しているオルソのセルサイズは1〜2m？程度（最高レベルが17）なのが残念なところです。
竹林をデジタイズしてポリゴンデータを作るには，もう少し解像度が欲しいところです。
おしまい。
赤木 実さんがプラグインを開発しておられました。
プラグインの説明サイトはここです↓。
http://space.geocities.jp/bischofia_vb/qgis-plugins/TileLayerPlugin/
説明の通りにやってみたらオルソ画像が表示できたので，その手順をメモしておきます。
まずQGISを開き，プラグイン/プラグインの管理とインストールをクリック↓
アルファベット順に並んでいるプラグインから，TileLayer Pluginを選んでインストールをクリック↓
これでプラグイン自体のインストールは終了です。
前述の赤木さんのプラグイン説明サイトから，以下の部分を，テキストエディタ（例えばメモ帳）にコピペする↓。
スペースを，キーボードの[Delete]キーで削除し，代わりに[Tab]キーでタブを挿入する↓。
編集後のテキストファイルを，文字コード＝UTF-8，拡張子＝tsvで，適当なフォルダに保存する。
私の場合は，D:\TileLayer_Setting\GSI_Tile.tsv↓
日本語のフォルダ名やファイル名は使わないことをお勧めします。
QGISで，プラグイン/タイルレイヤプラグイン/設定をクリック↓
開く窓の外部レイヤ情報ディレクトリ右側のアイコンをクリック↓
先ほどのフォルダ（ディレクトリ）を選び，[フォルダーの選択]をクリック↓
[OK]をクリック↓
これで設定は完了ですが，いきなり地理院地図やオルソを選んでも，どこを表示すればいいのか，QGISが判断できないので，まず最初に，別の地図を開くことにします。
ここでは，以前このブログで取り上げた（http://gpsrsgis.seesaa.net/article/386679119.html），立山町の森林計画図No.1を開くことにします。
地図画像を選ぶ↓
狭い範囲の地図なら何でもいいです。
地図画像オリジナルの，Tokyo測地系第VII系（EPSG＝30167）を選んでOK↓
地図画像が開いたところ↓
プロジェクト/プロジェクトのプロパティをクリック↓
オンザフライCRS変換を有効にし，地理院タイル画像のCRSである，EPSG＝3857を選んでOK↓
QGIS画面右下のEPSGが3857になった↓
プラグイン/タイルレイヤプラグイン/タイルレイヤを追加するをクリック↓
とりあえず，標準地図2012を反転表示して，追加をクリック↓
画面左下の青いプログレスバー？が左右に動く（サーバにアクセス中）↓
地理院の地図画像が表示される↓
再び，タイルレイヤを追加するから，今度は電子国土オルソを選んで追加をクリック↓
地図画面右上に何やら警告が出る↓
オルソ画像のレベルは15〜17なのに，現在のズームレベルが14なので，表示できないという意味だと思われる。
＋アイコンをクリックして倍率を上げる（ズームレベルの数字を大きくする）とオルソが見える↓
もう少し倍率を上げたところ↓
はい。めでたしめでたし。
欲を言えば，オリジナルのオルソ画像は地上分解能（セルサイズ）が約20cmなのに，地理院サーバで提供しているオルソのセルサイズは1〜2m？程度（最高レベルが17）なのが残念なところです。
竹林をデジタイズしてポリゴンデータを作るには，もう少し解像度が欲しいところです。
おしまい。
2014年02月08日
森林計画図のbmp出力
QGIS（v2.0.1）でマージして地理座標系に変換した立山町の計画図を，カシミール３Ｄが読める，ビットマップ（bitmap）形式に変換（画像フォーマットの変換）する手順を書いておきます。
QGISを起動して，ラスタ/変換/変換（形式変換）をクリック↓
開く変換（形式変換）ウィンドウで，入力レイヤ右側の[選択]をクリック↓
前の記事で再投影した，"ALL_TateKei_LL".tifを選んで[開く]をクリック↓
変換（形式変換）ウィンドウで，出力ファイル右側の[選択]をクリック↓
開くウィンドウ下のファイルの種類で，「MS Windows･･（*.bmp，*.BMP）」を選ぶ↓。
探しにくいですが，赤丸を付けたスライダを参考に発見してください。
入力ファイルと同じ，"ALL_TateKei_LL"と命名する（拡張子はbmpが付く）↓
出力ファイル指定後の変換ウィンドウ↓
下のボックスに，"･･-of GTiff･･"と表示されていることに注目！
QGISはうましかなので，ビットマップを指定しても，ジオティフで出力するようにコマンドを誤訳します。
このバグに気付くまで，ずいぶんと時間を無駄にしました。
タダより高いものはない。
参考までに，画像フォーマットとコマンドの解説サイト
GDAL Raster Formatsはここ → http://www.gdal.org/formats_list.html
gdal_translateはここ → http://www.gdal.org/gdal_translate.html
変換ウィンドウの出力サイズ左側にチェックを付け，32%と入力する↓
32%の根拠は後で説明します。
変換ウィンドウ下のボックス右側の，ペン型アイコン（編集）をクリックし，編集モードにする↓。
"GTiff"を”BMP”に書き換えて，[OK]をクリック↓
こんな画面のまま↓，処理が進み，ディスクが回る音がする。
しばらく待つとbmp出力が終わり，こんな画面になる↓。
[OK]と[閉じる]を押す。
QGISは終了してもよい。
形式変換後のエクスプローラの画面はこんな感じ↓
縦1/3，横1/3に間引いたので，ファイルサイズは約1/10になった。
32%の根拠を説明しておく。
前提）
・森林計画図（基本図も）は縮尺1/5,000で作成されている。
・1インチは2.54cmである。
・計画図を400dpi（dot/inch）でスキャンすれば，
1ドット（ピクセル）あたりの地上距離（ピクセルサイズ）は，
2.54(cm) x 5,000 / 400 = 31.75(cm) = 0.3175(m) となる。
地理座標系に変換する前の，平面直角座標系のときの"ALL_TateKei.tif"のプロパティを見ると，pixel sizeは約0.317x0.317mになっていることがわかる。
私の経験上，カシミール３Ｄで計画図（基本図）を表示する際のピクセルサイズ（セルサイズ）は1m程度がちょうど良いと考えており，約0.32mのセルサイズを1mに落とす必要があります。
そこで，sを縮小率（％）とすると，
0.32(m) x 100/s = 1(m) という式が立ち，
32/s=1となり，
s=32（％）と答えが出るわけです。
今回の事例では，400dpiでスキャンした計画図画像を使いましたが，
300dpi（ピクセルサイズは0.423m）でも
200dpi（ピクセルサイズは0.635ｍ）でも
最終的にカシミール３Ｄ用の，ピクセルサイズ1.0m程度の画像を作るには十分な解像度です。
300dpiでスキャンした画像は，bmp化のときに縮小率を42％に，200dpiでは64％に指定すれば，ピクセルサイズ約1mのビットマップ画像ができる。
計画図のbmp化おわり
_
QGISを起動して，ラスタ/変換/変換（形式変換）をクリック↓
開く変換（形式変換）ウィンドウで，入力レイヤ右側の[選択]をクリック↓
前の記事で再投影した，"ALL_TateKei_LL".tifを選んで[開く]をクリック↓
変換（形式変換）ウィンドウで，出力ファイル右側の[選択]をクリック↓
開くウィンドウ下のファイルの種類で，「MS Windows･･（*.bmp，*.BMP）」を選ぶ↓。
探しにくいですが，赤丸を付けたスライダを参考に発見してください。
入力ファイルと同じ，"ALL_TateKei_LL"と命名する（拡張子はbmpが付く）↓
出力ファイル指定後の変換ウィンドウ↓
下のボックスに，"･･-of GTiff･･"と表示されていることに注目！
QGISはうましかなので，ビットマップを指定しても，ジオティフで出力するようにコマンドを誤訳します。
このバグに気付くまで，ずいぶんと時間を無駄にしました。
タダより高いものはない。
参考までに，画像フォーマットとコマンドの解説サイト
GDAL Raster Formatsはここ → http://www.gdal.org/formats_list.html
gdal_translateはここ → http://www.gdal.org/gdal_translate.html
変換ウィンドウの出力サイズ左側にチェックを付け，32%と入力する↓
32%の根拠は後で説明します。
変換ウィンドウ下のボックス右側の，ペン型アイコン（編集）をクリックし，編集モードにする↓。
"GTiff"を”BMP”に書き換えて，[OK]をクリック↓
こんな画面のまま↓，処理が進み，ディスクが回る音がする。
しばらく待つとbmp出力が終わり，こんな画面になる↓。
[OK]と[閉じる]を押す。
QGISは終了してもよい。
形式変換後のエクスプローラの画面はこんな感じ↓
縦1/3，横1/3に間引いたので，ファイルサイズは約1/10になった。
32%の根拠を説明しておく。
前提）
・森林計画図（基本図も）は縮尺1/5,000で作成されている。
・1インチは2.54cmである。
・計画図を400dpi（dot/inch）でスキャンすれば，
1ドット（ピクセル）あたりの地上距離（ピクセルサイズ）は，
2.54(cm) x 5,000 / 400 = 31.75(cm) = 0.3175(m) となる。
地理座標系に変換する前の，平面直角座標系のときの"ALL_TateKei.tif"のプロパティを見ると，pixel sizeは約0.317x0.317mになっていることがわかる。
私の経験上，カシミール３Ｄで計画図（基本図）を表示する際のピクセルサイズ（セルサイズ）は1m程度がちょうど良いと考えており，約0.32mのセルサイズを1mに落とす必要があります。
そこで，sを縮小率（％）とすると，
0.32(m) x 100/s = 1(m) という式が立ち，
32/s=1となり，
s=32（％）と答えが出るわけです。
今回の事例では，400dpiでスキャンした計画図画像を使いましたが，
300dpi（ピクセルサイズは0.423m）でも
200dpi（ピクセルサイズは0.635ｍ）でも
最終的にカシミール３Ｄ用の，ピクセルサイズ1.0m程度の画像を作るには十分な解像度です。
300dpiでスキャンした画像は，bmp化のときに縮小率を42％に，200dpiでは64％に指定すれば，ピクセルサイズ約1mのビットマップ画像ができる。
計画図のbmp化おわり
_
