R语言创建陆地掩模：代码解析

以下代码用于创建陆地掩模：

# 创建陆地掩模
mask <- r <- raster(dir$mask)
mask[] <- 1
mask[is.na(r[])] <- 0
mask <- rotate(mask)
plot(mask)

代码解释：

1. mask <- r <- raster(dir$mask): 读取栅格文件 (dir$mask) 并创建一个新的栅格对象 (mask)。
2. mask[] <- 1: 将所有像元值设置为1。
3. mask[is.na(r[])] <- 0: 将原始栅格对象 (r) 中的缺失值 (NA) 的像元值设置为0。
4. mask <- rotate(mask): 对栅格对象进行旋转。
5. plot(mask): 绘制掩模图像。

通过以上代码，您将创建一个代表陆地的掩模，其中陆地部分的像元值为1，非陆地部分的像元值为0。该掩模可用于其他地理空间分析任务，例如计算陆地覆盖面积或进行空间数据过滤。