Dear all, I am trying to simulate Tomosynthesis reconstruction with stationary flat panel detector where source rotates around Y axis to acquire 21 projections (system geometry diagram is attached). Geometry file geometry.xml is generated with Reg23ProjectionGeometry. I generated projection images of a cube volume and then I tried to reconstruct with FDK algorithm. Even though I am able to generate projections, FDK reconstruction yields the volume with all voxel values zero. The commands I used for generating projections and reconstructing volume are
rtkforwardprojections -g geometry.xml -o projections.mhd -i volume.mhd --spacing 0.1518,0.1518,1 --dimension 1420,1416,21 --fp=CudaRayCast rtkfdk -p . -r projections.mhd -o fdk.mhd -g geometry.xml --spacing=0.25,0.25,0.25 --origin=-31.875,-31.875,25 --dimension=256,256,256 --hardware=cuda I am also attaching geometry xml file, mhd header files of the original volume, projections and reconstructed volume. With regards, Shiras
<?xml version="1.0"?>
<!DOCTYPE RTKGEOMETRY>
<RTKThreeDCircularGeometry version="2">
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>367.054653749843</SourceToIsocenterDistance>
<SourceToDetectorDistance>367.054653749843</SourceToDetectorDistance>
<SourceOffsetX>171.160396004983</SourceOffsetX>
<Matrix>
-367.054653749843 0 171.160396004983 0
0 -367.054653749843 0 0
0 0 1 -367.054653749843
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>374.171210667071</SourceToIsocenterDistance>
<SourceToDetectorDistance>374.171210667071</SourceToDetectorDistance>
<SourceOffsetX>154.986790107861</SourceOffsetX>
<Matrix>
-374.171210667071 0 154.986790107861 0
0 -374.171210667071 0 0
0 0 1 -374.171210667071
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>380.575511418293</SourceToIsocenterDistance>
<SourceToDetectorDistance>380.575511418293</SourceToDetectorDistance>
<SourceOffsetX>138.518158046896</SourceOffsetX>
<Matrix>
-380.575511418293 0 138.518158046896 0
0 -380.575511418293 0 0
0 0 1 -380.575511418293
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>386.255365053032</SourceToIsocenterDistance>
<SourceToDetectorDistance>386.255365053032</SourceToDetectorDistance>
<SourceOffsetX>121.785848799231</SourceOffsetX>
<Matrix>
-386.255365053032 0 121.785848799231 0
0 -386.255365053032 0 0
0 0 1 -386.255365053032
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>391.199959647073</SourceToIsocenterDistance>
<SourceToDetectorDistance>391.199959647073</SourceToDetectorDistance>
<SourceOffsetX>104.821713266521</SourceOffsetX>
<Matrix>
-391.199959647073 0 104.821713266521 0
0 -391.199959647073 0 0
0 0 1 -391.199959647073
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>395.399882883573</SourceToIsocenterDistance>
<SourceToDetectorDistance>395.399882883573</SourceToDetectorDistance>
<SourceOffsetX>87.6580436449317</SourceOffsetX>
<Matrix>
-395.399882883573 0 87.6580436449317 0
0 -395.399882883573 0 0
0 0 1 -395.399882883573
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>398.847139969944</SourceToIsocenterDistance>
<SourceToDetectorDistance>398.847139969944</SourceToDetectorDistance>
<SourceOffsetX>70.3275119551068</SourceOffsetX>
<Matrix>
-398.847139969944 0 70.3275119551068 0
0 -398.847139969944 0 0
0 0 1 -398.847139969944
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>401.535168856393</SourceToIsocenterDistance>
<SourceToDetectorDistance>401.535168856393</SourceToDetectorDistance>
<SourceOffsetX>52.8631078491209</SourceOffsetX>
<Matrix>
-401.535168856393 0 52.8631078491209 0
0 -401.535168856393 0 0
0 0 1 -401.535168856393
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>403.458852727157</SourceToIsocenterDistance>
<SourceToDetectorDistance>403.458852727157</SourceToDetectorDistance>
<SourceOffsetX>35.2980758128016</SourceOffsetX>
<Matrix>
-403.458852727157 0 35.2980758128016 0
0 -403.458852727157 0 0
0 0 1 -403.458852727157
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>404.614529740652</SourceToIsocenterDistance>
<SourceToDetectorDistance>404.614529740652</SourceToDetectorDistance>
<SourceOffsetX>17.6658518829611</SourceOffsetX>
<Matrix>
-404.614529740652 0 17.6658518829611 0
0 -404.614529740652 0 0
0 0 1 -404.614529740652
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>405</SourceToIsocenterDistance>
<SourceToDetectorDistance>405</SourceToDetectorDistance>
<SourceOffsetX>0</SourceOffsetX>
<Matrix>
-405 0 0 0
0 -405 0 0
0 0 1 -405
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>404.614529740652</SourceToIsocenterDistance>
<SourceToDetectorDistance>404.614529740652</SourceToDetectorDistance>
<SourceOffsetX>-17.6658518829611</SourceOffsetX>
<Matrix>
-404.614529740652 0 -17.6658518829611 0
0 -404.614529740652 0 0
0 0 1 -404.614529740652
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>403.458852727157</SourceToIsocenterDistance>
<SourceToDetectorDistance>403.458852727157</SourceToDetectorDistance>
<SourceOffsetX>-35.2980758128016</SourceOffsetX>
<Matrix>
-403.458852727157 0 -35.2980758128016 0
0 -403.458852727157 0 0
0 0 1 -403.458852727157
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>401.535168856393</SourceToIsocenterDistance>
<SourceToDetectorDistance>401.535168856393</SourceToDetectorDistance>
<SourceOffsetX>-52.8631078491209</SourceOffsetX>
<Matrix>
-401.535168856393 0 -52.8631078491209 0
0 -401.535168856393 0 0
0 0 1 -401.535168856393
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>398.847139969944</SourceToIsocenterDistance>
<SourceToDetectorDistance>398.847139969944</SourceToDetectorDistance>
<SourceOffsetX>-70.3275119551068</SourceOffsetX>
<Matrix>
-398.847139969944 0 -70.3275119551068 0
0 -398.847139969944 0 0
0 0 1 -398.847139969944
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>395.399882883573</SourceToIsocenterDistance>
<SourceToDetectorDistance>395.399882883573</SourceToDetectorDistance>
<SourceOffsetX>-87.6580436449317</SourceOffsetX>
<Matrix>
-395.399882883573 0 -87.6580436449317 0
0 -395.399882883573 0 0
0 0 1 -395.399882883573
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>391.199959647073</SourceToIsocenterDistance>
<SourceToDetectorDistance>391.199959647073</SourceToDetectorDistance>
<SourceOffsetX>-104.821713266521</SourceOffsetX>
<Matrix>
-391.199959647073 0 -104.821713266521 0
0 -391.199959647073 0 0
0 0 1 -391.199959647073
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>386.255365053032</SourceToIsocenterDistance>
<SourceToDetectorDistance>386.255365053032</SourceToDetectorDistance>
<SourceOffsetX>-121.785848799231</SourceOffsetX>
<Matrix>
-386.255365053032 0 -121.785848799231 0
0 -386.255365053032 0 0
0 0 1 -386.255365053032
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>380.575511418293</SourceToIsocenterDistance>
<SourceToDetectorDistance>380.575511418293</SourceToDetectorDistance>
<SourceOffsetX>-138.518158046896</SourceOffsetX>
<Matrix>
-380.575511418293 0 -138.518158046896 0
0 -380.575511418293 0 0
0 0 1 -380.575511418293
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>374.171210667071</SourceToIsocenterDistance>
<SourceToDetectorDistance>374.171210667071</SourceToDetectorDistance>
<SourceOffsetX>-154.986790107861</SourceOffsetX>
<Matrix>
-374.171210667071 0 -154.986790107861 0
0 -374.171210667071 0 0
0 0 1 -374.171210667071
</Matrix>
</Projection>
<Projection>
<GantryAngle>0</GantryAngle>
<SourceToIsocenterDistance>367.054653749843</SourceToIsocenterDistance>
<SourceToDetectorDistance>367.054653749843</SourceToDetectorDistance>
<SourceOffsetX>-171.160396004983</SourceOffsetX>
<Matrix>
-367.054653749843 0 -171.160396004983 0
0 -367.054653749843 0 0
0 0 1 -367.054653749843
</Matrix>
</Projection>
</RTKThreeDCircularGeometry>
fdk.mhd
Description: Binary data
projections.mhd
Description: Binary data
volume.mhd
Description: Binary data
_______________________________________________ Rtk-users mailing list [email protected] http://public.kitware.com/mailman/listinfo/rtk-users
