Matlab ERA5数据nc格式读取及可视化

最新推荐文章于 2025-03-24 20:02:50 发布

原创最新推荐文章于 2025-03-24 20:02:50 发布 · 3k 阅读

36 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #开发语言

背景：

前两天帮同门画2021年4月19日晚上，渤海湾的风速及温度图。
在ERA5网站上选择需要的数据类型及时间

数据种类：

在这里插入图片描述

Matlab读取nc文件数据并可视化：

clc;
clear all;
close all
%% 读取数据并出图，风速
figure
    set(gcf,'position',[0 40 750 750]);%设置图形窗口位置和大小

path1='渤海湾风温.nc'; % 请自行填补
% path1='C:\Users\小炎\Desktop\try.nc'; % 请自行填补
ncdisp('渤海湾风温.nc')
wind=double(ncread(path1,'si9')); %   是风速
nlon=double(ncread(path1,'longitude')); %读取经度变量
nlat=double(ncread(path1,'latitude')); %读取纬度变量
%ncdisp(path1);%显示nc文件的所有结构，以便大概了解里面的内容
%ncdisp(path1,'siconc');%显示指定变量的内容，注意一定要是变量variables才可以
%ncdisp(path1,'/','min');%简单显示结构以及定义
%ncdisp(path1,'/','full');%全部显示所有结构和定义信息
asadasd = wind(:,:,1)
title({'渤海湾风速示意图',''},'fontsize',14,'fontweight','bold');
[mlat,mlon]=meshgrid(nlat,nlon);

time=19; %提取这个时间点的数据
wind_time = wind(:,:,time);
m_proj('mercator','lon',[116+30/60 124],'lat',[36 41+50/60]);
%     m_gshhs_f('save','mycoast')  % 海岸线精度，我记得有这几个等级i、h、f、c、l
%     m_grid('linewidth',0.1,'tickdir','out','backcolor',[1 1 1]); %边框及网格设定
    m_scale('auto');hold on; 
    hold on
%     m_usercoast('mycoast','patch',[0.74 0.74 0.74],'edgecolor',[0.5 0.5 0.5]);hold on  % 加载海岸线，设陆地颜色 and 海岸线颜色
h = colorbar;
    set(get(h,'label'),'string','风速(m/s)','Fontsize',25); % 给颜色栏命名
%m_pcolor(mlon,mlat,ice_month); # 显示原像素（分辨率较低）
m_contourf(mlon,mlat,wind_time,'linestyle','none'); %等值线 ---类似于插值后的效果
%m_contourf(mlon,mlat,ice_month); %等值线 ---类似于插值后的效果
%colormap(flipud(m_colmap('jet')));%设置颜色legend
colormap('jet');%设置颜色legend
hold on;
shading flat;
    m_gshhs_f('save','mycoast')  % 海岸线精度，我记得有这几个等级i、h、f、c、l

% m_coast('patch',[0.7451 0.7451 0.7451],'edgecolor','none'); % ;  
    m_usercoast('mycoast','patch',[0.74 0.74 0.74],'edgecolor',[0.5 0.5 0.5]);hold on  % 加载海岸线，设陆地颜色 and 海岸线颜色

% 格网格式设置
m_grid('xtick',6,'linestyle',':','tickdir','out','LineWidth',0.5,'fontsize',15);
m_plot(122+13.7/60,37+27.8/60,'*k','markersize',20)
exportgraphics(gcf,['G:\小炎（博）\课题组大论文\渤海湾风速示意图.tif'], 'Resolution', 600)