MATLAB 如何绘制GIF动图，以全球SST展示厄尔尼诺现象为例

Hua_Zhang

matlab绘制GIF主要分为俩步：
1、批量导出图片
SST（longitude,latitude,time）就要批量导出绘制时间序列的月平均SST（本人在1979-2017年数据中提取07-16年的十年数据绘制）

2、GIF动图生成
由批量导出图制作GIF(本人提取了97-98年的月平均数据，就是合成厄尔尼诺年的SST动图)

动图不知道为啥上传不了，URL也不让放，SST数据太大了，你懂的。。。
所以自己下载吧，只为赚点积分和给大家提供一点帮助而已


代码

1. %% load data
   
2. clc; close all; clear ;
   
3. input_path = 'D:\file_data\课件_秋季\气候统计方法和应用\气候统计方法与应用_homework\气候统计3\';
   
4. out_path = 'D:\file_data\课件_秋季\气候统计方法和应用\气候统计方法与应用_homework\气候统计3\out_fig\';
   
5. lon = ncread([input_path,'1979-2017.nc'],'longitude');
   
6. lat = ncread([input_path,'1979-2017.nc'],'latitude');
   
7. time = ncread([input_path,'1979-2017.nc'],'time');
   
8. sst = ncread([input_path,'1979-2017.nc'],'sst');
   
9. 
   
10. %% data process
    
11. mon_s = 1;mon_e = 456;
    
12. SST_07_16 = sst(:,:,mon_s:mon_e);
    
13. SST_07_16(SST_07_16 > 35) = nan;
    
14. SST_07_16(SST_07_16 < -4) = nan;
    
15. nt = length(time(mon_s:mon_e));
    
16. yr = nt/12;
    
17. [sst_ano] = ano(SST_07_16,yr);
    
18. [xc,yc] = meshgrid(lon,lat);
    
19. 
    
20. for mon = 1:nt
    
21.     %yr = (mon -1) / 12 + 1;
    
22.     figure
    
23.     contourf(xc,yc,sst_ano(:,:,mon)');colorbar;caxis([-1.5 1.5]);
    
24.     colormap;
    
25.     grid on;
    
26.     xlabel('longitude(E)','FontSize',12);ylabel('latitude(S)','FontSize',12) ;
    
27.     set(gca,'fontsize',12);
    
28.     title(['SST ano TIME=',num2str(mon),' mon '],'fontsize',12);
    
29.     %gtext({'(Pa)'},'fontsize',12);
    
30.     set(gcf,'position',[200,300,800,600]);
    
31.     saveas(gcf,[out_path,num2str(mon),'mon_ano.jpg']);
    
32.     close
    
33. end
    
34. 
    
35. disp('okk');
    
36. 
    
37. %% clc ; clear ;close all;
    
38. filename = 'SST_eino_1.gif';
    
39. key_s = 215;
    
40. key_e = 240;
    
41. for idx = key_s:key_e%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
42.    str = strcat(num2str(idx),'mon_ano.jpg');
    
43.    K = imread(str);
    
44.    [A,map] = rgb2ind(K,256);
    
45. if idx == key_s%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
46.     imwrite(A,map,filename,'gif','LoopCount',inf,'DelayTime',0.1);
    
47. else
    
48.     imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0.1);
    
49. end
    
50. end
    
51. 
    
52. disp('finish!!');
    
53. 
    
54. 
    
55. function [A_anomaly] = ano(A,n)
    
56. %% Description
    
57. % ano is the code to get a muti-year monthly average data that removes monthly averages
    
58. % separately.
    
59. % A is the original data . It is designed specifically for 3D matricies of data such as  
    
60. % sea surface temperatures where dimensions 1 and 2 are spatial dimensions (e.g., lat
    
61. % and lon; lon and lat; x and y, etc.), and the third dimension represents different
    
62. % slices or snapshots of data in time.  
    
63. % n is the times of years about original data.
    
64. % A_anomaly is a muti-year monthly average data that removes monthly averages separately.
    
65. 
    
66. %% Author Info
    
67. % huazhang
    
68. % 27/11/2018
    
69. 
    
70. for i1=1:12
    
71.     a=A(:,:,i1);
    
72.    for j1=1:n
    
73.     i2=i1+12*(j1-1);
    
74.     a=a+A(:,:,i2);
    
75.    end
    
76.    a_all(:,:,i1)=a./j1;
    
77. end
    
78. 
    
79. for i1=1:12
    
80.    for j1=1:n
    
81.     i2=i1+12*(j1-1);
    
82.     a=A(:,:,i2)-a_all(:,:,i1);
    
83.     A_anomaly(:,:,i2)=a;
    
84.    end
    
85. end
    
86. end
    

复制代码