用MATLAB读取和编写nc文件

是兜兜啊

clear all;close all;clc;
%% 创建nc文件
ncid = netcdf.create('D:\BOB\marine\uv_cmems_lev1.nc','CLOBBER');    %创建一个存放数据的nc文件

%% 定义维度参数
    % dimid = netcdf.defDim(ncid,dimname,dimlen)   
    dimidx = netcdf.defDim(ncid,'lon',121);
    dimidy = netcdf.defDim(ncid,'lat',109);
    dimidz = netcdf.defDim(ncid,'time',2048);
   
%% 定义新的变量
     varid = netcdf.defVar(ncid,'uvel','double',[dimidx dimidy dimidz]);
    varid2 = netcdf.defVar(ncid,'vvel','double',[dimidx dimidy dimidz]);%%不知道为什么之前写为nc_short怎么也没有办赋值，很难受
    % varid = netcdf.defVar(ncid,varname,xtype,dimids)
    varid3 = netcdf.defVar(ncid,'lon','NC_FLOAT',[dimidx]);
    varid4 = netcdf.defVar(ncid,'lat','NC_FLOAT',[dimidy]);
    varid5 = netcdf.defVar(ncid,'time','double',[dimidz]);
   
%% 新变量的属性
     netcdf.putAtt(ncid,varid,'units','m s-1');                                                     
     netcdf.putAtt(ncid,varid,'long_name','Eastward velocity');
     netcdf.putAtt(ncid,varid2,'units','m s-1');
     netcdf.putAtt(ncid,varid2,'long_name','Northward velocity');
    netcdf.putAtt(ncid,varid3,'units','degress_east');
    netcdf.putAtt(ncid,varid3,'long_name','Longitude');
    netcdf.putAtt(ncid,varid3,'valid_min','= 80');
    netcdf.putAtt(ncid,varid3,'valid_max','= 90');
    netcdf.putAtt(ncid,varid4,'units','degress_north');
    netcdf.putAtt(ncid,varid4,'long_name','Latitude');
    netcdf.putAtt(ncid,varid4,'valid_min','= 4');
    netcdf.putAtt(ncid,varid4,'valid_max','= 13');
    netcdf.putAtt(ncid,varid4,'step','= 0.083336');
    netcdf.putAtt(ncid,varid5,'units','hours since 1950-01-01 00:00:00');
    netcdf.endDef(ncid);
   
%% 给新变量赋值
%读取数据
    file1='D:/BOB/data/uv/u0/u0_cmems_lev1.nc';
    file2='D:/BOB/data/uv/v0/v0_cmems_lev1.nc';
    u0=ncread(file1,'uo');
    v0=ncread(file2,'vo');
    lon0=ncread(file1,'longitude');
    lat0=ncread(file1,'latitude');
    t=ncread(file1,'time');
    u0=squeeze(u0);v0=squeeze(v0);
    u1=u0(:,:,1:2048);v1=v0(:,:,1:2048);
    t1=t(1:2048,1);
   
%赋值
    netcdf.putVar(ncid,varid,u1);
    netcdf.putVar(ncid,varid2,v1);
    netcdf.putVar(ncid,varid3,lon0);
    netcdf.putVar(ncid,varid4,lat0);
    netcdf.putVar(ncid,varid5,t1);
    netcdf.close(ncid);
   
%% 检查新的nc文件
file3='D:\BOB\marine\uv_cmems_lev1.nc';
% ncdisp(file3);
u3=ncread(file3,'uvel');
u31=squeeze(u3(1,1,:));
v3=ncread(file3,'vvel');
v31=squeeze(v3(1,1,:));
lon1=ncread(file3,'lon');
lat1=ncread(file3,'lat');
u_t=ncread(file3,'time');

是兜兜啊

忘了，给数据去除异常值，无语   

ssh(ssh==6.103701889514923e-04)=nan;
ssh(ssh==-6.103701889514923e-04)=nan;

wangzhaoying

谢谢楼主！嘻嘻

香草拿铁

嚯
这个好
matlab处理nc格式数据简单易学

lliwanling

楼主你好，请问如果一个四维数组（经度，纬度，时间，高度），对时间维循环写入数据，但是每个时间维对应的高度维数据已知，该怎么写入呀~

是兜兜啊

lliwanling 发表于 2020-3-1 00:08
楼主你好，请问如果一个四维数组（经度，纬度，时间，高度），对时间维循环写入数据，但是每个时间维对应的 ...

你用ncdisp读取原数据每个变量的特征，具体的高度你可以作为ssh或者其他变量的属性写进去就好了吧，我也是个瓜皮，不知道建议是否有用了

海那边的大脸猫

楼主您好，我有一个500M的nc文件，我把数据读入Matlab，再编写nc文件把数据保存出来，结果新的nc文件变成了2G，xtype也没有改，不知是哪里有问题，请问您有遇到过这种情况吗？

是兜兜啊

海那边的大脸猫 发表于 2020-3-11 10:34
楼主您好，我有一个500M的nc文件，我把数据读入Matlab，再编写nc文件把数据保存出来，结果新的nc文件变成了 ...

没有哎，你是按照我的程序写的吗？要不你程序贴出来看一下？菜鸟互相讨论一下，哈哈哈哈哈

海那边的大脸猫

是兜兜啊 发表于 2020-3-11 13:47
没有哎，你是按照我的程序写的吗？要不你程序贴出来看一下？菜鸟互相讨论一下，哈哈哈哈哈

哈哈，谢谢你啊，我后来发现程序没有问题，我发现还是数据取的范围大了，把范围缩小后，就变为几百兆了。没想到就扩大了一点范围，所占的内存大了这么多

肖申克的梦

楼主你好请问如果变量time的数据是字符串的话该怎么设置维度？