Changeset 1118 for trunk/src/series


Ignore:
Timestamp:
Jun 2, 2022, 12:22:11 PM (2 years ago)
Author:
sommeria
Message:

civdata from get_field repaired, detail data output activated in stereo_civ

File:
1 edited

Legend:

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

    r1107 r1118  
    149149j2_series_Civ2=j2_series_Civ1;
    150150
    151 % if isempty(j1_series_Civ1)
    152 %     FrameIndex_A_Civ1=i1_series_Civ1;
    153 %     FrameIndex_B_Civ1=i2_series_Civ1;
    154 %     j1_series_Civ1=ones(size(i1_series_Civ1));
    155 %     j2_series_Civ1=ones(size(i1_series_Civ1));
    156 % else
    157 %     FrameIndex_A_Civ1=j1_series_Civ1;
    158 %     FrameIndex_B_Civ1=j2_series_Civ1;
    159 % end
    160151if isempty(PairCiv2)
    161152    FrameIndex_A_Civ2=FrameIndex_A_Civ1;
     
    289280       
    290281    Civ1Dir=OutputDir;
    291 
    292 %         ncfile=fullfile_uvmat(RootPath_A,Civ1Dir,[RootFile_A,'_All'],'.nc',NomTypeNc,i2_series_Civ1(ifield),[],...
    293 %             j1_series_Civ1(ifield),j2_series_Civ1(ifield));
    294        
    295282       
    296283        ncfile2=fullfile_uvmat(RootPath_A,Civ1Dir,RootFile_A,'.nc',NomTypeNc,i2_series_Civ1(ifield),[],...
     
    527514            mask=imread(par_civ2.Mask);
    528515        end
    529 %         ibx2=ceil(par_civ2.CorrBoxSize(1)/2);
    530 %         iby2=ceil(par_civ2.CorrBoxSize(2)/2);
    531516        par_civ2.SearchBoxShift=[Shiftx(nbval>=1)./nbval(nbval>=1) Shifty(nbval>=1)./nbval(nbval>=1)];
    532517        par_civ2.Grid=[par_civ2.Grid(nbval>=1,1)-par_civ2.SearchBoxShift(:,1)/2 par_civ2.Grid(nbval>=1,2)-par_civ2.SearchBoxShift(:,2)/2];% grid taken at the extrapolated origin of the displacement vectors
     
    685670            mask=imread(par_civ3.Mask);
    686671        end
    687 %         ibx2=ceil(par_civ3.CorrBoxSize(1)/2);
    688 %         iby2=ceil(par_civ3.CorrBoxSize(2)/2);
    689672        par_civ3.SearchBoxShift=[Shiftx(nbval>=1)./nbval(nbval>=1) Shifty(nbval>=1)./nbval(nbval>=1)];
    690673        par_civ3.Grid=[par_civ3.Grid(nbval>=1,1)-par_civ3.SearchBoxShift(:,1)/2 par_civ3.Grid(nbval>=1,2)-par_civ3.SearchBoxShift(:,2)/2];% grid taken at the extrapolated origin of the displacement vectors
     
    755738    end
    756739
    757    
    758740     %% Patch3
    759741    if isfield (Param.ActionInput,'Patch3')
     
    785767        Data.Civ3_FF(ind_good)=FFres;
    786768        Data.CivStage=Data.CivStage+1;
    787        
    788            
    789          % get z from u and v (displacements)
    790        
     769           
     770         % get z from u and v (displacements)     
    791771        Data.Xmid=Rangx(1)+(Rangx(2)-Rangx(1))*(Data.Civ3_X-0.5)/(Npx-1);%temporary coordinate (velocity taken at the point middle from imgae 1 and 2)
    792772        Data.Ymid=Rangy(2)+(Rangy(1)-Rangy(2))*(Data.Civ3_Y-0.5)/(Npy-1);%temporary coordinate (velocity taken at the point middle from imgae 1 and 2)
     
    803783   
    804784    %% write result in a netcdf file if requested
    805 %     if LSM ~= 1 % store all data
    806 %         if exist('ncfile','var')
    807 %             errormsg=struct2nc(ncfile,Data);
    808 %             if isempty(errormsg)
    809 %                 disp([ncfile ' written'])
    810 %             else
    811 %                 disp(errormsg)
    812 %             end
    813 %         end
    814 %     else
     785    if LSM ~= 1 % store all data
     786        if exist('ncfile','var')
     787            errormsg=struct2nc(ncfile,Data);
     788            if isempty(errormsg)
     789                disp([ncfile ' written'])
     790            else
     791                disp(errormsg)
     792            end
     793        end
     794    else
    815795       % store only phys data
    816        % Data_light.ListVarName={'Xphys','Yphys','Zphys','Civ3_C','Xmid','Ymid','Uphys','Vphys','Error'};
    817        % Data_light.VarDimName={'nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3'};
    818796        Data_light.ListVarName={'Xphys','Yphys','Zphys','Civ3_C','DX','DY','Error'};
    819797        Data_light.VarDimName={'nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3'};
     
    828806        Data_light.Xphys=Data.Xphys(ind_good);
    829807        Data_light.Civ3_C=Data.Civ3_C(ind_good);
    830 %         Data_light.Xmid=Data.Xmid(ind_good);
    831 %         Data_light.Ymid=Data.Ymid(ind_good);
    832808        Data_light.DX=Data.Uphys(ind_good);
    833809        Data_light.DY=Data.Vphys(ind_good);
     
    842818       end
    843819       
    844 %     end
     820    end
    845821end
    846822disp(['ellapsed time for the loop ' num2str(toc) ' s'])
Note: See TracChangeset for help on using the changeset viewer.