下载地址
https://marine.copernicus.eu
数据在这里
链接:https://pan.baidu.com/s/1RR_i6LtJqutaQrXzXV0tzg 密码:536g
clear all; close all; clc
%% 读取数据
file='northwest_pacific.nc';
sla=ncread(file,'sla'); u=ncread(file,'ugos'); v=ncread(file,'vgos');
lon=ncread(file,'longitude'); lat=ncread(file,'latitude');
time=ncread(file,'time'); [lat,lon]=meshgrid(lat,lon);
%将time数组中的数字转化为字符串日期
date=datestr(double(time)+datenum('19500101','yyyymmdd'),'yyyymmdd');
k=1;
%% 画中尺度涡运动图
for i=1:length(time)
figure(1)
set(gcf,'visible','off')
set(gcf,'units','centimeters','position',[0 0 20 15])
set(gca,'position',[0.05 0.05 0.95 0.9])
m_proj('Mercator','lon',[110 160],'lat',[5 40]);
m_pcolor(lon,lat,squeeze(sla(:,:,i))); shading interp;%画海表高度异常场
colormap([m_colmap('diverging',256)]); caxis([-0.4 0.4]);
colorbar('fontsize',10);
m_gshhs_l('patch',[.5 .8 .0],'edgecolor','none');
m_grid('box','none','xtick',10,'ytick',7,'fontsize',10);
hold on
h=m_quiver(lon,lat,squeeze(u(:,:,i)),squeeze(v(:,:,i)),2,'color','k',...
'maxheadsize',2,'linewidth',1);%画地转流场
m_text(117,28,date(i,5:8),'color','k','fontsize',30,'FontWeight','bold')%标上日期
frame=getframe(gcf);
im=frame2im(frame);%制作gif文件,图像必须是index索引图像
[I,map]=rgb2ind(im,256);
if k==1
imwrite(I,map,'eddy_motion_new.gif','gif','Loopcount',inf,...
'DelayTime',0.2);%loopcount只是在i==1的时候才有用
else
imwrite(I,map,'eddy_motion_new.gif','gif','WriteMode','append',...
'DelayTime',0.2);%DelayTime用于设置gif文件的播放快慢
end
clf;k=k+1;
end