%多年月平均 无等值线 视频gif图片一起生成 colorbar不均但是是一定范围内的数据取一个色块
clear all
lon=ncread('tos_OISST_L4_AVHRR-only-v2_200601-201012.nc','lon');
lat=ncread('tos_OISST_L4_AVHRR-only-v2_200601-201012.nc','lat');
a=find(lon>=99 & lon<=120);
b=find(lat>=0 & lat<=24);
lon_num=length(a);
lat_num=length(b);
lon_1=lon(a);
lat_1=lat(b);
sst1=ncread('tos_OISST_L4_AVHRR-only-v2_200101-200512.nc','tos',[397 361 1],[lon_num lat_num 60]);
sst2=ncread('tos_OISST_L4_AVHRR-only-v2_200601-201012.nc','tos',[397 361 1],[lon_num lat_num 60]);
sst3=ncread('tos_OISST_L4_AVHRR-only-v2_201101-201412.nc','tos',[397 361 1],[lon_num lat_num 48]);
sst=cat(3,sst1,sst2,sst3);
sst=reshape(sst,[84,96,12,14]);
sst(sst==1.000000020040877e+20)=nan;
mov=moviein(12);
for i=1:12
figure(i);
ssti=mean(sst(:,:,i,:),4);
ssti=squeeze(ssti);
ssti=ssti-273.13;
dlevels=[5,10,15,20,22,24,26,27,27.5,28,28.5,29,30,31,32,34,36,38];
temp=ssti;
for k=1:length(dlevels)-1
ssti(find(temp>dlevels(k) & temp<=dlevels(k+1)))=k;
end
ssti(find(temp==dlevels(1)))=1;
colormap jet;
caxis([0 length(dlevels)-1]) ;
hold on
m_proj('miller','lon',[99,120],'lat',[0,24]);
m_pcolor(lon_1,lat_1,ssti');
hold on
m_coast('patch',[0.6 0.6 0.6],'edgecolor',[0 0 0],'linewidth',0.5);
m_grid('box','fancy','tickdir','out','linestyle','none','fontsize',13);
h=colorbar('fontsize',10,'fontname','Times New Roman');
set(h,'Ticks',[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18],'TickLabels',dlevels) ;
set(get(h,'title'),'string','Temp/℃')
set (gcf,'color','w')
A = ['南海',int2str(i),'月的逐年月平均海表面温度'];
title(A,'fontsize',13)
print(gcf,A,'-dpng','-r1000')
mov(:,i)=getframe(gcf);
if i==1
[I,map]=rgb2ind(mov(:,i).cdata,256);
imwrite(I,map,'南海2001-2014年的逐年月平均海表面温度.gif','DelayTime',.5)
else
imwrite(rgb2ind(mov(:,i).cdata,map),map,'南海2001-2014年的逐年月平均海表面温度.gif','WriteMode','append','DelayTime',.5)
end
end
fprintf('播放电影中...\n');
movie2avi(mov, '南海2001-2014年的逐年月平均海表面温度.avi', 'fps', 1);%输出动画
SST的相关处理.rar_ZZ9_sst_sst的相关处理_sst趋势_取年份
版权申诉
5星 · 超过95%的资源 43 浏览量
2022-07-15
20:47:20
上传
评论
收藏 5KB RAR 举报
小波思基
- 粉丝: 72
- 资源: 1万+
最新资源
- Yolov8改进---注意力机制:Polarized Self-Attention,效果秒杀CBAM、SE.html
- 人才网站设计-asp.net+sql-(系统源码)
- asp.net+sql人才网站设计-含系统源码
- C#应用的用户配置窗体方案
- python实现绘制爱心图形的代码
- JAVAWEB项目-校园订餐系统项目源码.zip
- flink-1.19.0-bin-scala-2.12.tgz flink-1.16.3-bin-scala-2.12.tgz
- javaWeb项目-物资管理系统项目源码.zip
- javaweb项目-物流配货项目源码.zip
- 使用C++基于颜色纹理特征的人脸活体检测实现-附项目源码.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈