micaps第11类数据怎么读取？

solitaire

方便streamslice()\quiver(x,y,u,v)等函数用

solitaire

雪お轩漠 发表于 2015-3-27 20:19
我觉得也是- -，本来觉得给太明显伸手的程序不太好的。。不过鉴于是以前那个一部分。。。直接抽出来粘贴 ...

主要之前没接触过，不过昨天看懂了  其他类型的数据也差不多能都出来了

雪お轩漠

function dia=diamond11(filename)
%此函数用于读取diamond11类型数据
%filename为要读取文件的路径，由主程序按照循环生成
%*********本程序内有部分数据未使用，但是读取过程已经包含，具体对应请查看
%Micaps数据结构，如有需要可自行设置接口使用
  fid=eval(['fopen(''' filename ''')']);
  Type=fscanf(fid,'%s%s%s',3);
  Time=fscanf(fid,'%d%d%d%d%d',6);
  Pointpara=fscanf(fid,'%f%f%f%f%f%f',6);
  Linepara=fscanf(fid,'%f%f',2);
  Temp=struct('year',[],'month',[],'day',[],'hour',[],'presc',[],...
      'lev',[],'londis',[],'latdis',[],'lonbeg',[],'lonend',[],...
      'latbeg',[],'latend',[],'latgrinum',[],'longrinum',[],...
      'datau',[],'datav',[]);
  Temp.year=Time(1);
  Temp.month=Time(2);
  Temp.day=Time(3);
  Temp.hour=Time(4);
  Temp.presc=Time(5);
  Temp.lev=Time(6);
  Temp.londis=Pointpara(1);
  Temp.latdis=Pointpara(2);
  Temp.lonbeg=Pointpara(3);
  Temp.lonend=Pointpara(4);
  Temp.latbeg=Pointpara(5);
  Temp.latend=Pointpara(6);
  Temp.latgrinum=Linepara(1);
  Temp.longrinum=Linepara(2);
  n=1;
  while n<=Linepara(2)
    Temp.datau(n,:)=fscanf(fid,'%f',Linepara(1));
    Temp.datav(n,:)=fscanf(fid,'%f',Linepara(1));
    n=n+1;
  end
  dia=Temp;
  fclose(fid);

供参考

雪お轩漠

另外我觉得，像这种数据格式的数据，应当在知道结构的基础上读取并不难。我觉得这种程序都求现成，怕是不太好呢。

rjys1219

学习了，很好哈

二爷名声在外

雪お轩漠 发表于 2015-3-27 14:16
function dia=diamond11(filename)
%此函数用于读取diamond11类型数据
%filename为要读取文件的路径，由 ...

绝对是做毕业设计的！

solitaire

雪お轩漠 发表于 2015-3-27 14:16
function dia=diamond11(filename)
%此函数用于读取diamond11类型数据
%filename为要读取文件的路径，由 ...

十分感谢  我学习一下

雪お轩漠

二爷名声在外 发表于 2015-3-27 17:07
绝对是做毕业设计的！

我觉得也是- -，本来觉得给太明显伸手的程序不太好的。。不过鉴于是以前那个一部分。。。直接抽出来粘贴了- -！。。不过都懒得解释了囧。。最近毕设真是翻天了要

amoschx

终于知道所谓的第*类型数据是怎么看出来的了，是不是diamond n,就是第n类数据啊？

dtczqb

第五类数据格式：用于TLOGP和剖面图  diamond  5  屏幕上需显示的内容  年  月  日  时次  总站点数 区站号  经度  纬度  跋海高度  单站内容长度  百帕  高度  温度  露点 风向  风速  百帕