R

ランドサット8号のデータ2

rLandsat8 , raster , sp , colorRamps パッケージ

記事「ランドサット8号のデータ」の続き
1
2
3
4
5
#データは/home/user/landsat8/LC81110352015214LGN00 にあるとする
#library(rLandsat8)
#setwd("/home/user/landsat8")
#product  <- "LC81110352015214LGN00"
#l <- ReadLandsat8(product)

記事「ランドサット7号のデータ」でやったように鳥取市周辺を切り出す。

可視光
1
2
3
4
5
6
7
8
9
10
11
12
13
#library(raster)
l8=stack(l[[2]]$blue,l[[2]]$green,l[[2]]$red)
Crop <- c(420000,430000,3926000,3934600)
tottoriRGB <- crop(l8,Crop)
#作成したデータを保存
#writeRaster(tottoriRGB, filename="tottori_RGB.tif")
#このデータを使うときには
#library("raster")
#band blue,green,red
tottoriRGB<- stack("tottori_RGB.tif")
#png("landsat8_004.png",width=1000,height=1000)
plotRGB(tottoriRGB,3,2,1,stretch="lin")
#dev.off()

brightness temperature,NDVI,LSWI
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
#計算して
btemp <- ToAtSatelliteBrightnessTemperature(l, band="tirs1")
NDVI <- ToNDVI(l)
LSWI <- ToLSWI(l)
l8=stack(btemp,NDVI,LSWI)
#Crop <- c(420000,430000,3926000,3934600)
tottoriDATA <- crop(l8,Crop)
#作成したデータを保存
#writeRaster(tottoriDATA, filename="tottori_DATA.tif")
#このデータを使うときには
#library("raster")
tottoriDATA<- stack("tottori_DATA.tif")
# 絶対温度を摂氏に変換してからプロット
library(colorRamps)
#png("landsat8_005.png",width=1000,height=1000)
plot(tottoriDATA[[1]]- 273.15,col=blue2green2red(30))
#dev.off()

1
2
3
4
#Normalized Difference Vegetation Index(NDVI)
#png("landsat8_006.png",width=1000,height=1000)
plot(tottoriDATA[[2]])
#dev.off()

1
2
3
4
#Land Surface Water Index(LSWI)
#png("landsat8_007.png",width=1000,height=1000)
plot(tottoriDATA[[3]])
#dev.off()