Changeset 176 for trunk

Timestamp:

Jan 7, 2011, 1:50:08 AM (14 years ago)

Author:

sommeria

Message:

avi2png cleaned and generalised

File:

• trunk/src/series/avi2png.m (modified) (8 diffs)

trunk/src/series/avi2png.m

@@ -2 +2 @@
 %----------------------------------------------------------------------
 function GUI_input=avi2png(num_i1,num_i2,num_j1,num_j2,Series)
+%% INPUT PARAMETERS (to edit)
+increment=4% frame increment: the frequency of the png images will be (initial frequency)/increment.
+colorweight=[0 0.5 0.5]; % relative weight of color components [r g b] for the resulting B/W image
+colorweight=colorweight/sum(colorweight)
 
+%% default output (set the input options in the GUI series)
 if ~exist('num_i1','var')
     GUI_input={};
@@ -8 +13 @@
 end
 hh=guidata(Series.hseries);
-increment=4;
-set(hh.incr_i,'String',num2str(increment))
+set(hh.incr_i,'String',num2str(increment))% preset the increment in the GUI
+
+%% set file names
 nbfield=length(num_i1);
-% basename=[fullfile(Series.RootPath,Series.RootFile) ];
 aviname=fullfile(Series.RootPath,[Series.RootFile Series.FileExt]);
 rootname=fullfile(Series.RootPath,Series.RootFile);
@@ -19 +24 @@
 basename=fullfile(rootname,'frame');
 
-%enable waitbar and stop buttons on the series interface:
+%% enable waitbar and stop buttons on the series interface:
 hRUN=findobj(Series.hseries,'Tag','RUN');% handle of the RUN button
 hwaitbar=findobj(Series.hseries,'Tag','waitbar');%handle of the waitbar
@@ -25 +30 @@
 waitbarpos(3)=Series.WaitbarPos(3);% width of the waitbar
 
+%% read the movie object
 if isequal(lower(Series.FileExt),'.avi')
+    display('opening the avi movie ...')
     hhh=which('mmreader');%look for the existence of 'mmreader'for movie reading
     if ~isequal(hhh,'')&& mmreader.isPlatformSupported()
@@ -34 +41 @@
     end
 end
-%main loop
+
+%% main loop on frames
 for ifile=1:nbfield
      stopstate=get(hRUN,'BusyAction');
@@ -42 +50 @@
         set(hwaitbar,'Position',waitbarpos)%update waitbar on the series interface
         drawnow
-        D=read_image(aviname,'movie',num_i1(ifile),MovieObject);
-        C=uint8(D);% transform to 8 bit integers
+        A=read_image(aviname,'movie',num_i1(ifile),MovieObject);
+        if ndims(A)==3% convert color image to B/W
+            A=double(A);
+            A=colorweight(1)*A(:,:,1)+colorweight(2)*A(:,:,2)+colorweight(3)*A(:,:,3);
+            A=uint8(A);% transform to 8 bit integers
+        end
         new_index=1+floor((num_i1(ifile)-num_i1(1))/increment);
         filename=[basename '_' num2str(new_index) '.png'];%create image name
-        imwrite(C,filename,'BitDepth',8);%write image
+        imwrite(A,filename,'BitDepth',8);%write image
+        display(['new frame '  num2str(new_index) ' written as png image'])
      end
 end
 
-%create xml file with timing: 
+%% create xml file with timing: 
 info=aviinfo(aviname);
 t=xmltree;
@@ -57 +70 @@
 % A AJOUTER
 % Heading.Project='';
-Heading.ImageName=[Series.RootFile '_' num2str(num_i1(1)) '.png'];
+Heading.ImageName='frame_1.png';
 t=struct2xml(Heading,t,uid);
 [t,uid]=add(t,1,'element','Camera');
 Camera.TimeUnit='s';
-Camera.BurstTiming.FrameFrequency=info.FramesPerSecond/increment;
-Camera.BurstTiming.Dti=1;
+% Camera.BurstTiming.FrameFrequency=info.FramesPerSecond/increment;
+Camera.BurstTiming.Dti=increment/info.FramesPerSecond;
 Camera.BurstTiming.NbDti=numel(num_i1)-1;
 Camera.BurstTiming.Time=(num_i1(1)-1)/info.FramesPerSecond;%time of the first frame of the avi movie
@@ -84 +97 @@
         A=imread(filename);
 end
-siz=size(A);
-if length(siz)==3;%color images
-    A=sum(double(A(:,:,2:3)),3)/2;% average green and blue components
-end
+