尘埃cd 2024-03-03 22:54 采纳率: 25%
浏览 16
已结题

R语言读取nc按月均值转为tif

现在拿到一个土壤nc文件,一年有365个波段,波段名字就是1-365
使用brick方法读取的时候,crop方法一直报错,按raster方法读取就可以执行crop。
现在裁剪完的clipdata,想问怎么按月均值提取出来,导出成每个月的tif文件?
或者也可以转化为我怎么从一个多图层raster中取子集出来做mean统计


```r
library(ncdf4) # 加载ncdf4包
library(raster)
library(terra)

###裁剪矢量
cmr = vect("F:/extratAreaW84jwd.shp")

path <- "d:/datadownload/SMCI_1km_10cm.nc"
soilwater<-rast(path )
clipdata <- soilwater %>% crop(cmr) %>% mask(cmr) #裁剪范围

```

  • 写回答

25条回答 默认 最新

  • 喝茶品人生 2024-03-04 10:08
    关注
    获得1.00元问题酬金
    monthly_means <- lapply(1:12, function(month) {
      month_indices <- seq(month, 365, by = 12)
      monthly_mean <- calc(clipdata, fun = function(x) mean(x[month_indices]))
      return(monthly_mean)
    })
    
    for (i in 1:12) {
      writeRaster(monthly_means[[i]], filename = paste0("month_", i, "_mean.tif"), format = "GTiff")
    }
    
    
    评论

报告相同问题?

问题事件

  • 系统已结题 3月11日
  • 修改了问题 3月4日
  • 修改了问题 3月4日
  • 创建了问题 3月3日

悬赏问题

  • ¥15 Yolov5-7➕mobilebetv3,报错list index out of range
  • ¥15 有没有人会打学生成绩管理系统呀
  • ¥15 在使用Fiddler和夜神模拟器抓包的时候一直出现443该怎么办啊QAQ搜了好几个笔记都没有解决
  • ¥15 3x7的二维数组A、B、C,A中的任意1个数组元素与B的任意1个数组元素、同时又与C的任意1个数组元素比较,把不同位置出现相同数的比较称为无意义,反之称为有意义,把有意义的比较打印输出。
  • ¥20 预测模型怎么处理原始数据(随机森林)
  • ¥20 请问discuz3.5如何实现插入ckplayer全能播放器功能呢?
  • ¥15 thingsboard代码编译出错误
  • ¥15 博途v18仿真报错怎么解决
  • ¥15 欧姆龙plc枕式包装机 ST编程
  • ¥15 为啥快手广告联盟的广告这么难出来