Changeset 1197 for trunk/src/series

Timestamp:

Mar 11, 2026, 5:31:25 PM (4 weeks ago)

Author:

sommeria

Message:

bugs repaired

Location:

trunk/src/series

Files:

• civ_series.m (modified) (1 diff)
• merge_proj.m (modified) (1 diff)
• turb_correlation_x.m (modified) (9 diffs)

trunk/src/series/civ_series.m

@@ -276 +276 @@
     time_patch2=0;
     if ~isempty(RUNHandle)% update the waitbar in interactive mode with GUI series  (checkrun=1)
-        update_waitbar(WaitbarHandle,ifield/NbField)
         if  checkrun && ~strcmp(get(RUNHandle,'BusyAction'),'queue')
             disp('program stopped by user')

trunk/src/series/merge_proj.m

@@ -354 +354 @@
             end
             [MaskData,~,errormsg] = read_field(maskname,'image');
+            if ~isempty(errormsg)
+                disp(['error reading ' maskname ' :' errormsg])
+                return
+            end
             if ~isempty(NbSlice_calib)
                 MaskData.ZIndex=mod(i1_series{iview}(index)-1,NbSlice_calib{iview})+1;%Zindex for phys transform

trunk/src/series/turb_correlation_x.m

@@ -94 +94 @@
     
 %% root input file(s) name, type and index series
-RootPath=Param.InputTable(:,1);
-RootFile=Param.InputTable(:,3);
-SubDir=Param.InputTable(:,2);
-NomType=Param.InputTable(:,4);
-FileExt=Param.InputTable(:,5);
+RootPath=Param.InputTable{:,1};
+RootFile=Param.InputTable{:,3};
+%SubDir=Param.InputTable(:,2);
+NomType=Param.InputTable{1,4};
+%FileExt=Param.InputTable(:,5);
 [filecell,i1_series,i2_series,j1_series,j2_series]=get_file_series(Param);
 %%%%%%%%%%%%
@@ -108 +108 @@
 % i1_series(iview,fileindex) expresses the same indices as a 1D array in file indices
 %%%%%%%%%%%% NbView=1 : a single input series
-NbView=numel(i1_series);%number of input file series (lines in InputTable)
+%NbView=numel(i1_series);%number of input file series (lines in InputTable)
 NbField_j=size(i1_series{1},1); %nb of fields for the j index (bursts or volume slices)
 NbField_i=size(i1_series{1},2); %nb of fields for the i index
 NbField=NbField_j*NbField_i; %total number of fields
 
-%% determine the file type on each line from the first input file 
-ImageTypeOptions={'image','multimage','mmreader','video','cine_phantom'};
-NcTypeOptions={'netcdf','civx','civdata'};
-for iview=1:NbView
-    if ~exist(filecell{iview,1}','file')
-        msgbox_uvmat('ERROR',['the first input file ' filecell{iview,1} ' does not exist'])
-        return
-    end
-    [FileInfo{iview},MovieObject{iview}]=get_file_info(filecell{iview,1});
-    FileType{iview}=FileInfo{iview}.FileType;
-    CheckImage{iview}=~isempty(find(strcmp(FileType{iview},ImageTypeOptions)));% =1 for images
-    CheckNc{iview}=~isempty(find(strcmp(FileType{iview},NcTypeOptions)));% =1 for netcdf files
-    if ~isempty(j1_series{iview})
-        frame_index{iview}=j1_series{iview};
-    else
-        frame_index{iview}=i1_series{iview};
-    end
-end
-
-%% calibration data and timing: read the ImaDoc files
-XmlData=[];
-[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
-        msgbox_uvmat('WARNING',['times of series differ by (max) ' num2str(diff_time)])
-    end   
-end
-
-%% coordinate transform or other user defined transform
-transform_fct='';%default
-if isfield(Param,'FieldTransform')&&~isempty(Param.FieldTransform.TransformName)
-    addpath(Param.FieldTransform.TransformPath)
-    transform_fct=str2func(Param.FieldTransform.TransformName);
-    rmpath(Param.FieldTransform.TransformPath)
-end
-
-%%%%%%%%%%%% END STANDARD PART  %%%%%%%%%%%%
- % EDIT FROM HERE
-
-%% check the validity of  input file types
-if CheckImage{1}
-    FileExtOut='.png'; % write result as .png images for image inputs
-elseif CheckNc{1}
-    FileExtOut='.nc';% write result as .nc files for netcdf inputs
+%% determine the file type on each line from the first input file
+
+[FileInfo,MovieObject]=get_file_info(filecell{1,1});
+FileType=FileInfo.FileType;
+
+if ~isempty(j1_series{1})
+    frame_index=j1_series{1};
 else
-    msgbox_uvmat('ERROR',['invalid file type input ' FileType{1}])
-    return
-end
-
+    frame_index=i1_series{1};
+end
 
 %% settings for the output file
-NomTypeOut=nomtype2pair(NomType{1});% determine the index nomenclature type for the output file
+FileExtOut='.nc';% write result as .nc files for netcdf inputs
+NomTypeOut=nomtype2pair(NomType);% determine the index nomenclature type for the output file
 first_i=i1_series{1}(1);
 last_i=i1_series{1}(end);
@@ -177 +138 @@
 %% Set field names and velocity types
 InputFields{1}=[];%default (case of images)
-if isfield(Param,'InputFields')
-    InputFields{1}=Param.InputFields;
-end
-
+%%%%%%%%%%%%%%%%%%%%%
+if isfield(Param.IndexRange,'TimeSource') && find(regexp(Param.IndexRange.TimeSource,'^var:'))==1
+    Param.InputFields.TimeVarName=regexprep(Param.IndexRange.TimeSource,'^var:','');
+end
+%%%%%%%%%%%%%%%%%%%%
+% if isfield(Param,'InputFields')
+%     InputFields{1}=Param.InputFields;
+% end
+Param.InputFields.FieldName={'Uprime','Vprime'};
 nbfiles=0;
 nbmissing=0;
@@ -197 +163 @@
 disp('loop for mean started')
 for index=1:NbField
-    update_waitbar(WaitbarHandle,index/NbField)
-    if ~isempty(RUNHandle)&& ~strcmp(get(RUNHandle,'BusyAction'),'queue')
-        disp('program stopped by user')
-        break
-    end
-    [Field,tild,errormsg] = read_field(filecell{1,index},FileType{iview},InputFields{iview},frame_index{iview}(index));
+    %     if ~isempty(RUNHandle)&& ~strcmp(get(RUNHandle,'BusyAction'),'queue')
+    %         disp('program stopped by user')
+    %         break
+    %     end
+    [Field,~,errormsg] = read_field(filecell{1,index},FileType,Param.InputFields,frame_index(index));
+    
+    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+    Field.U=Field.Uprime;
+    Field.V=Field.Vprime;
+    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     if index==1 %first field
         if ~isfield(Field,'U')||~isfield(Field,'V')
@@ -208 +178 @@
             return
         end
+        
         [npy,npx]=size(Field.U);
         UMean=zeros(npy,npx);
         VMean=zeros(npy,npx);
         Counter=zeros(npy,npx);
-                % transcripts the global attributes
+        % transcripts the global attributes
         if isfield(Field,'ListGlobalAttribute')
             DataOut.ListGlobalAttribute= Field.ListGlobalAttribute;
@@ -236 +207 @@
 disp('loop for correlation started')
 for index=1:NbField
-    update_waitbar(WaitbarHandle,index/NbField)
     if ~isempty(RUNHandle)&& ~strcmp(get(RUNHandle,'BusyAction'),'queue')
         disp('program stopped by user')
         break
     end
-    [Field,tild,errormsg] = read_field(filecell{1,index},FileType{iview},InputFields{iview},frame_index{iview}(index));
-
+    [Field,~,errormsg] = read_field(filecell{1,index},FileType,Param.InputFields,frame_index(index));
+    Field.U=Field.Uprime;
+    Field.V=Field.Vprime;
     %%%%%%%%%%%% MAIN RUNNING OPERATIONS  %%%%%%%%%%%%
     if index==1 %first field
@@ -281 +252 @@
         VVCorr(ishift+npcorr+1,:,:)=Field.V.*V_shift;
         UVCorr(ishift+npcorr+1,:,:)=Field.U.*V_shift;
-        FFCorr(ishift+npcorr+1,:,:)=FF | FF_shift;      
+        FFCorr(ishift+npcorr+1,:,:)=FF | FF_shift;
     end
     DataOut.UUCorr=DataOut.UUCorr+UUCorr;
@@ -322 +293 @@
 
 %% writing the result file as netcdf file
-OutputFile=fullfile_uvmat(RootPath{1},OutputDir,RootFile{1},FileExtOut,NomTypeOut,first_i,last_i,first_j,last_j);
+OutputFile=fullfile_uvmat(RootPath,OutputDir,RootFile,FileExtOut,NomTypeOut,first_i,last_i,first_j,last_j);
  %case of netcdf input file , determine global attributes
  errormsg=struct2nc(OutputFile,DataOut); %save result file
@@ -331 +302 @@
  end
 
-
-%% open the result file with uvmat (in RUN mode)
-if checkrun
-    uvmat(OutputFile)% open the last result file with uvmat
-end