1、　全球/地区温度图

（1）读取数据

clear all

setup_nctoolbox %调用工具包

tic %计时

%%

nc=ncgeodataset('tmpsfc.gdas.199401.grb2');   %读文件

tem_1=nc.variables   %浏览数据类型

%%

a1=nc.geovariable(tem_1(1));%取得数据类型为Temperature_surface的数据

b1=a1.data(1,:,:); %第一个时间点温度数据

c1=squeeze(b1)-273.16;%删除单一维度，换为摄氏温度

%%

a2=nc.geovariable(tem_1(2));%取得数据类型为lat的数据，纬度

b2=a2.data(:,1)%提取数据

%%

a3=nc.geovariable(tem_1(3));%经度

b3=a3.data(:,1)%

%%

a4=nc.geovariable(tem_1(4));%取得数据类型为time的数据

b4=a4.data(:,1)%

%%

save A b2 b3 c1

toc

读取的是NCEP/CFSR数据，1994年1月的温度数据。

该数据一共四项。温度，纬度，经度，时间。温度中是3维数据组织形式

时间数据

经度数据

纬度数据

可以看出该数据集的数据组织形式，经度纬度构成世界地图，记录了744个时间点的温度数据。1小时采集一次数据。

温度数据

保存为mat格式数据，为画图做准备。

（2）画图

clear all

load A

[Plg,Plt]=meshgrid(b3',b2');%形成网格

 m_proj('hammer-aitoff','clongitude',-150);%投影模式

m_pcolor(Plg,Plt,c1);

shading flat;

colormap('jet');%颜色选择

hold on;

m_pcolor(Plg-360,Plt,c1);

shading flat; %着色模式

colormap('jet');

m_coast();

m_grid('xaxis','middle');

% h=colorbar('h');

% set(get(h,'title'),'string','1991年1月全球温度');

c=colorbar('southoutside','fontsize',12)

c.Label.String = '1994年1月全球温度';

c.Label.FontSize = 15;

（3）找出最大、最小温度的经纬度

clear all

load C

Max_col=max(c1);%列最大值

Max_row=max(c1,[],2)%行最大值

Max=max(max(c1));

[x1,y1]=find(c1==max(max(c1)));%x 行，y 列

T_1=Plt(x1,y1)%纬度

T_2=Plg(x1,y1)%经度

Min_col=min(c1);%列最大值

Min_row=min(c1,[],2)%行最大值

Min=min(min(c1));

[x,y]=find(c1==min(min(c1)));%x 行，y 列

T_x=Plt(x,y)%纬度

T_y=Plg(x,y)%经度

　

可以看出最热52度，在澳大利亚那块（142.8123E，23.2610S）；最冷-62度，在北极圈那块（89.9999E，66.3486N）。

（4）中国（地区）温度图

clear all

load A

LATLIMS=[3 54];

LONLIMS=[72 134];%选定边界范围

[Plg,Plt]=meshgrid(b3',b2');%形成网格

 m_proj('lambert','lon',LONLIMS,'lat',LATLIMS);%投影模式

m_pcolor(Plg,Plt,c1);

shading flat;

colormap('jet');%颜色选择

hold on;

m_pcolor(Plg-360,Plt,c1);

shading flat; %着色模式

colormap('jet');

m_coast();

m_grid('box','fancy','tickdir','in');

% h=colorbar('h');

% set(get(h,'title'),'string','1991年1月全球温度');

c=colorbar('southoutside','fontsize',12)

c.Label.String = '1994年1月中国温度';

c.Label.FontSize = 15;

该方法是读取全球数据，只展示部分

（5）

改进的区域方法，读取该区域的数据数据，

clear all

setup_nctoolbox %调用工具包

tic %计时

%%

min_lat=115;

max_lat=279;

min_lon=231;

max_lon=430;  %区域经纬度范围，在数据中的位置

%%

nc=ncgeodataset('tmp2m.gdas.199401.grb2');   %读文件

tem_1=nc.variables   %浏览数据类型

%%

N1=nc.size(tem_1(1));%读取数据大小,可以看出数据的组织形式

a1=nc.geovariable(tem_1(1));%取得数据类型为Temperature_surface的数据

b1=a1.data(1,1,min_lat:max_lat,min_lon:max_lon); %第一个时间点温度数据

c1=squeeze(b1)-273.16;%删除单一维度，换为摄氏温度

%%

N2=nc.size(tem_1(2));

a2=nc.geovariable(tem_1(2));%取得数据类型为lat的数据，纬度

b2=a2.data(min_lat:max_lat,1);%提取数据

%%

N3=nc.size(tem_1(3))

a3=nc.geovariable(tem_1(3));%经度

b3=a3.data(min_lon:max_lon,1);%

%%

N4=nc.size(tem_1(4))

a4=nc.geovariable(tem_1(4));%取得数据类型为time的数据

b4=a4.data(:,1);%

%%

N5=nc.size(tem_1(5));%读取数据大小

a5=nc.geovariable(tem_1(5));%取得数据类型为time的数据

b5=a5.data(:,1);%

%%

save tem b2 b3 c1

toc


clear all

load tem

LATLIMS=[3 54];

LONLIMS=[72 134];%选定边界范围

[Plg,Plt]=meshgrid(b3',b2');%形成网格

 m_proj('lambert','lon',LONLIMS,'lat',LATLIMS);%投影模式

m_pcolor(Plg,Plt,c1);

shading flat;

colormap('jet');%颜色选择

hold on;

m_pcolor(Plg-360,Plt,c1);

shading flat; %着色模式

colormap('jet');

m_coast();

m_grid('box','fancy','tickdir','in');

% h=colorbar('h');

% set(get(h,'title'),'string','1991年1月全球温度');

c=colorbar('southoutside','fontsize',12)

c.Label.String = '1994年1月中国温度';

c.Label.FontSize = 15;

（6）读取方式的改变

clear all

setup_nctoolbox %调用工具包

tic %计时

%%

min_lat=115;

max_lat=279;

min_lon=231;

max_lon=430;  %区域经纬度范围，在数据中的位置

%%

nc=ncgeodataset('tmp2m.gdas.199401.grb2');   %读文件

tem_1=nc.variables   %浏览数据类型

%%

N1=nc.size(tem_1(1));%读取数据大小

b1=nc.data(tem_1(1),[1,1,min_lat,min_lon],[1,1,max_lat,max_lon]);%初始的读取位置，最终的位置

c1=squeeze(b1)-273.16;

%%

N2=nc.size(tem_1(2));%读取数据大小

b2=nc.data(tem_1(2),[min_lat],[max_lat]);

%%

N3=nc.size(tem_1(3));%读取数据大小

b3=nc.data(tem_1(3),[min_lon],[max_lon]);

save tem b2 b3 c1

toc

可以看出海面2米的温度组织形式，4维数据，包含了744个时间点，1个位置，576个纬度点，1152个经度点。

读取方式是初始位置用一个数组表示，终止位置用一个数组表示

clear all

load tem

LATLIMS=[3 54];

LONLIMS=[72 134];%选定边界范围

[Plg,Plt]=meshgrid(b3',b2');%形成网格

 m_proj('lambert','lon',LONLIMS,'lat',LATLIMS);%投影模式

m_pcolor(Plg,Plt,c1);

shading flat;

colormap('jet');%颜色选择

hold on;

m_pcolor(Plg-360,Plt,c1);

shading flat; %着色模式

colormap('jet');

m_coast();

m_grid('box','fancy','tickdir','in');

% h=colorbar('h');

% set(get(h,'title'),'string','1991年1月全球温度');

c=colorbar('southoutside','fontsize',12)

c.Label.String = '1994年1月中国温度';

c.Label.FontSize = 15;

(7)多个时间的平均值

clear all

setup_nctoolbox %调用工具包

tic %计时

%%

min_lat=115;

max_lat=279;

min_lon=231;

max_lon=430;  %区域经纬度范围，在数据中的位置

%%

nc=ncgeodataset('tmp2m.gdas.199401.grb2');   %读文件

tem_1=nc.variables   %浏览数据类型

N1=nc.size(tem_1(1));%读取数据大小

d=zeros(1,165,200);%预定义最后的数值存放空间

f=0%验证预留数

%%

for n=1:10  %选取10个时间点

    b1=nc.data(tem_1(1),[n,1,min_lat,min_lon],[n,1,max_lat,max_lon]);%初始的读取位置，最终的位置

    c(n,:,:)=squeeze(b1)-273.16;

    d=d+c(n,:,:); %最终结果

end

for n=1:10

    f=f+c(n,1,1);

end

e=squeeze(d);

% save tem

toc

f值等于e的第一个值，说明计算正确.最后计算平均值即可。　

2、风向

（1）数据读取

clear all

setup_nctoolbox

tic

%% 读取数据文件

wind= ncgeodataset('wnd10m.gdas.199401.grb2');

wind_list = wind.variables;%文件的列表情况

%%

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);

%%

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);

%%

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米

%%

wind_speed=sqrt(u_2.^2+v_2.^2);  %矢量合成

save wind u_2 v_2 wind_speed

toc

（2）展示

clear all

load lon_lat %载入坐标数据

load wind  %载入风速数据

LATLIMS=[3 15];

LONLIMS=[115 134];%选定边界范围

m_proj('lambert','lon',LONLIMS,'lat',LATLIMS);%投影模式

min_lat=115;

max_lat=279;

min_lon=231;

max_lon=430; 

m_coast;

m_grid('box','fancy','tickdir','in');%边缘经纬度宽度

hold on;%图层合并

m_quiver(Plg(min_lat:max_lat,min_lon:max_lon),Plt(min_lat:max_lat,min_lon:max_lon),...

    u_2(min_lat:max_lat,min_lon:max_lon),v_2(min_lat:max_lat,min_lon:max_lon)...

    ,7.5,'r');%经度，纬度，向东，向北的速度分量，转化为矢量箭头

xlabel('global surface winds 1994/01');

（3）

clear all

load lon_lat

load wind

m_proj('hammer-aitoff','clongitude',-150);%投影模式

m_quiver(Plg,Plt,u_2,v_2,15,'r');

shading flat;

colormap('jet');%颜色选择

hold on;

m_quiver(Plg-360,Plt,u_2,v_2,15,'r');

shading flat; %着色模式

colormap('jet');

m_coast();

m_grid('xaxis','middle');

xlabel('global surface winds 1994/01');

　　

3、气压

（1）一个时间点的气压

clear all

clear all

setup_nctoolbox

tic

%% 读取数据文件

pre= ncgeodataset('pressfc.gdas.199401.grb2');

pre_list = pre.variables;%文件的列表情况

%%

size_of_pre = pre.size(pre_list(1));%620小时,576纬度，1152经度

pre_1=pre.data(pre_list(1),[1,1,1],[1,576,1152]); %

pre_2=squeeze(pre_1);

%%

size_of_2 = pre.size(pre_list(2));%v分量的数据尺寸,777小时，1个高度，经纬度数据，4D数据

size_of_3 = pre.size(pre_list(3));

size_of_4 = pre.size(pre_list(4));

size_of_5 = pre.size(pre_list(5));

size_of_8 = pre.size(pre_list(8));

t_4=pre.data(pre_list(4),1,620);

t_8=pre.data(pre_list(8),1:124);

toc

可以看出时间点交错开了，所以最终处理数据时还要按顺序排列起来。

clear all

load lon_lat

load pre

 m_proj('hammer-aitoff','clongitude',-150);%投影模式

m_pcolor(Plg,Plt,pre_2);

shading flat;

colormap('jet');%颜色选择

hold on;

m_pcolor(Plg-360,Plt,pre_2);

shading flat; %着色模式

colormap('jet');

m_coast();

m_grid('xaxis','middle');

% h=colorbar('h');

% set(get(h,'title'),'string','1991年1月全球温度');

c=colorbar('southoutside','fontsize',12)

c.Label.String = '1994年1月全球气压';

c.Label.FontSize = 15;

（2）

clear all

setup_nctoolbox

tic

%% 读取数据文件

pre= ncgeodataset('pressfc.gdas.199401.grb2');

pre_list = pre.variables;%文件的列表情况

size_of_1 = pre.size(pre_list(1));%620小时,576纬度，1152经度

size_of_2 = pre.size(pre_list(2));%v分量的数据尺寸,777小时，1个高度，经纬度数据，4D数据

size_of_3 = pre.size(pre_list(3));

size_of_4 = pre.size(pre_list(4));

size_of_5 = pre.size(pre_list(5));

size_of_8 = pre.size(pre_list(8));

t_4=pre.data(pre_list(4),1,620);

t_8=pre.data(pre_list(8),1:124);

pre_sum=zeros(744,576,1152);

%%

for n=1:744

    if mod(n-1,6)==0  %取余数

        a=floor(n/6);

        n_time=n-a*5;%下标位置

        pre_sum(n,:,:)=pre.data(pre_list(5),[n_time,1,1],[n_time,576,1152]); %

    else

        a=floor((n-1)/6);

        n_time=n-a-1;%下标位置

        pre_sum(n,:,:)=pre.data(pre_list(1),[n_time,1,1],[n_time,576,1152]);

    end

end

%% end

pre_data=squeeze(mean(pre_sum))%求平均值

save pre_data_199401 pre_data

toc

综合运用余数和向下取整的方法，每隔五个从另外的数据中取一个数。最后用mean（）函数求平均值。

4、湿度

　　

　　

　　

　　

　　

　　

　　

　　

　　

　

Matlab——m_map指南（4）——实例的更多相关文章

  1. Matlab——m_map指南（3）——实例

    m_map 实例 1. clear all m_proj('ortho','lat', 48,'long',-123');%投影方式,范围 m_coast('patch','r');%红色填充 m_g ...

  2. Matlab——m_map指南（2）

    3.海岸线和深度测量 3.1.1 海岸线选项 m_coast('line', ...optional line arguments ); m_coast('line', ...optional lin ...

  3. MATLAB——m_map指南（1）

    1.例图 (1) clear all m_proj('oblique mercator');%确定投影方式和绘图界线 m_coast;%画出海岸线 m_grid;%添加格网 第一行代码初始化投影,对于 ...

  4. MATLAB神经网络原理与实例精解视频教程

    教程内容:<MATLAB神经网络原理与实例精解>随书附带源程序.rar9.随机神经网络.rar8.反馈神经网络.rar7.自组织竞争神经网络.rar6.径向基函数网络.rar5.BP神经网 ...

  5. C语言与MATLAB接口 编程与实例 李传军编着

    罗列一下以前自己学习C语言与MATLAB混编的笔记,顺便复习一遍. <C语言与MATLAB接口 编程与实例 李传军编着>(未看完,目前看到P106) 目录P4-8 ************ ...

  6. matlab load函数用法 实例

    一 语法: load(filename) load(filename,variables) load(filename,'-ascii') load(filename,'-mat') load(fil ...

  7. 【翻译】EXTJS 编码风格指南与实例

    原文:EXTJS Code Style Guide with examples Ext JS风格指南: 熟知的且易于学习 快速开发,易于调试,轻松部署 组织良好.可扩展和可维护 Ext JS应用程序的 ...

  8. matlab文件读写处理实例（三）——读取文件特定行

    (1) 读取文件特定行 CODE: ; ;        if nline==line               fprintf(fidout,'%s\n',tline);         data ...

  9. matlab文件读写处理实例（二）——textread批量读取文件

    问题:对文件夹下所有文件进行批量读取,跳过文件头部分,读取每个文件数据部分的7,8,9列,保存到变量并且输出到文件. 数据: 文件夹11m\

随机推荐

  1. 简单说 通过CSS实现 文字渐变色 的两种方式

    说明 这次的重点就在于两个属性, background 属性 mask 属性 这两个属性分别是两种实现方式的关键. 解释 方式一 效果图 代码 <!DOCTYPE html> <ht ...

  2. 前端每日实战：55# 视频演示如何用纯 CSS 创作一个太阳、地球、月亮的运转模型

    效果预览 按下右侧的"点击预览"按钮可以在当前页面预览,点击链接可以全屏预览. https://codepen.io/comehope/pen/RJjQYY 可交互视频 此视频是可 ...

  3. 利用virtualenvwrapper创建虚拟环境出现错误“/usr/bin/python: No module named virtualenvwrapper”

    Linux:CentOS7 python: 系统默认python版本2.7,利用python启动 自己安装python版本3.8,利用python3启动 问题描述: 在上述环境中利用virtualen ...

  4. GitOps入门与实践：如何集成Git和K8S？

    也许你之前听说过GitOps,但是对其并不了解.在本文中,我将对其进行简单介绍,它其实是一个应用程序开发和管理中的一个术语,其核心思想是将应用系统的声明性基础架构和应用程序存放在Git的版本控制库中. ...

  5. angular root在css和less的写法

    背景:angular7项目,angular cli创建了一个模板 新建了一个common.css 本身modle会有一个专属的less文件 写样式 父组件下的子组件样式修改 什么情况下用:root  ...

  6. 第八章、小节二vuex

    a.用vuex首先先安装vuex npm install vuex --save b.在src目录下创建store文件夹,在store中创建index.js存放各个状态 c.在一个模块化的打包系统中, ...

  7. js中所有函数的参数（按值和按引用）都是按值传递的,怎么理解？

    我觉着我可能对这块有点误解,所以单独开个博说下自己的理解,当然是研究后的正解了. 1,参数传递是基本类型,看个例子: function addTen(num){ num += 10; return n ...

  8. 数据结构 4 时间复杂度、B-树 B+树 具体应用与理解

    前言 面试中,经常会问到有关于MYSQL 索引的相关概念,我们之前也都学过有关树的概念.以及二叉树.二叉查找树.红黑树等.这一节,来关注经常是数据库索引中使用的B-树 在说这些之前,我们需要了解时间复 ...

  9. MySQL的万字总结（缓存，索引，Explain，事务，redo日志等）

    hello,小伙伴们,好久不见,MySQL系列停更了差不多两个月了,也有小伙伴问我为啥不更了呢?其实我去看了MySQL的全集,准备憋个大招,更新篇长文(我不会告诉你是因为我懒的). 好了,话不多说,直 ...

  10. 手机抓包app在python中使用

    使用python+airtesr+无线模式控制手机 官方文档中,在airtest.readthedocs.io/zh_CN/lates…有一段介绍如何连接安卓手机的例子: 但是这个线接模板,无线模式的 ...