File: TestOrientation.py

package info (click to toggle)
gdcm 3.0.21-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 26,880 kB
  • sloc: cpp: 203,477; ansic: 78,582; xml: 48,129; python: 3,459; cs: 2,308; java: 1,629; lex: 1,290; sh: 334; php: 128; makefile: 117
file content (56 lines) | stat: -rw-r--r-- 1,633 bytes parent folder | download | duplicates (5)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
############################################################################
#
#  Program: GDCM (Grassroots DICOM). A DICOM library
#
#  Copyright (c) 2006-2011 Mathieu Malaterre
#  All rights reserved.
#  See Copyright.txt or http://gdcm.sourceforge.net/Copyright.html for details.
#
#     This software is distributed WITHOUT ANY WARRANTY; without even
#     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
#     PURPOSE.  See the above copyright notice for more information.
#
############################################################################

import gdcm
import os,sys

if __name__ == "__main__":
  success = True

  ori1 = (1,0,0,0,1,0)
  ori2 = (1,0,0,0,0,1)
  ori3 = (0,1,0,0,0,1)

  label1 = gdcm.Orientation.GetLabel( gdcm.Orientation.GetType( ori1 ) )
  if label1 != 'AXIAL':
    print("Found:",label1)
    success = False
  label2 = gdcm.Orientation.GetLabel( gdcm.Orientation.GetType( ori2 ) )
  if label2 != 'CORONAL':
    print("Found:",label2)
    success = False
  label3 = gdcm.Orientation.GetLabel( gdcm.Orientation.GetType( ori3 ) )
  if label3 != 'SAGITTAL':
    print("Found:",label3)
    success = False

  image = gdcm.Image()
  image.SetNumberOfDimensions(2)
  print(image)
  print(image.GetDimensions())
  print(image.GetOrigin())
  print(image.GetSpacing())
  print(image.GetDirectionCosines())

  image.SetNumberOfDimensions(3)
  image.SetDimensions( (512,256,128) )
  print(image.GetDimensions())

  #dircos = gdcm.DirectionCosines( (1,0,0,0,0,-1) )
  dircos = gdcm.DirectionCosines()
  print(dircos)
  #print dircos.Cross()

  # Test succeed ?
  sys.exit(success == False)