||
# zhufei 2019-3-24 rm(list = ls()) library(ncdf4) setwd('C:/Users/zhufe/Desktop/R_spatial') # 设置nc文件所在的路径 nc <- nc_open("E4E81D21036EAC6D84FEA2521AA84E3D_ferret_listing.nc") pcp <- ncvar_get( nc = nc, varid = 'pcp') lon <- ncvar_get( nc = nc, varid = 'lon') lat <- ncvar_get( nc = nc, varid = 'lat') # 建立索引 id <- rep(1:10,each = 8) # 根据pcp第三纬的长度,然后除以8,确定1:10这个变量 #id <- id[1:60623] # 如果不能整除,变成整除然后取子集 e <- data.frame(id) # 建立数据框,第一列是id for (i in 1:80) { # 数据的重新整理 变成一列 pre <- as.matrix(pcp[,,i]) # 第三纬每一纬变成矩阵 dim(pre) <- c(100*70,1) # 按照经纬度重塑数据 e[i,2:7001] <- pre # 保存数据 } # 每八列求和 a <- data.frame() for (j in 1:10) { # pcp第三维的长度确定j/8的长度 a[j,1:7001] <- colSums(subset(e,id == j), na.rm = T,dims = 1) # 每八行计算一天的降水量 }
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 08:15
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社