爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3130|回复: 3

[程序设计] ncread读取数据在纬向上格点的排序问题

[复制链接]

新浪微博达人勋

发表于 2016-6-19 13:55:34 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
初学matlab,请问一下大家,这样读出来的变量var([1:lon_num],[1:lat_num],[1:time_num]),第二维[1:lat_num]是从南向北排序的还是反过来排的呢?
  1. lon1=70;%longitude of the west boundry.
  2. lon2=150;%longitude of the east boundry.
  3. lat1=0;%latitude of the south boundry.
  4. lat2=60;%latitude of the north boundry.

  5. year1=1961;%year begin.
  6. year2=2014;%year end.
  7. %-------------------------------
  8. lonw=lon1/2.5+1;
  9. lone=lon2/2.5+1;
  10. lats=(90-lat1)/2.5+1;
  11. latn=(90-lat2)/2.5+1;
  12. %--------------------------------
  13. lon=ncread('E:\Extreme_preci\data\ncep\mon\slp.mon.mean.nc','lon');
  14. lat=ncread('E:\Extreme_preci\data\ncep\mon\slp.mon.mean.nc','lat');
  15. time=ncread('E:\Extreme_preci\data\ncep\mon\slp.mon.mean.nc','time');
  16. a=find(lon>=lon1 & lon<=lon2);
  17. b=find(lat>=lat1 & lat<=lat2); % a和b作为经度和纬度的指针
  18. d=find(time>=1411296 & time<=1883904);
  19. lon_num=length(a);
  20. lat_num=length(b);
  21. time_num=length(d);
  22. var=ncread('E:\Extreme_preci\data\ncep\mon\slp.mon.mean.nc','slp',[a(1) b(1) d(1)],[lon_num lat_num time_num]);
复制代码
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-6-19 21:27:58 | 显示全部楼层
求大神指点
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-6-20 09:11:13 | 显示全部楼层
一般是从南向北,你可以ncdisp一下该nc的文件头说明,或者直接输出lat维看看不就知道了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-6-20 11:37:29 | 显示全部楼层
读出来查看一下不就知道了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表