将WRF文件转换为光栅堆栈！投影开发者社区

相关文章推荐

一身肌肉的板栗 · 极米h3设置教程 - 抖音· 5 月前 ·

玩手机的饭卡 · 科技新闻· 6 月前 ·

威武的大象 · 迦勒底神谕译文：残篇一 | 机核 GCORES· 6 月前 ·

帅呆的炒粉 · 最“软骨头”的世界500强，东风汽车竟靠捡洋 ...· 1 年前 ·

苦恼的冲锋衣 · 漫画全集在线观看 - 哔哩哔哩漫画· 1 年前 ·

腾讯云

开发者社区

文档意见反馈控制台

首页 TVP

文章/答案/技术大牛

发布

ro_rast <- nc_open("wrf_CAM0_daily_pre.nc") 
pre <- ncvar_get(ro_rast, "pre") ro_rast$dim$lon$vals -> lon ro_rast$dim$lat$vals -> lat ro_rast$dim$ncl2$vals -> time rm(ro_rast) 
r1_brick <- brick(pre, xmn=min(lat), xmx=max(lat), ymn=min(lon), ymx=max(lon), crs=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs+ towgs84=0,0,0")) 
names(r1_brick)<- seq(as.Date('2018-06-01'), as.Date('2018-08-31'), 'days') 
# convert names of layer into date par(mar = c(2, 2, 2, 2)) 
cam1_mean <- t(calc(r1_brick, sum)) 
# seasonal sum precipitation
 cam1 <- flip(cam1_mean, direction = 2)
library(akima)# intepolation 
lonlat_reg <- expand.grid(lon = seq(min(lon), max(lon), length.out = 1499), 
                          lat = seq(min(lat), max(lat), length.out = 749))
test <- interp(x = as.vector(lon), y = as.vector(lat), z = as.vector(pre), 
               xo = unique(lonlat_reg[,"lon"]), yo = unique(lonlat_reg[,"lat"]), 
               duplicate = "error", linear = FALSE, extrap = FALSE)
test <- interp(x = as.vector(lon), y = as.vector(lat), z = as.vector(pre),