とらりもんHOME  Index  Search  Changes  Login

大気地形効果補正プログラム

2016/06/17 Jin Katagi 書きかけです

AtmosphericCorrection

弘前大学の飯倉先生が公開した、オープンソースの大気地形効果補正プログラム。

大気地形補正には衛星画像のほか、「土地被覆図」「空間基盤地図」が必要。

ツールは6SV(radiative transfer code), Python 2.X, Pythonのライブラリはnumpy,cv2,scipy,osgeo,gdalが必要。

インストール

gitから必要なものをダウンロード

git clone https://github.com/y-iikura/AtmosphericCorrection.git

データの準備

以下、ディレクトリ構造は以下のようになっているとする

work/6SV
    /GSI
    /Landsat
    /AtmosphericCorrection

衛星画像の準備

Earth ExplorerでLandsat7 ETM+のデータをダウンロード。

とりあえずつくば周辺のデータをダウンロード。 Entity ID:LE71070352016005ASN00 Acquisition Date:05-JAN-16 Path:107 Row:35

mv work/Landsat
cd work/Landsat
# 解凍
tar xzvf LE71070352016005ASN00.tar.gz

土地被覆図の準備

JAXAの高解像度土地利用土地被覆図ホームページからダウンロード。

衛星画像に合わせて、つくば周辺(N36E140)のタイル(LC_N36E140.tif)をダウンロード。

# 必要な箇所だけ取り出す
gdal_translate -projwin 140.0000000 36.2500000 140.2500000 36.0833333 LC_N36E140.tif dummy.tif
# UTMに変換
gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3100 dummy.tif LC_N36E140_UTM.tif

                    

空間基盤図の準備

国土地理院の基盤地図情報サイトからダウンロード。

データのダウンロードには登録が必要なので、事前に登録しておく。

基盤地図情報のダウンロード -> 基盤地図情報数値標高モデル -> 10 m メッシュ -> 5446 -> 544010,544011,544020,544021を選択

すべてチェック -> まとめてダウンロード

ダウンロードしたものはunzipしておく。

tifへの変換

ダウンロードしたものはxmlファイルなので、Geotiffに変換する。

mkdir xml2dem && cd xml2dem
wget http://space.geocities.jp/bischofia_vb/python/fgddem/fgddem-20140128.zip
unzip fgddem-20140128.zip
rm fgddem-20140128.zip
# DLしたツールをxmlデータのあるディレクトリにコピー
cp ./fgddem/fgddem.py ../GSI/
cd ../GSI/
# まとめてtifに変換
ls -1 FG-GML*.xml | xargs -I{} python2 fgddem.py {}
# UTMに変換 (EPSG:4612に注意)
for dem_file in `ls -1 FG-GML*.tif`; do
  gdalwarp  -s_srs EPSG:4612 -t_srs EPSG:3100 ${dem_file} ${dem_file%.tif}_UTM.tif
done

# 4つのtifを結合する (output:dem.tif)
gdalbuildvrt merge.vrt *_UTM.tif
gdal_translate -a_srs epsg:3100 merge.vrt dem.tif
rm merge.vrt

6SVのダウンロード

サイト

mkdir 6SV && cd 6SV
wget http://6s.ltdri.org/files/update/6SV2.1.tar
tar xvf 6SV2.1.tar
# インストール
make

プログラムの実行 (書きかけです)。

衛星画像に合わせてパラメータを設定していく。

Last modified:2016/06/17 21:09:33
Keyword(s):
References:[memo_Katagi]