ERA5-land数据NC转tif,然后按照经纬度提取数据——基于R

最新推荐文章于 2024-08-31 11:30:20 发布

Fei Zhu

最新推荐文章于 2024-08-31 11:30:20 发布

阅读量1.7k

点赞数

分类专栏： R 文章标签： java 数据库开发语言

本文链接：https://blog.csdn.net/qq_29665133/article/details/126716770

版权

R 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

这段代码展示了如何使用R语言的raster和ncdf4库读取.nc文件，提取特定变量，并将数据提取到指定的经纬度。代码中提到的变量名(varname)需要用户自行确定，最终目的是将数据转换为TIFF格式，但实际转换部分未在代码中实现。

摘要由CSDN通过智能技术生成

抛砖引玉，for循环版的需要用的人写吧！

才开始我想转为tif的，所以有savename这个变量。

  library(raster)
  library(ncdf4)
  # NC文件所在额路径
  path <- "/Users/feizhu/Desktop/PRMS/era"
  filename <- paste0(path,"/",list.files(path,pattern = ".nc$"))
  
  #savename <- gsub(".nc",".tif",filename)
  #varname 变量名，可以通过nc_open()函数确定这个变量
  
  raster = stack(filename[1], varname = varname)
  
  
  aa = t(extract(raster,c(lon,lat)))
  # 数据第一列是时间，第二列是这个经纬度的对应的数据