Matlab处理后的数据写为二进制用于GrADS画图——fwrite函数使用

daqixuezi

感谢楼主，小白一枚受益匪浅

哈密小冬瓜

请教楼主 ，144*73*365这种三维的数据该怎么导出成dat数据，其中的循环改怎么写

Aires

哈密小冬瓜 发表于 2016-10-14 17:43
请教楼主 ，144*73*365这种三维的数据该怎么导出成dat数据，其中的循环改怎么写

对时间循环就可以了。
for i=1:365
     fwrite...
end

哈密小冬瓜

Aires 发表于 2016-10-14 19:24
对时间循环就可以了。
for i=1:365
     fwrite...

之前看到说是要先把三维的扩展成四维的格式，然后再进行写入吗？批量写入37年的做时间循环，不知道哪里错了
for i=1:37
    data=[];
    data=B{1,i};
    day=size(data,3);
    %C=reshape(data,[144 73 1 day]);
     str=strcat('e:\north\preciwaterindex\npwi.',num2str(i+1978),'.dat');
       fid=fopen(str,'wb');
          for m=1:day
           count=fwrite(fid,data(:,:,m),'float32');
          end
       fclose(fid);
end

Aires

哈密小冬瓜 发表于 2016-10-14 19:40
之前看到说是要先把三维的扩展成四维的格式，然后再进行写入吗？批量写入37年的做时间循环，不知道哪里错 ...

看起来没啥问题啊，哪里报错了？

哈密小冬瓜

Aires 发表于 2016-10-14 20:12
看起来没啥问题啊，哪里报错了？

错误使用 fwrite
文件标识符无效。使用 fopen 生成有效的文件标识符。

Aires

哈密小冬瓜 发表于 2016-10-14 20:19
错误使用 fwrite
文件标识符无效。使用 fopen 生成有效的文件标识符。

str=strcat('e:\north\preciwaterindex\npwi.',num2str(i+1978),'.dat');
改为：
str=['e:\north\preciwaterindex\npwi.',num2str(i+1978),'.dat'];
试一下

哈密小冬瓜

Aires 发表于 2016-10-14 22:47
str=strcat('e:\north\preciwaterindex\npwi.',num2str(i+1978),'.dat');
改为：
str=['e:\north\preci ...

谢楼主 还有一个小问题 对于二维的经度*纬度的mat数据要转置再处理 请问这种三维数据需要转置吗

Aires

哈密小冬瓜 发表于 2016-10-15 08:25
谢楼主 还有一个小问题 对于二维的经度*纬度的mat数据要转置再处理 请问这种三维数据需要转置吗

这要看你二维数据是怎么排列的了。fwrite是按列写入，然后你还要根据GrADS的二进制文件的存储格式来看你的矩阵怎么处理。你自己查一下fwrite的使用方式，以及GrADS存二进制文件的格式。

dovely

ocean university of china。。。