码迷,mamicode.com
首页 > 其他好文 > 详细

NCEP CFSR数据读取

时间:2018-09-08 20:03:39      阅读:1254      评论:0      收藏:0      [点我收藏+]

标签:png   bsp   div   cfs   1.2   size   data   经纬度   win   

一、

NCEP CFSR再分析数据,时间分辨率是1小时,空间分辨率是 0.312 * 0.312。

1、整体读取数据情况

clear all
setup_nctoolbox
tic
%% 读取数据文件
wind= ncgeodataset(‘wnd10m.cdas1.201612.grb2‘);
wind_list = wind.variables;%文件的列表情况

技术分享图片

wind_list包含了

技术分享图片

2、u分量

size_of_u = wind.size(wind_list(1));%u分量的数据尺寸,777小时,1个高度,经纬度数据,4D数据
data_u=wind.geovariable(wind_list(1));%取得数据类型为风速u的数据
u_1=data_u.data(1,1,:,:); %
u_2=squeeze(u_1);

技术分享图片

3、v 分量

%% v 分量
size_of_v = wind.size(wind_list(2));%v分量的数据尺寸,777小时,1个高度,经纬度数据,4D数据
data_v=wind.geovariable(wind_list(2));%取得数据类型为风速v的数据
v_1=data_v.data(1,1,:,:); %
v_2=squeeze(v_1);

技术分享图片

4、

高度只有一个值,10米

%% 高度分量
size_of_h= wind.size(wind_list(5));%v分量的数据尺寸,777小时,1个高度,经纬度数据,4D数据
data_h=wind.geovariable(wind_list(5));%取得数据类型为风速v的数据
v_1=data_h.data(1); %高度10米

技术分享图片

5、经度

 

6、纬度

  

7、时间分量

NCEP 一小时采集一次数据,所以 12月份 共 31*24=744 个数据

%% 时间数据
size_of_time = wind.size(wind_list(6));%u分量的数据尺寸,777小时,1个高度,经纬度数据,4D数据
data_time=wind.geovariable(wind_list(6));%取得数据类型为风速u的数据
time=data_time(1:size_of_time);

技术分享图片

技术分享图片

 

  

  

 

 

  

NCEP CFSR数据读取

标签:png   bsp   div   cfs   1.2   size   data   经纬度   win   

原文地址:https://www.cnblogs.com/ruo-li-suo-yi/p/9610160.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!