@@ -51,19 +51,23 @@ def main(args)-> None:
51
51
Lower = np .array (img .TransformPhysicalPointToContinuousIndex (Lower )).astype (int )
52
52
Upper = np .array (img .TransformPhysicalPointToContinuousIndex (Upper )).astype (int )
53
53
54
+ for i in range (3 ):
55
+ if Lower [i ] > Upper [i ]:
56
+ Lower [i ], Upper [i ] = Upper [i ], Lower [i ]
54
57
# Bounds checking
55
58
img_size = img .GetSize ()
56
59
Lower = [max (0 , l ) for l in Lower ]
57
60
Upper = [min (img_size [i ], u ) for i , u in enumerate (Upper )]
58
-
59
61
60
- # Ensure non-zero size for all dimensions
61
- for i in range (3 ):
62
- if Lower [i ] == Upper [i ]:
63
- if Upper [i ] < img_size [i ] - 1 :
64
- Upper [i ] += 1
65
- elif Lower [i ] > 0 :
66
- Lower [i ] -= 1
62
+
63
+ # # Ensure non-zero size for all dimensions and that lower < upper
64
+
65
+ # for i in range(3):
66
+ # if Lower[i] == Upper[i]:
67
+ # if Upper[i] < img_size[i] - 1:
68
+ # Upper[i] += 1
69
+ # elif Lower[i] > 0:
70
+ # Lower[i] -= 1
67
71
68
72
# Crop the image
69
73
crop_image = img [Lower [0 ]:Upper [0 ],
@@ -84,14 +88,18 @@ def main(args)-> None:
84
88
85
89
os .makedirs (os .path .dirname (ScanOutPath ), exist_ok = True )
86
90
87
- # try:
91
+ try :
88
92
89
- sitk .WriteImage (crop_image ,ScanOutPath )
93
+ sitk .WriteImage (crop_image ,ScanOutPath )
90
94
91
- # except:
92
- # import sys
93
- # print("Error for patient: ",patient)
94
- # print('The error says: ',sys.exc_info()[0])
95
+ except :
96
+ import sys
97
+ print ("Error for patient: " ,patient )
98
+ print ('The error says: ' ,sys .exc_info ()[0 ])
99
+ print ('Lower: ' ,Lower )
100
+ print ('Upper: ' ,Upper )
101
+ print ('Lower[2]:' ,Lower [2 ])
102
+ print ('Upper[2]:' ,Upper [2 ])
95
103
96
104
with open (args .logPath ,'r+' ) as log_f :
97
105
log_f .write (str (index ))
0 commit comments