Changeset 847 for trunk/src/proj_field.m
- Timestamp:
- Jan 15, 2015, 6:05:57 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/proj_field.m
r813 r847 1329 1329 XIndexMin=1; 1330 1330 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 % end1336 % if testXMin%test_direct(NbDim)==11337 % XIndexMin=(XMin-Coord{NbDim}(1))/DX+1;% matrix index corresponding to the min x value for the new field1338 % XIndexMax=(XMax-Coord{NbDim}(1))/DX+1;% matrix index corresponding to the max x value for the new field1339 % Xbound(1)=Coord{NbDim}(1)+DX*(XIndexMin-1);% x value corresponding to XIndexMin1340 % Xbound(2)=Coord{NbDim}(1)+DX*(XIndexMax-1);% x value corresponding to XIndexMax1341 % else1342 % 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 % end1347 1331 YIndexRange(1)=ceil(min(YIndexMin,YIndexMax));%first y index to select from the previous field 1348 1332 YIndexRange(1)=max(YIndexRange(1),1);% avoid bound lower than the first index … … 1461 1445 VarName=FieldData.ListVarName{ivar}; 1462 1446 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 1464 1448 else 1465 1449 ProjData.(VarName)=interp2(X,Y,double(FieldData.(VarName)(:,:,1)),XI,YI,'*linear'); 1466 1450 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 TEST1451 ProjData.(VarName)=cat(3,ProjData.(VarName),interp2(X,Y,double(FieldData.(VarName)(:,:,icolor)),XI,YI,'*linear')); 1468 1452 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) 1469 1458 end 1470 1459 ListVarName=[ListVarName VarName]; … … 1540 1529 eval(['ProjData.' WName '=NormVec_X*ProjData.' UName '+ NormVec_Y*ProjData.' VName '+ NormVec_Z* ProjData.' WName ';']); 1541 1530 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 % end1546 1531 end 1547 1532 end
Note: See TracChangeset
for help on using the changeset viewer.