<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Interesting, then maybe it is possible to standardise it.<br>I would prefer to have the normalisation in RTK if possible.<br>A value named KVNormChamber is read on this line: <a href="https://github.com/SimonRit/RTK/blob/master/src/rtkXimImageIO.cxx#L73">https://github.com/SimonRit/RTK/blob/master/src/rtkXimImageIO.cxx#L73</a></div><div dir="ltr">But if I'm reading the section from the Thesis you linked right, I guess we'd also need the I_0 value.<br>I think this must come from a file in the Calibrations folder, so it'll have to be an optional feature and require a CLI flag pointing to the right file (depending on the bow-tie type used if any).</div><div>I'm guessing the "KVNormChamber" value corresponds to NC_I, but is NC_I_0 then "KVNormChamber" in the calibration projection, or is it maybe the kV value from the scan.xml file ?</div><div><br></div><div><br>P.S. Shameless promotion of my own project:</div><div dir="ltr">I wouldn't recommend getting CbctRecon just to normalise projections. However, for the a priori scatter correction it will make your life easier:<br><ul><li>Take a look at Testing/test_end_to_end.cpp for a somewhat short script doing the full correction (If you want to script it) <a href="https://gitlab.com/agravgaard/cbctrecon/-/blob/Develop/Testing/test_end_to_end.cpp#L176">https://gitlab.com/agravgaard/cbctrecon/-/blob/Develop/Testing/test_end_to_end.cpp#L176</a></li><li>Configure and build the whole thing if you want a user interface (UI works best on Windows) for the scatter correction: <a href="https://gitlab.com/agravgaard/cbctrecon/-/tree/Develop/#cbctrecon">https://gitlab.com/agravgaard/cbctrecon/-/tree/Develop/#cbctrecon</a></li><li>It's optimised a lot and takes about 70 seconds for a full scatter correction (also including file-reading) on my Linux machine.</li><li>Although expanded from the original implementation by Yang et al., it avoids the "Correction Factor"</li><li>TBI: Two ongoing projects to automise the process, a minimal UI project and a clinical integration project (reconstructing as projections are captured with a frame-grabber).</li></ul></div><div dir="ltr"><br></div><div dir="ltr"><br>P.P.S. There's quite a different between C and C++, in particular C++ 17 (or 20), below is a few of my favourite ressources for learning C++ (although I'm mostly self-taught by looking at other peoples code)<br><ul><li>Jonathan Boccara's FluentC++ blog: <a href="https://www.fluentcpp.com/posts/">https://www.fluentcpp.com/posts/</a> For writing readable (expressive) C++ and understanding what's in the standard library</li><li>Jason Turner's Weekly C++ series on youtube: <a href="https://www.youtube.com/playlist?list=PLs3KjaCtOwSZ2tbuV1hx8Xz-rFZTan2J1">https://www.youtube.com/playlist?list=PLs3KjaCtOwSZ2tbuV1hx8Xz-rFZTan2J1</a></li><li>Bartek Filipek's blog: <a href="https://www.bfilipek.com/">https://www.bfilipek.com/</a> For C++17, lambdas, and in general modern C++</li><li>Craig Scott's book: <a href="https://crascit.com/professional-cmake/">https://crascit.com/professional-cmake/</a> The single best ressource for all your CMake needs.</li><li>Rainer Grimm's blog: <a href="http://www.modernescpp.com/index.php">http://www.modernescpp.com/index.php</a> For understanding threading, atomics, parallelism and functional programming in C++</li><li>Matt Godbolt's: <a href="https://godbolt.org/">https://godbolt.org/</a> Online compiler with all the compilers for quick and easy test of code snippets and compiler flags (you might also learn some assembly if you use it as often as I do)</li></ul></div><div dir="ltr"><br>/Andreas<br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><p style="margin-bottom:0.0001pt"><span lang="EN-US" style="color:rgb(31,73,125)"><span lang="EN-US" style="font-size:11pt;line-height:16.8667px;font-family:Calibri,sans-serif">__________________________________</span></span></p><p style="margin-bottom:0.0001pt"><span lang="EN-US" style="color:rgb(31,73,125)"><span lang="EN-US" style="font-size:11pt;line-height:16.8667px;font-family:Calibri,sans-serif">Andreas Gravgaard Andersen</span></span></p><p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Danish Center for Particle Therapy, </span></p><p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Aarhus University Hospital<u></u><u></u></span></p><p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Palle Juul-Jensens Blvd. 99,</span></p><p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">8200, Aarhus</span></p><div><p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Mail: <a href="mailto:agravgaard@protonmail.com" target="_blank">agravgaard@protonmail.com</a></span></p><p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt">Cell: </span><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"> </span><a value="+4523382411" style="font-family:Calibri,sans-serif;font-size:11pt;color:rgb(17,85,204)">+45 3165 8140</a></p></div></div></div></div></div></div></div></div></div></div><br></div></div></div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 27 May 2020 at 06:00, Brani Rusanov <<a href="mailto:brani_rusanov@hotmail.com">brani_rusanov@hotmail.com</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 style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div>
<div id="gmail-m_5529535081606338027appendonsend"></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi Andreas,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Your input is invaluable, thanks for your suggestions! You are correct, I am attempting to apply the a priori scatter correction. Is the normalisation procedure you mentioned regarding the calibration-chamber similar to what has been mentioned on page 63-64
of this thesis (<a href="https://pdfs.semanticscholar.org/d90c/d77b61d10d2c24df0686bc0be1f38424d67d.pdf" id="gmail-m_5529535081606338027LPlnk674520" target="_blank">https://pdfs.semanticscholar.org/d90c/d77b61d10d2c24df0686bc0be1f38424d67d.pdf</a>), specifically Norm Chamber (NC) values located in
raw projection header?<br>
<br>
I will admit, I am a complete novice when it comes to C coding and using RTK. For correctly normalising the Varian projections do you recommend installing your "cbctrecon" package? If so, which function performs the Varian normalisation? Or can the normalisation
be done by scaling the projection pixel values manually? My plan is to manually scale the CB projections and perform a reconstruction until the HU looks right. A more robust method is preferable, however. </div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Your advice is greatly appreciated,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Brani</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5529535081606338027divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Andreas Andersen <<a href="mailto:andreasga22@gmail.com" target="_blank">andreasga22@gmail.com</a>><br>
<b>Sent:</b> Tuesday, 26 May 2020 7:15 PM<br>
<b>To:</b> Brani Rusanov <<a href="mailto:brani_rusanov@hotmail.com" target="_blank">brani_rusanov@hotmail.com</a>><br>
<b>Cc:</b> <a href="mailto:rtk-users@public.kitware.com" target="_blank">rtk-users@public.kitware.com</a> <<a href="mailto:rtk-users@public.kitware.com" target="_blank">rtk-users@public.kitware.com</a>><br>
<b>Subject:</b> Re: [Rtk-users] Unable to Forward Project CT Volume With Correct Orientation</font>
<div> </div>
</div>
<div>
<div dir="ltr">I'm happy to hear it worked!
<div><br>
I don't know the complete Varian pipeline for correction, there should in addition to subtracting the bow-tie be something like a calibration-chamber-based correction, but I have not figured out how to do this (let me know if you do)
<div>RTK does not have filters to take care of this (yet). It would be possible to do, but you'd have to know which image in the calibration folder is the correct one. (There may be several and I think they are named by the local Varian engineers - not standardised)</div>
<div>I do a simple subtraction in my own OpenCL implementation, but looping through projections, extracting 2D from 3D, and using the ITK subtraction image filter would probably be simpler.</div>
<div>
<div>(<i>If I should find documentation from Varian describing this process with some standard, I'll write it into the reader for Xim images. However, I'm afraid their own solution might specific to each gantry or engineering team</i>)</div>
<div><br>
</div>
<div>It could sound like you want to do the "a priori scatter correction", I've implemented this for Varian based projections (hnd and xim), see:
<a href="https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fgitlab.com%2Fagravgaard%2Fcbctrecon&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580829813&sdata=Hiyv%2BwFoqjeHcqAH3jIG%2B5aiyqKnurLJ0qQ6ax4N3ug%3D&reserved=0" target="_blank">
gitlab.com/agravgaard/cbctrecon</a> (The Develop branch should be preferred at the moment)</div>
<div>The correction factor from the original Elekta implementation doesn't work without a normalisation as mentioned above, so I'm doing some tricks to ensure correctness of projection intensity for the subtraction, but I hope to implement a solution/normalisation
based on numbers from the machine instead.</div>
<div><br>
</div>
<div>/Andreas<br clear="all">
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<p style="margin-bottom:0.0001pt"><span lang="EN-US" style="color:rgb(31,73,125)"><span lang="EN-US" style="font-size:11pt;line-height:16.8667px;font-family:Calibri,sans-serif">__________________________________</span></span></p>
<p style="margin-bottom:0.0001pt"><span lang="EN-US" style="color:rgb(31,73,125)"><span lang="EN-US" style="font-size:11pt;line-height:16.8667px;font-family:Calibri,sans-serif">Andreas Gravgaard Andersen</span></span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Danish Center for Particle Therapy, </span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Aarhus University Hospital<u></u><u></u></span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Palle Juul-Jensens Blvd. 99,</span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">8200, Aarhus</span></p>
<div>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Mail: <a href="mailto:agravgaard@protonmail.com" target="_blank">agravgaard@protonmail.com</a></span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt">Cell: </span><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"> </span><a value="+4523382411" style="font-family:Calibri,sans-serif;font-size:11pt;color:rgb(17,85,204)">+45
3165 8140</a></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
</div>
</div>
<br>
<div>
<div dir="ltr">On Tue, 26 May 2020 at 12:37, Brani Rusanov <<a href="mailto:brani_rusanov@hotmail.com" target="_blank">brani_rusanov@hotmail.com</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi Andreas,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
You are correct, I had failed to flip the y and z coordinates of the offset parameter in the header. I had thought that the transformmatrix would take care of this. Now the simulated projection looks as desired. Thank you!</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
You mentioned the possibility of applying pre-processing on the CB projections for a bow-tie filter. Do you know what kind of pre-processing is available within RTK for bow-tie filters? I ask because I plan on subtracting the measured CBCT projections from
the simulated CT-registered projections. However, the measured CBCT projections have the characteristic attenuation pattern on the left and right side of the image corresponding to a bow-tie filter. Can RTK remove these two regions of greater attenuation? </div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Also, would you happen to know if rtkprojections, when converting the raw projection float values to attenuation, is compatible with Varian raw projection data? I ask this because the studies I have read that perform this procedure use Elekta projection data.
The output attenuation values for my measured projection data in air region is about ~2 mm-1 after conversion with rtkprojections. Whereas, I previously converted the HU of the registered CT volume to attenuation using the formula: mu = (HU + 1024)/2^16 which
shows air to be ~0.019mm-1 once forward projected.</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you kindly for your advice,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Brani</div>
<div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5529535081606338027x_gmail-m_7439177778697677005divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Andreas Andersen <<a href="mailto:andreasga22@gmail.com" target="_blank">andreasga22@gmail.com</a>><br>
<b>Sent:</b> Tuesday, 26 May 2020 4:51 PM<br>
<b>To:</b> Brani Rusanov <<a href="mailto:Brani_Rusanov@hotmail.com" target="_blank">Brani_Rusanov@hotmail.com</a>><br>
<b>Cc:</b> <a href="mailto:rtk-users@public.kitware.com" target="_blank">rtk-users@public.kitware.com</a> <<a href="mailto:rtk-users@public.kitware.com" target="_blank">rtk-users@public.kitware.com</a>><br>
<b>Subject:</b> Re: [Rtk-users] Unable to Forward Project CT Volume With Correct Orientation</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>What do you more precisely mean by "The output projections do not look like the raw projections I have from the CBCT"</div>
<div>Are they flipped, are they all black, are they much brighter or darker, or do they not align?<br>
<br>
</div>
<div>If alignment: Consider if you need to sort the projections by gantry angle before visual comparison, the Xim reader just reads sequentially by the given regex, so a ProBeam CB (with two imagers) may interleave projections from different angles. This is
unlikely to be a problem for a gantry with only one imager.</div>
<div><br>
</div>
<div>If flipped or all black: Your transformation matrix looks correct to me, so I think the problem may be elsewhere.</div>
<div dir="ltr">
<div>It's always a mess to figure out. Be careful, that when you flip axes manually like that that you also remember to flip offsets (origin), spacing, etc. accordingly.</div>
<div>In particular orientation (directions), I can see it says RAI for "AnatomicalOrientation", but I think you'll need the "Orientation" tag with direction cosines if you want to specify the directions manually.</div>
<div><a href="https://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fitk.org%2FWiki%2FITK%2FMetaIO%2FDocumentation%23MetaObject_Tags&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580839808&sdata=v2SPprbq57AoHf3VyamdMKStgCNsVK0HozMmmJbI3dM%3D&reserved=0" target="_blank">https://itk.org/Wiki/ITK/MetaIO/Documentation#MetaObject_Tags</a></div>
<div><br>
</div>
<div>All this being said, I would suggest to instead use ITK's readers and the built-in filters for image transformation if that's even necessary, to avoid the hell of trying to keep track of all these values manually.</div>
<div><br>
</div>
<div>If intensity is darker or brighter: Check if you have preprocessed the CB projections correctly for bow-tie filter and the like.</div>
<div><br>
</div>
<div>I hope this gives a hint to where the problem might be.</div>
<div><br>
</div>
<div>/Andreas<br clear="all">
<div>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<p style="margin-bottom:0.0001pt"><span lang="EN-US" style="color:rgb(31,73,125)"><span lang="EN-US" style="font-size:11pt;line-height:16.8667px;font-family:Calibri,sans-serif">__________________________________</span></span></p>
<p style="margin-bottom:0.0001pt"><span lang="EN-US" style="color:rgb(31,73,125)"><span lang="EN-US" style="font-size:11pt;line-height:16.8667px;font-family:Calibri,sans-serif">Andreas Gravgaard Andersen</span></span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Danish Center for Particle Therapy, </span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Aarhus University Hospital<u></u><u></u></span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Palle Juul-Jensens Blvd. 99,</span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">8200, Aarhus</span></p>
<div>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span lang="DA" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Mail: <a href="mailto:agravgaard@protonmail.com" target="_blank">agravgaard@protonmail.com</a></span></p>
<p style="margin:0cm 0cm 0.0001pt;font-size:12pt;font-family:"Times New Roman",serif;color:black">
<span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt">Cell: </span><span style="color:rgb(31,73,125);font-family:Calibri,sans-serif;font-size:11pt"> </span><a value="+4523382411" style="font-family:Calibri,sans-serif;font-size:11pt;color:rgb(17,85,204)">+45
3165 8140</a></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
</div>
</div>
<br>
<div>
<div dir="ltr">On Tue, 26 May 2020 at 05:08, Brani Rusanov <<a href="mailto:Brani_Rusanov@hotmail.com" target="_blank">Brani_Rusanov@hotmail.com</a>> wrote:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
Hi all,</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
I am attempting to forward project a CT volume (previously registered to a CBCT volume) using the geometry of a CBCT volume to generate raw projections. I also have access to the raw CBCT projections. I have used rktprojections to stack the raw CBCT projections
and generate .mhd/raw files. </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
Next, I used rtkvariabprobeamgeometry on the scan.xml file that came with the raw CBCT projections to generate a RTK compatible geometry.xml file. </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
Next, I have stacked the registered CT dicoms into a .mhd/raw format. The trouble I'm having is creating the projections using rtkforwardprojections. Specifically, I cannot find the correct transformmatrix in the CT header to tell rtkforwardprojetions the conversion
from DICOM to IEC geometry. The output projections do not look like the raw projections I have from the CBCT. This may be a stupid question, but I cannot figure it out. From my own estimation, the transformation should be 1 0 0 0 0 -1 0 1 0 based on this diagram
(<a href="https://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.semanticscholar.org%2Fpaper%2Fnuclear-science-and-technology-Radiation-Quantities-Brink-McNitt-Gray%2F1ec27be3617dfc4e03bba86dcd92db8f3323de21%2Ffigure%2F3&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580849808&sdata=%2BcHqKs1UZFalgC6Aug0OUp8qmGor5gm2ubuS5FmwoJg%3D&reserved=0" rel="noopener noreferrer" style="margin:0px" target="_blank">https://www.semanticscholar.org/paper/nuclear-science-and-technology-Radiation-Quantities-Brink-McNitt-Gray/1ec27be3617dfc4e03bba86dcd92db8f3323de21/figure/3</a>
<div style="margin:16px 0px;max-width:800px;min-width:424px">
<table style="padding:12px;width:800px;border-width:1px;border-style:solid;border-color:rgb(200,200,200);border-radius:2px">
<tbody>
<tr style="border-spacing:0px">
<td>
<div style="margin:0px 12px 0px 0px;height:104.547px;overflow:hidden"><a href="https://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.semanticscholar.org%2Fpaper%2Fnuclear-science-and-technology-Radiation-Quantities-Brink-McNitt-Gray%2F1ec27be3617dfc4e03bba86dcd92db8f3323de21%2Ffigure%2F3&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580859815&sdata=E3D3%2FE%2FFGybdE%2FFPZGFIuN8nGTwrd8NNiLP%2BhtXthdE%3D&reserved=0" style="margin:0px" target="_blank"><img alt="" height="104" width="240" style="margin: 0px; display: block;" src="https://ai2-s2-public.s3.amazonaws.com/figures/2017-08-08/1ec27be3617dfc4e03bba86dcd92db8f3323de21/14-Figure1.4-1.png"></a></div>
</td>
<td style="width:512.667px">
<div style="margin:0px 8px 12px 0px;font-weight:300;font-size:21px;font-family:wf_segoe-ui_light,"Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif">
<a href="https://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.semanticscholar.org%2Fpaper%2Fnuclear-science-and-technology-Radiation-Quantities-Brink-McNitt-Gray%2F1ec27be3617dfc4e03bba86dcd92db8f3323de21%2Ffigure%2F3&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580869808&sdata=jvzkLIg7%2FGfktUYgbxwd7Um%2BIAG3jnivsoo0sitiFEw%3D&reserved=0" style="margin:0px" target="_blank">Figure
1.4 from nuclear science and technology Radiation Quantities and Units , Dose to the Patients , and Image Quality in Computed Tomography ( CT ) ( RAD UNITS ) Author : | Semantic Scholar</a></div>
<div style="margin:0px 8px 12px 0px;font-size:14px;font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;color:rgb(102,102,102);max-height:100px;overflow:hidden">
Figure 1.4 (left) The IEC patient coordinate system (the patient is lying in supine position); (right) transformation of coordinates between the IEC and DICOM coordinate systems [IEC, 2000] (cited in ICRU Report 71, 2004, page 35) - "nuclear science and technology
Radiation Quantities and Units , Dose to the Patients , and Image Quality in Computed Tomography ( CT ) ( RAD UNITS ) Author :"</div>
<div style="margin:0px;font-size:14px;font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif;color:rgb(166,166,166)">
<a href="https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.semanticscholar.org%2F&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580879808&sdata=rVbKGVfyqhrgyhcDV%2BON%2F%2FavYdAuXYUzWsKQgLpps%2FA%3D&reserved=0" target="_blank">www.semanticscholar.org</a></div>
</td>
</tr>
</tbody>
</table>
<span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt">)</span></div>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
My questions are: Is this the correct transformmatrix to apply?<br>
Do I need to also change the offset values in order to get the correct projections? <br>
Could this have something to do with the fact that I am using variant data?<br>
<br>
The .mhd for the registered CT is:<br>
<span style="margin:0px">ObjectType = Image<br>
</span>
<div style="margin:0px">NDims = 3<br>
</div>
<div style="margin:0px">BinaryData = True<br>
</div>
<div style="margin:0px">BinaryDataByteOrderMSB = False<br>
</div>
<div style="margin:0px">CompressedData = False<br>
</div>
<div style="margin:0px">TransformMatrix = 1 0 0 0 0 -1 0 1 0<br>
</div>
<div style="margin:0px">Offset = -131.03451538085938 -131.03451538085938 -91.809638977050781 <br>
</div>
<div style="margin:0px">CenterOfRotation = 0 0 0<br>
</div>
<div style="margin:0px">AnatomicalOrientation = RAI<br>
</div>
<div style="margin:0px">ElementSpacing = 0.51285523409720002 0.51285523409720002 1.9958648681640625<br>
</div>
<div style="margin:0px">DimSize = 512 512 93<br>
</div>
<div style="margin:0px">ElementType = MET_DOUBLE<br>
</div>
<span style="margin:0px">ElementDataFile = CIRS_6A.raw<br>
</span><br>
Similarly, the .mhd for the raw CBCT projections is:<br>
<span style="margin:0px">ObjectType = Image<br>
</span>
<div style="margin:0px">NDims = 3<br>
</div>
<div style="margin:0px">BinaryData = True<br>
</div>
<div style="margin:0px">BinaryDataByteOrderMSB = False<br>
</div>
<div style="margin:0px">CompressedData = False<br>
</div>
<div style="margin:0px">TransformMatrix = 1 0 0 0 1 0 0 0 1<br>
</div>
<div style="margin:0px">Offset = -198.46200000000002 -148.798 0<br>
</div>
<div style="margin:0px">CenterOfRotation = 0 0 0<br>
</div>
<div style="margin:0px">AnatomicalOrientation = RAI<br>
</div>
<div style="margin:0px">ElementSpacing = 0.38800000000000001 0.38800000000000001 1<br>
</div>
<div style="margin:0px">DimSize = 1024 768 501<br>
</div>
<div style="margin:0px">ElementType = MET_FLOAT<br>
</div>
<span style="margin:0px">ElementDataFile = CIRS_6A_SCGH_P.raw</span><br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="margin:0px"><br>
</span></div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="margin:0px">Thank you kindly for your response,</span></div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="margin:0px">Brani</span></div>
<br>
</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://nam05.safelinks.protection.outlook.com/?url=https%3A%2F%2Fpublic.kitware.com%2Fmailman%2Flistinfo%2Frtk-users&data=02%7C01%7C%7C63220d57e95844eb634508d8016629c0%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C637260885580889790&sdata=P7oUcTzE8acsSOlau9pS9WtPtRqDU1MyE9vQ%2BxPxRnk%3D&reserved=0" rel="noreferrer" target="_blank">https://public.kitware.com/mailman/listinfo/rtk-users</a><br>
</blockquote>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</blockquote></div>