<div dir="ltr"><div><div><div><div>Hi,<br></div>Ok. Can you try rtkramp with only one projection? Then you shouldn't encounter this regions issue (that should be solved one day but that's a different story). <br></div>Can you also indicate which version of CUDA you're using?<br></div>FYI, we use ITK for the FFT of the kernel of the ramp filter and cufft for the FFT of the projections. So if they are not in the same format, you're in trouble and that might be the issue here... The weird part is that I didn't reproduce this problem on linux (with and without FFTW) but it might be a CUDA version difference. I'll give it another try tomorrow when you provide the CUDA version.<br></div>Thanks again for reporting the issue,<br>Simon<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 12, 2015 at 11:39 PM, Yang K Park <span dir="ltr"><<a href="mailto:theday79@gmail.com" target="_blank">theday79@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 link="#0563C1" vlink="#954F72" lang="EN-US"><div><p>Hi Simon,<u></u><u></u></p><p><u></u> <u></u></p><p>Thanks for the advice. The current RTK was compiled with FFTW to use the phase extraction feature from Amsterdam shroud.<u></u><u></u></p><p><u></u> <u></u></p><p>Regarding your suggestions,<u></u><u></u></p><p><u></u> <u></u></p><p style="margin-left:.5in"><u></u><span>1)<span style="font:7.0pt "Times New Roman"">      </span></span><u></u>rtkramp w/ or w/o CUDA:<u></u><u></u></p><p>(command line: “ rtkramp --path projPath --regexp .*.his --hardware cuda --output Ramp_wCUDA.mha” )<u></u><u></u></p><p><u></u> <u></u></p><p>Then, I have encountered following error.<u></u><u></u></p><p>“CudaFFTRampImageFilter<..> CudaFFTConvolutionImageFilter assumes that input requested and buffered regions are equal.”<u></u><u></u></p><p><u></u> <u></u></p><p><img src="cid:image001.png@01D11D71.1A70E010" height="88" width="646"><u></u><u></u></p><p><u></u> <u></u></p><p>Without CUDA, it was successfully done.<u></u><u></u></p><p><u></u> <u></u></p><p style="margin-left:.5in"><u></u><span>2)<span style="font:7.0pt "Times New Roman"">      </span></span><u></u>I’ve re-compiled <u>ITK w/o FFTW</u> and built RTK again. r<u>tkfdk worked fine.<u></u><u></u></u></p><p style="margin-left:.75in"><u></u><span>-<span style="font:7.0pt "Times New Roman"">          </span></span><u></u>I’ve also repeated “rtkramp+CUDA” with this rtk version and the error was same as above.<u></u><u></u></p><p style="margin-left:.75in"><u></u> <u></u></p><p>In conclusion, my <u>FFTW seems to affect CUDA-based fdk reconstruction only</u>.<u></u><u></u></p><p>But I do need FFTW for the phase extraction. What should I do?<u></u><u></u></p><p><u></u> <u></u></p><p>Thanks.<span class="HOEnZb"><font color="#888888"><u></u><u></u></font></span></p><span class="HOEnZb"><font color="#888888"><p><u></u> <u></u></p><p>Yang<u></u><u></u></p></font></span><div><div class="h5"><p><u></u> <u></u></p><p><u></u> <u></u></p><p>-----Original Message-----<br>From: <a href="mailto:simon.rit@gmail.com" target="_blank">simon.rit@gmail.com</a> [mailto:<a href="mailto:simon.rit@gmail.com" target="_blank">simon.rit@gmail.com</a>] On Behalf Of Simon Rit<br>Sent: Wednesday, November 11, 2015 4:16 PM<br>To: Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>><br>Cc: <a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a><br>Subject: Re: [Rtk-users] rtkfdk doesn't work correctly after a recent update</p><p><u></u> <u></u></p><p>Hi,<u></u><u></u></p><p>I didn't manage to reproduce the two problems on linux with the same ITK version.<u></u><u></u></p><p>For point 1, did you copy itkMacro.h by hand? Because I don't see how replacing <itkMacro.h> by "itkMacro.h" can help if itkMacro.h is not in the same folder as rtkCudaUtilities.hcu For point 2, no clue. Do you use FFTW? It looks like a ramp filter issue. Can you try rtkramp --path projPath --regexp .*.his, with and without --hardware cuda? Maybe you can provide your geometry file in case it's a specific geometric issue.<u></u><u></u></p><p>Simon<u></u><u></u></p><p><u></u> <u></u></p><p>On Wed, Nov 11, 2015 at 7:50 PM, Yang K Park <<a href="mailto:theday79@gmail.com" target="_blank"><span style="color:windowtext;text-decoration:none">theday79@gmail.com</span></a>> wrote:<u></u><u></u></p><p>> Hi all,<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> I’ve updated rtk with the up-to-date revision and have encountered a <u></u><u></u></p><p>> couple of problems.<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> 1.       Compile error: RTK couldn’t find itkMacro.h while it compiled<u></u><u></u></p><p>> rtkCudaUtilities.hcu. It was solved by a minor fetching:<u></u><u></u></p><p>><u></u> <u></u></p><p>> #include <itkMacro.h>   à #include "itkMacro.h" in rtkCudaUtilities.hcu<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> 2.       Output of rtkfdk with CUDA was not the same as before. It was very<u></u><u></u></p><p>> blurry as seen in the attached image.<u></u><u></u></p><p>><u></u> <u></u></p><p>> a.       With “FDK with CPU” or “SART with CUDA”, there was no problem.<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> I’m using VS2013 with win7, itk4.8.1.<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> The command line I used was :<u></u><u></u></p><p>><u></u> <u></u></p><p>> rtkfdk --geometry 1111.xml --path projPath --verbos --regexp .*.his <u></u><u></u></p><p>> --output output.mha --spacing 1,1,1 --dimension 400,200,400 --hardware <u></u><u></u></p><p>> cuda --pad 1 --hann 5 --hannY 5<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> Does anyone have any idea about this?<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> Thanks.<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> Yang<u></u><u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>><u></u> <u></u></p><p>> _______________________________________________<u></u><u></u></p><p>> Rtk-users mailing list<u></u><u></u></p><p>> <a href="mailto:Rtk-users@public.kitware.com" target="_blank"><span style="color:windowtext;text-decoration:none">Rtk-users@public.kitware.com</span></a><u></u><u></u></p><p>> <a href="http://public.kitware.com/mailman/listinfo/rtk-users" target="_blank"><span style="color:windowtext;text-decoration:none">http://public.kitware.com/mailman/listinfo/rtk-users</span></a><u></u><u></u></p><p>><u></u> <u></u></p></div></div></div></div></blockquote></div><br></div>