Hi,
I am trying to understand the RTK geometry, but I am a bit confused. I
wrote the following:
#include "rtkThreeDCircularProjectionGeometry.h"
typedef rtk::ThreeDCircularProjectionGeometry GeometryType;
int main(int argc, char *argv[])
{
int argOffset = 1;
double gantryAngle = atof(argv[argOffset++]);
double outOfPlaneAngle = atof(argv[argOffset++]);
double inPlaneAngle = atof(argv[argOffset++]);
double sid = atof(argv[argOffset++]);
double sdd = atof(argv[argOffset++]);
GeometryType::Pointer rtkGeo = GeometryType::New();
rtkGeo->AddProjection(sid, sdd, gantryAngle, 0.0, 0.0, outOfPlaneAngle,
inPlaneAngle, 0.0, 0.0);
std::cout << rtkGeo->GetSourcePosition(0) << std::endl;
return EXIT_SUCCESS;
}
With inputs 0 0 0 400 1000, I expect to get the source at [0,0,-400 1] but
I get [0,0,400,1]. The reason for my expectation is this sentence from
geometry description file:
"With all parameters set to 0, the detector is normal to the z direction of
the fixed coordinate system". I drew a sketch (please see attached) of what
I understood from this sentence, and was wondering if I am correct?
It looks like I am missing something trivial, but couldn't figure out what.
Many thanks in advance.
Kind regards,
Serkan
_______________________________________________
Rtk-users mailing list
[email protected]
http://public.kitware.com/mailman/listinfo/rtk-users
