人気ブログランキング | 話題のタグを見る

引越

諸事情により引っ越します。
http://forestry-gis.blogspot.jp/


# by thorb38 | 2017-01-22 19:29 | 独り言

高度なデジタイジングツール解説

(QGIS2.14)

FOSS4G Advent Calendar 2016 6日目の記事です。
http://qiita.com/advent-calendar/2016/foss4g
ここで一気にハードルを下げますw
2つめに投稿した方が良かったと思いつつ、ちょいと先走って登録してしまったので…


さて2年ほど前に、『諸君、私はポリゴンが好きだ』と仰った方がいましたが
(僕もポリゴン大好きです)
http://d.hatena.ne.jp/wata909/20141222
大した技術もない僕にとっては、QGISの機能が底上げされていくということが非常にありがたいわけで。

このところQGISのデジタイジングツールも飛躍的に向上してきましたが、あまり解説しているサイトを見かけないので まとめておこうと思った次第。
(まぁ今さらまとめるまでもない、というのが多くの認識かもしれませんが…)


まず、デジタイズ時に欠かせいないスナップオプションですが、QGISの設定から、
高度なデジタイジングツール解説_a0122776_18473943.png
『ドックウィンドウ内でスナップオプションを開く』にチェックを入れておくと
高度なデジタイジングツール解説_a0122776_18491079.png
スナップオプションのウィンドウを独立させて常に表示しておけるので便利です。


では、改めてデジタイジングツールの説明に入ります。


高度なデジタイジングツール解説_a0122776_18495809.png

高度なデジタイズパネル

高度なデジタイジングツール解説_a0122776_18504480.png
このアイコンをクリックするとパネルが開きます。
高度なデジタイジングツール解説_a0122776_18595863.png
d【距離】、a【角度】、x【x座標】、y【y座標】を直接入力することができます。
また、値を相対値で表示したり、値をロックすることができます。

右上のオプションから
高度なデジタイジングツール解説_a0122776_19012715.png
作図時にスナップする角度や、スナップオプションの適用方法を設定することができます。

正確な作図をする際に役立ちます。
高度なデジタイジングツール解説_a0122776_19034516.png

トレース

高度なデジタイジングツール解説_a0122776_21125920.png
既存のレイヤをトレースして作図することができます。
始点と終点をクリックすれば、
高度なデジタイジングツール解説_a0122776_21084718.png
高度なデジタイジングツール解説_a0122776_21094293.png
ラインやポリゴンの縁をトレースしてくれますが、
高度なデジタイジングツール解説_a0122776_21105147.png
思った方にトレースしてくれない場合は経由点をクリックしましょう。


Undo Redo

高度なデジタイジングツール解説_a0122776_21141846.png
取消と再実行。説明いりませんね。


地物の回転

高度なデジタイジングツール解説_a0122776_21152034.png
回転させたい地物をクリックすると、入力パネルが表示されます。
マウスを動かして回転させても良いし、パネルに直接入力でも可。
高度なデジタイジングツール解説_a0122776_21170169.png
スナップ先を指定することで、指定した角度毎に回転させられます。


地物の簡素化

高度なデジタイジングツール解説_a0122776_21181157.png

簡素化させたい地物をクリックすると、簡素化ツールが開きます。
高度なデジタイジングツール解説_a0122776_21190028.png
許容範囲を設定しながら、
高度なデジタイジングツール解説_a0122776_21194786.png
地物の頂点数をどのくらい減らすか調整しましょう。
高度なデジタイジングツール解説_a0122776_21201796.png
高度なデジタイジングツール解説_a0122776_21205995.png
個人的に使ったことがない機能ですが、海岸線を簡素化して地図をイラスト化したり、ハンディGPSのトラックログを間引きするのに便利かもしれませんね。


リングの追加

高度なデジタイジングツール解説_a0122776_21220410.pngいわゆるドーナツ状のポリゴンを作るための機能です。


高度なデジタイジングツール解説_a0122776_21225313.png
高度なデジタイジングツール解説_a0122776_21225324.png
高度なデジタイジングツール解説_a0122776_21225421.png

部分の追加

高度なデジタイジングツール解説_a0122776_21235456.png元となる地物を選択した後、
高度なデジタイジングツール解説_a0122776_21244460.png
部分を作図して追加します。
高度なデジタイジングツール解説_a0122776_21244569.png
高度なデジタイジングツール解説_a0122776_21244524.png
見かけ上 複数のポリゴンですが、1つの地物として扱うことができます(マルチポリゴン)。


リングの充填

高度なデジタイジングツール解説_a0122776_21260432.pngドーナツの穴を作りつつ、


高度なデジタイジングツール解説_a0122776_21263587.png
高度なデジタイジングツール解説_a0122776_21263568.png
高度なデジタイジングツール解説_a0122776_21263552.png
新しい地物で埋めてくれます(ポリゴンが重なっているわけではない)。


リングの削除

高度なデジタイジングツール解説_a0122776_21283453.png穴を埋めます。


高度なデジタイジングツール解説_a0122776_21225421.png
高度なデジタイジングツール解説_a0122776_21225313.png

部分の削除

高度なデジタイジングツール解説_a0122776_21333329.pngマルチポリゴンの一部を削除できます。


高度なデジタイジングツール解説_a0122776_21343226.png
高度なデジタイジングツール解説_a0122776_21343291.png
地物の数は変わりません。


地物の変形

高度なデジタイジングツール解説_a0122776_21353777.pngラインを引くことで、その形に地物を変形できます。


高度なデジタイジングツール解説_a0122776_21363356.png
高度なデジタイジングツール解説_a0122776_21363487.png
高度なデジタイジングツール解説_a0122776_21364328.png
高度なデジタイジングツール解説_a0122776_21364412.png
隣接したポリゴンを変形した場合、
高度なデジタイジングツール解説_a0122776_21375302.png
高度なデジタイジングツール解説_a0122776_21375468.png
高度なデジタイジングツール解説_a0122776_21375401.png
トポロジ編集してくれます。


曲線のオフセット

高度なデジタイジングツール解説_a0122776_21400975.png個人的に使い道が思い浮かばないので、何が出来るのかは、スクショを参考に。


高度なデジタイジングツール解説_a0122776_21482191.png

地物の分割

高度なデジタイジングツール解説_a0122776_21425366.png
地物をぶった斬ることが出来ます。


高度なデジタイジングツール解説_a0122776_21494517.png
高度なデジタイジングツール解説_a0122776_21495426.png

他にラインレイヤがあったりした場合は、上記のトレース機能を併用することができます。
高度なデジタイジングツール解説_a0122776_21495720.png
高度なデジタイジングツール解説_a0122776_21495884.png


部分の分割

高度なデジタイジングツール解説_a0122776_21525099.png『地物の分割』との違いが分かりにくいですが、『部分の分割』の場合は地物は1つのままです(マルチポリゴン化)。


高度なデジタイジングツール解説_a0122776_21535524.png
高度なデジタイジングツール解説_a0122776_21535541.png
高度なデジタイジングツール解説_a0122776_21544148.png
『地物の分割』を使った場合は、スクショの例では地物が4つになります。
高度なデジタイジングツール解説_a0122776_21551178.png
高度なデジタイジングツール解説_a0122776_21551183.png
個人的には、「部分に分割」なのでは、と思ってます。


地物の結合

高度なデジタイジングツール解説_a0122776_20495045.png結合させたい地物を選択して、


高度なデジタイジングツール解説_a0122776_20513778.png
高度なデジタイジングツール解説_a0122776_20513820.png
引き継ぐ属性データを設定すれば地物が1つになります。
高度なデジタイジングツール解説_a0122776_21244460.png
例では属性フィールドが1つだけですが、フィールド毎に引き継ぐデータを設定することができます。
また、隣接していないポリゴンを選択して結合した場合は、
高度なデジタイジングツール解説_a0122776_20572840.png
高度なデジタイジングツール解説_a0122776_20513869.png
マルチポリゴンになります。


地物の属性結合

高度なデジタイジングツール解説_a0122776_20572756.png地物の属性のみを結合します。フィールド毎に引き継ぐデータを設定できるのは『地物の結合』と同じです。


高度なデジタイジングツール解説_a0122776_20572840.png
今回は、属性データBを適用してみます。
高度なデジタイジングツール解説_a0122776_21004590.png
高度なデジタイジングツール解説_a0122776_20584387.png
『地物の結合』とは違い、地物は2つのままです。
高度なデジタイジングツール解説_a0122776_20584318.png
高度なデジタイジングツール解説_a0122776_20584475.png

ポイントシンボルの回転

高度なデジタイジングツール解説_a0122776_21031858.pngシンボルを回転させるには、少し下準備が必要です。

高度なデジタイジングツール解説_a0122776_21032067.png
まず、ポイントデータの属性テーブルを開いて、新しいフィールドを追加します。
高度なデジタイジングツール解説_a0122776_21032139.png
名前は何でも良いです。タイプは整数値、長さは最低でも3に。
高度なデジタイジングツール解説_a0122776_21032117.png
データは、NULLのままで良いです。

次に、レイヤのプロパティを開いてスタイル設定をします。
高度なデジタイジングツール解説_a0122776_21044477.png
回転の右にあるボタンをクリックして、追加したフィールドで設定します。

これで、地物をクリックすればシンボルが回転できるようになっています。
高度なデジタイジングツール解説_a0122776_21044407.png
属性データに直接入力でも良いです。
高度なデジタイジングツール解説_a0122776_21044473.png

さて、如何だったでしょうか。
徹底解説というほどに丁寧に説明したわけでもありませんが(スミマセン…)、多少なりとも参考になったと思ってくれる方がいれば幸いであります。



以下、ダラダラと雑談

今回、デジタイズツールを取り上げたのは、自分が最も使っているツールだからです。
せっかくのGISなんだから、解析して云々だとか意思決定支援ツールとして、だとか色々とあるんでしょうが、個人的には その前段となるデジタイズが最も楽しい。
そんな僕がOpenStreetMapに手を出すと
http://dohokugeo.blogspot.jp/2015/12/blog-post.html
こんなことになります。それはさておき…


うちの職場、北海道に307万ha(北海道の約4割)の土地を持っています。
この広大な面積の図面も管理しているわけですが、未だに修正は手書きです。
手書きで修正した図面を(コピーしたりスキャンしたり諸々歪ませて)毎年、上局に提出しています。それがさらに上局にある図面に(恐らく)手作業で書き写され、その後 外部にGISデータの作成が発注され、デジタル化されて1年後くらいに帰ってきます(詳しくは知らないけど、おおよその流れはこんな感じのはず…)。
常に、1~2年遅れのデータをGISで使っているというわけです。
さらに言えば、データ作成を請け負った事業者が(多分)林業に精通しているわけもなく、業界内の人間から見たら「それはないだろ」と言いたくなるような間違いがあったりします。
あと、ドーナツ状であるべきポリゴンも全て紐付だったり……
高度なデジタイジングツール解説_a0122776_21110795.png
こいつが原因で、解析やら何やらでエラー吐きまくりです。まぁ簡単に修正できるので良いんですが、修正方法を知るまでは泣きそうでした。
(いまだに紐付ポリゴンとか隙間のあるポリゴンを見ると殺意が湧きますが)

結局、1年も待っていられないし、待っていても間違ったデータが来るのが分かっているので、僕は常に自分で修正・最新化しているわけです。その上で、QGISを使える他の職員と共有したり、誤りを上局に報告しています。
ある意味で、図面修正しているだけで幸せな僕には良い状況なのかも……いや良くないだろ。


まぁそんなこんなで、
高度なデジタイズツールが実装された時の僕の喜びようを分かっていただけるのでは。


Let’s digitizing!


# by thorb38 | 2016-12-06 21:15 | GPS・GIS

Freehand raster georeferencer

(QGIS2.14)
さて、前回ラスタをドラッグで移動させる方法を紹介しましたが、
http://denthor.exblog.jp/22581100/
あくまでも位置情報のついたラスタを移動させるものでした。

今回のFreehand raster georeferencerプラグインは、位置情報のないラスタをQGISに表示し、マウスドラッグで操作して位置情報を与えることができます。

操作は非常に簡単です。
【AD】をクリックし、取り込むラスタデータを選択します(jpg、bmp、png)。
Freehand raster georeferencer_a0122776_18213487.png




【MO】をクリックしてから、マウスドラッグでラスタを移動できます。
Freehand raster georeferencer_a0122776_18262144.png
Freehand raster georeferencer_a0122776_18275318.jpg

【RO】で回転。
Freehand raster georeferencer_a0122776_18305917.png
Freehand raster georeferencer_a0122776_18315656.jpg


【SC】【AD】で変形。
Freehand raster georeferencer_a0122776_18325973.png
Freehand raster georeferencer_a0122776_1833232.png
Freehand raster georeferencer_a0122776_1835555.jpg

【T-】【T+】透明度の変更。
Freehand raster georeferencer_a0122776_18371590.png
Freehand raster georeferencer_a0122776_18371899.png






【II】でエクスポート(jpg、bmp、png)。ワールドファイルも作成されます。
Freehand raster georeferencer_a0122776_18394943.png

Freehand raster georeferencer_a0122776_18431845.gif

# by thorb38 | 2016-07-14 18:48 | GPS・GIS

ラスタデータをちょっとズラす

(QGIS2.8)

例えば、紙図面をデジタル化してGISに表示させたい、というような時。
スキャナでPCに取り込んで、ジオリファレンサで位置情報を与えるわけですが、ピッタリ合うことなんて ほとんどないわけです。

こんな感じ(わざと大きめにズラしてます)
ラスタデータをちょっとズラす_a0122776_20432264.png

さてそんなとき、もう1度ジオリファレンスし直すか否か。
誤差の大きさにもよるのでしょうが、一時的にしか使わないようなデータだったりすると、わざわざやり直すのも面倒な訳で…

「マウスでドラッグして、ちょいとズラしたいなぁ」と思ったこと ありませんか?

というわけで、プラグイン【rasmover】の紹介です。

ラスタレイヤを選択した後に、
ラスタデータをちょっとズラす_a0122776_20451355.pngrasmoverのアイコンをクリック。




その後に、画像をドラッグしてズラしましょう。
ラスタデータをちょっとズラす_a0122776_20465923.png
画面上では、赤い始点と移動先を示す線が表示され、ドラッグ後に画像が移動します。
※ラスタデータのファイル名などに日本語が混じっているとエラーが出ます。

いかがでしょうか
ラスタデータをちょっとズラす_a0122776_2049086.png
だいぶマシになりましたね。
移動させたラスタは【moved】という名前で一時レイヤとして追加されます。
ラスタデータをちょっとズラす_a0122776_20491254.png
【名前をつけて保存】しておかないと、QGISを終了させた時に消えるので注意してください。


ジオリファレンサについては、↓ こちらが大変参考になります。
http://www.geopacific.org/opensourcegis/gcngisbook/QGIS_book/7b2c17ae0/7b2c27ae07b2c17bc0

# by thorb38 | 2016-03-09 20:54 | GPS・GIS

Cartogramプラグインを使ってみた

統計地図(カルトグラム)とは…
主題地図の一種であり、人口や移動時間、国民総生産などを土地面積と距離に置き換えた地図である。面積統計地図と距離統計地図の2種類あり、それぞれ数値を地図化することにより地形が歪んでいる。
Wikipediaより引用

ということで、QGISにCartogramのプラグインがあったので試してみました。
(QGIS2.8)

テーマは、
北海道はどれくらい札幌なのか?

「北海道といえば札幌」とは道北在住者として言いたくはないのですが、人口は200万人を突破しようか、という勢いであり、野球もサッカーもバスケットボールも、プロチームは全て札幌が本拠地です。
『北海道日本ハムファイターズ』と名乗ったところで、札幌以外の人が日ハムに愛着が持てるかと言ったら、かなり微妙な気がします。
(地方から熱心に応援に行っている人たちがいることも承知していますが…)

まぁ、それはともかく
では北海道が実際のところ、どれくらい札幌なのか視覚化してみよう、というのが今回の趣旨です。

まず、データを集めます。
・北海道の地図データ(Shapefile)
  国土数値情報の行政区域
  http://nlftp.mlit.go.jp/ksj/
・市町村別人口データ(csv)
  e-Stat 平成22年度国勢調査
  http://www.e-stat.go.jp/SG1/estat/eStatTopPortal.do
・観光客数
  北海道観光入込客数調査報告書
  http://www.pref.hokkaido.lg.jp/kz/kkd/irikomi.htm


QGISに、北海道の地図を表示します。
(統計データのない北方領土を削除したり、小さい島を削除したり、飛び地になっている市町村をマルチポリゴン化したり、多少編集しています)
Cartogramプラグインを使ってみた_a0122776_19302081.png


統計データを、属性データに追加します。
(csvファイルを若干、編集してます)
Cartogramプラグインを使ってみた_a0122776_19331299.jpg

属性データにcsvファイルを追加する方法の詳細は、こちら↓
https://sites.google.com/site/qgismemo/040-bian-ji/shu-xingdetanicsvmatahadbffairuwo-jie-heshitai


あとは、Cartogramプラグインを起動して、
Cartogramプラグインを使ってみた_a0122776_19344888.jpg
レイヤ(北海道の地図)と、変換元にする属性データ(人口データ等)を選択するだけです。
※Cartogramプラグインは、変換元にマイナスのデータや【0】【null】が含まれていると動きません。


人口(2010年)によるカルトグラム↓
Cartogramプラグインを使ってみた_a0122776_19355818.png

観光客数↓
Cartogramプラグインを使ってみた_a0122776_19363611.png
後志地方が盛り返している印象。


北海道でも、オープンデータの公開が徐々に広まってきており
室蘭市 http://www.city.muroran.lg.jp/main/org2260/odlib.php
森町 http://www.town.hokkaido-mori.lg.jp/docs/2014101000041/
札幌市 http://www4.city.sapporo.jp/opendata/
北海道庁 http://www.pref.hokkaido.lg.jp/ss/jsk/opendata/opendata.htm

様々なデータが扱いやすい状態で公開されるケースが増えてきています。
これらのデータを使って面白い地図が作れるかもしれませんね。



余談:
このネタは本来、もうちょい色々とデータを使ったりして煮詰めた上でFOSS4G Advent Calendar 2015に投稿するつもりでしたが、色々あって間に合いませんでした。
今年の年末まで温めておいても仕方がないので、公開しておきます。

今から年末のネタ探しが心配ですw
# by thorb38 | 2016-01-23 19:49 | GPS・GIS