0

** ムービーの作成方法を理解したので、コードは将来誰にとっても役立つ場合に備えて、正しいものを反映するように変更されています。このスクリプトは、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); 
4

2 に答える 2