Changeset 847 for trunk/src/proj_field.m


Ignore:
Timestamp:
Jan 15, 2015, 6:05:57 PM (6 years ago)
Author:
sommeria
Message:

various

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/proj_field.m

    r813 r847  
    13291329                        XIndexMin=1;
    13301330                    end
    1331 %                         YIndexMin=(Coord{NbDim-1}(1)-YMax)/DY+1;
    1332 % %                         YIndexMax=(Coord{NbDim-1}(1)-YMin)/DY+1;
    1333 %                         Ybound(2)=Coord{NbDim-1}(1)-DY*(YIndexMax-1);
    1334 %                         Ybound(1)=Coord{NbDim-1}(1)-DY*(YIndexMin-1);
    1335 %                     end
    1336 %                     if testXMin%test_direct(NbDim)==1
    1337 %                         XIndexMin=(XMin-Coord{NbDim}(1))/DX+1;% matrix index corresponding to the min x value for the new field
    1338 %                         XIndexMax=(XMax-Coord{NbDim}(1))/DX+1;% matrix index corresponding to the max x value for the new field
    1339 %                         Xbound(1)=Coord{NbDim}(1)+DX*(XIndexMin-1);%  x value corresponding to XIndexMin
    1340 %                         Xbound(2)=Coord{NbDim}(1)+DX*(XIndexMax-1);%  x value corresponding to XIndexMax
    1341 %                     else
    1342 %                         XIndexMin=(Coord{NbDim}(1)-XMax)/DX+1;
    1343 %                         XIndexMax=(Coord{NbDim}(1)-XMin)/DX+1;
    1344 %                         Xbound(2)=Coord{NbDim}(1)+DX*(XIndexMax-1);
    1345 %                         Xbound(1)=Coord{NbDim}(1)+DX*(XIndexMin-1);
    1346 %                     end
    13471331                    YIndexRange(1)=ceil(min(YIndexMin,YIndexMax));%first y index to select from the previous field
    13481332                    YIndexRange(1)=max(YIndexRange(1),1);% avoid bound lower than the first index
     
    14611445                        VarName=FieldData.ListVarName{ivar};
    14621446                        if size(FieldData.(VarName),3)==1
    1463                             ProjData.(VarName)=interp2(X,Y,double(FieldData.(VarName)),XI,YI,'*linear');
     1447                            ProjData.(VarName)=interp2(X,Y,double(FieldData.(VarName)),XI,YI,'*linear');%interpolation fct
    14641448                        else
    14651449                            ProjData.(VarName)=interp2(X,Y,double(FieldData.(VarName)(:,:,1)),XI,YI,'*linear');
    14661450                            for icolor=2:size(FieldData.(VarName),3)% project 'color' components
    1467                                 ProjData.(VarName)=cat(3,ProjData.(VarName),interp2(X,Y,double(FieldData.(VarName)(:,:,icolor)),XI,YI,'*linear')); %TO TEST
     1451                                ProjData.(VarName)=cat(3,ProjData.(VarName),interp2(X,Y,double(FieldData.(VarName)(:,:,icolor)),XI,YI,'*linear'));
    14681452                            end
     1453                        end
     1454                        if isa(FieldData.(VarName),'uint8')
     1455                            ProjData.(VarName)=uint8(ProjData.(VarName));%put result to integer 8 bits if the initial field is integer (image)
     1456                        elseif isa(FieldData.(VarName),'uint16')
     1457                            ProjData.(VarName)=uint16(ProjData.(VarName));%put result to integer 16 bits if the initial field is integer (image)
    14691458                        end
    14701459                        ListVarName=[ListVarName VarName];
     
    15401529                    eval(['ProjData.' WName '=NormVec_X*ProjData.' UName '+ NormVec_Y*ProjData.' VName '+ NormVec_Z* ProjData.' WName ';']);
    15411530                end
    1542                 %                 if ~isequal(Psi,0)
    1543                 %                     eval(['ProjData.' UName '=cos(Psi)* ProjData.' UName '- sin(Psi)*ProjData.' VName ';']);
    1544                 %                     eval(['ProjData.' VName '=sin(Psi)* ProjData.' UName '+ cos(Psi)*ProjData.' VName ';']);
    1545                 %                 end
    15461531            end
    15471532        end
Note: See TracChangeset for help on using the changeset viewer.