<div dir="ltr"><div>Hi,</div><div>There is an <a href="https://github.com/SimonRit/RTK/blob/master/include/rtkAdditiveGaussianNoiseImageFilter.h" target="_blank">AdditiveGaussianNoiseImageFilter</a>. In practice, I personally do it in numpy with the python wrappings, see <a href="https://doi.org/10.1118/1.4945418">this article</a> where I used this code:</div><div><br></div>I0=1e7<br><div>dH2O=0.01879 #mm^-1 at 75 keV</div>if I0!=0:<br>    slarray = rtk.GetArrayFromImage(sl)<br> 
   slarray = I0*np.exp(-1.*dH2O*slarray)                                
                                                         slarray = 
np.maximum(np.random.<span class="gmail-il">poisson</span>(slarray), 1)<br>    slarray = np.log(I0/slarray)/dH2O<br>    slarray = rtk.GetImageFromArray(slarray.astype(np.float32))<br>    slarray.CopyInformation(sl)<br><div>    sl = slarray</div><div><br></div><div>and wrote "The same simulations were repeated with Poisson noise.<br>The Shepp Logan densities were weighted by 0.01879 mm −1 ,<br>i.e., the linear attenuation coefficient of water at 75 keV. The<br>number of photons received per detector pixel without object<br>in the beam was constant for all pixels in both geometries and<br>equal to 10^7 ."</div><div>Simon<br></div><div class="gmail-adL"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, May 15, 2020 at 10:04 AM <<a href="mailto:gabriele.belotti.bergamo@gmail.com" target="_blank">gabriele.belotti.bergamo@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 lang="IT"><div><p class="MsoNormal"><span lang="EN-US">Dear RTK users and developers,<br><br>I’m creating a CBCT dataset from existing CTs and would like to introduce noise into the generated projections before reconstruction.<br>Is there an RTK implementation for this feature? Or is there any software, to you knowledge, that would efficiently add this noise (in terms of computation time)?<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">Thanks in advance,<br>Gabriele Belotti<u></u><u></u></span></p></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="https://public.kitware.com/mailman/listinfo/rtk-users" rel="noreferrer" target="_blank">https://public.kitware.com/mailman/listinfo/rtk-users</a><br>
</blockquote></div>