[Insight-developers] GDCM Dicom Series Reader

[Vincent Magnotta]

I have been looking at the output of the GDCM series reader. I believe 
that there is a bug in the setting of the origin in the resulting image. 
I collected a set of images that were prescribed as Axial, Sagittal, and 
Coronal. Here is the result of the origin setting from the image series 
reader:

0002: (2D Axial T1)

============ ITK Output ============
Origin: [-120, -128.475, -82.7712]
Image Position Patient = -120.00000002461\-128.47457597179\50.228813171387
Image Orientation (Patient) = 1\-2.0510349e-010\0\2.051034897e-010\1\0

===== Image Positions from DICOM Header ====
-120.00000002461 -128.47457597179 50.228813171387 
-120.00000002461 -128.47457597179 43.228813171387 
-120.00000002461 -128.47457597179 36.228813171387 
-120.00000002461 -128.47457597179 29.228813171387 
-120.00000002461 -128.47457597179 22.228813171387 
-120.00000002461 -128.47457597179 15.228813171387 
-120.00000002461 -128.47457597179 8.2288131713867 
-120.00000002461 -128.47457597179 1.2288135290146 
-120.00000002461 -128.47457597179 -5.7711863517761 
-120.00000002461 -128.47457597179 -12.771186828613 
-120.00000002461 -128.47457597179 -19.771186828613 
-120.00000002461 -128.47457597179 -26.771186828613 
-120.00000002461 -128.47457597179 -33.771186828613 
-120.00000002461 -128.47457597179 -40.771186828613 
-120.00000002461 -128.47457597179 -47.771186828613 
-120.00000002461 -128.47457597179 -54.771186828613 
-120.00000002461 -128.47457597179 -61.771186828613 
-120.00000002461 -128.47457597179 -68.771186828613 
-120.00000002461 -128.47457597179 -75.771186828613 
-120.00000002461 -128.47457597179 -82.771186828613 




0003: (2D Coronal T1)
============ ITK Output ============
Origin: [-120, -68.0847, 104.407]
Image Position Patient = -120\57.915252685547\104.40677928925
Image Orientation (Patient) = 1\0\0\0\0\-1

===== Image Positions from DICOM Header ====
-120 57.915252685547 104.40677928925 
-120 50.915252685547 104.40677928925 
-120 43.915252685547 104.40677928925 
-120 36.915252685547 104.40677928925 
-120 29.915254592896 104.40677928925 
-120 22.915254592896 104.40677928925 
-120 15.915254592896 104.40677928925 
-120 8.9152545928955 104.40677928925 
-120 1.9152542352676 104.40677928925 
-120 -5.0847458839417 104.40677928925 
-120 -12.084745407104 104.40677928925 
-120 -19.084745407104 104.40677928925 
-120 -26.084745407104 104.40677928925 
-120 -33.084747314453 104.40677928925 
-120 -40.084747314453 104.40677928925 
-120 -47.084747314453 104.40677928925 
-120 -54.084747314453 104.40677928925 
-120 -61.084747314453 104.40677928925 
-120 -68.084747314453 104.40677928925 

0004: (2D Sagittal T1)
============ ITK Output ============
Origin: [47.0763, -122.373, 104.407]
Image Position Patient = -57.923728942871\-122.37288141251\104.40677928925
Image Orientation (Patient) = 0\1\0\0\0\-1


===== Image Positions from DICOM Header ====

-57.923728942871 -122.37288141251 104.40677928925 
-50.923728942871 -122.37288141251 104.40677928925 
-43.923728942871 -122.37288141251 104.40677928925 
-36.923728942871 -122.37288141251 104.40677928925 
-29.923728942871 -122.37288141251 104.40677928925 
-22.923728942871 -122.37288141251 104.40677928925 
-15.923728942871 -122.37288141251 104.40677928925 
-8.9237289428711 -122.37288141251 104.40677928925 
-1.9237288236618 -122.37288141251 104.40677928925 
5.0762710571289 -122.37288141251 104.40677928925 
12.076271057129 -122.37288141251 104.40677928925 
19.076271057129 -122.37288141251 104.40677928925 
26.076271057129 -122.37288141251 104.40677928925 
33.076271057129 -122.37288141251 104.40677928925 
40.076271057129 -122.37288141251 104.40677928925 
47.076271057129 -122.37288141251 104.40677928925 



 Based on this data, it suggests the Axial image is correct. The 
Sagittal and Coronal images are not correct. The origin is always being 
set based on the scanner coordinate systm (x,y,z) and not that of the 
images. I am going to try to make the Tcon to discuss this further if 
there is time.

Vince



-- 
Assistant Professor
Department of Radiology
0453-D JCP
200 Hawkins Drive
Iowa City, IA 52242
E-mail: vincent-magnotta at uiowa.edu
Phone: 319-356-8255
Fax: 319-353-6275
Website: http://www.radiology.uiowa.edu