"""
Usage:
python GetPortionCSAHeader.py input.dcm
Footnote:
SIEMENS is not publishing any information on the CSA header. So any info extracted
is at your own risk.
"""
import sys
import gdcm
if __name__ == "__main__":
file = sys.argv[1]
r = gdcm.Reader()
r.SetFileName( file )
if not r.Read():
sys.exit(1)
ds = r.GetFile().GetDataSet()
csa_t1 = gdcm.CSAHeader()
csa_t2 = gdcm.CSAHeader()
t1 = csa_t1.GetCSAImageHeaderInfoTag();
print t1
t2 = csa_t2.GetCSASeriesHeaderInfoTag();
print t2
if ds.FindDataElement( t1 ):
csa_t1.LoadFromDataElement( ds.GetDataElement( t1 ) )
print csa_t1
bvalues = csa_t1.GetCSAElementByName( "B_value" )
print bvalues
diffgraddir = csa_t1.GetCSAElementByName( "DiffusionGradientDirection" )
print diffgraddir
if ds.FindDataElement( t2 ):
csa_t2.LoadFromDataElement( ds.GetDataElement( t2 ) )
gdt = csa_t2.GetCSAElementByName( "GradientDelayTime" )
print gdt
bv = gdt.GetByteValue();
str = bv.GetPointer()
print str.split("\\")