src/Main.hs

   1 module Main
   2 where
   3
   4 import Data.Array.Repa (
   5   Array,
   6   DIM1,
   7   DIM3,
   8   Shape,
   9   Z(..),
  10   (:.)(..),
  11   index,
  12   reshape,
  13   toList
  14   )
  15
  16 import Data.Array.Repa.IO.Vector (readVectorFromTextFile)
  17
  18 --import Grid(make_grid, zoom)
  19
  20
  21 read_mri_data :: IO (Array DIM1 Int)
  22 read_mri_data = readVectorFromTextFile "./data/mridata.txt"
  23
  24 mri_shape :: DIM3
  25 mri_shape = (Z :. 256 :. 256 :. 109)
  26
  27 main :: IO ()
  28 main = do
  29   mridata <- read_mri_data
  30   let mridata2 = reshape mri_shape mridata
  31   let tmp = Data.Array.Repa.toList mridata2
  32   print tmp
  33 --  let g = make_grid 1 (Data.Array.Repa.toList mridata2)
  34 --  let output = zoom g 2
  35 --  print "Hello, world."