[Rtk-users] Digital Breast Tomosynthesis - Iterative reconstruction

[Simon Rit]

Please answer to the mailing list.
This makes sense then, by doing your own code you did not include the
displaced detector filter. You have to shunt it in the iterative algorithms
if you want to do the same.
Thanks for sharing and we'll let you know if we implement something that
could fix your issue.
Simon

On Thu, Jun 11, 2015 at 3:46 PM, Padraig Looney <padraig.looney at gmail.com>
wrote:

> I've not used rtkfdk. The code for the back projection that works for me
> is
>
>         typedef rtk::ProjectionGeometry<3> ProjectionGeometry;
>         ProjectionGeometry::Pointer baseGeom = geometry.GetPointer();
>         typedef rtk::BackProjectionImageFilter< FloatImageType
> ,FloatImageType> FDKCPUType;
>         FDKCPUType::Pointer feldkamp = FDKCPUType::New();
>         feldkamp->SetInput( 0, constantImageSource->GetOutput() );
>         feldkamp->SetInput( 1, projectionStack);
>         feldkamp->SetGeometry( baseGeom );
>         feldkamp->Update();
>
> We did the same with the rtk::FDKConeBeamReconstructionFilter but as I
> mention we have to adapt the ramp filter.
>
> On Thu, Jun 11, 2015 at 2:43 PM, Simon Rit <simon.rit at creatis.insa-lyon.fr
> > wrote:
>
>> Can you explain why I get an error that you don't get? Did you write your
>> own c++ software or did you use rtkfdk?
>> I think it would be interesting for everyone if you could share this.
>> Thanks,
>> Simon
>>
>> On Thu, Jun 11, 2015 at 3:22 PM, Padraig Looney <padraig.looney at gmail.com
>> > wrote:
>>
>>> Thanks for looking into this.
>>>
>>> Yes. We can use the back projection filters fine. Although, we have yet
>>> to modify the ramp filters to suit our usecase.
>>>
>>>
>>>
>>> On Thu, Jun 11, 2015 at 2:12 PM, Simon Rit <
>>> simon.rit at creatis.insa-lyon.fr> wrote:
>>>
>>>> Hi,
>>>> I have received your data and run the attached test. I also get the
>>>> error
>>>> "Description: itk::ERROR: Cannot account for too large detector
>>>> displacements, a part of space must be covered by all projections."
>>>> wtih rtkfdk. Did you use rtkfdk for testing FDK? Your email suggested
>>>> that it worked for FDK but it doesn't work for me...
>>>> Given your geometry, this is the expected behavior. If you still want
>>>> to try the iterative reconstruction, you'll have to shunt the displaced
>>>> detector filter in the pipelines using the documentation, see for example this
>>>> image <http://www.openrtk.org/Doxygen/dot_inline_dotgraph_20.png> for
>>>> sart.
>>>> Simon
>>>>
>>>> PS: maybe we will consider changing the error to a simple warning and
>>>> doing nothing in such cases, I'll discuss it with the local iterative pope
>>>>
>>>> On Wed, Jun 10, 2015 at 1:40 PM, Simon Rit <
>>>> simon.rit at creatis.insa-lyon.fr> wrote:
>>>>
>>>>> Hi,
>>>>> Can you send the mhd header and the geometry file to reproduce the
>>>>> bug? You don't have to send the data, we should be able to simulate it.
>>>>> Thanks,
>>>>> Simon
>>>>>
>>>>> On Wed, Jun 10, 2015 at 1:35 PM, Padraig Looney <
>>>>> padraig.looney at gmail.com> wrote:
>>>>>
>>>>>> Dear All,
>>>>>>
>>>>>>
>>>>>> In my group we have been using RTK for digital breast tomosynthesis.
>>>>>> We have successfully used the back projection FDK algorithms to reconstruct
>>>>>> our stack of planes. I have been trying to use the SART filter instead to
>>>>>> compare. The geometry has remained the same as for the case of the back
>>>>>> projection but I now get an ITK error at runtime
>>>>>>
>>>>>> itk::ERROR: Cannot account for detector displacement larger than 50%
>>>>>> of panel size. Corner inf=-0 and corner sup=-163.422
>>>>>>
>>>>>> Any help would be appreciated.
>>>>>>
>>>>>> Thanks
>>>>>> Pádraig
>>>>>>
>>>>>> _______________________________________________
>>>>>> Rtk-users mailing list
>>>>>> Rtk-users at public.kitware.com
>>>>>> http://public.kitware.com/mailman/listinfo/rtk-users
>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/rtk-users/attachments/20150611/06f242a6/attachment-0002.html>