Hi Louie, Thanks for the report. Could you make a quick test case that does all this so that I test (and try to solve it) locally? You can maybe modify rtkdrawgeometricphantomtest.cxx with your test. Thanks, Simon
On Mon, Jul 18, 2016 at 5:00 PM, louie <[email protected]> wrote: > Dear rtk users, > > I ran into a small problem with the geometric phantom filter > (rtk::DrawGeometricPhantomImageFilter). > > Here comes the description: > > Input: > Geometry: Magnification = 1 (virtual detector at the image center) > > GeometricFile: [Box] A=50. B=50. C=50. x=0. > y=0. z=0. beta=0 gray=5. > Spacing: 1 > Dimension: 150 > centered. > > Projection stack: > Spacing: 1 > Dimension: 1500 (much larger than the projection) > > Projector: > CudaForwardProjectionImageFilter (step=0.1) > > Problem: > If I connect the output of drawGeometricPhantom filter to the input of the > projector (input 1), the result has non-zero pixels on the projection image > which is not part of the projection. > If I write and read the output drawGeometricPhantom filter to and from the > disc first, the result is correct, i.e. there is nothing but the projection > of the object. > > Question: > Can anyone reproduce this? Why is there such a difference? > > Best regards, > Louie > > _______________________________________________ > Rtk-users mailing list > [email protected] > http://public.kitware.com/mailman/listinfo/rtk-users > _______________________________________________ Rtk-users mailing list [email protected] http://public.kitware.com/mailman/listinfo/rtk-users
