<div dir="ltr"><div>Hi,</div><div>Please keep the mailing list in the loop.</div><div>Are the distances in pixel units? If not, you need to set the spacing of the projections according to the pixel spacing of your projections.</div><div>Simon<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 19, 2024 at 8:56 AM yx deng <<a href="mailto:yxd15227813801@gmail.com">yxd15227813801@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi,Simon<br>I had red the doc of 

<span style="color:rgb(0,0,0);font-family:Roboto,sans-serif;background-color:rgb(249,250,252)"><b>RTK 3D circular projection geometry</b></span>   carefully,    <span style="color:rgb(0,0,0);font-family:Roboto,sans-serif;font-size:14px">GantryAngle is the rotation around y, and the code of setting geometry  as follows:<br></span>using GeometryType = rtk::ThreeDCircularProjectionGeometry;<br>  GeometryType::Pointer geometry = GeometryType::New();<br>  unsigned int          numberOfProjections = 1201;<br>  double                firstAngle = 0;<br>  double                angularArc = 360;<br>  unsigned int          sid = 66.;// source to isocenter distance<br>  unsigned int          sdd = 199.;  // source to detector distance<br>  for (unsigned int noProj = 0; noProj < numberOfProjections; noProj++)<br>  {<br><br>    double angle = firstAngle + noProj * angularArc / numberOfProjections;<br>    double ProjOffsetX = 0.;<br>    double ProjOffsetY = 0.;<br>    double OutoffPlanAngle = 0;<br>    double InPlanAngle = 0; <br>    double SourceOffSetX = 0.;<br>    double SourceOffSetY = 0.;<br>    geometry->AddProjection(<br>      sid, sdd, angle, ProjOffsetX, ProjOffsetY, OutoffPlanAngle, InPlanAngle, SourceOffSetX, SourceOffSetY);<br>  }<br><span style="color:rgb(0,0,0);font-family:Roboto,sans-serif;font-size:14px">I also tried to change the value of </span>

OutoffPlanAngle and InPlanAngle, but  it didn't work.Looking forward to your answer.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Simon Rit <<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>> 于2024年6月18日周二 20:24写道:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi,</div><div>It seems to be a geometry problem. I would suggest to check if your geometry object is adequately set, see doc <a href="http://www.openrtk.org/Doxygen/DocGeo3D.html" target="_blank">here</a>.</div><div>Simon<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 18, 2024 at 1:32 PM yx deng <<a href="mailto:yxd15227813801@gmail.com" target="_blank">yxd15227813801@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><pre style="color:rgb(0,0,0)"><i>Hi Simon,</i>
<pre><i>I'm a image processing engineer doing project in image</i><i> reconstruction using rtk.I met a question that the shape of the region corresponding to the reconstructed target changes from a rectangle to a diamond in the slice.
As the pictures show:The result should be a circle, not a diamond.Please give me some advices to correct the reconstruction ,thanks a lot.
</i></pre><img src="cid:ii_lxkahm6p0" alt="企业微信截图_1718708067845.png" width="132" height="132" style="margin-right: 0px;"> <img src="cid:ii_lxkahm7a2" alt="企业微信截图_17187080755468.png" width="133" height="133" style="margin-right: 0px;"> <img src="cid:ii_lxkahm7i3" alt="企业微信截图_17187080843559.png" width="134" height="134" style="margin-right: 0px;"></pre><pre style="color:rgb(0,0,0)">the projection images'size is 972*768*1200,
origin is (-485.5,-383.5,0),spacing is (1.0,1.0,1.0)</pre><pre style="color:rgb(0,0,0)">the volume's size is 768*768*128</pre><pre style="color:rgb(0,0,0)"><pre>origin is (-383.5,-383.5,-63.5),</pre><pre>spacing is (1.0,1.0,1.0).</pre><pre>the class of CudaFDKConeBeamReconstructionFilter was applied.</pre>


<pre><br></pre>



</pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)"></pre><pre style="color:rgb(0,0,0)"><i>
</i></pre><br></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>