Imagesc movie with multiple axes/matrices

1 次查看（过去 30 天）
Evan Kias
Evan Kias 2013-1-4
I would like to do the same as this person:
http://www.mathworks.com/matlabcentral/answers/17727-plotting-multiple-matrices-in-the-same-plot-so-i-can-seee-how-the-matrices-develop-trough-time
Although, I want multiple axes/matrices in the same movie. I wrote the following code to achieve this:
for i = 1:l
% Switch to first axes
axes(syy_ax)
imagesc(syy(:,:,i))
colorbar
% Switch to second axes
axes(wd_ax)
imagesc(wd(:,:,i))
colorbar
all_mov(i) = getframe;
end
where syy_ax and wd_ax are axis handles, and syy and wd are 3D matrices. I realize how slow this is, but don't know an alternative. I tried specifying the axis in the plot function, but I don't think imagesc supports this. Could anyone suggest another way, please?
  1 个评论
Jurgen
Jurgen 2013-1-4
Have you considered concatenating the matrices into 1 matrix if possible?
Also the post you refer to is trying to display something over time, while your use of getframe implies you want to make a movie object?

请先登录，再进行评论。

请先登录，再回答此问题。

回答（3 个）

Image Analyst
Image Analyst 2013-1-4
编辑：Image Analyst 2013-1-4
Well how slow is it? I would think it would be too fast, if anything, and you'd want to put in pauses so the person viewing it will see more than a blip as the movie goes whizzing by. Just how many frames per second is it doing? Or is it so slow we're talking seconds per frame?
Also I'm puzzled by you using imagesc() rather than imshow()? So, you're okay with that bizarre, arbitrary, default colormap that imagesc() applies? Just by the luck of the draw, that happens to be the best colormap for your images? I'm surprised - you're lucky. For me, almost always I don't want the default colormap that imagesc chooses.
Also, please see the FAQ: http://matlab.wikia.com/wiki/FAQ#How_can_I_create_a_movie_from_my_MATLAB_figures.3F You can use export_fig() to save your figure, with both axes on it, out to files that can then be made into frames of your movie.
Evan Kias
Evan Kias 2013-1-4
Its slow, just under two frames per second. Here's a copy of the complete code.
fig1 = figure('color','white');
maximize(fig1)
syy_ax = axes;
set(syy_ax,'Position',[0.05 0.05 0.4 0.95])
set(syy_ax,'XTick',[],'YTick',[],'XColor','w')
hold(syy_ax)
%
wd_ax = axes;
set(wd_ax,'Position',[0.55 0.05 0.4 0.95])
set(wd_ax,'XTick',[],'YTick',[],'XColor','w')
hold(wd_ax)
max_syy = max_val(syy);
max_wd = max_val(wd);
syylims = [0 ceil(max_syy)];
wdlims = [0 ceil(max_wd)];
s = size(syy);
l = s(3);
for i = 1:l
axes(syy_ax)
imagesc(syy(:,:,i),syylims)
colorbar
%
axes(wd_ax)
imagesc(wd(:,:,i),wdlims)
colorbar
all_mov(i) = getframe;
end
maximize and max_val are custom functions that only run once outside the loop. My matrices (syy and wd) are both size = [20 10 3200]. I am eventually going to need larger matrices and 3-5 more axes in the figure. I just thought because of how slow it is processing I am doing something terribly wrong.
In regards to the unattractive coloring, I don't mind so much right now. The resulting images are color contours of stress and work values. I plan on using interp2 to smooth it out and can change the colormap then if necessary. Any thoughts? Thanks so much!
  8 个评论
显示 6更早的评论
Evan Kias
Evan Kias 2013-1-10
I had some success using the following code. It implements, imshow, the 'CData' property and uses export_fig instead of getframe. I altered my strategy for now to include a plot instead of another imshow on the second axis. I think the plot command within the loop is a slow operation (although not slower than export_fig). Is there an alternative that is similar to the 'CData' property? Also, would you please recommend a freeware that can string together my png's into a movie? Apparently video encoding is a huge area, and I hope to stay out of that one! Thanks!
fig1 = figure('color','white');
maximize(fig1)
pix_array = syy_pix_half;
plot_array = -syy_mcs_half;
img_ax = axes;
set(img_ax,'Position',[0.45 0.15 0.6 0.7])
set(img_ax,'XTick',[],'YTick',[],'XColor','w')
set(img_ax,'DataAspectRatio',[1 1 1])
hold(img_ax)
%
pl_ax = axes;
set(pl_ax,'Position',[0.10 0.15 0.4 0.75])
set(pl_ax,'XLim',[0 12000],'YLim',[0 8e6])
xlabel('Steps')
ylabel('Stress (Pa)')
hold(pl_ax)
max_arr = max_val(pix_array);
arrlims = [0 ceil(max_arr)];
s = size(pix_array);
l = s(3);
h1 = imshow(pix_array(:,:,1), arrlims, 'Parent', img_ax);
colorbar('peer',img_ax,'SouthOutside');
for i = 6000:7000
set(h1, 'CData', pix_array(:,:,i));
plot(plot_array(1:i),'Parent',pl_ax);
drawnow
filename = ['img_' num2str(i)];
export_fig(filename)
end
Jurgen
Jurgen 2013-1-11
You can do it in MATLAB with avifile or videowriter, getframe should be faster than export_fig if you intialize a struct array.
That aside, freeware options that come to mind are imagemagick & ffmpeg (both commandline interfaces I believe).
If you insist on saving images, saveas() might be faster than exportfig.
Lastly, either when dumping images or grabbing frames, you may suppress the visible output to speed things up.

请先登录，再进行评论。

Jurgen
Jurgen 2013-1-5
I dont use colorbars often, so this was my initial idea without them:
syy = reshape(syy,[size(syy,1) size(syy,2) 1 size(syy,3)]);
wd = reshape(wd,[size(wd,1) size(wd,2) 1 size(wd,3)]);
DATA = [syy wd]; clear syy wd %data is now formatted (n-m-1-p) and concatenated
MOV = immovie(DATA,jet(256)); %or some other map (< 256).
close;
Note you may need to convert DATA to indexed format (i.e. integers). A quick fix for this would be real2rgb (from the FEX) then rgb2ind.
fps = 5; %framerate
N = 1; %iterations
movie(MOV,N,fps) %Plays movie MOV, N times at fps frames per second.

类别

在 Help Center 和 File Exchange 中查找有关 Animation 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by