回答此问题可获得 20 贡献值,回答如果被采纳可获得 50 分。
<p>将此作为源:<a href="https://stackoverflow.com/questions/39881171/how-to-concatenate-monthly-trmm-netcdf-files-into-a-single-netcdf-file-using-nco">How to concatenate monthly TRMM netCDF files into a single netCDF file using NCO or R on windows 7?</a></p>
<pre><code>install.packages("ncdf4")
library(ncdf4)
install.packages("abind")
library(abind)
install.packages("RNetCDF")
library(RNetCDF)
install.packages("ncdf.tools")
library(ncdf.tools)
filenames=read.csv('TRMM.filenames.csv',head=F)
filenames=as.character(filenames[,1])
n.lon=4
n.lat=7
NA.matrix=matrix(rep(NA,n.lon*n.lat),nrow=n.lon)
prcp=array(NA.matrix,c(n.lon,n.lat,1))
for (i in 1:length(filenames)){ncdata=nc_open(filenames[i])
+ nc=ncvar_get(ncdata,"precipitation") prcp=abind(prcp,nc)}
prcp=prcp[,,-1]
dim(prcp)
saveRDS(prcp,'TRMM.all.rds')
</code></pre>
<p>我可以创建一个rds文件。但是,我很想把它保存为nc文件。我尝试创建一个新的netCDF文件,其中包含12个步骤(每个月一个)时间维度,方法是:</p>
^{pr2}$
<p>现在每个月的降水量都是挑战
在第一个脚本之后,我尝试使用<code>ncvar_put</code>函数,但没有成功。在</p>
<pre><code>filenames1=read.csv('TRMM.filenames.csv',head=F)
filenames1=as.character(filenames1[,1])
for (i in 1:length(filenames1)){ncdata1=nc_open(filenames1[i])
nc1=ncvar_get(ncdata1,"precipitation")
prcp1=abind(prcp1,nc1)}
n.lon1=4
n.lat1=7
data2d<-(4*7)
for (i in 1:length(filenames1))
ncvar_put( precip.nccreate, precip.ncvar, data2d, start=c(1), count=c(1) )
precip.nccreate<- nc_create( "precip.nccreate.nc", precip.ncvar, force_v4=FALSE, verbose=FALSE )
</code></pre>
<p>我得到了</p>
<blockquote>
<p>Error in ncvar_put(precip.nccreate, precip.ncvar, data2d, start = c(1), :
object 'precip.nccreate' not found</p>
<p>Error in nc_create("precip.nccreate.nc", precip.ncvar, force_v4 = FALSE, : object 'precip.ncvar' not found</p>
</blockquote>
<p>总之,我想我只是想找到一种简单的方法,将多个netcdf文件连接到一个netcdf中。在</p>
<p>谢谢</p>