经验 - R - Raster Package

R中的Raster Package提供了对栅格数据类型的大量操作， 包括读取、输出、栅格计算等。

1. writeRaster

此函数将Raster对象输出为外部文件，可以为binary、ASCII、GeoTiff等多种格式。

注意输出binary文件时，为 行优先模式 输出。

2. 绘制raster

Raster Package示例代码：

# 基本绘制命令
r <- raster(nrows=10, ncols=10)
r <- setValues(r, 1:ncell(r))
plot(r)

# 在一个图片里绘制两个栅格，不同的数值范围，同样的尺度。
r[] <- runif(ncell(r))
r2 <- r/2
brks <- seq(0, 1, by=0.1) 
nb <- length(brks)-1 
cols <- rev(terrain.colors(nb))
par(mfrow=c(1,2))
plot(r, breaks=brks, col=cols, lab.breaks=brks, zlim=c(0,1), main='first') 
plot(r2, breaks=brks, col=cols, lab.breaks=brks, zlim=c(0,1), main='second')