我有WRF输出的netCDF文件与1499_749_91尺寸生产的“墨卡托”投影在非洲之角。我想把netCDF文件转换成光栅堆栈来进行进一步的分析。我一直在尝试不同的选择，但对我来说没有用。我在错误的地点得到了价值。我在这方面需要帮助，任何帮助都是非常感谢的。

以下是代码：

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),