[Rtk-users] geometryReader->GetOutputObject() number of input files

Marc Vila Oliva marc.vila-oliva at creatis.insa-lyon.fr
Mon Jun 10 12:59:42 EDT 2013


Oks! I have MVSC 2010 installed. I will give it a try when I have your
file. 

On Mon, 2013-06-10 at 18:50 +0200, Laura Trubuil wrote:
> If you work on visual studio 2010, as me, I can provide you the 
> procedure I have followed.
> Just let me time to write an understandable document.
> 
> Le 10/06/2013 18:34, Simon Rit a écrit :
> > I tried to reproduce what you mention but couldn't. Here is the patch
> > I have applied :
> > diff --git a/applications/rtkbackprojections/rtkbackprojections.cxx
> > b/applications/rtkbackprojections/rtkbackprojections.cxx
> > index 6949e04..1b772bc 100644
> > --- a/applications/rtkbackprojections/rtkbackprojections.cxx
> > +++ b/applications/rtkbackprojections/rtkbackprojections.cxx
> > @@ -56,6 +56,7 @@ int main(int argc, char * argv[])
> >     TRY_AND_EXIT_ON_ITK_EXCEPTION( geometryReader->GenerateOutputInformation() )
> >     if(args_info.verbose_flag)
> >       std::cout << " done." << std::endl;
> > +  std::cout <<
> > geometryReader->GetOutputObject()->GetGantryAngles().size() <<
> > std::endl;
> >
> >     // Create an empty volume
> >     typedef rtk::ConstantImageSource< OutputImageType > ConstantImageSourceType;
> >
> > And this gives me 15 on the standard output, as expected. Can you send
> > us something reproductible?
> > Simon
> >
> > On Mon, Jun 10, 2013 at 6:25 PM, Laura Trubuil
> > <laura.trubuil at mail.medecom.fr> wrote:
> >> Hi Simon,
> >> What I am referring to is not really an example sorry.
> >> I am referring to the directory applications/rtkbackprojections. Using the
> >> rtksimulatedgeometry.exe and then the rtkbackprojections.exe (giving them
> >> the right parameters -o output.mha  -p ... etc), I perform reconstruction.
> >>
> >> I am studying the rtkbackprojections.cxx file in debug mode. I wanted to
> >> know exactly what was generated by the   geometryReader->GetOutputObject()
> >> in the line
> >> bp->SetGeometry( geometryReader->GetOutputObject() );
> >>
> >> So I typed :
> >> rtk::BackProjectionImageFilter<OutputImageType,
> >> OutputImageType>::GeometryPointer geometryreaderOutputObject =
> >> geometryReader->GetOutputObject();
> >>
> >> and looked inside my geometryreaderOutputObject.
> >>
> >> You can find attached my geometry.xml file.
> >>
> >> Regards,
> >>
> >> Laura
> >>
> >>
> >> Le 10/06/2013 18:11, Simon Rit a écrit :
> >>
> >> Hi Laura,
> >> Can you be a bit more specific on which example you are referring to?
> >> You seem to be using the geometry reader, maybe you could provide the
> >> geometry xml file.
> >> Thanks,
> >> Simon
> >>
> >> On Mon, Jun 10, 2013 at 5:44 PM, Laura Trubuil
> >> <laura.trubuil at mail.medecom.fr> wrote:
> >>
> >> Dear RTK users,
> >> I am studiying the rtkbackprojections application example. While looking at
> >> it in debug mode, I found something I cannot explain.
> >> The data I study, generated by geometryReader->GetOutputObject(), should
> >> tell me I have 15 projections.
> >> The gantry angle should be between 0 and 360°, and each step 24°.
> >> But,  if I look at the object generated, the gantry angle starts at 0, then
> >> goes until the 14th projection to 360°... And then do it again from the 15th
> >> to 29th. I finally have 2 series of projections (30 instead of 15). Why?
> >>
> >> Regards,
> >>
> >> Laura
> >>
> >> _______________________________________________
> >> Rtk-users mailing list
> >> Rtk-users at openrtk.org
> >> http://public.kitware.com/cgi-bin/mailman/listinfo/rtk-users
> >>
> >>
> 
> 
> _______________________________________________
> Rtk-users mailing list
> Rtk-users at openrtk.org
> http://public.kitware.com/cgi-bin/mailman/listinfo/rtk-users





More information about the Rtk-users mailing list