[Rtk-users] Source position - ThreeDCircularProjectionGeometry

[Serkan Cimen]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20161021/5e31f524/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Geometry.png
Type: image/png
Size: 53721 bytes
Desc: not available
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20161021/5e31f524/attachment.png>