Matlab拼接（merge）两个nc文件

最新推荐文章于 2023-03-25 19:01:13 发布
最新推荐文章于 2023-03-25 19:01:13 发布
阅读量3.2k 收藏 16
点赞数 4
分类专栏： MATLAB+NC
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jerry003/article/details/110295117
版权
创建 nc 文件主要分为以下几个步骤
  1. define input and output directories
  2. define dimensions
  3. define variables
  4. read nc files and merge variables（实际中可跳过）
  5. put variables into nc
  6. put attribute into nc
  7. define global attribution

拼接前的 ncdisp 结果为：

>> ncdisp('domain.lnd.140x80pt_ECMWF_noocean.201129_b22.nc')
Source:
           I:\merge_China_srfdata\domain.lnd.140x80pt_ECMWF_noocean.201129_b22.nc
Format:
           classic
Global Attributes:
           title              = 'CESM domain data:'
           Conventions        = 'CF-1.0'
           source_code        = 'SVN $Id: gen_domain.F90 65202 2014-11-06 21:07:45Z mlevy@ucar.edu $'
           SVN_url            = ' $URL: https://svn-ccsm-models.cgd.ucar.edu/tools/mapping/gen_domain/trunk/src/gen_domain.F90 $'
           Compiler_Optimized = 'TRUE'
           hostname           = ''
           history            = 'created by +1.0000000000000000'
           source             = 'map_140x80pt_ECMWF_noocean_to_140x80pt_ECMWF_nomask_aave_da_201129.nc'
           map_domain_a       = '/datahuge/wdy_cesmfiles/cesm_2_2_0/components/clm/tools/mkmapgrids/SCRIPgrid_140x80pt_ECMWF_noocean_c201129.nc'
           map_domain_b       = '/datahuge/wdy_cesmfiles/cesm_2_2_0/components/clm/tools/mkmapgrids/SCRIPgrid_140x80pt_ECMWF_nomask_c201129.nc'
           map_grid_file_ocn  = '/datahuge/wdy_cesmfiles/cesm_2_2_0/components/clm/tools/mkmapgrids/SCRIPgrid_140x80pt_ECMWF_noocean_c201129.nc'
           map_grid_file_atm  = '/datahuge/wdy_cesmfiles/cesm_2_2_0/components/clm/tools/mkmapgrids/SCRIPgrid_140x80pt_ECMWF_nomask_c201129.nc'
Dimensions:
           n  = 11200
           ni = 140
           nj = 80
           nv = 4
Variables:
    xc  
           Size:       140x80
           Dimensions: ni,nj
           Datatype:   double
           Attributes:
                       long_name = 'longitude of grid cell center'
                       units     = 'degrees_east'
                       bounds    = 'xv'
    yc  
           Size:       140x80
           Dimensions: ni,nj
           Datatype:   double
           Attributes:
                       long_name = 'latitude of grid cell center'
                       units     = 'degrees_north'
                       bounds    = 'yv'
    xv  
           Size:       4x140x80
           Dimensions: nv,ni,nj
           Datatype:   double
           Attributes:
                       long_name = 'longitude of grid cell verticies'
                       units     = 'degrees_east'
    yv  
           Size:       4x140x80
           Dimensions: nv,ni,nj
           Datatype:   double
           Attributes:
                       long_name = 'latitude of grid cell verticies'
                       units     = 'degrees_north'
    mask
           Size:       140x80
           Dimensions: ni,nj
           Datatype:   int32
           Attributes:
                       long_name   = 'domain mask'
                       note        = 'unitless'
                       coordinates = 'xc yc'
                       comment     = '0 value indicates cell is not active'
    area
           Size:       140x80
           Dimensions: ni,nj
           Datatype:   double
           Attributes:
                       long_name   = 'area of grid cell in radians squared'
                       coordinates = 'xc yc'
                       units       = 'radian2'
    frac
           Size:       140x80
           Dimensions: ni,nj
           Datatype:   double
           Attributes:
                       long_name   = 'fraction of grid cell that is active'
                       coordinates = 'xc yc'
                       note        = 'unitless'
                       filter1     = 'error if frac> 1.0+eps or frac < 0.0-eps; eps = 0.1000000E-11'
                       filter2     = 'limit frac to [fminval,fmaxval]; fminval= 0.1000000E-02 fmaxval=  1.000000'

拼接要求：将两个140x80的网格，拼成一个280x80的网格区域，原来文件中的所有变量都进行合并。

clear all; clc
% define input and output directories
inDir   = 'I:\merge_China_srfdata';
outDir  = 'I:\merge_China_srfdata';

inFile_b21  = 'domain.lnd.140x80pt_ECMWF_noocean.201128_b21.nc';
inDirFile_b21 = [inDir, '\', inFile_b21];
inFile_b22  = 'domain.lnd.140x80pt_ECMWF_noocean.201129_b22.nc';
inDirFile_b22 = [inDir, '\', inFile_b22];

outFile = 'domain.lnd.280x80pt_South_China_ECMWF_noocean.201201.nc';
outDirFile = [outDir, '\', outFile];
outid = netcdf.create(outFile,'64BIT_OFFSET');

% define dimensions
ndimID  = netcdf.defDim(outid, 'n', 11200);
nidimID = netcdf.defDim(outid, 'ni', 280);
njdimID = netcdf.defDim(outid, 'nj', 80);
nvdimID = netcdf.defDim(outid, 'nv', 4);

% define variables
xc_id   = netcdf.defVar(outid, 'xc', 'NC_DOUBLE', [nidimID, njdimID]);
yc_id   = netcdf.defVar(outid, 'yc', 'NC_DOUBLE', [nidimID, njdimID]);
xv_id   = netcdf.defVar(outid, 'xv', 'NC_DOUBLE', [nvdimID, nidimID, njdimID]);
yv_id   = netcdf.defVar(outid, 'yv', 'NC_DOUBLE', [nvdimID, nidimID, njdimID]);
mask_id = netcdf.defVar(outid, 'mask', 'NC_INT', [nidimID, njdimID]);
area_id = netcdf.defVar(outid, 'area', 'NC_DOUBLE', [nidimID, njdimID]);
frac_id = netcdf.defVar(outid, 'frac', 'NC_DOUBLE', [nidimID, njdimID]);

netcdf.endDef(outid);  % 重要，结束定义模式！

% read nc files and merge variables
% dayon comment:  cat(dimensionID, f01, f02), dimensionID 表示拼接的维度，顺序是从左至右，f01在前，f02在后。
xc_b21 = ncread(inDirFile_b21, 'xc');
xc_b22 = ncread(inDirFile_b22, 'xc');
xc_cmb = cat(1, xc_b21, xc_b22);  

yc_b21 = ncread(inDirFile_b21, 'yc');
yc_b22 = ncread(inDirFile_b22, 'yc');
yc_cmb = cat(1, yc_b21, yc_b22);

xv_b21 = ncread(inDirFile_b21, 'xv');
xv_b22 = ncread(inDirFile_b22, 'xv');
xv_cmb = cat(2, xv_b21, xv_b22);

yv_b21 = ncread(inDirFile_b21, 'yv');
yv_b22 = ncread(inDirFile_b22, 'yv');
yv_cmb = cat(2, yv_b21, yv_b22);

mask_b21 = ncread(inDirFile_b21, 'mask');
mask_b22 = ncread(inDirFile_b22, 'mask');
mask_cmb = cat(1, mask_b21, mask_b22);

area_b21 = ncread(inDirFile_b21, 'area');
area_b22 = ncread(inDirFile_b22, 'area');
area_cmb = cat(1, area_b21, area_b22);

frac_b21 = ncread(inDirFile_b21, 'frac');
frac_b22 = ncread(inDirFile_b22, 'frac');
frac_cmb = cat(1, frac_b21, frac_b22);

% put variables into nc
netcdf.putVar(outid, xc_id, xc_cmb);
netcdf.putVar(outid, yc_id, yc_cmb);
netcdf.putVar(outid, xv_id, xv_cmb);
netcdf.putVar(outid, yv_id, yv_cmb);
netcdf.putVar(outid, mask_id, mask_cmb);
netcdf.putVar(outid, area_id, area_cmb);
netcdf.putVar(outid, frac_id, frac_cmb);

netcdf.reDef(outid); % 重要，重新进入定义模式！

% put attribute into nc
netcdf.putAtt(outid, xc_id, 'longname', 'longitude of grid cell center');
netcdf.putAtt(outid, xc_id, 'units', 'degrees_east');
netcdf.putAtt(outid, xc_id, 'bounds', 'xv');

netcdf.putAtt(outid, yc_id, 'longname', 'latitude of grid cell center');
netcdf.putAtt(outid, yc_id, 'units', 'degrees_north');
netcdf.putAtt(outid, yc_id, 'bounds', 'yv');

netcdf.putAtt(outid, xv_id, 'longname', 'longitude of grid cell verticies');
netcdf.putAtt(outid, xv_id, 'units', 'degrees_east');

netcdf.putAtt(outid, yv_id, 'longname', 'latitude of grid cell verticies');
netcdf.putAtt(outid, yv_id, 'units', 'degrees_north');

netcdf.putAtt(outid, mask_id, 'longname', 'domain mask');
netcdf.putAtt(outid, mask_id, 'note', 'unitless');
netcdf.putAtt(outid, mask_id, 'coordinates', 'xc yc');
netcdf.putAtt(outid, mask_id, 'comment', '0 value indicates cell is not active');

netcdf.putAtt(outid, area_id, 'longname', 'area of grid cell in radians squared');
netcdf.putAtt(outid, area_id, 'coordinates', 'xc yc');
netcdf.putAtt(outid, area_id, 'units', 'radian2');

netcdf.putAtt(outid, frac_id, 'longname', 'fraction of grid cell that is active');
netcdf.putAtt(outid, frac_id, 'coordinates', 'xc yc');
netcdf.putAtt(outid, frac_id, 'note', 'unitless');
netcdf.putAtt(outid, frac_id, 'filter1', 'error if frac> 1.0+eps or frac < 0.0-eps; eps = 0.1000000E-11');
netcdf.putAtt(outid, frac_id, 'filter2', 'limit frac to [fminval,fmaxval]; fminval= 0.1000000E-02 fmaxval=  1.000000');

% define global attributions
netcdf.putAtt(outid, netcdf.getConstant('NC_GLOBAL'), 'title', 'CESM domain data');
netcdf.putAtt(outid, netcdf.getConstant('NC_GLOBAL'), 'Conventions', 'CF-1.0');
netcdf.putAtt(outid, netcdf.getConstant('NC_GLOBAL'), 'source_code', 'SVN $Id: gen_domain.F90 65202 2014-11-06 21:07:45Z mlevy@ucar.edu $');
netcdf.putAtt(outid, netcdf.getConstant('NC_GLOBAL'), 'SVN_url', '$URL: https://svn-ccsm-models.cgd.ucar.edu/tools/mapping/gen_domain/trunk/src/gen_domain.F90 $');
netcdf.putAtt(outid, netcdf.geetConstant('NC_GLOBAL'), 'Compiler_Optimized', 'TRUE');
netcdf.putAtt(outid, netcdf.getConstant('NC_GLOBAL'), 'hostname', '');

netcdf.close(outid);
大雨海深
关注
专栏目录
Matlab批量裁剪nc文件并合并
qq_41651338的博客
09-29 238
【代码】Matlab批量裁剪nc文件并合并。
将多个月的nc数据文件合并成一个（月平均）
weixin_43259541的博客
02-04 9900
将多个月的nc数据文件合成一个（月平均） 最近在处理nc数据，希望可以获得多年的月平均数据。但是ERA这个产品从网站上下载到的是每个月的日数据，所以我希望可以把他合成一个，以每个月月平均数据存储的nc文件。 首先是引入要用的库： from netCDF4 import Dataset #nc库功能不用说，读nc文件肯定要用 import numpy as np from pandas import Series import netCDF4 as nc import tkinter as tk from
matlab函数merge_MATLAB数据合并方法
weixin_39777404的博客
12-19 3780
在数据处理当中，我们经常碰到数据集合并的需要。在一些语言如SAS，SQL中，数据集合并是一个非常简单的问题，但是在Matlab中，数据集合并则显得稍微复杂了一点。在Matlab中，要想合并数据，有两个办法，第一个办法就是使用dataset这种数据类型，这个statistictoolbox中的数据类型，专门为统计分析开发的。mat,cell,等数据类型转换成dataset数据类型可以查看mat2d...
matlab将多个nc文件进行合成TIF
PeanutbutterBoh的博客
04-07 1547
多个nc转TIF
matlab函数merge_如何在MATLAB中合并这些数据？
weixin_39638057的博客
12-19 894
In the sample text file below, if column 3 contains a 1 then the corresponding data of column 2 should be merged with the data of the previous row in column 2. For example, the 40 in row 2 should be a...
多个netcdf文件的合并（matlab）
schumacher2016的博客
09-26 8586
说明：下载CCMP再分析风场，一个nc文件存储一天的数据，而ECMWF-Interim再分析风场则是一个nc文件存储一个月的数据，有时候我们需要的是一个nc文件里头包含很多个月的数据，这时候就需要把每天或者每月的数据合并到一个nc文件里头了。 工具：matlab 2010b netcdf包 代码如下： %% % desciption: merge multiple netcdf file...
MATLAB 矩阵合并函数
08-20
本代码主要利用MATLAB工具实现矩阵合并函数，简单明了，易于理解
Merge PDF-Documents：将多个pdf文档合并为一个pdf文档-matlab开发
05-31
此函数使用 apache pdfbox（包含在 Matlab 中）将多个 pdf 文件合并为一个文档。 用法： inputFiles = {'foo.pdf', 'bar.pdf}; outputFileName = 'foobar.pdf'; mergePdfs(inputFiles, outputFileName);
matlab两表格合并,在matlab中将两个单元合并到一个单个单元中(merge two cell in one single cell in matlab)...
weixin_29205061的博客
03-16 4434
在matlab中将两个单元合并到一个单个单元中(merge two cell in one single cell in matlab)我有两个牢房。 其中一个是Cell 1='1007_s_at' 780 'DDR1''1053_at' 5982 'RFC2''117_at' 3310 'HSPA6''121_at' 7849 'PAX8''1255_g_at'...
Matlab合并排序法程序
04-19
使用matlab编写的合并排序法，merge_main为示例主函数，merge_func为调用排序的求解函数
ncmerge:将多个不规则的NetCDF网格合并到等矩形网格中。-开源
04-29
ncmerge可以非常快速地将基于单元格的不规则网格合并到单个等矩形（平板Carree）网格中，以简化可视化工作流程。 它首先从坐标网格创建地图文件，然后在重新捕获数据值时执行快速查找。 支持3D时变变量。
matlab函数merge_Merge into函数用法详解及示例
weixin_39929813的博客
12-29 1163
Oracle在9i中引入了merge命令，通过这个 merge into语句，能够在一个SQL语句中对一个表同时执行update和insert操作。当然是update还是insert得依据于你的指定的条件判断的，mergeinto 可以实现用B表来更新A表数据，如果A表中有匹配的记录，则更新数据；没有匹配的记录，则把B表的记录插入A表。merge命令可以从一个或多个数据源中选择行来update或i...
合并nc文件在同一时间内的不同变量（CDO）
weixin_44119544的博客
09-22 2009
使用语句：cdo merge ifile1 ifile2 ifile3 ofile 备注：ifile 为输入文件 应该具有不同的变量，但是这三个文件应处在同一时间（时间步长）;ofile 为输出文件名（可自定义）
Python与Matlab的数据合并
xiaotian的博客
03-25 1072
Python与Matlab数据合并
Python不同文件夹下nc文件的整合
lulingxiaoc的博客
03-13 1997
在一个主文件夹中有许许多多的子文件夹（文件夹1，2，3，4 , , , n），每个子文件夹下都会存储一个nc文件（或多个），每天凌晨12：00呢又会更新子文件夹（子文件夹n+1）。能否每天都能将最后倒数七个子文件夹中的nc文件整合成一个nc（通过时间序列）？
MERGE 用法
u013000942的博客
10-28 370
MERGE INTO SYS_USER_POS PS   USING (SELECT U.USERID, M.ORG_ID, M.POS_ID            FROM SYS_USER U           INNER JOIN MID_SYS_USER_PMP_SYNC M ON U.ACCOUNT = M.ACCOUNT WHERE M.CHECK_STATUS = '0')
多个NC时间序列文件合并为一个nc文件
qq_41057480的博客
10-13 6423
【代码】多个NC时间序列文件合并为一个nc文件。
matlab合并多个nc文件
最新发布
04-07
要合并多个nc文件，可以使用MATLAB中的ncread和ncwrite函数。以下是一个示例代码，演示如何读取多个nc文件，并将它们合并成一个单独的nc文件。 ```matlab % 输入文件名 filename1 = 'file1.nc'; filename2 = 'file2.nc'; filename3 = 'file3.nc'; output_filename = 'merged_file.nc'; % 读取第一个文件的变量名和维度信息 vars1 = ncread(filename1, 'var'); dims1 = ncread(filename1, 'dim'); % 创建新文件并写入第一个文件的变量和维度信息 ncid = netcdf.create(output_filename, 'NC_WRITE'); for i = 1:length(dims1) dimid(i) = netcdf.defDim(ncid, dims1{i}, length(ncread(filename1, dims1{i}))); end for i = 1:length(vars1) varid(i) = netcdf.defVar(ncid, vars1{i}, 'double', dimid); end netcdf.endDef(ncid); for i = 1:length(vars1) netcdf.putVar(ncid, varid(i), ncread(filename1, vars1{i})); end % 读取并追加第二个和第三个文件的变量和维度信息 vars2 = ncread(filename2, 'var'); dims2 = ncread(filename2, 'dim'); vars3 = ncread(filename3, 'var'); dims3 = ncread(filename3, 'dim'); for i = 1:length(dims2) if ~ismember(dims2{i}, dims1) dimid(length(dims1)+i) = netcdf.defDim(ncid, dims2{i}, length(ncread(filename2, dims2{i}))); end end for i = 1:length(vars2) if ~ismember(vars2{i}, vars1) varid(length(vars1)+i) = netcdf.defVar(ncid, vars2{i}, 'double', dimid(length(dims1)+find(strcmp(dims2, ncread(filename2, vars2{i}))))); end end for i = 1:length(dims3) if ~ismember(dims3{i}, [dims1 dims2]) dimid(length(dims1)+length(dims2)+i) = netcdf.defDim(ncid, dims3{i}, length(ncread(filename3, dims3{i}))); end end for i = 1:length(vars3) if ~ismember(vars3{i}, [vars1 vars2]) varid(length(vars1)+length(vars2)+i) = netcdf.defVar(ncid, vars3{i}, 'double', dimid(length(dims1)+length(dims2)+find(strcmp(dims3, ncread(filename3, vars3{i}))))); end end netcdf.endDef(ncid); for i = 1:length(vars2) netcdf.putVar(ncid, varid(length(vars1)+i), ncread(filename2, vars2{i})); end for i = 1:length(vars3) netcdf.putVar(ncid, varid(length(vars1)+length(vars2)+i), ncread(filename3, vars3{i})); end % 关闭文件 netcdf.close(ncid); ``` 在这个示例代码中，我们首先读取第一个nc文件的变量名和维度信息，并创建一个新的nc文件。然后，我们读取第二个和第三个nc文件的变量和维度信息，并将它们追加到新文件中。注意，我们需要检查变量和维度是否已经存在于新文件中，如果存在，就不需要重复定义。最后，我们将所有变量写入新文件中，并关闭文件。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

大雨海深

感谢您的支持和鼓励

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值