X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FMRI.hs;h=5f2ceea483437bc7a9216f066e86e1d32867288d;hb=15ccf02b51c72a7bb6eb098c157005a78cd3f84d;hp=2dd6f6592e6e3bb7ac8d1eefbe73a4764e3cf787;hpb=65eed42ae2b9ca947bd7c08c2a102c1cf5dce8cc;p=spline3.git

diff --git a/src/MRI.hs b/src/MRI.hs
index 2dd6f65..5f2ceea 100644
--- 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
 
@@ -139,3 +140,43 @@ blue_dbl_data =
 z_slice :: Elt a => Int -> Array DIM3 a -> Array DIM2 a
 z_slice n arr =
   slice arr (Any :. n :. All :. All)
+
+
+transpose_zx :: Elt a => Array DIM3 a -> Array DIM3 a
+transpose_zx arr =
+  traverse arr
+           (\(Z :. zdim :. ydim :. xdim) -> (Z :. xdim :. ydim :. zdim))
+           (\_ -> (\(Z :. z :. y :. x) -> arr ! (Z :. x :. y :. z)))
+
+
+z_slice3 :: Elt a => Int -> Array DIM3 a -> Array DIM3 a
+z_slice3 n arr
+  | n == 0 = transpose_zx $ current R.++ next
+  | n == zdim-1 = transpose_zx $ previous R.++ current
+  | otherwise = transpose_zx $ previous R.++ current R.++ next
+  where
+    (Z :. zdim :. _ :. _) = extent 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
+
+
+write_values_chunk_to_bitmap :: Values3D -> FilePath -> IO ()
+write_values_chunk_to_bitmap v3d path
+  | zdim /= 3 = error "not a 3xMxN chunk"
+  | otherwise = write_values_slice_to_bitmap target_slice path
+  where
+    (Z :. zdim :. _ :. _) = extent v3d
+    target_slice = z_slice 1 v3d
+