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

[Simon Rit]

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
>
>