全部版块 我的主页
论坛 站务区 十一区 新手入门区
945 1
2018-05-29
如下为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

有没有大神知道为什么啊 ~~~
看了好久了还是不明白为什么

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2018-5-29 11:35:27
为什么我看到内容被删了啊   是不是没有人能看到啊
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群