Changeset 1143 for trunk/src/series

Timestamp:

May 9, 2024, 6:51:30 PM (5 months ago)

Author:

sommeria

Message:

new display of false vectors

Location:

trunk/src/series

Files:

• civ_input.fig (modified) (previous)
• civ_input.m (modified) (2 diffs)
• civ_series.m (modified) (12 diffs)
• extract_rdvision.m (modified) (2 diffs)
• test_patch_tps.m (modified) (2 diffs)

trunk/src/series/civ_input.m

@@ -81 +81 @@
     set(handles.CheckThreshold,'Visible','on')
     set(handles.CheckDeformation,'Value',0)% desactivate 
-    set(handles.num_SubDomainSize(1),'String','250')
-    set(handles.num_SubDomainSize(2),'String','500')
+%     set(handles.num_SubDomainSize(1),'String','250')
+%     set(handles.num_SubDomainSize(2),'String','500')
 end
 switch Param.Action.ActionName
@@ -97 +97 @@
 FileType='image';%fdefault
 FileInfo=[];
-if isfield(SeriesData,'FileInfo')...
+if isfield(SeriesData,'FileInfo')
     FileType=SeriesData.FileInfo{1}.FileType;% info on the first input file series
     FieldType=SeriesData.FileInfo{1}.FieldType;% info on the first input file series

trunk/src/series/civ_series.m

@@ -2 +2 @@
 % --- call the sub-functions:
 %   civ: PIV function itself
-%   fix: removes false vectors after detection by various criteria
+%   detect_false: put a flag to false vectors after detection by various criteria
 %   filter_tps: make interpolation-smoothing 
 %------------------------------------------------------------------------
@@ -22 +22 @@
 %     Param.ActionInput: substructure with the parameters provided by the GUI civ_input
 %                      .Civ1: parameters for civ1cc
-%                      .Fix1: parameters for fix1
+%                      .Fix1: parameters for detect_false1
 %                      .Patch1:
 %                      .Civ2: for civ2
@@ -131 +131 @@
             iview_A=2;% the second line is used for the input images of Civ2
         end
-%         if strcmp(Param.ActionInput.ListCompareMode,'shift')
-%             iview_B=iview_A+1; % the second image series is on the next line of the input table
-%         end
         if iview_A~=0
             RootPath_A=Param.InputTable{iview_A,1};
@@ -237 +234 @@
             NbField=numel(i1_series_Civ2);
         else
-            NbField=numel(i1_series_Civ1);% no image used (only fix or patch) TO CHECK
+            NbField=numel(i1_series_Civ1);% no image used (only detect_false or patch) TO CHECK
         end
 
@@ -347 +344 @@
     ImageName_A='';ImageName_B='';%default
     VideoObject_A=[];VideoObject_B=[];
+    
     %% Civ1
     % if Civ1 computation is requested
@@ -544 +542 @@
     %% Fix1
     if isfield (Param.ActionInput,'Fix1')
-        disp('fix1 started')
+        disp('detect_false1 started')
         if ~isfield (Param.ActionInput,'Civ1')% if we use existing Civ1, remove previous data beyond Civ1
             Fix1_attr=find(strcmp('Fix1',Data.ListGlobalAttribute));
@@ -563 +561 @@
         nbvar=length(Data.ListVarName);
         Data.VarAttribute{nbvar}.Role='errorflag';
-        Data.Civ1_FF=int8(fix(Param.ActionInput.Fix1,Data.Civ1_F,Data.Civ1_C,Data.Civ1_U,Data.Civ1_V));
+        Data.Civ1_FF=int8(detect_false(Param.ActionInput.Fix1,Data.Civ1_F,Data.Civ1_C,Data.Civ1_U,Data.Civ1_V));
         Data.CivStage=2;
     end
@@ -823 +821 @@
         if exist('ncfile','var')
             CivFile=ncfile;
-            [Data,tild,tild,errormsg]=nc2struct(CivFile);%read civ1 and fix1 data in the existing netcdf file
+            [Data,tild,tild,errormsg]=nc2struct(CivFile);%read civ1 and detect_false1 data in the existing netcdf file
             if ~isempty(errormsg)
                 disp_uvmat('ERROR',errormsg,checkrun)
@@ -842 +840 @@
     %% Fix2
     if isfield (Param.ActionInput,'Fix2')
-        disp('fix2 started')
+        disp('detect_false2 started')
         list_param=fieldnames(Param.ActionInput.Fix2)';
         Fix2_param=regexprep(list_param,'^.+','Fix2_$0');% insert 'Fix1_' before  each string in ListFixParam
@@ -857 +855 @@
                 Data.VarDimName=[Data.VarDimName {'nb_vectors2'}];
             end
-            Data.vec_FixFlag=fix(Param.Fix2,Data.vec2_F,Data.vec2_C,Data.vec2_U,Data.vec2_V,Data.vec2_X,Data.vec2_Y);
+            Data.vec_FixFlag=detect_false(Param.Fix2,Data.vec2_F,Data.vec2_C,Data.vec2_U,Data.vec2_V);
         else
             Data.ListVarName=[Data.ListVarName {'Civ2_FF'}];
@@ -863 +861 @@
             nbvar=length(Data.ListVarName);
             Data.VarAttribute{nbvar}.Role='errorflag';
-            Data.Civ2_FF=double(fix(Param.ActionInput.Fix2,Data.Civ2_F,Data.Civ2_C,Data.Civ2_U,Data.Civ2_V));
+            Data.Civ2_FF=double(detect_false(Param.ActionInput.Fix2,Data.Civ2_F,Data.Civ2_C,Data.Civ2_U,Data.Civ2_V));
             Data.CivStage=Data.CivStage+1;
         end
@@ -1293 +1291 @@
 end
 
-%'RUN_FIX': function for fixing velocity fields:
-%-----------------------------------------------
-% RUN_FIX(filename,field,flagindex,thresh_vecC,thresh_vel,iter,flag_mask,maskname,fileref,fieldref)
-%
-%filename: name of the netcdf file (used as input and output)
-%field: structure specifying the names of the fields to fix (depending on civ1 or civ2)
-%.vel_type='civ1' or 'civ2';
-%.nb=name of the dimension common to the field to fix ('nb_vectors' for civ1);
-%.fixflag=name of fix flag variable ('vec_FixFlag' for civ1)
-%flagindex: flag specifying which values of vec_f are removed:
-% if flagindex(1)=1: vec_f=-2 vectors are removed
-% if flagindex(2)=1: vec_f=3 vectors are removed
-% if flagindex(3)=1: vec_f=2 vectors are removed (if iter=1) or vec_f=4 vectors are removed (if iter=2)
-%iter=1 for civ1 fields and iter=2 for civ2 fields
-%thresh_vecC: threshold in the image correlation vec_C
-%flag_mask: =1 mask used to remove vectors (0 else)
-%maskname: name of the mask image file for fix
-%thresh_vel: threshold on velocity, or on the difference with the reference file fileref if exists
-%inf_sup=1: remove values smaller than threshold thresh_vel, =2, larger than threshold
-%fileref: .nc file name for a reference velocity (='': refrence 0 used)
-%fieldref: 'civ1','filter1'...feld used in fileref
-
-function FF=fix(Param,F,C,U,V,X,Y)
+
+function FF=detect_false(Param,F,C,U,V)
 FF=zeros(size(F));%default
-
-%criterium on warn flags
-FlagName={'CheckFmin2','CheckF2','CheckF3','CheckF4'};
-FlagVal=[-2 2 3 4];
-for iflag=1:numel(FlagName)
-    if isfield(Param,FlagName{iflag}) && Param.(FlagName{iflag})
-        FF=(FF==1| F==FlagVal(iflag));
-    end
-end
-%criterium on correlation values
+% FF=-2, for correlation max at edge
+% FF=-1, for too small correlation
+% FF=1, for velocity outside bounds
+% FF=2 for exclusion by difference with the smoothed field
+FF(F==-2)=-2;
 if isfield (Param,'MinCorr')
-    FF=FF==1 | C<Param.MinCorr;
+     FF(C<Param.MinCorr)=-1;
 end
 if (isfield(Param,'MinVel')&&~isempty(Param.MinVel))||(isfield (Param,'MaxVel')&&~isempty(Param.MaxVel))
     Umod= U.*U+V.*V;
     if isfield (Param,'MinVel')&&~isempty(Param.MinVel)
-        FF=FF==1 | Umod<(Param.MinVel*Param.MinVel);
+        U2Min=Param.MinVel*Param.MinVel;
+        FF(Umod<U2Min)=1;
     end
     if isfield (Param,'MaxVel')&&~isempty(Param.MaxVel)
-        FF=FF==1 | Umod>(Param.MaxVel*Param.MaxVel);
-    end
-end
+         U2Max=Param.MaxVel*Param.MaxVel;
+        FF(Umod>U2Max)=1;
+    end
+end
+
+
+%criterium on warn flags
+% FlagName={'CheckFmin2','CheckF2','CheckF3','CheckF4'};
+% FlagVal=[-2 2 3 4];
+% for iflag=1:numel(FlagName)
+%     if isfield(Param,FlagName{iflag}) && Param.(FlagName{iflag})
+%         FF=(FF==1| F==FlagVal(iflag));
+%     end
+% end
+% %criterium on correlation values
+% if isfield (Param,'MinCorr')
+%     FF=FF==1 | C<Param.MinCorr;
+% end
+% if (isfield(Param,'MinVel')&&~isempty(Param.MinVel))||(isfield (Param,'MaxVel')&&~isempty(Param.MaxVel))
+%     Umod= U.*U+V.*V;
+%     if isfield (Param,'MinVel')&&~isempty(Param.MinVel)
+%         FF=FF==1 | Umod<(Param.MinVel*Param.MinVel);
+%     end
+%     if isfield (Param,'MaxVel')&&~isempty(Param.MaxVel)
+%         FF=FF==1 | Umod>(Param.MaxVel*Param.MaxVel);
+%     end
+% end
 
 

trunk/src/series/extract_rdvision.m

@@ -308 +308 @@
 %             m.Data=data;
     %%%%%%%
-        timestamp=zeros(1,numel(m.Data));
-        for ii=1: numel(m.Data)
-            timestamp(ii)=m.Data(ii).timestamp;
-        end
-        if isequal(Param.IndexRange.first_i,1)
+    timestamp=zeros(1,numel(m.Data));
+    for ii=1: numel(m.Data)
+        timestamp(ii)=m.Data(ii).timestamp;
+    end
+    if isequal(Param.IndexRange.first_i,1)
         [nbfield1,nbfield2,msg]=copyfile_modif(filexml,timestamp,newxml); %copy the xml file in the upper folder
         [XmlData,errormsg]=imadoc2struct(newxml);% check reading of the new xml file
@@ -319 +319 @@
             return
         end
-        timestamp=timestamp(1:nbfield1*nbfield2);
-        timestamp=reshape(timestamp,nbfield2,nbfield1);
-        difftime=XmlData.Time(2:end,2:end)'-timestamp;
-        disp(['time from xml and timestamp differ by ' num2str(max(max(abs(difftime))))])
-        if max(abs(difftime))>0.01
-            checkpreserve=1;% will not erase the initial files, possibility of error
-        end      
-        % checking consistency with the xml file
-        %     if ~isequal(SeqData.nb_frames,numel(timestamp))
-        %         disp_uvmat('ERRROR',['inconsistent number of images ' num2str(SeqData.nb_frames) ' with respect to the xml file: ' num2str(numel(timestamp))] ,checkrun);
-        %         return
-        %     end
+        if numel(timestamp)~=nbfield1*nbfield2
+            disp('WARNING: total image number defined by the xml file differs from  the number of frames ')
+        else
+            timestamp=reshape(timestamp,nbfield2,nbfield1);
+            difftime=XmlData.Time(2:end,2:end)'-timestamp;
+            disp(['time from xml and timestamp differ by ' num2str(max(max(abs(difftime))))])
+            if max(abs(difftime))>0.01
+                checkpreserve=1;% will not erase the initial files, possibility of error
+            end
+        end
     else
        [nbfield1,nbfield2,msg]=copyfile_modif(filexml,timestamp,''); 

trunk/src/series/test_patch_tps.m

@@ -220 +220 @@
 for index=1:1%numel(filecell)
     Data=read_field(filecell{1,index},'civdata',Param.InputFields);
-            ind_good=find(Data.FF==0|Data.FF==20);%keep good civ data, and also the ones excluded by the criterium of discrepancy betwween smoothed and raw fields (FF=20)
+            ind_good=find(Data.FF==0|Data.FF==20|Data.FF==2);%keep good civ data, and also the ones excluded by the criterium of discrepancy betwween smoothed and raw fields (FF=2 or 20 (old convention))
         NbGood=numel(ind_good);
         Xin=Data.X(ind_good);
@@ -252 +252 @@
         DataOut.(['U_' str_i{irho}])=U_smooth;
         DataOut.(['V_' str_i{irho}])=V_smooth;
-        % DataOut.FF(ind_false)=FieldSmooth(irho);
     end
     time=toc