Add the round_array function and use it in main3d.

[spline3.git] / src / MRI.hs

diff --git a/src/MRI.hs b/src/MRI.hs
index 52be3f7b1cba92c96d9bb64142ed7b3dd4acb546..510e4c011191f8cc78c99fac464a04175a8ed9af 100644 (file)
--- a/src/MRI.hs
+++ b/src/MRI.hs
@@ -6,6 +6,7 @@ import Data.Bits
 import Data.Array.Repa                  as R
 import Data.Array.Repa.IO.Binary        as R
 import Data.Array.Repa.IO.ColorRamp     as R
+import Data.Array.Repa.IO.BMP           as R (writeComponentsToBMP)
 
 import Values
 
@@ -69,11 +70,15 @@ swap_bytes :: (Shape sh) => (RawData sh) -> (RawData sh)
 swap_bytes arr =
     R.force $ R.map flip16 arr
 
+
 bracket_array :: (Shape sh) => (Values sh) -> (RawData sh)
 bracket_array arr =
-  R.force $ R.map f arr
-  where
-    f = bracket
+  R.force $ R.map bracket arr
+
+
+round_array :: (Shape sh) => (Values sh) -> (RawData sh)
+round_array arr =
+  R.force $ R.map round arr
 
 
 flip_y :: RawData3D -> RawData3D
@@ -158,3 +163,14 @@ z_slice3 n arr
     previous = transpose_zx $ reshape mri_slice3d (z_slice (n-1) arr)
     current = transpose_zx $ reshape mri_slice3d (z_slice n arr)
     next = transpose_zx $ reshape mri_slice3d (z_slice (n+1) arr)
+
+
+write_values_slice_to_bitmap :: Values2D -> FilePath -> IO ()
+write_values_slice_to_bitmap v3d path =
+  R.writeComponentsToBMP path routput goutput boutput
+  where
+    arr_bracketed = bracket_array v3d
+    colors = values_to_colors $ R.map fromIntegral arr_bracketed
+    routput = R.map (\(red, _,     _)    -> red)   colors
+    goutput = R.map (\(_,   green, _)    -> green) colors
+    boutput = R.map (\(_,   _,     blue) -> blue)  colors