Fanbeam command projection directions and attributions clarification
8 views (last 30 days)
Show older comments
oumi k
on 13 Apr 2023
Commented: Bjorn Gustavsson
on 15 Apr 2023
In the fan beam projection in matlab, what is the starting point of the projection acquisitions? Is the source considered placed on the right side on the positive x axis and the detectors on the other side or the opposite?
Also does the source with detectors rotate clockwise or anticlockwise for 360 degrees?
Regarding the detector positions in the sinogram which range from for example -40 to 40 degrees, when the source is at 0 degrees for the fan beam rotation assuming the source is placed horizontally at the beginning (along positive x axis), is -40 attributed to the uppermost detector or the lowest?
0 Comments
Accepted Answer
Bjorn Gustavsson
on 14 Apr 2023
The way I interpret the help and documentation the vertex with the source are located at negative y for rotation-angle zero. For the x-prime coordinate when converting the arc-geometry for the sensors to linear coordinates it is stated that:
This linear spacing is measured on the x-prime axis. The x-prime axis for each column, COL, of F is oriented at FAN_ROTATION_ANGLES(COL) degrees counterclockwise from the x-axis. The origin of both axes is the center pixel of the image.
As to which sensor is to the right on the sensor-array I think this part of one of the examples in the help:
% Convert angular positions to linear distance along x-prime axis
FposArc = D*tan(FposArcDeg*pi/180);
shows that the positive FposArcDeg is to the right in the detector array.
You could always try with a simpler test-image:
TestImg = zeros(128);
TestImg(32:48,32:48) = 1;
TestImg(78:89,64:85) = 5;
TestImg(20:25,89:100) = 3;
That should make it more straightforward to untangle what is where looking in this or that direction.
HTH
2 Comments
More Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!