<div dir="ltr"><div><div><div>Hi,<br></div>Attached is the python script using <a href="http://wiki.openrtk.org/index.php?title=SimpleRTK">SimpleRTK</a> that I have used to do the simulation of figure 6 of <a href="http://www.creatis.insa-lyon.fr/~srit/biblio/rit2016.pdf">this publication</a>. The paragraph where I describe what I'm doing starts with "The same simulations were repeated with Poisson noise". The interesting part of the script for you is line 158 where if I'm given a number of photons I0, I add Poisson noise.<br></div>Simon<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 29, 2016 at 12:09 AM, Solomon Tang <span dir="ltr"><<a href="mailto:solomoncztang@gmail.com" target="_blank">solomoncztang@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thanks for the feedback Chao and Simon,<div><br></div><div>My geometry was using default RTK sdd/sid settings. I have now changed it to match the DICOM header from the original images (1085.6 SDD, 595 SID), but not much has qualitatively changed.</div><div><br></div><div>How do you suggest adding photon noise? I have discovered an itkShotNoiseImageFilter but I'm not sure what is an acceptable scaling level. I'm assuming the reconstructed image should be passed through the filter, and not the projection.</div><div><div class="h5"><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 28, 2016 at 1:34 AM, Simon Rit <span dir="ltr"><<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<br>
I don't expect a drastic change, only a slight loss of spatial 
resolution if the ray distance at the isocenter (I agree with Chao that 
it plays an important role) is larger than the original voxel size. 
Maybe it's there but you would need to zoom more to see it.<br>
You would see a more realistic difference if you were adding photon noise to your data.<span><font color="#888888"><br>
Simon</font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 28, 2016 at 9:14 AM, Chao Wu <span dir="ltr"><<a href="mailto:wuchao04@gmail.com" target="_blank">wuchao04@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">What is the magnification factor of your geometry?</div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>2016-06-27 23:56 GMT+02:00 Solomon Tang <span dir="ltr"><<a href="mailto:solomoncztang@gmail.com" target="_blank">solomoncztang@gmail.com</a>></span>:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr">Hi Simon,
<div><br></div><div>I am using RTK to simulate CT acquisitions using different detector sizes to see how this impact on image quality might change some of our in-house metrics. </div><div><br></div><div>The images I have linked to below have been created using rtkforwardprojections with different projection spacings (0.3 isometric and 0.75 isometric) reconstructed with rtkfdk with the same pixel spacing and image dimensions (0.4688x0.4688x0.6 | 512x512x225). The CUDA projection stepsize is equal to the projection spacing. The window levels between images of their respective rows are the same. </div><div><br></div><div>I am simply wondering if the differences between these images are realistic. I would expect the image with a detector size than is more than twice as large as the original would be drastically different when in fact they turn out to be incredibly similar. Are the assumptions made about projection spacing == cuda stepsize == simulated hardware detector size incorrect?</div><div><a href="http://goog_1486088111" target="_blank"><br></a></div><div><a href="https://gyazo.com/e86436826f687a2db4b234699d050450" target="_blank">https://gyazo.com/e86436826f687a2db4b234699d050450<br></a></div><div><br></div><div><a href="https://gyazo.com/ca9612218f082e78ba3082950a27fa4c" target="_blank">https://gyazo.com/ca9612218f082e78ba3082950a27fa4c</a><span><font color="#888888"><br></font></span></div><span><font color="#888888"><div><br></div><div>Solomon</div></font></span></div>
<br></div></div>_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@public.kitware.com" target="_blank">Rtk-users@public.kitware.com</a><br>
<a href="http://public.kitware.com/mailman/listinfo/rtk-users" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/rtk-users</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@public.kitware.com" target="_blank">Rtk-users@public.kitware.com</a><br>
<a href="http://public.kitware.com/mailman/listinfo/rtk-users" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/rtk-users</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div><br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"></div></div></div></div>
</div></div></div></div>
</blockquote></div><br></div>