将.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 (将#修改为@)