在R中使用if语句从栅格创建新栅格

library(raster)
raster_data <- raster("path_to_raster_file")

new_raster <- raster(nrows=ncell(raster_data), ncols=ncell(raster_data), 
                     extent(raster_data), res(raster_data))

for (i in 1:ncell(raster_data)) {
  if (raster_data[i] > threshold) {
    new_raster[i] <- 1  # 设置新栅格像素值为1
  } else {
    new_raster[i] <- 0  # 设置新栅格像素值为0
  }
}

writeRaster(new_raster, "path_to_output_raster_file", format="GTiff")