<div dir="ltr"><div>Hi,<br>So you checked that projections are read but did you check what's read with a viewer, e.g. <a href="http://vv.creatis.insa-lyon.fr">vv</a>? Can you also tell us what's the offset and the spacing in the reader.mha file?<br></div>Simon<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 25, 2015 at 7:27 AM, safdary mohsen <span dir="ltr"><<a href="mailto:ieee.safdary@gmail.com" target="_blank">ieee.safdary@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>Hi </div><div><br></div><div> what is genuine rtkfdk.exe ?</div><div><br></div><div>source code is attached to mail and command line is:</div><div>rtkfdk --geometry geometry.xml --regexp .*.hnd --path Scan0 --output slice30.mha --verbose --spacing 0.25,0.25,0.25 --dimension 1024,1,1024 --origin -127.875,30,-127.875 --hardware cuda<br></div><div><br></div><div><font face="Consolas">
<font color="#0000ff" size="2"><font color="#0000ff" size="2"><font color="#0000ff" size="2"></font></font></font></font><p><br></p><font face="Consolas" size="2"><p><br></p><font face="Consolas" size="2"><p><br></p><p><br></p>
</font></font></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 24, 2015 at 5:51 PM, Yang-Kyun 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"><p dir="ltr">Hi Safdary,</p>
<p dir="ltr">So you might have tried the genuine rtkfdk.exe for your recon. Would you give me your command line you used?</p><span><font color="#888888">
<p dir="ltr">Yang</p></font></span><div><div>
<div class="gmail_quote">On Jun 24, 2015 3:48 AM, "safdary mohsen" <<a href="mailto:ieee.safdary@gmail.com" target="_blank">ieee.safdary@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><div dir="ltr"><div>Hi </div><div> </div><div>I cheched projectionreader and i abtained that all projection are readed.</div><div>( by using:</div><div><p><font face="Consolas" color="#0000ff" size="2"><font face="Consolas" color="#0000ff" size="2"><font face="Consolas" color="#0000ff" size="2">typedef</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"> itk::</font></font><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2">ImageFileWriter</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"><</font></font><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2">OutputImageType</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">> </font></font><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2">WriterType</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">;</font></font></p><font face="Consolas" size="2"><font face="Consolas" size="2">
</font></font><p><font face="Consolas" size="2"><font face="Consolas" size="2"> </font></font><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2">WriterType</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">::</font></font><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2">Pointer</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"> writer = </font></font><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2"><font face="Consolas" color="#2b91af" size="2">WriterType</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">::New();</font></font></p><font face="Consolas" size="2"><font face="Consolas" size="2">
</font></font><p><font face="Consolas" size="2"><font face="Consolas" size="2"> writer->SetFileName( </font></font><font face="Consolas" color="#a31515" size="2"><font face="Consolas" color="#a31515" size="2"><font face="Consolas" color="#a31515" size="2">"reader.mha"</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"> );</font></font></p></div><font face="Consolas" size="2"><font face="Consolas" size="2"><p> writer->SetInput( reader->GetOutput() );</p><div>
</div><p> writer->Update();</p></font></font><div><font face="Consolas" size="2"><font face="Consolas" size="2"></font></font>)</div><div><br></div><div>and</div><div><br></div><div>geometry is work for cpu mode means outputs are ok .</div><div><br></div><div><br></div><div>I am confused for solving this problem!!!!!!!!!!!!!!!!!</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 16, 2015 at 11:09 AM, Cyril Mory <span dir="ltr"><<a href="mailto:cyril.mory@creatis.insa-lyon.fr" target="_blank">cyril.mory@creatis.insa-lyon.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div bgcolor="#FFFFFF" text="#000000">
Hi Safdary,<br>
<br>
Just to be sure: have you tried reconstructing in a larger volume
(something really huge), just to make sure you are not trying to
reconstruct some portion of space where, indeed, there is nothing ?
<br>
You can try that easily, just replace your spacing, dimension and
origin with<br>
--spacing 16,16,16 --dimension 256,256,256 --origin
-2040,-2040,-2040<br>
in your rtkfdk command line.<span><font color="#888888"><br>
<br>
Cyril</font></span><div><div><br>
<br>
<div>On 06/15/2015 05:19 PM, safdary mohsen
wrote:<br>
</div>
</div></div><blockquote type="cite"><div><div>
<div dir="ltr">
<div class="MsoNormal"><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt">Hi
Yang</span></div>
<div class="MsoNormal"><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"><br>
</span></div>
<div class="MsoNormal"><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"><br>
</span></div>
<div class="MsoNormal"><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.66px">
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">I
am use rtkfdk.exe but the results are same as my
code, </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">means
all pixels are zero!!!</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"><u><br>
</u></span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">regards</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">safdari</span></p>
</span></span></div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Mon, Jun 15, 2015 at 7:10 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:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div vlink="purple" link="blue" lang="EN-US">
<div>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Hi
Safdary, </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Again,
please cc <a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a> when you
send your email, since other users (much better
experts than me!) can give us any suggestions.</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Thanks.</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Yang</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<div>
<div style="border-width:1pt medium medium;border-style:solid none none;border-color:rgb(225,225,225) currentColor currentColor;padding:3pt 0in 0in">
<p class="MsoNormal"><b><span style="font-family:"Calibri",sans-serif;font-size:11pt">From:</span></b><span style="font-family:"Calibri",sans-serif;font-size:11pt"> Yang
K Park [mailto:<a href="mailto:theday79@gmail.com" target="_blank">theday79@gmail.com</a>] <br>
<b>Sent:</b> Monday, June 15, 2015 10:37 AM<br>
<b>To:</b> 'safdary mohsen'<br>
<b>Cc:</b> <a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a><br>
<b>Subject:</b> RE: [Rtk-users] error on build
rtk in visual studio 2012 when rtk use cuda</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Hi
Safdary,</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">I
recommend you to try rtkfdk.exe first (with cuda
option) to see if it works for you.</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">The
example of the command line would be:</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:black;font-family:"Calibri",sans-serif;font-size:11pt">rtkfdk
--geometry geom.xml --path proj_dir --regexp
.*.his --output test_output.mha --verbos
--spacing 1,1,1 --dimension 512,10,512 --origin
-127.875,29.5,-127.875 </span><span style="color:red;font-family:"Calibri",sans-serif;font-size:11pt">--hardware
cuda</span></p>
<p class="MsoNormal"><span style="color:red;font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">If
it works for your image, then you can mimic the
original rtkfdk code to make your own one.</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Good
luck!</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Yang</span><span style="color:red;font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p class="MsoNormal"><span style="color:black;font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:black;font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><b><span style="font-family:"Calibri",sans-serif;font-size:11pt">From:</span></b><span style="font-family:"Calibri",sans-serif;font-size:11pt">
safdary mohsen [<a href="mailto:ieee.safdary@gmail.com" target="_blank">mailto:ieee.safdary@gmail.com</a>]
<br>
<b>Sent:</b> Saturday, June 13, 2015 11:58 PM<br>
<b>To:</b> Yang K Park<br>
<b>Subject:</b> Re: [Rtk-users] error on build
rtk in visual studio 2012 when rtk use cuda</span></p>
<p class="MsoNormal"> </p>
<div>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">Hi
Yang</span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"> </span></p>
<p class="MsoNormal"><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">I
have tried to adjust the window/level to see
something inside of the image but all pixels
set to zero.</span></p>
<div>
<p class="MsoNormal"><u><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt">My
code is:</span></u></p>
</div>
<div>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#define</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:rgb(111,0,138);font-family:Consolas;font-size:10pt">USE_CUDA</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> main(</span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
, </span><span style="color:blue;font-family:Consolas;font-size:10pt">char</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> **)</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">{</span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">////
variable</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">const</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
Dimension = 3;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">float</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">OutputPixelType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#ifdef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:rgb(111,0,138);font-family:Consolas;font-size:10pt">USE_CUDA</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> itk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">CudaImage</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">< </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">OutputPixelType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">,
Dimension > </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">OutputImageType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#else</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
itk::Image< OutputPixelType, Dimension
> OutputImageType;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#endif</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">////////////////////////////////////////////////////////////////////////////////////////////////////</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
Generate file names of projections</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">itk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">RegularExpressionSeriesFileNames</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Pointer</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> names
= itk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">RegularExpressionSeriesFileNames</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::New();</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">names->SetDirectory(std::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">string</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">(</span><span style="color:rgb(163,21,21);font-family:Consolas;font-size:10pt">"d://RTK/build/bin/Release/Scan1/"</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">));</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">names->SetRegularExpression(std::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">string</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">(</span><span style="color:rgb(163,21,21);font-family:Consolas;font-size:10pt">"Proj_.*.hnd"</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">));</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:green;font-family:Consolas;font-size:10pt">////////////////////////////////////////////////////////////////////////////////////////////////////</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">////
geometry</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
Defines the RTK geometry object</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> rtk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ThreeDCircularProjectionGeometry</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">GeometryType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">GeometryType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Pointer</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
geometry = </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">GeometryType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::New();</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
Projection matrices</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
numberOfProjections = 86;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
firstAngle = 0;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
angularArc = 360;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
sid = 1000; </span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
sdd = 1200; </span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> isox =
0; </span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> isoy =
0 ;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">for</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">(</span><span style="color:blue;font-family:Consolas;font-size:10pt">unsigned</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:blue;font-family:Consolas;font-size:10pt">int</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
noProj=0; noProj<numberOfProjections;
noProj++)</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">{</span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">double</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> angle
= (</span><span style="color:blue;font-family:Consolas;font-size:10pt">float</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">)firstAngle
+ (</span><span style="color:blue;font-family:Consolas;font-size:10pt">float</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">)noProj
* angularArc / (</span><span style="color:blue;font-family:Consolas;font-size:10pt">float</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">)numberOfProjections;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">geometry->AddProjection(sid,</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">sdd,</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">angle,</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">isox,</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">isoy);</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">}</span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">/////////////////////////////////////////////////////////////////////////////////////////////////////</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
Projections reader</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> rtk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ProjectionsReader</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">< </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">OutputImageType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> > </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ReaderType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ReaderType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Pointer</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> reader
= </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ReaderType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::New();</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">reader->SetFileNames(names->GetFileNames());</span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
Create reconstructed image
////////////////////////////////////</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> rtk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSource</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">< </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">OutputImageType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> > </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSourceType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSourceType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Pointer</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
ImageSource = </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSourceType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::New();</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSourceType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">PointType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
origin;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSourceType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">SizeType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
sizeOutput;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ConstantImageSourceType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">SpacingType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
spacing;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">origin[0]
= -127.75;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">origin[1]
= -127.75;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">origin[2]
= -127.75;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">sizeOutput[0]
= 512;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">sizeOutput[1]
= 512;</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">sizeOutput[2]
= 512;</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">spacing[0]
= .5;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">spacing[1]
= .5;</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">spacing[2]
= .5;</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">ImageSource->SetOrigin(
origin );</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">ImageSource->SetSpacing(
spacing );</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">ImageSource->SetSize(
sizeOutput );</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">ImageSource->SetConstant(
0. ); </span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
getchar();
//////////////////////////////////////</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">std::cout
<< </span><span style="color:rgb(163,21,21);font-family:Consolas;font-size:10pt">"Performing
reconstruction"</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
<< std::endl;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
FDK reconstruction filtering</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
FDK reconstruction filtering</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#ifdef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:rgb(111,0,138);font-family:Consolas;font-size:10pt">USE_CUDA</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> rtk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">CudaFDKConeBeamReconstructionFilter</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">FDKType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#elif</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
USE_OPENCL</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
rtk::OpenCLFDKConeBeamReconstructionFilter
FDKType;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#else</span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
rtk::FDKConeBeamReconstructionFilter<
OutputImageType > FDKType;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">#endif</span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">FDKType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Pointer</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
feldkamp = </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">FDKType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::New();</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">feldkamp->SetInput(
0, ImageSource->GetOutput() );</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">feldkamp->SetInput(
1, reader->GetOutput());</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">feldkamp->SetGeometry(
geometry);</span></p>
<p><span style="color:rgb(111,0,138);font-family:Consolas;font-size:10pt">TRY_AND_EXIT_ON_ITK_EXCEPTION</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">(
feldkamp->Update() );</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">std::cout
<< </span><span style="color:rgb(163,21,21);font-family:Consolas;font-size:10pt">"Writing
output image"</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">
<< std::endl;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:green;font-family:Consolas;font-size:10pt">//
Writer</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> itk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Image</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">< </span><span style="color:blue;font-family:Consolas;font-size:10pt">float</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">, 3
> </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ImageType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">typedef</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> itk::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ImageFileWriter</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"><</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">ImageType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">> </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">WriterType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">WriterType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::</span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">Pointer</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> writer
= </span><span style="color:rgb(43,145,175);font-family:Consolas;font-size:10pt">WriterType</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">::New();</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">writer->SetFileName(
</span><span style="color:rgb(163,21,21);font-family:Consolas;font-size:10pt">"output11.mha"</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> );</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">writer->SetInput(
feldkamp->GetOutput() );</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">writer->Update();</span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> </span></p>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">getchar();</span></p>
<p><span style="color:blue;font-family:Consolas;font-size:10pt">return</span><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt"> 0;</span><span style="color:rgb(31,73,125);font-family:"Calibri",sans-serif;font-size:11pt"></span></p>
</div>
</div>
<div>
<p><span style="color:rgb(31,73,125);font-family:Consolas;font-size:10pt">}</span></p>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
<br>
<fieldset></fieldset>
<br>
</div></div><span><pre>_______________________________________________
Rtk-users mailing list
<a href="mailto:Rtk-users@public.kitware.com" target="_blank">Rtk-users@public.kitware.com</a>
<a href="http://public.kitware.com/mailman/listinfo/rtk-users" target="_blank">http://public.kitware.com/mailman/listinfo/rtk-users</a>
</pre>
</span></blockquote>
<br>
</div>
</blockquote></div><br></div>
</blockquote></div>
</div></div></blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@public.kitware.com">Rtk-users@public.kitware.com</a><br>
<a href="http://public.kitware.com/mailman/listinfo/rtk-users" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/rtk-users</a><br>
<br></blockquote></div><br></div>