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

[Laura Trubuil]

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