DIM3,
Z(..),
(:.)(..),
- toList
)
-import Values
---import Grid(make_grid, zoom)
+import System.Environment (getArgs)
+
+import Values (read_values_3d, write_values_1d)
+import Grid (make_grid, zoom)
mri_shape :: DIM3
-mri_shape = (Z :. 256 :. 256 :. 109)
+mri_shape = (Z :. 256 :. 256 :. 1)
main :: IO ()
main = do
- mridata <- read_values_3d mri_shape "./data/mridata.txt"
- let tmp = Data.Array.Repa.toList mridata
- print tmp
--- let g = make_grid 1 (Data.Array.Repa.toList mridata2)
--- let output = zoom g 2
--- print "Hello, world."
+ args <- getArgs
+ let color = head args
+ let in_file = "./data/MRbrain.40." ++ color
+ let out_file = "MRbrain.40." ++ color ++ ".out"
+ mridata <- read_values_3d mri_shape in_file
+ let g = make_grid 1 mridata
+ let output = zoom g (8,8,1)
+ write_values_1d output out_file