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

Timestamp:

Mar 11, 2025, 9:51:39 PM (3 days ago)

Author:

sommeria

Message:

add loglin histogram in uvmat and image rescale in civ_series

File:

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

trunk/src/series/civ_series.m

@@ -51 +51 @@
 %% set the input elements needed on the GUI series when the action is selected in the menu ActionName or InputTable refreshed
 if isstruct(Param) && isequal(Param.Action.RUN,0)% function activated from the GUI series but not RUN
+    if 0==1 %never satisfied but trigger compilation with the appropriate transform functions ('eval' inactive for compilation)
+        ima_rescale
+    end
     path_series=fileparts(which('series'));
     addpath(fullfile(path_series,'series'))
@@ -66 +69 @@
     Data.VelType='off';% menu for selecting the velocity type (options 'off'/'one'/'two',  'off' by default)
     Data.FieldName='on';% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default)
-    Data.FieldTransform = 'off';%can use a transform function
+    Data.FieldTransform = 'on';%can use a transform function
     Data.ProjObject='off';%can use projection object(option 'off'/'on',
     Data.Mask='off';%can use mask option   (option 'off'/'on', 'off' by default)
@@ -99 +102 @@
 % CheckOutputFile=isfield(Param,'OutputSubDir');%= 1 in test use for TestPatch (no nc file produced)
 
-%% input files and indexing (skipped in Test mode)
+%% input files and indexing 
 hseries=findobj(allchild(0),'Tag','series');
 RUNHandle=findobj(hseries,'Tag','RUN');%handle of RUN button in GUI series
@@ -112 +115 @@
 if isfield(Param,'InputTable')
     [filecell,i1_series,i2_series,j1_series,j2_series]=get_file_series(Param);
-    %         iview_A=0;% series index (iview) for the first image series
     iview_B=0;% series index (iview) for the second image series (only non zero for option 'shift' comparing two image series )
     if Param.ActionInput.CheckCiv1
@@ -119 +121 @@
         iview_A=2;% the second line is used for the input images
     end
-    %         if iview_A~=0
     RootPath_A=Param.InputTable{iview_A,1};
     RootFile_A=Param.InputTable{iview_A,3};
@@ -253 +254 @@
     end
 end
-% end
+
+%% introduce input image transform
+if isfield(Param,'FieldTransform')&&~isempty(Param.FieldTransform.TransformName)
+        addpath(Param.FieldTransform.TransformPath)
+    transform_fct=str2func(Param.FieldTransform.TransformName);
+    rmpath(Param.FieldTransform.TransformPath)
+end
+
 
 %%%%% MAIN LOOP %%%%%%
@@ -276 +284 @@
         end
     end
-    %     %     if CheckInputFile
-    %     if CheckOutputFile
     OutputPath=fullfile(Param.OutputPath,Param.Experiment,Param.Device);
-    %         if iview_A==0 % no nc file has been entered
-    %             ncfile=fullfile_uvmat(OutputPath,Param.InputTable{1,2},Param.InputTable{1,3},Param.InputTable{1,5},...
-    %                 NomTypeNc,i1_series_Civ1(ifield),i2_series_Civ1(ifield),j1_series_Civ1(ifield),j2_series_Civ1(ifield));
-    %         else% an existing nc file has been entered
-    %             if iview_A==1% if Civ1 is performed
-    %                 Civ1Dir=OutputDir;
-    %             else
-    %                 Civ1Dir=Param.InputTable{1,2};
-    %             end
+
     if strcmp(Param.ActionInput.ListCompareMode,'PIV')
         ncfile=fullfile_uvmat(OutputPath,OutputDir,RootFile_A,'.nc',NomTypeNc,i1_series_Civ1(ifield),i2_series_Civ1(ifield),...
@@ -295 +293 @@
             j1_series_Civ1(ifield),j2_series_Civ1(ifield));
     end
-    %         end
     ncfile_out=ncfile;% by default
     
@@ -380 +377 @@
             end
             [par_civ1.ImageB,VideoObject_B] = read_image(ImageName_B,FileType_B,VideoObject_B,FrameIndex_B_Civ1(ifield));
+            
         catch ME % display errors in reading input images
             if ~isempty(ME.message)
@@ -385 +383 @@
                 continue
             end
+        end
+        %% user defined image transform
+        if ~isempty(transform_fct)
+               par_civ1 =transform_fct(par_civ1,Param);
         end
         
@@ -597 +599 @@
             par_civ2.Grid(:,2)=reshape(GridY,[],1);% increases with array index
         end
+        
+                %% user defined image transform
+        if ~isempty(transform_fct)
+               par_civ2 =transform_fct(par_civ2,Param);
+        end
+        
         
         % get the guess from patch1 or patch2 (case 'CheckCiv3')