<div dir="ltr">Thanks Simon! It works great!<div style>Merry Christmas and Happy New Year!</div><div style><br></div><div style>Andy</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2012/12/27 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>
Indeed, you also need to comment the update of the ramp filter:<br>
m_RampFilter->Update();<br>
This time, I tested it and it works!<br>
<span class="HOEnZb"><font color="#888888">Simon<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
On Thu, Dec 20, 2012 at 2:03 PM, Andy Shieh <<a href="mailto:hsieandy@gmail.com">hsieandy@gmail.com</a>> wrote:<br>
> Hi Simon,<br>
><br>
> Thanks. I tried disabling ramp filter in<br>
> rtkFDKConeBeamReconstructionFilter.txx with the new line, but then gets a<br>
> runtime error as this:<br>
><br>
> ExceptionObject caught with writer->Update()<br>
> itk::ExceptionObject (0000000000AFE978)<br>
> Location: "void __cdecl itk::DataObject::PropagateRequestedRegion(void)<br>
> throw(cl<br>
> ass itk::InvalidRequestedRegionError)"<br>
> File: ..\..\..\Code\Common\itkDataObject.cxx<br>
> Line: 397<br>
> Description: Requested region is (at least partially) outside the largest<br>
> possible region.<br>
><br>
> which I don't get if I use the original line. I tried lowering the<br>
> reconstructed region, but it didn't help.<br>
> My projection size is 1024*768, 0.388 mm in both directions. It's half beam<br>
> data, detector being displaced by 14.8 cm.<br>
> The reconstructed region I always use is 512*100*512, spacing being<br>
> 0.88,2,0.88 mm.<br>
><br>
> Cheers,<br>
> Andy.<br>
><br>
><br>
><br>
> 2012/12/17 Simon Rit <<a href="mailto:simon.rit@creatis.insa-lyon.fr">simon.rit@creatis.insa-lyon.fr</a>><br>
>><br>
>> 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>
>><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>
>><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<br>
>> > half-beam<br>
>> > data.<br>
>> > In some situation, I want to use rtkbackprojections (for cases I don't<br>
>> > want<br>
>> > to include ramp filtering). Does rtkbackprojections handle half-beam<br>
>> > data<br>
>> > correctly, too? I gave it a try and it seems like rtkbackprojections<br>
>> > 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)?<br>
>> > Thanks!<br>
>> ><br>
>> > Cheers,<br>
>> > Andy<br>
><br>
><br>
</div></div></blockquote></div><br></div>