X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=util%2Fview-mri-data.py;fp=util%2Fview-mri-data.py;h=01bddbdcfc8d91aa24f6ae3156a5d57ccb985a69;hb=b306f322ec05410a45155be3139d6824824aadca;hp=62ae9b77e8c8302f82240491509d59bcc4d6d3f4;hpb=42830c416e147bfe828863b2eb7fa12940ca2ad5;p=spline3.git

diff --git a/util/view-mri-data.py b/util/view-mri-data.py
index 62ae9b7..01bddbd 100755
--- a/util/view-mri-data.py
+++ b/util/view-mri-data.py
@@ -33,18 +33,19 @@ def cube_root(n):
 import numpy as np
 data = np.fromfile(sys.argv[1], dtype='>u2')
 
-original_data_length = 7143424
+original_data_length = 109*256*256
 multiplier = cube_root(len(data) // original_data_length)
 data.shape = (109*multiplier, 256*multiplier, 256*multiplier)
 data = data.T
 
 # Display the data
-mlab.figure(bgcolor=(0, 0, 0), size=(1000, 1000))
+mlab.figure( bgcolor=(0,0,0), size=(1000,1000) )
 
+# Our data is scalar (grayscale), but not equally spaced in all
+# directions.
 src = mlab.pipeline.scalar_field(data)
-
-# Our data is not equally spaced in all directions:
 src.spacing = [1, 1, 1.5]
 src.update_image_data = True
-v = mlab.pipeline.volume(src, vmax=2500, vmin=1400)
+
+mlab.pipeline.volume(src, vmin=1400, vmax=2500)
 mlab.show()