MATLAB绘制台风（Tropical cyclones)风场(Wind field)图、风场动画(Animati - ToB企服应用市场:ToB评测及商务社交产业平台

% 读取ERA5数据
% ncinfo("ERA5.nc") 用于了解NC文件结构
lon = ncread("ERA5.nc",'longitude'); %longitude
lat = ncread("ERA5.nc",'latitude'); %latitude
date = ncread("ERA5.nc",'time'); %时间
date = datetime(1900,1,1) + double(date(:))/24; %时间格式转化
u = ncread("ERA5.nc",'u10'); %速度u分量
v = ncread("ERA5.nc",'v10'); %速度v分量
speed = sqrt(u.^2+v.^2); %速度

复制代码

% 生成m_map project
m_proj('Equidistant Cylindrical','long',[110 160],'lat',[5 45]);
% 计算坐标对应地图位置
X=[];
Y=[];
for i=1:length(lon)
for j=1:length(lat)
[CX,CY]=m_ll2xy(lon(i),lat(j));
X(i,j)=CX;
Y(i,j)=CY;
end
end
for i=1:6:length(date)
figure;
% 绘制风场
hold on
contourf(X,Y,(speed(:,:,i)),15,'linestyle','none');
% 海岸线
m_coast('patch',[.8 .8 .8]);
% 绘制边框
m_grid('box','on','tickdir','in','backcolor',[1,1,1],'xtick',[],'ytick',[],'LineWidth',1.5);
hold off
% 生成colorbar
cb = colorbar;
cb.FontSize=14;
cb.Label.String='Wind speed (m/s)';
% 颜色设置
cr1 = flipud(slanCM(167));%彩色
cr2 = slanCM(25);%灰色
colormap([cr2(1:3:83,:);flipud(cr1(1:2:end-3,:))]);% 将低风速设置为灰色，突出高风速台风区域
% 日期
XLB = xlabel(string(date(i)),'FontSize',15,'FontWeight','bold');
set(XLB,'Units','normalized','Position',[0.5,0],'HorizontalAlignment','center','VerticalAlignment','top')
frame = getframe(gcf);
im = frame2im(frame);
[I,map] = rgb2ind(im,256);
if i == 1
imwrite(I,map,'windfield.gif','gif', 'Loopcount',inf,'DelayTime',1);
else
imwrite(I,map,'windfield.gif','gif','WriteMode','append','DelayTime',1);
end
close
end
winopen('windfield.gif')% 预览结果

复制代码

%% 提取1:12日每天0点的数据
index = hour(date) == 0& day(date) ~= 13;
date = date(index);
u = u(:,:,index);
v = v(:,:,index);
speed = speed(:,:,index);
%% 绘制台风风场图
% 分子图每幅图一个时刻
t = tiledlayout(3,4);
% 生成m_map project
m_proj('Equidistant Cylindrical','long',[110 160],'lat',[5 45]);
for i = 1:12
% 生成子图
nexttile(i);
hold on
% 绘制等高线图
contourf(X,Y,speed(:,:,i),15,'LineStyle','none')
caxis([0,32])% 颜色上下限
% 海岸线图
m_coast('patch',[.8 .8 .8]);
% m_gshhs_i('patch',[.8 .8 .8]);
% 绘制边框
m_grid('box','on','tickdir','in','backcolor',[1,1,1],'xtick',[],'ytick',[],'LineWidth',1.5);
hold off
% xlabel
XLB=xlabel(string(date(i)),'FontSize',15,'FontWeight','bold');
set(XLB,'Units','normalized','Position',[0.5,0],'HorizontalAlignment','center','VerticalAlignment','top')
end
% 设置间隙
t.Padding='none';
t.TileSpacing="none";
% 生成colorbar,实用Tiledlayout的设置
cb = colorbar;
cb.FontSize=14;
cb.Layout.Tile='east';
cb.Label.String='Wind speed (m/s)';
% 颜色设置
cr1 = flipud(slanCM(167));%彩色
cr2 = slanCM(25);%灰色
colormap([cr2(1:3:83,:);flipud(cr1(1:2:end-3,:))]);% 将低风速设置为灰色，突出高风速台风区域
% 输出图片设置
figWidth = 30; % 设置图片宽度
figHeight = 20; % 设置图片高度
set(gcf,'PaperUnits','centimeters'); % 图片尺寸所用单位
set(gcf,'PaperPosition',[0 0 figWidth figHeight]);
fileout = 'Windfield_Chan_hom'; % 输出图片的文件名
print(gcf,[fileout,'.tif'],'-r300','-dtiff'); % 设置图片格式、分辨率
winopen([fileout,'.tif'])

复制代码