Changeset 1146


Ignore:
Timestamp:
May 14, 2024, 7:54:48 PM (2 months ago)
Author:
sommeria
Message:

civ-series adjested for new false vector conventions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/series/civ_series.m

    r1144 r1146  
    817817                return
    818818            end
    819         elseif isfield(Param,'Civ2_X')% use Civ2 data as input in Param (test mode)
    820             Data.ListGlobalAttribute={};
    821             Data.ListVarName={};
    822             Data.VarDimName={};
    823             Data.Civ2_X=Param.Civ2_X;
    824             Data.Civ2_Y=Param.Civ2_Y;
    825             Data.Civ2_U=Param.Civ2_U;
    826             Data.Civ2_V=Param.Civ2_V;
    827             Data.Civ2_FF=Param.Civ2_FF;
     819%         elseif isfield(Param,'Civ2_X')% use Civ2 data as input in Param (test mode)
     820%             Data.ListGlobalAttribute={};
     821%             Data.ListVarName={};
     822%             Data.VarDimName={};
     823%             Data.Civ2_X=Param.Civ2_X;
     824%             Data.Civ2_Y=Param.Civ2_Y;
     825%             Data.Civ2_U=Param.Civ2_U;
     826%             Data.Civ2_V=Param.Civ2_V;
     827%             Data.Civ2_FF=Param.Civ2_FF;
    828828        end
    829829    end
     
    10031003
    10041004%% Apply mask
    1005 % Convention for mask IDEAS TO IMPLEMENT ?
     1005% Convention for mask, IDEAS NOT IMPLEMENTED
    10061006% mask >200 : velocity calculated
    10071007%  200 >=mask>150;velocity not calculated, interpolation allowed (bad spots)
     
    10201020    end
    10211021    check_undefined=(par_civ.Mask<200 & par_civ.Mask>=20 );
    1022     %     par_civ.ImageA(check_undefined)=0;% put image A to zero (i.e. the min image value) in the undefined  area
    1023     %     par_civ.ImageB(check_undefined)=0;% put image B to zero (i.e. the min image value) in the undefined  area
    10241022end
    10251023
     
    10581056            sizemask=sum(sum(mask1_crop))/(numel(subrange1_y)*numel(subrange1_x));%size of the masked part relative to the correlation sub-image
    10591057            if sizemask > 1/2% eliminate point if more than half of the correlation box is masked
    1060                 F(ivec)=3; %
     1058                F(ivec)=1; %
    10611059                utable(ivec)=NaN;
    10621060                vtable(ivec)=NaN;
     
    10721070        end
    10731071        %threshold on image minimum
    1074         if F(ivec)~=3
     1072        if F(ivec)~=1
    10751073            if check_MinIma && (image1_mean < par_civ.MinIma || image2_mean < par_civ.MinIma)
    1076                 F(ivec)=3;
     1074                F(ivec)=1;
    10771075                %threshold on image maximum
    10781076            elseif check_MaxIma && (image1_mean > par_civ.MaxIma || image2_mean > par_civ.MaxIma)
    1079                 F(ivec)=3;
    1080             end
    1081             if F(ivec)==3
     1077                F(ivec)=1;
     1078            end
     1079            if F(ivec)==1
    10821080                utable(ivec)=NaN;
    10831081                vtable(ivec)=NaN;
     
    11341132                            utable(ivec)=0;
    11351133                            vtable(ivec)=0;
    1136                             F(ivec)=3;
     1134                            F(ivec)=1;
    11371135                        end
    11381136                        ctable(ivec)=corrmax/sum_square;% correlation value
    11391137                    catch ME
    1140                         F(ivec)=3;
     1138                        F(ivec)=1;
    11411139                        disp(ME.message)
    11421140                    end
    11431141                else
    1144                     F(ivec)=3;
     1142                    F(ivec)=1;
    11451143                end
    11461144            end
Note: See TracChangeset for help on using the changeset viewer.