<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 color="#0000ff" face="Consolas" size="2"><font color="#0000ff" face="Consolas" size="2"><font color="#0000ff" face="Consolas" size="2">typedef</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"> itk::</font></font><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2">ImageFileWriter</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"><</font></font><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2">OutputImageType</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">> </font></font><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" 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 color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2">WriterType</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">::</font></font><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2">Pointer</font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2"> writer = </font></font><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" size="2"><font color="#2b91af" face="Consolas" 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 color="#a31515" face="Consolas" size="2"><font color="#a31515" face="Consolas" size="2"><font color="#a31515" face="Consolas" 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <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 class="HOEnZb"><font color="#888888"><br>
    <br>
    Cyril</font></span><div><div class="h5"><br>
    <br>
    <div>On 06/15/2015 05:19 PM, safdary mohsen
      wrote:<br>
    </div>
    </div></div><blockquote type="cite"><div><div class="h5">
      <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 lang="EN-US" vlink="purple" link="blue">
              <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>