NCL用存2维矩阵的方式存储3维矩阵

白兰

最近刚刚开始学习使用NCL，想用NCL存储3维矩阵，并写成txt文件输出，继而用matlab处理数据并画图。但存储数据量巨大的3D矩阵用write_table是不太好用的，我比较喜欢用write_matrix，用写二维矩阵的方式写三维矩阵，

最后用matlab编写ncl_reshape函数把维度重新提取出来，我认为这样比较方便，编辑好后以后一直可以用。
如果各位朋友有更好的方法，还请多多指教，谢谢~~^_^
重新提取维度的matlab的程序稍后上传，请在本站搜索“ncl_reshape提取ncl三维矩阵”。

兰溪之水

为啥不直接存成二进制格式或者NetCDF格式。。。

白兰

接上一贴“[经验总结] NCL用存2维矩阵的方式存储3维矩阵"
最近刚刚开始学习使用NCL，想用NCL存储3维矩阵，并写成txt文件输出，继而用matlab处理数据并画图。但存储数据量巨大的3D矩阵用write_table是不太好用的，我比较喜欢用write_matrix，用写二维矩阵的方式写三维矩阵，
最后用matlab编写ncl_reshape函数把维度重新提取出来，我认为这样比较方便，编辑好后以后一直可以用。
如果各位朋友有更好的方法，还请多多指教，谢谢~~^_^
重新提取维度的matlab的程序请见附件。（因为matlab内置的reshape函数与我NCL存储的矩阵规则不太一样，所以需要单独写此函数）

白兰

兰溪之水 发表于 2014-6-27 10:39
为啥不直接存成二进制格式或者NetCDF格式。。。

我刚刚开始学习使用NCL，以前都习惯用matlab，对您说的二进制格式和存netCDF格式的数据用的不是很顺，我还是个菜鸟，不过准备慢慢学习，之后想写成netCDF格式，毕竟matlab能直接读nc文件。受教了~

二爷名声在外

下次麻烦楼主能一个帖子说完的不要分成两个帖子，除非是必须分割不可的！多谢支持！

白兰

二爷名声在外 发表于 2014-6-27 10:49
下次麻烦楼主能一个帖子说完的不要分成两个帖子，除非是必须分割不可的！多谢支持！

好的，抱歉我第一次发帖，没意识到这个问题，我以为matlab和ncl的程序是要分版块发的。

白兰

function newdata=ncl_reshape(data,ny,nx,nlevel)
    newdata=zeros(nx,ny,nlevel);
    for i=1:nlevel
        num=ny*(i-1)+1;
        temp=data(num:num+ny-1,:);
        newdata(:,:,i)=temp';
    end

end

cyf287411599

这个真心不错哦