Changeset 609 for trunk/src/series/civ_series.m

Timestamp:

Apr 9, 2013, 8:20:00 PM (11 years ago)

Author:

sommeria

Message:

various bugs corrected after tests with Windows OS.

File:

• trunk/src/series/civ_series.m (modified) (5 diffs)

trunk/src/series/civ_series.m

@@ -58 +58 @@
     if ~exist(filecell{1,1},'file')
         msgbox_uvmat('WARNING','the first input file does not exist')
+    else
+        FileType=get_file_type(filecell{1,1});
+        if isempty(find(strcmp(FileType,{'civdata','image','multimage','mmreader','video'})));% =1 for images
+            msgbox_uvmat('ERROR',['bad input file type for ' mfilename ': an image or civdata file is needed'])
+        end
     end
     return
@@ -85 +90 @@
     MinIndex=cell2mat(Param.IndexRange.MinIndex);
     [filecell,i_series,tild,j_series]=get_file_series(Param);
+    [FileType_A,FileInfo,MovieObject_A]=get_file_type(filecell{1,1});
+    if strcmp(FileType_A,'civdata')% a civdata file has been introduced as input.
+         [FileType_A,FileInfo,MovieObject_A]=get_file_type(filecell{2,1});
+    end
     [i1_series_Civ1,i2_series_Civ1,j1_series_Civ1,j2_series_Civ1,check_bounds,NomTypeNc]=...
         find_pair_indices(PairCiv1,i_series{1},j_series{1},MinIndex,MaxIndex);
@@ -118 +127 @@
     
     NbField=numel(i1_series_Civ1);
-    [FileType_A,FileInfo,MovieObject_A]=get_file_type(filecell{1,1});
+
     FileType_B=FileType_A;
     MovieObject_B=MovieObject_A;
@@ -143 +152 @@
 
 %% get timing from the ImaDoc file or input video
-[XmlData,NbSlice_calib,time,errormsg]=read_multimadoc(RootPath,SubDir,RootFile,FileExt,i1_series,i2_series,j1_series,j2_series);
+XmlFileName=find_imadoc(RootPath,SubDir,RootFile,FileExt);
+time=[];
+if ~isempty(XmlFileName)
+    XmlData=imadoc2struct(XmlFileName);
+        if isfield(XmlData,'Time')
+            time=XmlData.Time;
+            TimeSource='xml';
+        end
+        if isfield(XmlData,'Camera')
+            if isfield(XmlData.Camera,'NbSlice')&& ~isempty(XmlData.Camera.NbSlice)
+                NbSlice_calib{iview}=XmlData.Camera.NbSlice;% Nbre of slices for Zindex in phys transform
+                if ~isequal(NbSlice_calib{iview},NbSlice_calib{1})
+                    msgbox_uvmat('WARNING','inconsistent number of Z indices for the two field series');
+                end
+            end
+            if isfield(XmlData.Camera,'TimeUnit')&& ~isempty(XmlData.Camera.TimeUnit)
+                TimeUnit=XmlData.Camera.TimeUnit;
+            end
+        end
+end
+if strcmp(InputTable{iview,4},'*')
+    if ~isempty(VideoObject)
+        imainfo=get(VideoObject);
+        time=zeros(imainfo.NumberOfFrames+1,2);
+        time(:,2)=(0:1/imainfo.FrameRate:(imainfo.NumberOfFrames)/imainfo.FrameRate)';
+        TimeSource='video';
+       % set(han:dles.Dt_txt,'String',['Dt=' num2str(1000/imainfo.FrameRate) 'ms']);%display the elementary time interval in millisec
+        ColorType='truecolor';
+    elseif ~isempty(imformats(regexprep(InputTable{iview,5},'^.',''))) || isequal(InputTable{iview,5},'.vol')%&& isequal(NomType,'*')% multi-frame image
+        if ~isempty(InputTable{iview,2})
+            imainfo=imfinfo(fullfile(InputTable{iview,1},InputTable{iview,2},[InputTable{iview,3} InputTable{iview,5}]));
+        else
+            imainfo=imfinfo([FileBase InputTable{iview,5}]);
+        end
+        ColorType=imainfo.ColorType;%='truecolor' for color images
+        if length(imainfo) >1 %case of image with multiple frames
+            nbfield=length(imainfo);
+            nbfield_j=1;
+        end
+    end
+end
+
 %TODO: get time_A and time_B
-% case of movies TODO TODO TODO
+
 if isempty(time) && (strcmp(FileType,'video') || strcmp(FileType,'mmreader'))
     set(handles.ListPairMode,'Value',1);
@@ -161 +211 @@
     end
     TimeUnit='s';
+end
 %%%%% MAIN LOOP %%%%%%