matlab2010b有关NCEP数据的个人经验与心得

sdzxlk

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

WEN

楼主，我在读取NETCDF文件过程中，想要求变量的和，prec2010 = prec2010_single(:,:,1:59），这条命令可以用来计算第1-59天的降水量， 那如何求不连续的数据之和，比如1月，6月和12月，这3个月的降水量之和，程序是要如何修改呢？

斥鷃

先把日数据处理成逐月的，并且扩维，前阵子刚好有个程序可以参考

1. yrstart=1979;
   
2. yrend=2014;
   
3. N=yrstart-yrend+1;
   
4. Q1_m=nan(134,54,27,12,N);
   
5. for k=yrstart:yrend
   
6.    data=ncread(['NCEP1.daily.Q1.' num2str(k) '.nc'],'Q1');
   
7.    m=size(data);
   
8.    if m(4)==365
   
9.       list=[31 28 31 30 31 30 31 31 30 31 30 31];
   
10.    else
    
11.       list=[31 29 31 30 31 30 31 31 30 31 30 31];
    
12.    end
    
13.       ind=1;
    
14.       for n=1:12
    
15.        Q1_m(:,:,:,n,k-yrstart+1)=mean(data(:,:,:,ind:ind+list(n)-1),4);
    
16.        ind=ind+list(n);
    
17.       end
    
18. end

复制代码

有逐月数据以后取特点位置数据sum求和即可，如
Q_used=nansum(Q_m(:,:,:,[1 6 12],:),4);%如果想消去一维，可以参考squeeze函数

mmdd

总结的好仔细

WEN

斥鷃 发表于 2016-1-9 13:02
先把日数据处理成逐月的，并且扩维，前阵子刚好有个程序可以参考

有逐月数据以后取特点位置数据sum求和 ...

谢谢楼主!楼主好强大！总结的很是全面！仔细呢！  
楼主，现在数据已经处理好了，下一步要作图，用MATLAB。画出了全国范围内的。不知道怎样可以截取成区域的？比如西藏省的那种？目前手头有各省市的白化BLN图，不知道有没有用。还有GIS里导出的SHP文件，貌似是可以添加进去的。但是不知道如何设置成白底，而且出图区域好像并没有变的很小，还是按数据源的区域范围出图的呢。捣鼓了好几天了，求楼主热心解答呀！跪谢！

斥鷃

WEN 发表于 2016-1-12 10:17
谢谢楼主!楼主好强大！总结的很是全面！仔细呢！  
楼主，现在数据已经处理好了，下一步要作图，用MATLA ...

matlab中地图边界与掩膜(去掉边界外区域)的实现（基于shape文件）
http://bbs.06climate.com/forum.p ... 40&fromuid=9098
(出处: 气象家园)

WEN

斥鷃 发表于 2016-1-13 18:16
matlab中地图边界与掩膜(去掉边界外区域)的实现（基于shape文件）
http://bbs.06climate.com/forum.php? ...

谢谢楼主啦。楼主附的帖子我也在坛子里看过，很浩大的一个工程呀。我再仔细琢磨琢磨哈！谢谢！

pr711f

感谢楼主分享！

shinexzz

很受用 谢谢~

shinexzz

很有用！！