module Main
where

import Data.Array.Repa (
  DIM3,
  Z(..),
  (:.)(..),
  )

import System.Environment (getArgs)

import Grid (make_grid, zoom)
import Values (read_values_3d, write_values_1d)

mri_shape :: DIM3
mri_shape = (Z :. 256 :. 256 :. 1)




main :: IO ()
main = do
  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 (4,4,1)
  write_values_1d output out_file