基于R语言读取NC数据并提取单个波段（netcdf转为tiff或其他格式）

ncfile = ncdf4::nc_open('E:\\GOSIF_GPP_2001.M01_Mean.nc')

names(ncfile$var)

#该方法生成的raster仅能包含一个波段，band表示读取的波段为nc文件中的第多少个波段

input_nc = 'E:\\GOSIF_GPP_2001.M01_Mean.nc'

varname = 'GPP'

nc2raster = raster(input_nc,varname = varname,band = 1)

#数据可视化检查，参考：【R语言】栅格数据快速制图-levelplot函数

png("E:\\plot1.png",

    height = 15,

    width = 20,

    units = 'cm',

    res = 1000)

print(levelplot(nc2raster))

dev.off()

#将整个nc数据所有波段导入stack

nc2raster = stack(input_nc,varname = varname)

#结果写入GeoTiff文件

output = 'E:\\1.tif'

writeRaster(nc2raster,output,format = 'GTiff',overwrite = TRUE)

#此处的overwrite必须为TRUE