[Insight-users] accessing the HEADER of a DICOM file

Evelyne Balteau e . balteau at ulg . ac . be
Wed, 10 Sep 2003 16:38:26 +0200


This is a multi-part message in MIME format.

------=_NextPart_000_003C_01C377B9.F51D2950
Content-Type: text/plain;
	charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

Hello !

I have written a code based on the GLSliceView.h header =
(InsightApplications) and=20
the itkDICOMImageSeriesTest, to view slices of a 3D volume from a series =
of DICOM files...
I successfully tested it on a wide range of (magnitude) images, but I =
had some problems with phase data.
I finally found the solution : the phase data are "rescaled" by =
substracting 4096 and multiplying by 2 (keeping unsigned short type) !!!
Well...
To write a code that could take this into account WHEN REQUIRED (for =
phase but not for magnitude images), I need to access the header of the =
DICOM files.
There are two fields (RescaleIntercept=3D-4096 and RescaleSlope=3D2) in =
this header giving informations on the rescaling factors.
I could read it with matlab, but I don't know how to do with ITK...
Is there some method to do so ??

Thanks for any help !!

Evelyne
------=_NextPart_000_003C_01C377B9.F51D2950
Content-Type: text/html;
	charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=3DContent-Type content=3D"text/html; =
charset=3Diso-8859-1">
<META content=3D"MSHTML 6.00.2800.1226" name=3DGENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=3D#ffffff>
<DIV><FONT face=3DArial size=3D2>Hello !</FONT></DIV>
<DIV><FONT face=3DArial size=3D2></FONT>&nbsp;</DIV>
<DIV><FONT face=3DArial size=3D2>I have written a code based on the =
GLSliceView.h=20
header (InsightApplications) and </FONT></DIV>
<DIV><FONT face=3DArial size=3D2>the itkDICOMImageSeriesTest, to view =
slices of a 3D=20
volume from a series of DICOM files...</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>I successfully tested it on a wide =
range of=20
(magnitude) images, but I had some problems with phase =
data.</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>I finally found the solution : the =
phase data are=20
"rescaled" by substracting 4096 and multiplying by 2 (keeping unsigned =
short=20
type) !!!</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>Well...</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>To write a code that could take this =
into account=20
WHEN REQUIRED (for phase but not for magnitude images), I need to access =
the=20
header of the DICOM files.</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>There are two fields =
(RescaleIntercept=3D-4096 and=20
RescaleSlope=3D2) in this header giving informations on the rescaling=20
factors.</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>I could read it with matlab, but I =
don't know how=20
to do with ITK...</FONT></DIV>
<DIV><FONT face=3DArial size=3D2>Is there some method to do so =
??</FONT></DIV>
<DIV><FONT face=3DArial size=3D2></FONT>&nbsp;</DIV>
<DIV><FONT face=3DArial size=3D2>Thanks for any help !!</FONT></DIV>
<DIV><FONT face=3DArial size=3D2></FONT>&nbsp;</DIV>
<DIV><FONT face=3DArial size=3D2>Evelyne</FONT></DIV></BODY></HTML>

------=_NextPart_000_003C_01C377B9.F51D2950--