<div dir="ltr"><div><div><div>Please answer to the mailing list.<br></div>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.<br></div>Thanks for sharing and we'll let you know if we implement something that could fix your issue.<br></div>Simon<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 11, 2015 at 3:46 PM, Padraig Looney <span dir="ltr"><<a href="mailto:padraig.looney@gmail.com" target="_blank">padraig.looney@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I've not used rtkfdk. The code for the back projection that works for me is <br><br> typedef rtk::ProjectionGeometry<3> ProjectionGeometry;<br> ProjectionGeometry::Pointer baseGeom = geometry.GetPointer();<br> typedef rtk::BackProjectionImageFilter< FloatImageType ,FloatImageType> FDKCPUType;<br> FDKCPUType::Pointer feldkamp = FDKCPUType::New();<br> feldkamp->SetInput( 0, constantImageSource->GetOutput() );<br> feldkamp->SetInput( 1, projectionStack);<br> feldkamp->SetGeometry( baseGeom );<br> feldkamp->Update();<br><br></div>We did the same with the rtk::FDKConeBeamReconstructionFilter but as I mention we have to adapt the ramp filter. <br></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 11, 2015 at 2:43 PM, Simon Rit <span dir="ltr"><<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>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?<br></div>I think it would be interesting for everyone if you could share this.<br></div>Thanks,<br></div>Simon<br></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 11, 2015 at 3:22 PM, Padraig Looney <span dir="ltr"><<a href="mailto:padraig.looney@gmail.com" target="_blank">padraig.looney@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Thanks for looking into this.<br><br></div>Yes. We can use the back projection filters fine. Although, we have yet to modify the ramp filters to suit our usecase.<br><br></div><br></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 11, 2015 at 2:12 PM, Simon Rit <span dir="ltr"><<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi,<br>I have received your data and run the attached test. I also get the error<br>"Description: itk::ERROR: Cannot account for too large detector displacements, a part of space must be covered by all projections."<br>wtih rtkfdk. Did you use rtkfdk for testing FDK? Your email suggested that it worked for FDK but it doesn't work for me...<br>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 <a href="http://www.openrtk.org/Doxygen/dot_inline_dotgraph_20.png" target="_blank">this image</a> for sart.<span><font color="#888888"><br>Simon<br></font></span><br></div>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<br></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 10, 2015 at 1:40 PM, Simon Rit <span dir="ltr"><<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi,<br></div>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.<br></div>Thanks,<br></div>Simon<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Jun 10, 2015 at 1:35 PM, Padraig Looney <span dir="ltr"><<a href="mailto:padraig.looney@gmail.com" target="_blank">padraig.looney@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div>Dear All,<br><br><br></div><div>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<br><br>itk::ERROR: Cannot account for detector displacement larger than 50% of panel size. Corner inf=-0 and corner sup=-163.422<br><br></div><div>Any help would be appreciated.<br><br></div><div>Thanks<br></div><div>Pádraig<br></div></div>
<br></div></div>_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@public.kitware.com" target="_blank">Rtk-users@public.kitware.com</a><br>
<a href="http://public.kitware.com/mailman/listinfo/rtk-users" target="_blank">http://public.kitware.com/mailman/listinfo/rtk-users</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>