2019年05月07日

QGISで画像解析

パロットBLUEGRASSに搭載されたSequoiaセンサで撮影した4バンド(G,R,RE,NIR)画像を,
Pix4Dmapperで処理すると,各バンドの反射率画像が作成されます。

Explorerの画面です↓。
0反射画像ALL.gif
波長の短いものから列挙すると,
Green(可視域緑,G)
Red(可視域赤,R )
Red_Edge{赤と近赤外域の間(勝手に「赤端」と命名する),RE}
NIR(近赤外域,NearInfraRed,NIR)

QGISで画像とヒストグラムを表示しました。
まずはGreen↓
1緑_画像.gif

2緑_noalpha_ヒスト.gif

次はRed↓
3赤_画像.gif

4赤_noalpha_ヒスト.gif

次はRE↓
5赤端_画像.gif

6赤端_noalpha_ヒスト.gif

最後にNIR↓
7近赤_画像.gif

8近赤_noalpha_ヒスト.gif


RとNIRを使ってNDVI(正規化植生指数)を計算してみます。
計算式は,
NDVI=(NIR−R)/(NIR+R)
です。
9NDVI計算1.gif

ラスタ/ラスタ計算機を開き,ラスタバンドや演算子をクリックして式を作ります↓
10NDVI計算2.gif
出力レイヤはNDVIです。

出力されたNDVIの画像です↓。
11NDVI計算3.gif
NDVI値の低い塊(黒っぽく見える)が3箇所確認できます。

NDVIのヒストグラムです↓。
12NDVI計算4.gif


さて,Pix4DmapperでSequoia画像を処理すると,自動でNDVIを計算してくれます。
その画像を呼び出して,ヒストグラムを表示したところです↓。
13NDVI計算5.gif

ふたつのヒストグラムはよく似ていますね。

というか,両者は完全に一致しています。
14NDVI計算6.gif
「完全に一致した。」は某国(亡国?)の首相の口癖ですね。

一致して当たり前で,違っていたら困ります。

このラスタ計算機を使えば,NDVI以外の植生指数,例えば,
---------------------------------------------------------------------
NDRE=(RE−R)/(NIR+RE)クロロフィル濃度
GNDVI=(NIR−G)/(NIR+G)光合成量,植物ストレス
---------------------------------------------------------------------
も計算できますね。

いずれ共同研究者である職場の同僚にやってもらいましょう。

最後に,SequoiaのRGBカメラ画像とヒストグラムも載せておきます↓。
15RGB_画像.gif

16RGB_ヒスト.gif

赤く見えるのは,NDVI値が小さかった,海岸クロマツ林のマツ枯れです。




posted by Dr.koba at 11:04| Comment(0) | QGIS(無料GISソフト) | このブログの読者になる | 更新情報をチェックする

2019年05月06日

Learn QGIS 第6章

頑張ったけど,途中で力尽きました。

Python Scriptうまく動きません。
タイプミスもあるけど,ほかにも不具合があるんじゃないかな?

パイソンアクション13枚
1パイソンaction.gif

2パイソンaction2.gif

3パイソンaction3.gif

4パイソンaction4.gif

5パイソンaction5.gif

6パイソンaction6.gif

7パイソンaction7.gif

8パイソンaction8.gif

9パイソンaction9.gif

10パイソンaction10.gif

11パイソンaction11.gif

12パイソンaction12.gif

8パイソンaction8.gif

パイソンベクタ14枚
14パイソンconsole_V.gif

15パイソンconsole_V2.gif

16パイソンconsole_V3.gif

17パイソンconsole_V4.gif

18パイソンconsole_V5.gif

19パイソンconsole_V6.gif

20パイソンconsole_V7.gif

21パイソンconsole_V8.gif

22パイソンconsole_V9.gif

23パイソンconsole_V10.gif

24パイソンconsole_V11.gif

25パイソンconsole_V12.gif

26パイソンconsole_V13.gif

27パイソンconsole_V14.gif

パイソンラスタ1枚
28パイソンconsole_R.gif

スタイル保存1枚
29スタイル保存.gif

スタイリングレイヤ1枚
30スタイリングレイヤー.gif

コードエディタ2枚
31コードエディタ.gif

32コードエディタ2.gif

p.219の上半分で挫折です。

いつか必要になったら,Pythonでがんがんスクリプトが書けるようにまた勉強することにします。

219/251ページ読了

2019/5/6(月)16:37記



posted by Dr.koba at 16:37| Comment(0) | QGIS(無料GISソフト) | このブログの読者になる | 更新情報をチェックする

2019年05月05日

Learn QGIS 第5章

この章も苦しかった!

まずはクリップ(切出し)ラスタを6枚
1クリップラスタ.gif

2クリップラスタ2.gif

3クリップラスタ3.gif

4クリップラスタ4.gif

5クリップラスタ5.gif

6クリップラスタ6.gif

DEMを13枚
7DEM.gif

8DEM2.gif

9DEM3.gif

10DEM4.gif

11DEM5.gif

12DEM6.gif

13DEM7.gif

14DEM8.gif

15DEM9.gif

16DEM10.gif

17DEM11.gif

18DEM12.gif

19DEM13.gif

土地被覆を2枚
20土地被覆reclass.gif

21土地被覆reclass2.gif

ラスタ2ベクタを2枚
22ラスタ2ベクタ.gif

23ラスタ2ベクタ2.gif

ベクタ2ラスタを3枚
24ベクタ2ラスタ.gif

25ベクタ2ラスタ2.gif

26ベクタ2ラスタ3.gif

フィールド統計を6枚
27フィールド統計.gif

28フィールド統計2.gif

29フィールド統計3.gif

30フィールド統計4.gif

31フィールド統計5.gif

32フィールド統計6.gif

ヒートマップを1枚
33ヒートマップ.gif

空港アクセス性を1枚
34空港アクセス性.gif

凸包を2枚
35凸包.gif

36凸包不可.gif

空港バッファを3枚
37空港バッファ.gif

38空港バッファx川.gif

39空港バッファx川2.gif

プロセス履歴を1枚
40プロセス履歴.gif

QGIS+GRASSを2枚
41QGIS+grass.gif

42QGIS+grass2.gif

六角ポリゴンを1枚
44六角ポリゴン.gif

ジオメトリの修正を1枚
45ジオメトリの修正.gif

川の選択を1枚
46川の選択.gif

川のバッファを2枚
47川のバッファ.gif

48川のバッファ2.gif

川沿いの森林を1枚
49川沿いの森林.gif

ディゾルブを2枚
50ディゾルブ.gif

51ディゾルブ2.gif

川沿いの森林率を1枚
52川沿いの森林率.gif

最後に,グラフィカルモデラーを15枚
53グラフィカルモデラー.gif

54グラフィカルモデラー2.gif

55グラフィカルモデラー3.gif

56グラフィカルモデラー4.gif

57グラフィカルモデラー5.gif

58グラフィカルモデラー6.gif

59グラフィカルモデラー7.gif

60グラフィカルモデラー8.gif

61グラフィカルモデラー9.gif

62グラフィカルモデラー10.gif

63グラフィカルモデラー11.gif

64グラフィカルモデラー12.gif

65グラフィカルモデラー13.gif

66グラフィカルモデラー14.gif

67グラフィカルモデラー15.gif






posted by Dr.koba at 11:03| Comment(0) | QGIS(無料GISソフト) | このブログの読者になる | 更新情報をチェックする

2019年05月03日

Learn QGIS 第4章

この章も手強かった。
サンプル通りの地図が作れなかった。

ラベリング
1ラベリング.gif

ラベリング_バッファなし
2ラベリング_バッファなし.gif

バッファあり
3ラベリング_バッファあり.gif

オフセットあり
4ラベリング_オフセットあり.gif

オフセットなし
5ラベリング_オフセットなし.gif

手動移動2枚
6ラベリング_手動移動.gif

7ラベリング_手動移動2.gif

ラベリング変更
8ラベリング_変更.gif

空港名ー標高
9ラベリング_空港名ー標高.gif

背景なし
10ラベリング_背景なし.gif

背景あり
11ラベリング_背景あり.gif

ラベリング_川_平行
12ラベリング_川_平行.gif

川_水平
13ラベリング_川_水平.gif

川_マージなし
14ラベリング_川_マージなし.gif

川_マージあり
15ラベリング_川_マージあり.gif

川_dissolve4枚
16川_dissovle.gif

17川_dissovle2.gif

18川_dissovle3.gif

19川_dissovle4.gif

川_ラップなし
20川_ラップなし.gif

川_ラップあり
21川_ラップあり.gif

川_ラベルフリー
22川_ラベルフリー.gif

地図作成11枚
23地図作成1.gif

24地図作成2.gif

25地図作成3.gif

26地図作成4.gif

27地図作成5.gif

28地図作成6.gif

29地図作成7.gif

30地図作成8.gif

31地図作成9.gif

32地図作成10.gif

33地図作成11.gif

スタイルの読み込み2枚
34スタイルの読み込み.gif

35スタイルの読み込み2.gif

グリッド追加2枚
36グリッド追加.gif

37グリッド追加2.gif

属性テーブル
38属性テーブル.gif

Atlas(地図帳)4枚
39Atlas.gif

40Atlas2.gif

41Atlas3.gif

42Atlas4.gif

webMAP3枚
43webmap.gif

44webmap2.gif

45webmap3.gif

3DwebMAP6枚
46_3dwebmap.gif

47_3dwebmap2.gif

48_3dwebmap3.gif

49_3dwebmap4.gif

50_3dwebmap5.gif

51_3dwebmap_EDGE.gif
グーグルchromeで3Dwebマップが表示できないのは残念。
Microsoftエッジは普段使ってないので残念。






posted by Dr.koba at 09:43| Comment(0) | QGIS(無料GISソフト) | このブログの読者になる | 更新情報をチェックする

2019年05月01日

Learn QGIS 第3章

平成から令和にまたがってしまった。

付録データのダウンロード4枚
1付録データDL.gif

2付録データDL2.gif

3付録データDL3.gif

4付録データDL4.gif

レイヤースタイル7枚
5レイヤースタイル1.gif

6レイヤースタイル2.gif

7レイヤースタイル3.gif

8レイヤースタイル4.gif

9レイヤースタイル5.gif

10レイヤースタイル6.gif

11レイヤースタイル7.gif

ストレッチ4枚
12ストレッチ1.gif

13ストレッチ2.gif

14ストレッチ3.gif

15ストレッチ4.gif

スタイリング
16スタイリング.gif

透明度100%と0%の2枚
17不透明度100%.gif

18不透明度0%.gif

乗算なしとありの2枚
19通常.gif

20乗算.gif

スタイル保存3枚
21スタイル保存.gif

22スタイル保存2.gif

23スタイル保存3.gif

ポイントスタイル4枚
24ポイントスタイル.gif

25ポイントスタイル2.gif

26ポイントスタイル3.gif

27ポイントスタイル4.gif

SVG5枚
28SVG.gif

29SVG2.gif

30SVG3.gif

31SVG4.gif

32SVG5.gif

シンボルレベル無効と有効の2枚
33シンボルレベル無効.gif

34シンボルレベル有効.gif

自作した川のシンボル
35シンボル_自作川.gif

透明な塗りつぶし
36透明な塗りつぶし.gif

湖と森
37LakesTrees.gif

ShapeBurstFill
38ShapeBurstFill.gif


ページ数は2章より少なかったけど手強かった!





posted by Dr.koba at 11:10| Comment(0) | QGIS(無料GISソフト) | このブログの読者になる | 更新情報をチェックする