<div dir="ltr">Hi Simon,<div><br></div><div style>Thanks. I tried disabling ramp filter in <span style="font-family:arial,sans-serif;font-size:14px">rtkFDKConeBeamReconstructionFi</span><span style="font-family:arial,sans-serif;font-size:14px">lter.txx with the new line, but then gets a runtime error as this:</span></div>
<div style><div><font face="arial, sans-serif"><span style="font-size:14px"><br></span></font></div><div><font face="arial, sans-serif"><span style="font-size:14px">ExceptionObject caught with writer->Update()</span></font></div>
<div><font face="arial, sans-serif"><span style="font-size:14px">itk::ExceptionObject (0000000000AFE978)</span></font></div><div><font face="arial, sans-serif"><span style="font-size:14px">Location: "void __cdecl itk::DataObject::PropagateRequestedRegion(void) throw(cl</span></font></div>
<div><font face="arial, sans-serif"><span style="font-size:14px">ass itk::InvalidRequestedRegionError)"</span></font></div><div><font face="arial, sans-serif"><span style="font-size:14px">File: ..\..\..\Code\Common\itkDataObject.cxx</span></font></div>
<div><font face="arial, sans-serif"><span style="font-size:14px">Line: 397</span></font></div><div><font face="arial, sans-serif"><span style="font-size:14px">Description: Requested region is (at least partially) outside the largest possib</span></font><span style="font-size:14px;font-family:arial,sans-serif">le region.</span></div>
<div><span style="font-size:14px;font-family:arial,sans-serif"><br></span></div><div style><span style="font-size:14px;font-family:arial,sans-serif">which I don't get if I use the original line. I tried lowering the reconstructed region, but it didn't help.</span></div>
<div style><font face="arial, sans-serif"><span style="font-size:14px">My projection size is 1024*768, 0.388 mm in both directions. It's half beam data, detector being displaced by 14.8 cm.</span></font></div><div style>
<font face="arial, sans-serif"><span style="font-size:14px">The reconstructed region I always use is 512*100*512, spacing being 0.88,2,0.88 mm. </span></font></div><div style><font face="arial, sans-serif"><span style="font-size:14px"><br>
</span></font></div><div style><font face="arial, sans-serif"><span style="font-size:14px">Cheers,</span></font></div><div style><font face="arial, sans-serif"><span style="font-size:14px">Andy.</span></font></div></div><div style>
<span style="font-family:arial,sans-serif;font-size:14px"><br></span></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2012/12/17 Simon Rit <span dir="ltr"><<a href="mailto:simon.rit@creatis.insa-lyon.fr" target="_blank">simon.rit@creatis.insa-lyon.fr</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Andy,<br>
Handling of the half-beam data (we call it displaced detector in RTK<br>
if I understand you correctly) is done on the projection level with<br>
the following filter :<br>
<a href="http://www.openrtk.org/Doxygen/classrtk_1_1DisplacedDetectorImageFilter.html" target="_blank">http://www.openrtk.org/Doxygen/classrtk_1_1DisplacedDetectorImageFilter.html</a><br>
<br>
You an check in rtkfdk.cxx how to use it.<br>
<br>
If you want to disable the ramp filter during reconstruction, a<br>
"simple" way would be to change the internal pipeline connections in<br>
rtk::FDKConeBeamReconstructionFilter. Replace in<br>
rtkFDKConeBeamReconstructionFilter.txx the line<br>
  m_BackProjectionFilter->SetInput( 1, m_RampFilter->GetOutput() );<br>
with<br>
  m_BackProjectionFilter->SetInput( 1, m_WeightFilter->GetOutput() );<br>
I hope this helps,<br>
Simon<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
On Sat, Dec 15, 2012 at 2:44 PM, Andy Shieh <<a href="mailto:hsieandy@gmail.com">hsieandy@gmail.com</a>> wrote:<br>
> Hi Simon,<br>
><br>
> I've been using rtkfdk to reconstruct and noticed that it handles half-beam<br>
> data.<br>
> In some situation, I want to use rtkbackprojections (for cases I don't want<br>
> to include ramp filtering). Does rtkbackprojections handle half-beam data<br>
> correctly, too? I gave it a try and it seems like rtkbackprojections would<br>
> result in some kind of voxel accumulation in the middle, which probably<br>
> means it doesn't handle half-beam data correctly like rtkfdk?<br>
> If that's the case, is it possible to make rtkbackproejctions support<br>
> half-beam data in the next patch (if it's not some massive work)? Thanks!<br>
><br>
> Cheers,<br>
> Andy<br>
</div></div></blockquote></div><br></div>