QGISによる衛星画像のカラー画像化．

[前回の続き](https://qiita.com/nigo1973/items/9bb6a11caac8e3e1e850)．

ダウンロードしたSentinel-2の衛星データからカラー画像を無料のGISソフトであるQGISを使って，カラー画像化します．

##1. 環境構築．

まずはQGISをダウンロードし環境を構築します．
私のPCはmacですので，OSX用のソフトをダウンロードします．

OS環境：mac OSX 10.12.6 macOS Sierra

QGISのホームページからダウンロードです．

[QGISホームページ](https://www.qgis.org/ja/site/)

ホームページのダウンロードをクリックし，自分の環境（Windows, MacOSXなど）にあわせたバージョンをダウンロードします．
私はmacなので，”QGIS macOS インストーラ Version 3.2”をクリックします．

ダウンロードし，パッケージを開くと，インストールの順番がついたファイルがあります．
QGISはPython環境で動作しているため，まずPython3.6をダウンロードしインストールします．

私はHomebrewでpython3.6環境を構築していましたが，これではQGISをインストールすることができませんでした．
インストールすると，”Python3.6が入っていません”とエラーメッセージがでてきます．

調べてみると，[Pythonのホームページからパッケージをダウンロードしインストールしないといけないみたい](https://qiita.com/ETA077/items/ec49b2e96a144be04f30)なので，
改めてPythonをインストールし直ししました．


そのあと，2.についているGDALをインストールし，3.QGISをインストールして環境構築終了です．

これで，QGIS　3.2.2-Bonnがインストールされました．

##2. QGISによる衛星画像のカラー化．
[QGISの使い方](https://qiita.com/tags/qgis)は他のページをご参考ください．
Googleで調べると多くの情報があります．

QGISのよいところは，多くのユーザがPluginという簡単な拡張機能を公開していること．
例えば，[GoogleMapを使いたいとき](http://d-and-i.hatenablog.com/entry/2016/10/04/090000)は，PluginでOpnelayersをインストールすれば使えます．

でも，QGISがVer.3 になってから，上記のPluginが見つかりませんでした．残念．
調べてみるとVe3になってから，今までのプラグインが使えなくなったみたいです．
これについては，また別記事で紹介します．

まずは，先程ダウンロードしたSentinel-2の画像をQGISで表示させます．
今回は，SentinelのData Hubよりダウロードした地震後の北海道厚真町の画像を表示させてみました．

![QGIS2.png](https://qiita-image-store.s3.amazonaws.com/0/187602/28f8bd69-8aaa-2626-93b4-5c62cb2249cc.png)

フォルダにあるB02, B03, B04のファイルを選んで，QGISのワークスペースにドラッグドロップします．

次に，この3枚の画像を合成しカラー化させます．これは，"Build Vertual Raster"を使います．

![QGIS4.png](https://qiita-image-store.s3.amazonaws.com/0/187602/4c5c0551-365e-53e0-6c3e-78ba050f118c.png)

上部のメニューから　ラスタ→その他→Build Vertual Rasterを選択します．

![QGIS5.png](https://qiita-image-store.s3.amazonaws.com/0/187602/f5c03433-f440-ef6d-babd-0f0adb9db1ad.png)

入力レイヤの右の選択ボタンから，ドラッグした3つのファイルを選択します．

![QGIS6.png](https://qiita-image-store.s3.amazonaws.com/0/187602/be2aec15-3f26-9373-f9ae-bf6bd706cf8f.png)

そして，下にある"バックグラウンドで実行"をクリックし，処理を実行します．
ここで，エラーが発生する場合があるます．どうも，２番めにインストールしたGDALのパスが設定されていないため，
その実行ファイルが見つからないことが原因みたいです．

このエラーが発生した場合は，[このページ](https://qiita.com/Yfuruchin/items/e54c41af09a83d406ac2)をご参考ください．
私はこの方法で解決できました（感謝）．

うまく実行されるとこのページが表示されます．
![QGIS7.png](https://qiita-image-store.s3.amazonaws.com/0/187602/ec1120b9-bee1-7732-2460-784c61e352c6.png)

ただ，みると分かる通り，青みがかった画像になっています．赤と青の画像の配分がうまくいっていないのが原因です．
そのため，赤色にB04，青色にB02を配分させます．

左下のウィンドウにできた”仮想”のレイヤで右クリックしてプロパティを選択します．

![QGIS8.png](https://qiita-image-store.s3.amazonaws.com/0/187602/a7026c7d-6565-6f6f-605f-4a4e747bbf09.png)

そして，プロパティウィンドウから，RGBのバンドの配分を変更します．

![QGIS9.png](https://qiita-image-store.s3.amazonaws.com/0/187602/44d9e8d2-c3f0-167c-8714-b76c94ed807a.png)


Redバンドにバンド３を，Blueバンドにバンド１を選び，”OK”をクリックします．
すると，先程の画像が以下になります．

![QGIS10.png](https://qiita-image-store.s3.amazonaws.com/0/187602/dff6c6f5-8e2d-ec67-676f-6b472b9db83f.png)

青みがかった画像から，いつもの衛星画像になったのが分かるかと思います．
北海道の厚真町の土砂崩れが報告されている箇所をみると，以下のようになっているのがわかります．
![QGIS11.png](https://qiita-image-store.s3.amazonaws.com/0/187602/2d5ac48f-9069-3b14-2fd2-edb6af049c0b.png)

厚真町の地域の山で，多くの土砂崩れが発生していたことがこの画像よりわかります．
心よりお見舞い申し上げますとともに，被害にあわれたみなさまの一日も早い再建をお祈り致します。

##3. おわりに
今回は，ダウンロードしたSentinel2の画像のカラー化まででした．

次回は，ここにGoogoleMapや国土地理院の地図情報のプラグインの方法について紹介します．

（追記：2019年3月26日）
　閲覧数が多く、みなさんの関心が高いことを知りました。
　近々、電波衛星（SAR satellite)であるSentinel−1のデータについても、紹介いたします。