<div dir="ltr"><div>Dear Chao,</div><div>Thanks for the question and sorry for the delay in answering it. The boolean allows to determine if the box goes from itk::Image::Origin to the coordinate of the last pixel or whether you want to add a half voxel border, see code <a href="https://github.com/SimonRit/RTK/blob/master/src/rtkBoxShape.cxx#L147-L168">here</a>. In SART, you want to divide the difference between the projection and the measure by the projection of a volume filled with one. Our projectors consider that the intensity is 0 before Origin and after the coordinate of the last pixel so I think the current version is the correct one. The <a href="https://github.com/SimonRit/RTK/blob/master/test/rtkforwardprojectiontest.cxx">rtkforwardprojectiontest.cxx</a> tests this and you can check in it that there is indeed no half pixel border. So I believe this is the right thing to do, yes.</div><div>Best regards,</div><div>Simon<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Aug 19, 2019 at 10:15 PM Chao Wu <<a href="mailto:wuchao04@gmail.com" target="_blank">wuchao04@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">L.S.,<div><br></div><div>I have noticed a behaviour change in SART after commit 5dc1d78.</div><div>In SARTConeBeamReconstructionFilter, m_RayBoxFilter->SetBoxFromImage(this->GetInput(0), false) attempts to pass false to bWithExternalHalfPixelBorder, although before this commit this variable was not used in RayBoxIntersectionImageFilter::SetBoxFromImage and then when rtk::BoxShape::SetBoxFromImage is called the default half pixel border flag was set to true. After this commit a false can finally be passed through. This gives slight different reconstruction results than before. I was just wondering which one is more correct, with or without the half pixel border for the volume used for the normalization. I tried to follow the external link of explanation given in the code of SetBoxFromImage but found it broken. Could anyone give a brief explanation? Thanks.</div><div><br></div><div>Best regards,</div><div>Chao</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></div>