とらりもんHOME  Index  Search  Changes  Login

フィルタ処理

2017/05/14 Jin Katagi (2017/07/11 石橋 書きかけです)

スペックルノイズの低減に戻る.

フィルタ処理(2017/07/11 石橋)

JAXA/ALOS-2のサイトから東京・高分解能3mモード (HH偏波)の画像をダウンロードする。

unzip 0000008656_001001_ALOS2014410740-140829.zip 

#画像をGRASSに取り込む。
grass #GRASS7を使用・UTM54
r.in.gdal in=IMG-HH-ALOS2014410740-140829-UBSL2.1GUA.tif out=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA

#とりあえず表示させてみる。
g.region rast=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA
r.colors map=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA color=grey.eq
d.mon wx0
d.rast IMG_HH_ALOS2014410740_140829_UBSL2.1GUA

#DN値を後方散乱係数に変換
r.mapcalc expression="IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_backscat = 10*log(IMG_HH_ALOS2014410740_140829_UBSL2.1GUA*IMG_HH_ALOS2014410740_140829_UBSL2.1GUA,10)-83"

#3×3ピクセルによる平滑化フィルター処理
r.mapcalc --overwrite "X=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_backscat"

r.mapcalc "IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_9=(X[0,1]+X[1,0]+X[0,-1]+X[-1,0]+X[1,1]+X[1,-1]+X[-1,1]+X[-1,-1]+X[0,0])/9"

#5×5ピクセルによる平滑化フィルター
r.mapcalc "IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_25=(X[-2,2]+X[-1,2]+X[0,2]+X[1,2]+X[2,2]+X[-2,1]+X[-1,1]+X[0,1]+X[1,1]+X[2,1]+X[-2,0]+X[-1,0]+X[0,0]+X[1,0]+X[2,0]+X[-2,-1]+X[-1,-1]+X[0,-1]+X[1,-1]+X[2,-1]+X[-2,-2]+X[-1,-2]+X[0,-2]+X[1,-2]+X[2,-2])/25"

#出力
r.out.gdal type=Float32 in=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_9 out=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_9.tif -f
r.out.gdal type=Float32 in=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_25 out=IMG_HH_ALOS2014410740_140829_UBSL2.1GUA_25.tif -f

PALSAR

PolSAR入門に戻る。

Last modified:2017/07/11 10:26:05
Keyword(s):
References:[スペックルノイズの低減]