<div dir="ltr"><div>Hi Stefanie,</div><div>In the example, the projections are simulated and not read in from a file. Anyhow, using the mha file format seems to be a good way too proceed and what you've done sounds good (the rei part).</div><div>The geometry is a bit tricky to master in RTK (as in any recon package I assume). The origin you've chosen does not seem right, if you check ITK's documentation, itk::Image::m_Origin is the coordinate of the first voxel and RTK uses the coordinate (0,0,0) (the origin of the coordinate system, not to be mixed with the origin of the image...). So if you want to center your image around the center of rotation (a typical choice), I would do here</div><div>origin.Fill( -0.5*(128-1)*0.2 ) such that the first pixel has coordinate -12.7 and the last pixel has coordinate 12.7 (in each dimension). </div><div>Something else is wrong looking at your image. If you share the projections and the code, we can probably help you to solve the other problem.</div><div>Cheers,</div><div>Simon<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 7, 2020 at 5:41 PM Kaser, Stefanie <<a href="mailto:Stefanie.Kaser@oeaw.ac.at">Stefanie.Kaser@oeaw.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">




<div dir="ltr">
<div id="gmail-m_-2066591298068327360divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p>Hi,</p>
<p><br>
</p>
<p>I am currently trying to get started with RTK by modifying the FirstReconstruction example. In fact, I tried to replace the given projections of a sphere by some projections I created myself with Geant4. However, the result is not really convincing and I
 am not sure if I messed with the geometry or if I got something else wrong.</p>
<p><br>
</p>
<p>I have 90 projections (128x128 pixels, 2 deg rotation each time); from source to detector it's 46cm and 30cm from source to the object (simple small cylinder) centre, which lies in the centre of rotation:</p>
<p><br>
</p>
<p></p>
<p>  unsigned int numberOfProjections = 90;</p>
<p>  double firstAngle = 0.; </p>
<p>  double angularArc = 178.; </p>
<p>  unsigned int sid = 300.;</p>
<p>  unsigned int sdd = 460.; </p>
<p><br>
</p>
<p>To read in the file, I used the ITK ImageFileReader, which works well. <span style="font-size:12pt">I am handing in my projections stacked in .mha format because it seemed that the sphere projections were generated and handed to the FDK reconstruction algorithm in
 a similar way </span><span style="font-size:12pt">(is that correct?)</span><span style="font-size:12pt">. My projection values are double type and are of a range between 0 (air) and 1.</span></p>
<p><br>
</p>
<p>For the FDK reconstruction I am choosing for my constant image source:</p>
<p>  sizeOutput.Fill(128); // output should be 3D image 128x128x128 pixels</p>
<p>  origin.Fill(0.); // I'm not entirely sure about this parameter? If my object is placed in the isocenter, this should work, right?</p>
<p>  spacing.Fill(0.2); // pixel length in each direction is 0.2mm</p>
<p><br>
</p>
<p>For the reconstruction itself, I didn't really change the original code but instead of 'rei', I handed my projections as input and then wrote the result to an output file. </p>
<p>I don't get any error when I run the program but the FDK gives some strange result (I am adding one screenshot of the first projection and the first slice of the reconstructed volume to my mail). </p>
<p>Maybe someone is familiar with this problem or sees what I am doing wrong? Thank you very much in advance!</p>
<p><br>
</p>
<p>Best regards,</p>
<p>Stefanie</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<p></p>
</div>
</div>

_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@public.kitware.com" target="_blank">Rtk-users@public.kitware.com</a><br>
<a href="https://public.kitware.com/mailman/listinfo/rtk-users" rel="noreferrer" target="_blank">https://public.kitware.com/mailman/listinfo/rtk-users</a><br>
</blockquote></div>