如下为R 程序:
library(ncdf4)
library(gstat)
library(sp)
nc = nc_open('g:/chla/modis_chl_a_4km_monthly/2016/A20163362016366.L3m_MO_CHL_chlor_a_4km.nc')
# chla1 = ncatt_get(nc,'chlor_a')
chla = ncvar_get(nc,'chlor_a',c(7129,1177),c(10,30))
lon = ncvar_get(nc,'lon',c(7129),c(10))
lat = ncvar_get(nc,'lat',c(1177),c(30))
nc_close(nc)
class(chla)
cc = data.frame(chla)
cc[is.na(cc)] = 0
data(cc)
class(cc)
names(cc)
coordinates(cc) <- c('x','y')
class(cc)
我想读取nc文件然后将缺测值插值,在网上搜的代码,但是在前面就卡住了,我想把矩阵chla转化呈一个"data.frame"形式,但是data(cc)时报错:
> data(cc)
Warning message:
In data(cc) : 没有‘cc’这个数据集
不知道为什么 而且设置坐标时也报错:
> coordinates(cc) <- c('x','y')
Error in `[.data.frame`(object, , value) : undefined columns selected
有没有大神知道为什么啊 ~~~
看了好久了还是不明白为什么