如何解决将.nc文件加载到R中,转换为保留波段名称的栅格
我对使用.nc文件是陌生的,并且正在努力获取我需要的GeoTIFF格式的输出以上传到GEE。
我已经从http://data.ceda.ac.uk/neodc/esacci/fire/data/burned_area/AVHRR-LTDR/grid/v1.0/Compressed的CCI档案中下载了火数据的2015年zip文件。并编写了一个循环以加载.nc文件并将其导出为GeoTIFF,这样我就可以在GEE中使用它们了。但是,当我在GEE中查看它们时,它们没有乐队名称。其次,我期望有11个波段,因为这里有11个变量详细说明:https://catalogue.ceda.ac.uk/uuid/4f377defc2454db9b2a6d032abfd0cbd。但是,GEE中有18个。所以我想知道我是否能够正确加载所有频段? (我需要全部使用)。最后,我不确定数据中是否包含时间变量。每个.nc文件代表一个月,我需要将它们全部合并到一个GEE图像集中,然后可以按日期进行过滤。
下面显示了我的基本代码,我将非常感谢任何人提供的任何帮助或指导!!谢谢!
### convert .nc files to .tif
### load libraries ###
library(ncdf4)
library(raster)
library(rgdal)
setwd("filepath\\Archive Fire CCI grid\\2015\\")`
### list all files in the folder
ncfiles <- list.files()
ncfileTest <- brick(ncfiles[1])
ncfileTest
##loop over every file in the folder,exporting tif with simplified name
for(i in ncfiles){
ncfile <- brick(i)
## file name doesnt agree with GEE so needs to be cleaned
cleanname <- gsub("-","_",i)
cleanname2 <- gsub("_AVHRR_LTDR_fv1.0.nc","",cleanname)
writeRaster(ncfile,file.path("Output file path\\ESA CCI Fire Archive",cleanname2),format="GTiff",overwrite=TRUE)
}
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。