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ソフト) | このブログの読者になる | 更新情報をチェックする