Hi Micah,<br>You did not do a mistake, you have successfully hit a bug. The truth is that I haven't used the divisions option in a very long time... This should now be fixed<br> <a href="https://github.com/SimonRit/RTK/commit/10f8e1b0a5f6b69c943d04513135856968e3e62b">https://github.com/SimonRit/RTK/commit/10f8e1b0a5f6b69c943d04513135856968e3e62b</a><br>
and I have added an automated test to avoid future inconveniences<br> <a href="https://github.com/SimonRit/RTK/commit/0932e38ad73d8e7214adda5103ced38c7b2754d0">https://github.com/SimonRit/RTK/commit/0932e38ad73d8e7214adda5103ced38c7b2754d0</a><br>
Thanks for the report and keep us a posted. By the way, I have not done a lot of work with tif inputs so you might have to modify the tif conversion from raw to attenuation<br> <a href="https://github.com/SimonRit/RTK/blob/master/code/rtkTiffLookupTableImageFilter.h">https://github.com/SimonRit/RTK/blob/master/code/rtkTiffLookupTableImageFilter.h</a><br>
Good luck!<br>Simon<br><br><div class="gmail_quote">On Wed, Oct 24, 2012 at 1:57 AM, M Miller <span dir="ltr"><<a href="mailto:croow@yahoo.com" target="_blank">croow@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font:inherit" valign="top"><div>I'm new to RTK and I'd like to reconstruct some projections into a volume. So far I've only been able to get an "InvalidRequestedRegionError", but I'd like to describe my steps with the hope that someone can point out what I've done wrong.</div>
<div><br>I built RTK using VS2010 x64, ITK 3.20.1, and CUDA 4.2. The sample programs seem to work fine as I can successfully walk through the FDK SheppLogan example.</div>
<div><br>My projections are 3143 16-bit tif files, each 2000x2000 pixels in size, covering 360deg of rotation. The center of rotation is 65.93mm from the xray source, and the detector is 870.86mm from the source. The detector is 2000x2000 pixels in size, with each pixel 0.2mmx0.2mm in size.</div>
<div> </div>
<div>I create a geometry file:<br>"rtksimulatedgeometry.exe --output=geo.xml --nproj=3143 --arc=360 --sdd=870.86 --sid=65.93"</div>
<div> </div>
<div>Attempt to reconstruct:<br>"rtkfdk.exe --verbose --lowmem --geometry d:\ct_projections\geo.xml --path=d:\ct_projections\test\ --regexp=tif$ --output=d:\ct_projections\rtkfdk.mha --divisions=2"</div>
<div> </div>
<div>Output:<br>>Regular expression matches 3143 file(s)<br>>Reading geometry information from d:\ct_projections\geo.xml...<br>>Reconstructing and writing... ExceptionObject caught with writer->Update()<br>><br>
>itk::InvalidRequestedRegionError (0000000000009FED10)<br>>Location: "void ___cdecl itk::DataObject::PropagateRequestedRegion(void) throw(class itk::InvalidRequestedRegionError)"<br>>File: ..\..\..\InsightToolkit-3.20.1\Code\Common\itkDataObject.cxx<br>
>Line: 397<br>>Description: Requested region is (at least partially) outside the largest possible region.</div>
<div> </div>
<div>I have replaced the original macro "TRY_AND_EXIT_ON_ITK_EXCEPTION( reader->Update() )" with just "reader->Update()", otherwise there is no debug information.</div>
<div><br>The "lowmem" switch is required to keep ITK from crashing, which I think is an unrelated problem. I use "divisions=2" because plastimatch would crash for me without specifying something similar. Setting the "hardware" switch to either cpu or cuda produces the same error. The client machine has 32 processors, 128GB RAM, Tesla C2070, and is running Win7 x64.</div>
<div> </div>
<div>Is there any extra information necessary or anything that needs clarification?<br>How can I use RTK to reconstruct a volume, if possible?</div>
<div> </div>
<div>Thanks<br>Micah</div></td></tr></tbody></table><br>_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@openrtk.org">Rtk-users@openrtk.org</a><br>
<a href="http://public.kitware.com/cgi-bin/mailman/listinfo/rtk-users" target="_blank">http://public.kitware.com/cgi-bin/mailman/listinfo/rtk-users</a><br>
<br></blockquote></div><br>