【R】GRASSGIS内のラスターデータの表示

めも。


GRASSGISのコマンドラインからRを起動。
library(spgrass6)
library(raster)


#GRASSから読み込み。
day_g=readRAST6("raster_name")
#rasterデータに変換
day=raster(day_g)

#余白の調整。タイトルをいれるので上側を少し空ける
par(oma=c(0,0,2,0),mfrow=c(1,1))
plot(day,axes=F);title("TITLE")

#ちなみにspplotを使っても作図可能。
#spplot(day_g,col.regions=terrain.colors(30))

#作図。pdfで保存
pdf("Filename.pdf")
plot(day,axes=F);title("TITLE")
dev.off()