Changeset 238 for trunk/src/geometry_calib.m
 Timestamp:
 Apr 18, 2011, 10:42:29 PM (10 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/geometry_calib.m
r231 r238 458 458 norm=abs(det(R)); 459 459 GeometryCalib.CalibrationType='linear'; 460 if (a_X1(2)/a_Y1(3))>0 461 epsilon=1; 460 462 GeometryCalib.fx_fy(1)=sqrt((a_X1(2)/a_Y1(3))*norm); 463 else 464 GeometryCalib.fx_fy(1)=sqrt((a_X1(2)/a_Y1(3))*norm); 465 epsilon=1; 466 end 461 467 GeometryCalib.fx_fy(2)=(a_Y1(3)/a_X1(2))*GeometryCalib.fx_fy(1); 462 468 GeometryCalib.CoordUnit=[];% default value, to be updated by the calling function 463 469 %GeometryCalib.Tx_Ty_Tz=[a_X1(1) a_Y1(1) 1]; 464 470 GeometryCalib.Tx_Ty_Tz=[a_X1(1)/GeometryCalib.fx_fy(1) a_Y1(1)/GeometryCalib.fx_fy(2) 1]; 465 R(1,:)= R(1,:)/GeometryCalib.fx_fy(1);471 R(1,:)=epsilon*R(1,:)/GeometryCalib.fx_fy(1); 466 472 R(2,:)=R(2,:)/GeometryCalib.fx_fy(2); 467 473 R=[R;[0 0]]; 468 GeometryCalib.R=[R [0;0; 1]];469 GeometryCalib.omc=(180/pi)*[acos(GeometryCalib.R(1,1)) 0 0] ;474 GeometryCalib.R=[R [0;0;1]]; 475 GeometryCalib.omc=(180/pi)*[acos(GeometryCalib.R(1,1)) 0 0] 470 476 % 471 477 % determine the tsai parameters for a view normal to the grid plane
Note: See TracChangeset
for help on using the changeset viewer.