<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
        {mso-style-name:msonormal;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
span.StileMessaggioDiPostaElettronica18
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;
        font-family:"Calibri",sans-serif;
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 2.0cm 2.0cm 2.0cm;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=IT link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span lang=EN-GB style='mso-fareast-language:EN-US'>Dear Simon and RTK users,<br><br>I’ve been experimenting on the generation of Half Fan CBCT images successfully from reprojections of CTs starting from Simon’s suggestions.<br>So far I was able to reconstruct images by displacing the detector in the X direction (+ or -) and completing a single rotation. Results were good and the FOV was of course larger than the one obtained from using the same virtual detector without displacement.<br><br>I’ve taken the simulation a step further and I’m currently creating a geometry which is similar to the combination of “</span><span lang=EN-GB>rtksimulatedgeometry -n 180 --proj_iso_x <displacement> -o g_1” and “rtksimulatedgeometry -n 180 --proj_iso_x <(-1)*displacement> -o g_2 -f 180” (I’m rotating first between 0° and 180° while displacing by half detector size on +X and then 180° and 360° while displacing by half detector size on -X).<br>With this single .xml I’m reprojecting a CT into a single .mha using rtkforwardprojections and then I’m using the output as input for rtkfdk.<br><br>My results however suffer from a centered artifact, of semi-cylindrical shape, in my opinion caused by the superimposition of rays from the two beams around the isocenter.<br>This is further supported by the fact that the more I displace the detector the smaller the artefact becomes (of course I can’t displace more than 50% of detector size).<br>I guess a possible solution would be to have a perfect half-cone x-ray beam by shaping it using a collimator, but I’m not sure how to proceed on this in the simulated environment.<br>Have you got any suggestions or observation on how to achieve a reconstruction based on this? (two rotations/acquistion given two opposite detector displacements)<br><br>Thanks in advance,<br>Gabriele<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB style='mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-GB style='mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><b><span lang=EN-GB>Da:</span></b><span lang=EN-GB> Simon Rit <simon.rit@creatis.insa-lyon.fr> <br><b>Inviato:</b> venerdì 11 ottobre 2019 13.10<br><b>A:</b> gabriele.belotti.bergamo@gmail.com<br><b>Cc:</b> rtk-users <rtk-users@public.kitware.com><br><b>Oggetto:</b> Re: [Rtk-users] Half Fan dataset<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><div><div><p class=MsoNormal><span lang=EN-GB>Hi,<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-GB>It's easy to generate, you need to offset your detector, either via the RTK geometry or by setting the first coordinate of the origin of your projection to something which makes the projection uncentered. For example, in the geometry :<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-GB>rtksimulatedgeometry -n 180 --proj_iso_x 100 -o g<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-GB>rtkprojectshepploganphantom -g g -o proj.mha<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-GB>rtkfdk -p . -g g -r proj.mha -o fdk.mha<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-GB>You can simulate from a CT image by following </span><a href="http://wiki.openrtk.org/index.php/RTK/Scripts/ForwardProjection"><span lang=EN-GB>this example</span></a><span lang=EN-GB>.<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-GB>Simon<o:p></o:p></span></p></div></div><p class=MsoNormal><span lang=EN-GB><o:p> </o:p></span></p><div><div><p class=MsoNormal><span lang=EN-GB>On Fri, Oct 11, 2019 at 9:58 AM <</span><a href="mailto:gabriele.belotti.bergamo@gmail.com" target="_blank"><span lang=EN-GB>gabriele.belotti.bergamo@gmail.com</span></a><span lang=EN-GB>> wrote:<o:p></o:p></span></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt'><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span lang=EN-GB>Dear RTK users and developers,<br><br>I’m currently experimenting with FDK reconstruction and I’m struggling to find a Half-Fan projection dataset to fiddle around.. Do you know where I can find one? I’ve taken into consideration generating a set of DRRs from an existing phantom. Any help or advice you can give me would be greatly appreciated, thanks!<br><br>Gabriele Belotti<o:p></o:p></span></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span lang=EN-GB> <o:p></o:p></span></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span lang=EN-GB> <o:p></o:p></span></p></div></div><p class=MsoNormal><span lang=EN-GB>_______________________________________________<br>Rtk-users mailing list<br></span><a href="mailto:Rtk-users@public.kitware.com" target="_blank"><span lang=EN-GB>Rtk-users@public.kitware.com</span></a><span lang=EN-GB><br></span><a href="https://public.kitware.com/mailman/listinfo/rtk-users" target="_blank"><span lang=EN-GB>https://public.kitware.com/mailman/listinfo/rtk-users</span></a><span lang=EN-GB><o:p></o:p></span></p></blockquote></div></div></body></html>