gdcmPython/testSuite.py

   1 import unittest
   2 import os
   3 from gdcmPython import *
   4
   5 class gdcmTestCase(unittest.TestCase):
   6    # The files whose name starts with a modality (e.g. CR-MONO1-10-chest.dcm)
   7    # come from Sebastien Barre's Dicom2 highly recommendable site
   8    # http://www.barre.nom.fr/medical/samples/index.html
   9    BarreFiles = [
  10       ####################################
  11       # CR modality examples:
  12       ####################################
  13       ["CR-MONO1-10-chest.dcm",
  14          [ ["Modality", "CR"],
  15            #"Transfer Syntax UID" is absent.
  16            ["Rows", "440"],
  17            ["Columns", "440"],
  18            ["Bits Stored", "10"],
  19            ["Bits Allocated", "16"],
  20            ["High Bit", "9"],
  21            ["Pixel Representation", "0"],
  22            ["Manufacturer", "FUJI PHOTO FILM CO. LTD."],
  23            ["Manufacturer Model Name", "9000"],
  24            ["Pixel Data", "gdcm::NotLoaded. Address:776 Length:387200"] ] ],
  25       ####################################
  26       # CT modality examples:
  27       ####################################
  28       ["CT-MONO2-16-brain.dcm",
  29          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"],  # Explicit VR, LE
  30            ["Modality", "CT"],
  31            ["Rows", "512"],
  32            ["Columns", "512"],
  33            ["Bits Stored", "16"],
  34            ["Bits Allocated", "16"],
  35            ["High Bit", "15"],
  36            ["Pixel Representation", "1"],
  37            ["Manufacturer", "Picker International, Inc."],
  38            ["Manufacturer Model Name", "PQ5000"],
  39            ["Pixel Data", "gdcm::NotLoaded. Address:1680 Length:524288"] ] ],
  40       ["CT-MONO2-16-ort.dcm",
  41          [ ["Transfer Syntax UID", "1.2.840.10008.1.2"],  # Implicit VR, LE
  42            ["Modality", "CT"],
  43            ["Rows", "512"],
  44            ["Columns", "512"],
  45            ["Bits Stored", "16"],
  46            ["Bits Allocated", "16"],
  47            ["High Bit", "15"],
  48            ["Pixel Representation", "1"],
  49            ["Manufacturer", "GE MEDICAL SYSTEMS"],
  50            ["Manufacturer Model Name", "HiSpeed CT/i"],
  51            ["Pixel Data", "gdcm::NotLoaded. Address:1674 Length:524288"] ] ],
  52       ["CT-MONO2-16-ankle.dcm",
  53          [ ["Transfer Syntax UID", "1.2.840.10008.1.2"],  # Implicit, little
  54            ["Modality", "CT"],
  55            ["Rows", "512"],
  56            ["Columns", "512"],
  57            ["Bits Stored", "16"],
  58            ["Bits Allocated", "16"],
  59            ["High Bit", "15"],
  60            ["Pixel Representation", "1"],
  61            ["Manufacturer", "GE MEDICAL SYSTEMS"],
  62            ["Manufacturer Model Name", "GENESIS_ZEUS"],
  63            ["Pixel Data", "gdcm::NotLoaded. Address:1148 Length:524288"] ] ],
  64       ["CT-MONO2-8-abdo.dcm",
  65          [ ["Transfer Syntax UID", "1.2.840.10008.1.2"],  # Implicit, little
  66            ["Modality", "CT"],
  67            ["Rows", "512"],
  68            ["Columns", "512"],
  69            ["Bits Stored", "8"],
  70            ["Bits Allocated", "8"],
  71            ["High Bit", "7"],
  72            ["Pixel Representation", "0"],
  73            ["Manufacturer", "FUJI"],
  74            ["Pixel Data", "gdcm::NotLoaded. Address:796 Length:262144"] ] ],
  75       ["CT-MONO2-12-lomb-an2.acr2",
  76          [ ["Modality", "CT"],
  77            #"Transfer Syntax UID" and "Photometric Interpretation" are absent.
  78            ["Rows", "512"],
  79            ["Columns", "512"],
  80            ["Bits Stored", "12"],
  81            ["Bits Allocated", "16"],
  82            ["High Bit", "11"],
  83            ["Pixel Representation", "0"],
  84            ["Manufacturer", "SIEMENS"],
  85            ["Pixel Data", "gdcm::NotLoaded. Address:1230 Length:524288"] ] ],
  86       ["CT-MONO2-16-chest.dcm",
  87          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.4.70"], # jpeg...
  88            ["Modality", "CT"],
  89            ["Photometric Interpretation", "MONOCHROME2"],
  90            ["Rows", "400"],
  91            ["Columns", "512"],
  92            ["Bits Stored", "16"],
  93            ["Bits Allocated", "16"],
  94            ["High Bit", "15"],
  95            ["Pixel Representation", "1"],
  96            ["Manufacturer", "GE MEDICAL SYSTEMS"],
  97            ["Pixel Data", "gdcm::NotLoaded. Address:1638 Length:143498"] ] ],
  98       ####################################
  99       ### MR modality examples:
 100       ####################################
 101       ["MR-MONO2-12-angio-an1.acr1",
 102          [ ["Recognition Code", "ACR-NEMA 1.0"],
 103            ["Modality", "MR"],
 104            #"Transfer Syntax UID" and "Photometric Interpretation" are absent.
 105            ["Rows", "256"],
 106            ["Columns", "256"],
 107            ["Bits Stored", "12"],
 108            ["Bits Allocated", "12"],
 109            ["High Bit", "11"],
 110            ["Pixel Representation", "0"],
 111            ["Manufacturer", "Philips"],
 112            # Note: 256*256*12/8 = 98304
 113            ["Pixel Data", "gdcm::NotLoaded. Address:650 Length:98304"] ] ],
 114       ["MR-MONO2-12-an2.acr2",
 115          [ ["Modality", "MR"],
 116            ["Photometric Interpretation", "MONOCHROME2"],
 117            ["Rows", "256"],
 118            ["Columns", "256"],
 119            ["Bits Stored", "12"],
 120            ["Bits Allocated", "16"],
 121            ["High Bit", "11"],
 122            ["Pixel Representation", "0"],
 123            ["Manufacturer", "SIEMENS"],
 124            ["Manufacturer Model Name", "MAGNETOM VISION"],
 125            ["Pixel Data", "gdcm::NotLoaded. Address:1866 Length:131072"] ] ],
 126       ["MR-MONO2-16-head.dcm",
 127          [ ["Transfer Syntax UID", "1.2.840.10008.1.2"], # Imp VR, LE
 128            ["Modality", "MR"],
 129            ["Photometric Interpretation", "MONOCHROME2"],
 130            ["Rows", "256"],
 131            ["Columns", "256"],
 132            ["Bits Stored", "16"],
 133            ["Bits Allocated", "16"],
 134            ["High Bit", "15"],
 135            ["Pixel Representation", "1"],
 136            ["Manufacturer", "GE MEDICAL SYSTEMS"],
 137            ["Manufacturer Model Name", "GENESIS_SIGNA"],
 138            ["Pixel Data", "gdcm::NotLoaded. Address:1804 Length:131072"] ] ],
 139       # MR-MONO2-12-shoulder example is present in the Jpeg section.
 140       # MR-MONO2-16-knee is not present in the test suite since it is too
 141       #    closely related to MR-MONO2-16-head.dcm to be of interest.
 142       ####################################
 143       # OT modality examples:
 144       ####################################
 145       # OT-MONO2-8-hip is not present in the test suite since it is too
 146       #    closely related to OT-MONO2-8-a7 to be of interest.
 147       ["OT-PAL-8-face.dcm",
 148          # Interest: has a color palette
 149          [#"Transfer Syntax UID", and "Manufacturer" are absent.
 150            ["Modality", "OT"],
 151            ["Photometric Interpretation", "PALETTE COLOR"],
 152            ["Rows", "480"],
 153            ["Columns", "640"],
 154            ["Bits Stored", "8"],
 155            ["Bits Allocated", "8"],
 156            ["High Bit", "7"],
 157            ["Pixel Representation", "0"],
 158            ["Pixel Data", "gdcm::NotLoaded. Address:1654 Length:307200"] ] ],
 159       # OT-MONO2-8-colon is not present in the test suite since it is too
 160       #    closely related to OT-MONO2-8-a7 to be of interest.
 161       ["OT-MONO2-8-a7.dcm",
 162          [#"Transfer Syntax UID" is absent.
 163            ["Modality", "OT"],
 164            ["Photometric Interpretation", "MONOCHROME2"],
 165            ["Rows", "512"],
 166            ["Columns", "512"],
 167            ["Bits Stored", "8"],
 168            ["Bits Allocated", "8"],
 169            ["High Bit", "7"],
 170            ["Pixel Representation", "0"],
 171            ["Manufacturer", "Philips Medical Systems"],
 172            ["Pixel Data", "gdcm::NotLoaded. Address:438 Length:262144"] ] ],
 173       ####################################
 174       # US (Ultra Sound) modality examples:
 175       ####################################
 176       ["US-RGB-8-esopecho.dcm",
 177          # Interest: RGB image
 178          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"],  # Exp VR, LE
 179            ["Modality", "US"],
 180            ["Photometric Interpretation", "RGB"],
 181            ["Rows", "120"],
 182            ["Columns", "256"],
 183            ["Bits Stored", "8"],
 184            ["Bits Allocated", "8"],
 185            ["High Bit", "7"],
 186            ["Pixel Representation", "0"],
 187            ["Manufacturer", "Acme Products"],
 188          # FIXME: 92160 / (120*256) = 3 bytes per pixel NOT 1. Maybe
 189          # it has something to do with [Samples Per Pixel] = [3] ???
 190            ["Pixel Data", "gdcm::NotLoaded. Address:904 Length:92160"] ] ],
 191       ["US-RGB-8-epicard.dcm",
 192          # Interest: Big endian transfert syntax
 193          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.2"],  # Big endian
 194            ["Modality", "US"],
 195            ["Photometric Interpretation", "RGB"],
 196            ["Rows", "480"],
 197            ["Columns", "640"],
 198            ["Bits Stored", "8"],
 199            ["Bits Allocated", "8"],
 200            ["High Bit", "7"],
 201            ["Pixel Representation", "0"],
 202            ["Manufacturer", "G.E. Medical Systems"],
 203            ["Manufacturer Model Name", "LOGIQ 700"],
 204          # FIXME: 921600/(480*640) = 3 bytes per pixel NOT 1. Maybe
 205          # it has something to do with [Samples Per Pixel] = [3] ???
 206            ["Implementation Version Name", "OFFIS-DCMTK-311"],
 207            ["Pixel Data", "gdcm::NotLoaded. Address:1012 Length:921600"] ] ],
 208    ]
 209
 210    MultiFrameFiles = [
 211    # The number of images is contained in the "Number of Frames" element
 212       ["MR-MONO2-8-16x-heart.dcm",
 213          # Interest: multi-frame
 214          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"], # ExpVR, LE
 215            ["Number of Frames", "16"],
 216            ["Modality", "MR"],
 217            ["Photometric Interpretation", "MONOCHROME2"],
 218            ["Rows", "256"],
 219            ["Columns", "256"],
 220            ["Bits Stored", "8"],
 221            ["Bits Allocated", "8"],
 222            ["High Bit", "7"],
 223            ["Pixel Representation", "0"],
 224            ["Manufacturer", "Philips"],
 225            # Note: 256*256*16=1048576 (16 is the number of Frames)
 226            ["Pixel Data", "gdcm::NotLoaded. Address:920 Length:1048576"] ] ],
 227       ["NM-MONO2-16-13x-heart.dcm",
 228          # Interest: NM modality, multi-frame
 229          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"], # ExpVR, LE
 230            ["Number of Frames", "13"],
 231            ["Modality", "NM"],
 232            ["Photometric Interpretation", "MONOCHROME2"],
 233            ["Rows", "64"],
 234            ["Columns", "64"],
 235            ["Bits Stored", "16"],
 236            ["Bits Allocated", "16"],
 237            ["High Bit", "15"],
 238            ["Pixel Representation", "0"],
 239            ["Manufacturer", "ACME Products"],
 240            # Note: 64*64*13*2=106496
 241            ["Pixel Data", "gdcm::NotLoaded. Address:1234 Length:106496"] ] ],
 242       ["US-MONO2-8-8x-execho.dcm",
 243          # Interest: multi-frame
 244          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"], # ExpVR, LE
 245            ["Number of Frames", "8"],
 246            ["Modality", "US"],
 247            ["Photometric Interpretation", "MONOCHROME2"],
 248            ["Rows", "120"],
 249            ["Columns", "128"],
 250            ["Bits Stored", "8"],
 251            ["Bits Allocated", "8"],
 252            ["High Bit", "7"],
 253            ["Pixel Representation", "0"],
 254            ["Manufacturer", "Acme Products"],
 255            ["Pixel Data", "gdcm::NotLoaded. Address:976 Length:122880"] ] ],
 256       ["US-PAL-8-10x-echo.dcm",
 257          # Interest: RLE (Run Length Encoded) compression, multiframe
 258          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.5"], # RLE
 259            ["Number of Frames", "10"],
 260            ["Modality", "US"],
 261            ["Photometric Interpretation", "PALETTE COLOR"],
 262            ["Rows", "430"],
 263            ["Columns", "600"],
 264            ["Bits Stored", "8"],
 265            ["Bits Allocated", "8"],
 266            ["High Bit", "7"],
 267            ["Pixel Representation", "0"],
 268            ["Manufacturer", "ACME Products"],
 269            ["Pixel Data", "gdcm::NotLoaded. Address:2428 Length:481182"] ] ],
 270       ["XA-MONO2-8-12x-catheter.dcm",
 271          # Interest: XA modality, multi-frame
 272          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.4.70"], # jpeg
 273            ["Number of Frames", "12"],  # 12 images
 274            ["Modality", "XA"],
 275            ["Photometric Interpretation", "MONOCHROME2"],
 276            ["Rows", "512"],
 277            ["Columns", "512"],
 278            ["Bits Stored", "8"],
 279            ["Bits Allocated", "8"],
 280            ["High Bit", "7"],
 281            ["Pixel Representation", "0"],
 282            ["Manufacturer", "Acme Products"],
 283            ["Pixel Data", "gdcm::NotLoaded. Address:1006 Length:920072"] ] ],
 284    ]
 285
 286    GdcmFiles = [
 287       ["gdcm-MR-SIEMENS-16.acr1",
 288          # Interest: good old Acr-Nema Version 1, see also Oddities below
 289          [ ["Recognition Code", "ACR-NEMA 1.0"],
 290            ["Modality", "MR"],
 291            ["Rows", "256"],
 292            ["Columns", "256"],
 293            ["Bits Stored", "12"],
 294            ["Bits Allocated", "16"],
 295            ["High Bit", "11"],
 296            ["Pixel Representation", "0"],
 297            ["Manufacturer", "SIEMENS"],
 298            ["Manufacturer Model Name", "GBS III"],
 299          # FIXME: this image looks padded at the end. The length of the file
 300          # is 140288. Notice that, 256*256*2 + 1024 = 131072 + 1024 = 132096
 301            ["Pixel Data", "gdcm::NotLoaded. Address:8192 Length:132096"],
 302            # Oddities: "Study ID" and "Series Number" are empty
 303            ["Study ID", ""],
 304            ["Series Number", ""] ] ],
 305            # Oddities: "Study" and "Serie Instance UID" are not present
 306       ["gdcm-MR-SIEMENS-16.acr2",
 307          # Interest: Acr-Nema Version 2
 308          [ ["Recognition Code", "ACR-NEMA 2.0"],
 309            ["Modality", "MR"],
 310            ["Rows", "512"],
 311            ["Columns", "512"],
 312            ["Bits Stored", "12"],
 313            ["Bits Allocated", "16"],
 314            ["High Bit", "11"],
 315            ["Pixel Representation", "0"],
 316            ["Manufacturer", "SIEMENS"],
 317            ["Manufacturer Model Name", "MAGNETOM VISION"],
 318            ["Study Instance UID",
 319             "1.3.12.2.1107.5.2.4.7630.20000918174641000"],
 320            ["Study ID", "1"],
 321            ["Series Instance UID",
 322             "1.3.12.2.1107.5.2.4.7630.20000918175714000007"],
 323            ["Series Number", "7"],
 324            ["Pixel Data", "gdcm::NotLoaded. Address:6052 Length:524288"] ] ],
 325       ["gdcm-US-ALOKA-16.dcm",
 326          # Interest: - possesses 3 LUTS: a Green (checked), a Red and BLue.
 327          #           - announced as implicit VR, but really explicit VR.
 328          [ ["Transfer Syntax UID", "1.2.840.10008.1.2"],  # Implicit VR, LE
 329            ["Modality", "US"],
 330            ["Photometric Interpretation", "PALETTE COLOR"],
 331            ["Segmented Green Palette Color LUT Data",
 332             "gdcm::NotLoaded. Address:89576 Length:113784"],
 333            ["Rows", "480"],
 334            ["Columns", "640"],
 335            ["Bits Stored", "16"],
 336            ["Bits Allocated", "16"],
 337            ["High Bit", "15"],
 338            ["Pixel Representation", "0"],
 339            ["Manufacturer", "ALOKA CO., LTD."],
 340            ["Manufacturer Model Name", "SSD-4000"],
 341            ["Pixel Data", "gdcm::NotLoaded. Address:258740 Length:614400"] ] ],
 342       ["gdcm-MR-PHILIPS-16.dcm",
 343          # Interest: - possesses a sequence
 344          #           - dicom file, with a recognition code of ACR-NEMA1
 345          [ ["Transfer Syntax UID", "1.2.840.10008.1.2"],  # Implicit VR, LE
 346            ["Recognition Code", "ACR-NEMA 1.0"],
 347            ["Modality", "MR"],
 348            ["Photometric Interpretation", "MONOCHROME2"],
 349            ["Rows", "256"],
 350            ["Columns", "256"],
 351            ["Bits Stored", "8"],
 352            ["Bits Allocated", "16"],
 353            ["High Bit", "7"],
 354            ["Pixel Representation", "0"],
 355            ["Manufacturer", "Philips Medical Systems"],
 356            ["Manufacturer Model Name", "Gyroscan Intera"],
 357            ["Sequence Variant", "OTHER"],
 358            ["Pixel Data", "gdcm::NotLoaded. Address:6584 Length:131072"] ] ],
 359       ["gdcm-MR-PHILIPS-16-Multi-Seq.dcm",
 360          # Interest: - possesses many sequences in group 0x0029
 361          #           - Big sequence 28808 bytes at (0x0029, 0x263d)
 362          #           - dicom file, with a recognition code of ACR-NEMA1
 363          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"],
 364            ["Recognition Code", "ACR-NEMA 1.0"],
 365            ["Modality", "MR"],
 366            ["Photometric Interpretation", "MONOCHROME2"],
 367            ["Rows", "128"],
 368            ["Columns", "128"],
 369            ["Bits Stored", "12"],
 370            ["Bits Allocated", "16"],
 371            ["High Bit", "11"],
 372            ["Pixel Representation", "0"],
 373            ["Manufacturer", "Philips Medical Systems"],
 374            ["Manufacturer Model Name", "Gyroscan Intera"],
 375            ["Sequence Variant", "OTHER"],
 376            ["Pixel Data", "gdcm::NotLoaded. Address:35846 Length:32768"] ] ],
 377       ["gdcm-MR-PHILIPS-16-NonRectPix.dcm",
 378          # Interest: - pixels are non rectangular
 379          #           - private elements are in explicit VR (which is normal)
 380          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"],  # Explicit VR, LE
 381            ["Modality", "MR"],
 382            ["Photometric Interpretation", "MONOCHROME2"],
 383            ["Rows", "64"],
 384            ["Columns", "160"],
 385            ["Bits Stored", "12"],
 386            ["Bits Allocated", "16"],
 387            ["High Bit", "11"],
 388            ["Pixel Representation", "0"],
 389            ["Manufacturer", "Philips Medical Systems"],
 390            ["Manufacturer Model Name", "Gyroscan Intera"],
 391            ["Pixel Spacing", "0.487416\\0.194966"],
 392            ["Pixel Data", "gdcm::NotLoaded. Address:5010 Length:20480"] ] ],
 393       ["gdcm-CR-DCMTK-16-NonSamplePerPix.dcm",
 394          # Interest: - Misses the "Samples Per Pixel" element which prevents
 395          #             e-film from reading it.
 396          #           - Misses the Manufacturer related info
 397          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.1"],  # Explicit VR, LE
 398            ["Modality", "CR"],
 399            ["Photometric Interpretation", "MONOCHROME2"],
 400            ["Rows", "750"],
 401            ["Columns", "750"],
 402            ["Bits Stored", "8"],
 403            ["Bits Allocated", "8"],
 404            ["High Bit", "7"],
 405            ["Pixel Representation", "0"],
 406            ["Implementation Version Name", "OFFIS_DCMTK_341"],
 407            ["Pixel Data", "gdcm::NotLoaded. Address:740 Length:562500"] ] ],
 408    ]
 409
 410    GdcmJpegFiles = [
 411       ["gdcm-JPEG-Extended.dcm",
 412          # Interest: Jpeg compression [Extended (2,4)]
 413          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.4.51"],
 414            ["Modality", "CT"],
 415            ["Photometric Interpretation", "MONOCHROME2"],
 416            ["Rows", "512"],
 417            ["Columns", "512"],
 418            ["Bits Stored", "12"],
 419            ["Bits Allocated", "16"],
 420            ["High Bit", "11"],
 421            ["Pixel Representation", "0"],
 422            ["Manufacturer", "SIEMENS"],
 423            ["Manufacturer Model Name", "Volume Zoom"],
 424            ["Pixel Data", "gdcm::NotLoaded. Address:2946 Length:192218"] ] ],
 425       ["MR-MONO2-12-shoulder.dcm",
 426          # Interest: Jpeg compression [Lossless, non-hierar. (14)]
 427          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.4.57"],
 428            ["Modality", "MR"],
 429            ["Photometric Interpretation", "MONOCHROME2"],
 430            ["Rows", "1024"],
 431            ["Columns", "1024"],
 432            ["Bits Stored", "12"],
 433            ["Bits Allocated", "16"],
 434            ["High Bit", "11"],
 435            ["Pixel Representation", "0"],
 436            ["Manufacturer", "Philips Medical Systems"],
 437            ["Manufacturer Model Name", "Gyroscan NT"],
 438            ["Pixel Data", "gdcm::NotLoaded. Address:1580 Length:718948"] ] ],
 439       ["gdcm-JPEG-LossLess3a.dcm",
 440          # Interest: - Jpeg compression [Lossless, hierar., first-order
 441          #             pred. 14, Select. Val. 1]
 442          #           - contains a sequence [circa (0008,2112)]
 443          #           - contains additional data after the pixels ???
 444          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.4.70"],
 445            ["Modality", "CT"],
 446            ["Photometric Interpretation", "MONOCHROME2"],
 447            ["Rows", "512"],
 448            ["Columns", "512"],
 449            ["Bits Stored", "12"],
 450            ["Bits Allocated", "16"],
 451            ["High Bit", "11"],
 452            ["Pixel Representation", "0"],
 453            ["Manufacturer", "SIEMENS"],
 454            ["Manufacturer Model Name", "Volume Zoom"] ] ],
 455          # FIXME: the pixel data looks like corrupted. As stated by "od -A d -c"
 456          #   0002528     à 177 020 \0  O  B \0 \0   ÿ  ÿ   ÿ   ÿ  þ  ÿ \0  à
 457          #   0002544    \0  \0  \0 \0  þ  ÿ \0  à 202  ? 003  \0  ÿ  Ø  ÿ  Ã
 458          # which we interpret as follows:
 459          #   - 2528: (à 177, 020 \0) = (0x7fe0, 0x0010) i.e. Pixel Data tag,
 460          #   - 2532: OB i.e. encapsulated encoded Pixel Data,
 461          #   - 2534: Two bytes reserved,
 462          #   - 2536: ÿ ÿ ÿ ÿ = 0xffffffff i.e. data element length is undefined.
 463          #   - 2540: (þ ÿ, \0 à) = (0xfffe, 0xe000) i.e. basic offset table
 464          #           item tag
 465          #   - 2544: \0 \0 \0 \0 i.e. offset table as length zero, i.e. no
 466          #           item value for offset table.
 467          #   - 2548: (þ ÿ, \0 à) = (0xfffe, 0xe000) i.e. item tag of first
 468          #           fragment of (encoded) pixel data.
 469          #   - 2552: 202 ? 003 \0 = 212866 bytes i.e. length of first fragment
 470          #           of (encoded) pixel data.
 471          # and so the next item tag should be found at the 2556+212866 th
 472          # octet i.e. at address 215422. But, as stated by od, we find:
 473          #   0215408   E   u Ö 026   Î   É 7 ¬ Ã ¸ ó ¿ ÿ Ù \f °
 474          #   0215424   ¶ 016 P   Ñ 002 016
 475          # and (\f °, ¶ 016) is NOT an OB item TAG which is required to be
 476          # present (at least there should be a sequence delimiter), refer to
 477          # PS 3.5-2001 page 50.
 478       ["gdcm-JPEG-LossLess3b.dcm",
 479          # Interest: - Jpeg compression [Lossless, hierar., first-order
 480          #             pred. 14, Select. Val. 1]
 481          #           - encoding is sligthly different from LossLess3a.dcm ???
 482          [ ["Transfer Syntax UID", "1.2.840.10008.1.2.4.70"],
 483            ["Modality", "CR"],
 484            ["Photometric Interpretation", "MONOCHROME2"],
 485            ["Rows", "2076"],
 486            ["Columns", "1876"],
 487            ["Bits Stored", "15"],
 488            ["Bits Allocated", "16"],
 489            ["High Bit", "14"],
 490            ["Pixel Representation", "0"],
 491            ["Manufacturer", "Philips Medical Systems"],
 492            ["Manufacturer Model Name", "Cassette Holder Type 9840 500 35201"],
 493            ["Pixel Data", "gdcm::NotLoaded. Address:3144 Length:4795668"] ] ],
 494    ]
 495
 496    def _BaseTest(self, FileSet):
 497       for entry in FileSet:
 498          fileName = os.path.join(GDCM_TEST_DATA_PATH, entry[0])
 499          toRead = gdcmHeader(fileName)
 500          valDict = toRead.GetPubElVal()
 501          for subEntry in entry[1]:
 502             element = subEntry[0]
 503             value   = subEntry[1]
 504             self.assertEqual(valDict[element], value,
 505                              ("Wrong %s for file %s (got %s, shoud be %s)"
 506                                % (element,fileName, valDict[element], value)) )
 507
 508    def testBarre(self):
 509       gdcmTestCase._BaseTest(self, gdcmTestCase.BarreFiles)
 510
 511    def testMultiFram(self):
 512       gdcmTestCase._BaseTest(self, gdcmTestCase.MultiFrameFiles)
 513
 514    def testFiles(self):
 515       gdcmTestCase._BaseTest(self, gdcmTestCase.GdcmFiles)
 516
 517    def testJpeg(self):
 518       gdcmTestCase._BaseTest(self, gdcmTestCase.GdcmJpegFiles)
 519
 520    def testWrite(self):
 521       import md5
 522       SourceFileName = os.path.join(GDCM_TEST_DATA_PATH,
 523                                     'gdcm-MR-PHILIPS-16-Multi-Seq.dcm')
 524       Source = gdcmFile(SourceFileName);
 525       Source.GetImageData()
 526       TargetFileName = "junk"
 527       Target = Source.WriteDcmImplVR(TargetFileName)
 528       Sign = '5af8739c15dd579dea223eb3930cacda'
 529       ComputeSign = md5.new(open(TargetFileName).read()).hexdigest()
 530       #print ComputeSign
 531       self.assertEqual(ComputeSign, Sign,
 532           ("Wrong signature for file %s (got %s, shoud be %s)"
 533            % (SourceFileName, ComputeSign, Sign)) )
 534
 535 if __name__ == '__main__':
 536    unittest.main()
 537