** ムービーの作成方法を理解したので、コードは将来誰にとっても役立つ場合に備えて、正しいものを反映するように変更されています。このスクリプトは、eqdconic マップのムービーを作成し、avi 形式で保存します。ムービーは 1255 フレームにわたって実行されます。また、画像の特定のポイントにドットをプロットし、映画に変化するタイトルを配置して、何月が実行されているかを示し、右側にカラーバーがあります.
使用される変数の一部は、別の場所で作成されました。それらを作成するためのコードは、コードを要約するために省略されています (そして、それらは私以外の誰にも役に立たないためです)。
% Create movie
nFrames = 34; % Number of frames
for k = 1:nFrames
% Eqdconic script
% Define figure and axes
fg1 = figure(1);
axesm('MapProjection','eqdconic', 'MapParallels', [], 'MapLatLimit',[-80 -59],'MapLonLimit',[190 251]) % 60-70S and 120-160W
framem on; gridm on; mlabel on; plabel on; hold all;
% Plot data
frame = dataPoint_movie(:,:,k);
image = contourfm(lat,lon,frame, 'LineStyle', 'none');
hold on
% Plot dot
plotm(-66.75,224,'k.','MarkerSize',30);
% Colorbar
caxis([0 100]);
h = colorbar;
ylabel(h,'Percent');
% Title: Days 1:1258 inclusive. 20100101 to 20130611
date = datenum(2009, 12, 31) + k; % Convert t into serial numbers
str = datestr(date, 'mmm yyyy'); % Show in the format mmm yyyy so title changes only once a month
title(str);
mov(k) = getframe(gcf); % gca would give only the image. gcf places the title and other attributes on the movie.
end
close(gcf)
% % Save as AVI file
movie2avi(mov, 'SeaIceConcentration.avi', 'compression', 'none', 'fps', 2);