R语言读取NC格式数据

数据准备

国家地球系统科学数据中心下载的全国逐月最低气温数据集，每3年数据共36个月的最低气温数据保存在一个nc文件中，现在想要将数据导出为逐月栅格数据，以15-17年数据为例。

R包加载

主要用到ncdf4包读取nc文件，raster包将读入数据转栅格数据并输出.

library(raster)
library(ncdf4)
library(rasterVis)
library(lattice)

NC数据导入

##读取单变量长时间序列nc数据
ncfile<-nc_open("E:/1data/tmn_2015_2017.nc")
names(ncfile$var)

names函数获取变量名tmn

导入为栅格数据

input_nc = 'E:/1data/tmn_2015_2017.nc'
varname = 'tmn'
##压缩呈栅格栈
raster1 = stack(input_nc, varname = varname)

设置输出路径名

outpath<-"E:/1data/月最低气温"
##为文件名设置做准备
x1<-c(1:12)
x2<-paste("2015",x1,sep = "")
x3<-paste("2016",x1,sep = "")
x4<-paste("2017",x1,sep = "")
x5<-c(x2,x3,x4)
out_filename<-sapply(x5, function(x)paste(x,".tif",sep = ""))      #输出tif格式
out_path<-sapply(out_filename, function(x){
  paste(outpath,x,sep = "/")})                #输出路径名

批量输出

for(i in 1:36){
  writeRaster(raster1[[i]],out_path[i], format = 'GTiff', overwrite = TRUE)
}

小结

学习R语言是一件激动人心的事情，在看了很多文章后终于可以站在巨人的肩膀上写一写自己的东西，自勉。
[文章参考]（https://www.jianshu.com/p/1f327f0c48e9）