<div dir="ltr"><div class="gmail_quote"><div dir="ltr">L.S.,<div><br></div><div>I was working on FDK in the past and interative reconstruction methods are still new to me.</div><div>I understand the concept of iteratvie methods but are not aware of technical details in implementation.</div><div><br></div><div>Recently I am trying SART but got streak artefacts in reconstructed slices, as well as dots with very high value (both negative and positive) at corners of slices.</div><div>When I checked intermediate images in the pipleline I found that those are introduced in itk::<wbr>DivideOrZeroOutImageFilter.</div><div>You can see from the attached picture: the left half shows the output of rtk::<wbr>RayBoxIntersectionImageFilter and the right half the output of itk::<wbr>DivideOrZeroOutImageFilter, both during processing of the first projection in the first iteration.</div><div>Apparently, although it contains the whole object, my volume is relatively small compared to the size of the detector images.</div><div>Then the rays intersecting the volume near corners and edges result in small values in the output of the raybox filter, and subsequently magnify the pixel values

<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">largely</span>
after division.</div><div>This may not be a problem if the detector images are noiseless, but in practice this will magnify the noise and they will stay as streaks and dots in slices.</div><div><br></div><div>To correct for this I have something in mind, such as making the volume bigger and cropping the detector images so that corners and edges of the volume do not project to the cropped detector; or increasing the threshold in the divide filter so that low values from edge/corner rays wll be zero out. Since I am lack of experiences in interative methods, my question is what the best or common practice will be to handle this? Thanks a lot.</div><div><br></div><div>Regards,</div><div>Chao</div><div><br></div></div>
</div><br></div>