とらりもん - フィルタ処理(2017/07/11 石橋) Diff

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

! !! [[スペックルノイズの低減]]を参照に戻る．

!! フィルタ処理(2017/07/11  石橋)
JAXA/ALOS-2の[[サイト|http://www.eorc.jaxa.jp/ALOS-2/doc/sam_jindex.htm]]から[[東京・高分解能3mモード (HH偏波)の画像|ftp://ftp.eorc.jaxa.jp/pub/ALOS-2/1501sample/020_tokyo/0000008656_001001_ALOS2014410740-140829.zip]]をダウンロードする。

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|http://pen.agbi.tsukuba.ac.jp//home/ryuiki_plus/public_html/sample_data/PALSAR2/IMG-HH-ALOS2014410740-140829-UBSL2.1GUA.tif.png]][[PALSAR|http://pen.envr.tsukuba.ac.jp//home/ryuiki_plus/public_html/sample_data/PALSAR2/IMG-HH-ALOS2014410740-140829-UBSL2.1GUA.tif.png]]

→ [[PolSAR入門]]に戻る。