Changeset 1134 for trunk/src/series

Timestamp:

Apr 19, 2024, 8:21:23 PM (6 months ago)

Author:

sommeria

Message:

adapted to reading hdf5 data from fluidimage, extract_rdvision improved

Location:

trunk/src/series

Files:

• extract_rdvision.m (modified) (10 diffs)
• merge_proj.m (modified) (9 diffs)

trunk/src/series/extract_rdvision.m

@@ -66 +66 @@
         ParamOut.WholeIndexRange='on';...% prescribes the file index ranges from min to max (options 'off'/'on', 'off' by default)
         ParamOut.NbSlice='off';%1; ...%nbre of slices, 1 prevents splitting in several processes, ('off' by default)
-        ParamOut.VelType='off';...% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
+    ParamOut.VelType='off';...% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
         ParamOut.FieldName='off';...% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
         ParamOut.FieldTransform = 'off';...%can use a transform function
@@ -72 +72 @@
         ParamOut.Mask='off';...%can use mask option   (option 'off'/'on', 'off' by default)
         ParamOut.CPUTime=0.1;% expected time for writting one image ( in minute)
-        ParamOut.OutputDirExt='.extract';%set the output dir extensionextract_rdvision.m
+    ParamOut.OutputDirExt='.extract';%set the output dir extension
     ParamOut.OutputSubDirMode='one'; %output folder given by the folder name of the first input line
     % detect the set of image folder
@@ -102 +102 @@
     set(hhseries.InputTable,'Data',InputTable)
     ParamOut.ActionInput.LogPath=RootPath;% indicate the path for the output info: 0_LOG ....
+
+    % check the names of .seq and .sqb files
+    iview=1;
+    switch Param.InputTable{iview,5}
+        case {'.seq','.sqb'}
+            filename_seq=fullfile(RootPath,Param.InputTable{iview,2},[Param.InputTable{iview,3} '.seq']);
+            filename_sqb=fullfile(RootPath,Param.InputTable{iview,2},[Param.InputTable{iview,3} '.sqb']);
+            if ~exist(filename_seq,'file')
+                msgbox_uvmat('ERROR',[filename_seq ' missing']);
+            end
+            if ~exist(filename_sqb,'file')
+                msgbox_uvmat('ERROR',[filename_sqb ' missing']);
+            end
+            filexml=[fullfile(RootPath,Param.InputTable{iview,2},Param.InputTable{iview,3}) '.xml'];%xml at the level of the image folder
+            if ~exist(filexml,'file')
+                msgbox_uvmat('ERROR',[filexml ' missing: needed to get the image organisation and timing ']);
+                return
+            end
+            [XmlData,errormsg]=imadoc2struct(filexml);
+            if ~isempty(errormsg)
+                msgbox_uvmat('ERROR',errormsg);
+                return
+            end
+            timexml=reshape(XmlData.Time(2:end,2:end)',1,[]);
+        otherwise
+            msgbox_uvmat('ERROR','bad input file : select .seq or .sqb for image extraction');
+            return
+    end
+    % get data from .seq file
+    s=ini2struct(filename_seq);
+    SeqData=s.sequenceSettings;
+    SeqData.nb_frames=str2double(s.sequenceSettings.numberoffiles);
+% reading the .sqb file
+    m = memmapfile(filename_sqb,'Format', { 'uint32' [1 1] 'offset'; ...
+        'uint32' [1 1] 'garbage1';...
+        'double' [1 1] 'timestamp';...
+        'uint32' [1 1] 'file_idx';...
+        'uint32' [1 1] 'garbage2' },'Repeat',SeqData.nb_frames);
+    
+        timestamp=zeros(1,numel(m.Data));
+        for ii=1: numel(m.Data)
+            timestamp(ii)=m.Data(ii).timestamp;
+        end
+        if numel(timestamp)<= numel(timexml)
+            timexml=timexml(1:numel(timestamp));
+        else
+            msgbox_uvmat('ERROR','time sequence defined by the xml file too small')
+            return
+        end
+        difftime=timestamp-timexml;
+        if max(difftime)>0.01
+        figure
+        plot(timestamp,difftime)
+        xlabel('timestamps(s)')
+        ylabel('time difference(s)')
+        title('discrepency timestamps-timexml') 
+        end
     return
 end
 
+
+%%%%%%%%%%%% STANDARD PART  %%%%%%%%%%%%
 ParamOut=[];
-%%%%%%%%%%%% STANDARD PART  %%%%%%%%%%%%
 %% read input parameters from an xml file if input is a file name (batch mode)
 
@@ -128 +186 @@
 % input file or frame indices i1_series,i2_series,j1_series,j2_series
 [filecell,i1_series,i2_series,j1_series,j2_series]=get_file_series(Param);
-
-%OutputDir=[Param.OutputSubDir Param.OutputDirExt];
  
 % numbers of slices and file indices
@@ -137 +193 @@
 
 %determine the file type on each line from the first input file 
-
 FileInfo=get_file_info(filecell{1,1});
 if strcmp(FileInfo.FileType,'rdvision')
-%     if ~isequal(FileInfo.NumberOfFrames,nbfield)
-%         disp_uvmat('WARNING',['the whole series of ' num2str(FileInfo.NumberOfextract_rdvision.mFrames) ' images must be extracted at once'],checkrun)
-%         %rmfield(OutputDir)
-% %         return
-%     end
+
     %% interactive input of specific parameters (for RDvision system)
     display('converting images from RDvision system...')
@@ -166 +217 @@
 %%%  loop on the cameras ( #iview)
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% RootPath=Param.InputTable(:,1);
-% RootFile=Param.InputTable(:,3);
-% SubDir=Param.InputTable(:,2);
-% NomType=Param.InputTable(:,4);
-% FileExt=Param.InputTable(:,5);
-
-% [XmlData,NbSlice_calib,time,errormsg]=read_multimadoc(RootPath,SubDir,RootFile,FileExt,i1_series,i2_series,j1_series,j2_series);
-% if size(time,1)>1
-%     diff_time=max(max(diff(time)));
-%     if diff_time>0
-%         disp_uvmat('WARNING',['times of series differ by (max) ' num2str(diff_time)],checkrun)
-%     end
-% end
-%
-%      nbfield2=size(time,1);
-checkpreserve=0;% if =1, will npreserve the original images, else it erases them at the end
 
 for iview=1:size(Param.InputTable,1)
@@ -196 +231 @@
     newxml=[newxml '.xml'];
     
-    %% get the names of .seq and .sqb files
+    %% get the names of .seq and .sqb files and save them to the log output folder .extract
     switch Param.InputTable{iview,5}
         case {'.seq','.sqb'}
@@ -202 +237 @@
             filename_sqb=fullfile(RootPath,Param.InputTable{iview,2},[Param.InputTable{iview,3} '.sqb']);
             errormsg='';
-            if isequal(Param.IndexRange.first_i,1)
-                
+            % backup of the seq, sqb and xml files for the first frame index
+            if isequal(Param.IndexRange.first_i,1)% backup of the seq, sqb and xml files for the first frame index
                 logdir=[Param.OutputSubDir Param.OutputDirExt];
                 [success,errormsg] = copyfile(filename_seq,[fullfile(RootPath,logdir,Param.InputTable{iview,3}) '.seq']); %copy the seq file in the upper folder
@@ -215 +250 @@
                 if check_xml
                     [success,errormsg] = copyfile(filexml,[fullfile(RootPath,logdir,Param.InputTable{iview,3}) '.xml']); %copy the original xml file in the upper folder
-                if ~success
-                    disp(errormsg)
-                end
+                    if ~success
+                        disp(errormsg)
+                    end
                 else
                     disp(['error:' filexml ' missing']);
@@ -310 +345 @@
         return
     end
-    
-    % check the existence of the expected output image files (from the xml)
-    
-    FileDir=SeqData.sequencename;
-     FileDir=regexprep(FileDir,'_Master_Dalsa_4M180$','');%suppress '_Master_Dalsa_4M180'
-%     for i1=1:numel(timestamp)/nbfield2
-%         for j1=1:nbfield2
-%             OutputFile=fullfile_uvmat(RootPath,FileDir,'img','.png',NomTypeNew,i1,[],j1);% TODO: set NomTypeNew from SeqData.mode
-%             try 
-%             A=imread(OutputFile);% check image reading (stop if error)
-%             catch ME
-%                 disp(['checking ' OutputFile])
-%                 disp(ME.message)
-%             end
-%         end
-%     end
+
 end
 

trunk/src/series/merge_proj.m

@@ -67 +67 @@
     ParamOut.FieldName='one';% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
     ParamOut.FieldTransform = 'on';%can use a transform function
-    ParamOut.TransformPath=fullfile(fileparts(which('uvmat')),'transform_field');% path to transform functions (needed for compilation only)
+    %%%%% list of possible transform functions (needed only for compilation)
+        ListTransform={'phys','phys_polar','sub_field'};%list of possible transform functions (needed only for compilation)
+        if 0==1 %never satisfied but trigger compilation with the appropriate transform functions
+            for ilist=1:numel(ListTransform)
+                eval(ListTransform)
+            end
+        end
+    ParamOut.TransformPath=fullfile(fileparts(which('uvmat')),'transform_field');% path to transform functions 
+    %%%%%%%%
     ParamOut.ProjObject='on';%can use projection object(option 'off'/'on',
     ParamOut.Mask='on';%can use mask option   (option 'off'/'on', 'off' by default)
@@ -86 +94 @@
     return
 end
-if 0==1
-    phys; % used to include phys when compiling is done
-end
+
 %%%%%%%%%%%% STANDARD PART (DO NOT EDIT) %%%%%%%%%%%%
 ParamOut=[]; %default output
@@ -115 +121 @@
 RootFile=Param.InputTable(:,3);
 SubDir=Param.InputTable(:,2);
-% NomType=Param.InputTable(:,4);
 FileExt=Param.InputTable(:,5);
 
@@ -126 +131 @@
 % i1_series(iview,ref_j,ref_i)... are the corresponding arrays of indices i1,i2,j1,j2, depending on the input line iview and the two reference indices ref_i,ref_j 
 % i1_series(iview,fileindex) expresses the same indices as a 1D array in file indices
-if ~isempty(hdisp),delete(hdisp),end;%end the waiting display
+if ~isempty(hdisp),delete(hdisp),end %end the waiting display
 
 NbView=numel(i1_series);%number of input file series (lines in InputTable)
@@ -134 +139 @@
 
 %% determine the file type on each line from the first input file 
-NcTypeOptions={'netcdf','civx','civdata'};
+NcTypeOptions={'netcdf','civx','civdata','pivdata_fluidimage'};
 for iview=1:NbView
     if ~exist(filecell{iview,1}','file')
@@ -178 +183 @@
 %% coordinate transform or other user defined transform
 transform_fct='';%default fct handle
+ checksub=0;
 if isfield(Param,'FieldTransform')&&~isempty(Param.FieldTransform.TransformName)
     currentdir=pwd;
@@ -187 +193 @@
             XmlData{iview}.TransformInput=Param.TransformInput;
         end
-    end       
+    end 
+    checksub=nargin(transform_fct);% number of input arguments for the selected transform fct
+    if checksub>2 && NbView>2
+        disp_uvmat('WARNING',['only the two first input file series will be combined by ' Param.FieldTransform.TransformName],checkrun)
+    end
 end
 %%%%%%%%%%%% END STANDARD PART  %%%%%%%%%%%%
@@ -305 +315 @@
         end
         
-        %% transform the input field (e.g; phys) if requested (no transform involving two input fields)
+        %% transform the input field iview (e.g; phys) if requested (no transform involving two input fields at this stage)
+        checksub=0;
         if ~isempty(transform_fct)
-            if nargin(transform_fct)>=2
+            checksub=nargin(transform_fct);
+            if checksub==2
                 Data{iview}=transform_fct(Data{iview},XmlData{iview});
-            else
+            elseif checksub==1
                 Data{iview}=transform_fct(Data{iview});
             end
@@ -338 +350 @@
     %%%%%%%%%%%%%%%% END LOOP ON VIEWS %%%%%%%%%%%%%%%%
 
-    %% merge the NbView fields
-    [MergeData,errormsg]=merge_field(Data);
+    %% merge the NbView fields 
+    if checksub<=2
+        [MergeData,errormsg]=merge_field(Data);%concatene all the input field series by fct merge_data
+    elseif checksub==3
+        MergeData=transform_fct(Data{1},XmlData{1},Data{2}); %combine the two input file series
+    else
+        MergeData=transform_fct(Data{1},XmlData{1},Data{2},XmlData{2});%combine the two input file series with calibration parameters
+    end
     if ~isempty(errormsg)
         disp_uvmat('ERROR',errormsg,checkrun);