Ignore:
Timestamp:
Jul 22, 2017, 7:20:02 PM (4 years ago)
Author:
sommeria
Message:

various bugs fixed

File:
1 edited

Legend:

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

    r985 r1009  
    453453        Data.Civ1_FF(ind_good)=FFres;
    454454        Data.CivStage=3;
    455              
    456        
    457        
    458        
    459 %               
    460 %          % get z from u and v (displacements)
    461 %       
    462 %         tempXmid=Rangx(1)+(Rangx(2)-Rangx(1))*(Data.Civ1_X-0.5)/(Npx-1);%temporary coordinate (velocity taken at the point middle from imgae 1 and 2)
    463 %         tempYmid=Rangy(1)+(Rangy(2)-Rangy(1))*(Data.Civ1_Y-0.5)/(Npy-1);%temporary coordinate (velocity taken at the point middle from imgae 1 and 2)
    464 %         tempUphys=Data.Civ1_U_smooth*(Rangx(2)-Rangx(1))/(Npx-1);
    465 %         tempVphys=Data.Civ1_V_smooth*(Rangy(2)-Rangy(1))/(Npy-1);
    466 %         [tempZphys,tempXphys,tempYphys,tempCiv3_E]=shift2z(tempXmid,tempYmid,tempUphys,tempVphys,XmlData); %Data.Xphys and Data.Xphys are real coordinate (geometric correction more accurate than xtemp/ytemp)
    467 %         temp=find(Data.Civ1_FF~=0);
    468 %         tempXmid(temp)=[];
    469 %         tempYmid(temp)=[];
    470 %         tempZphys(temp)=[];
    471 %             
     455               
    472456    end
    473457   
     
    651635        par_civ3.ImageA=par_civ1.ImageA;
    652636        par_civ3.ImageB=par_civ1.ImageB;
    653         %         if ~isfield(Param.Civ1,'ImageA')
    654 %         i1=i1_series_Civ3(ifield);
    655 %         i2=i1;
    656 %         if ~isempty(i2_series_Civ3)
    657 %             i2=i2_series_Civ3(ifield);
    658 %         end
    659 %         j1=1;
    660 %         if ~isempty(j1_series_Civ3)
    661 %             j1=j1_series_Civ3(ifield);
    662 %         end
    663 %         j2=j1;
    664 %         if ~isempty(j2_series_Civ3)
    665 %             j2=j2_series_Civ3(ifield);
    666 %         end
    667637        par_civ3.ImageWidth=size(par_civ3.ImageA,2);
    668638        par_civ3.ImageHeight=size(par_civ3.ImageA,1);
     
    844814%     else
    845815       % store only phys data
    846         Data_light.ListVarName={'Xphys','Yphys','Zphys','Civ3_C','Xmid','Ymid','Uphys','Vphys','Error'};
    847         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'};
    848         temp=find(Data.Civ3_FF==0);
    849         Data_light.Zphys=Data.Zphys(temp);
    850         Data_light.Yphys=Data.Yphys(temp);
    851         Data_light.Xphys=Data.Xphys(temp);
    852         Data_light.Civ3_C=Data.Civ3_C(temp);
    853         Data_light.Xmid=Data.Xmid(temp);
    854         Data_light.Ymid=Data.Ymid(temp);
    855         Data_light.Uphys=Data.Uphys(temp);
    856         Data_light.Vphys=Data.Vphys(temp);
    857         Data_light.Error=Data.Error(temp);
     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'};
     818        Data_light.ListVarName={'Xphys','Yphys','Zphys','Civ3_C','DX','DY','Error'};
     819        Data_light.VarDimName={'nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3','nb_vec_3'};
     820        Data_light.VarAttribute{1}.Role='coord_x';
     821         Data_light.VarAttribute{2}.Role='coord_y';
     822         Data_light.VarAttribute{3}.Role='scalar';
     823         Data_light.VarAttribute{5}.Role='vector_x';
     824         Data_light.VarAttribute{6}.Role='vector_y';
     825        ind_good=find(Data.Civ3_FF==0);
     826        Data_light.Zphys=Data.Zphys(ind_good);
     827        Data_light.Yphys=Data.Yphys(ind_good);
     828        Data_light.Xphys=Data.Xphys(ind_good);
     829        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);
     832        Data_light.DX=Data.Uphys(ind_good);
     833        Data_light.DY=Data.Vphys(ind_good);
     834        Data_light.Error=Data.Error(ind_good);
    858835       if exist('ncfile2','var')
    859836            errormsg=struct2nc(ncfile2,Data_light);
Note: See TracChangeset for help on using the changeset viewer.