Skip to main content
1 of 5
s0rce
  • 9.7k
  • 4
  • 47
  • 79

Here is my extremely slow volume rendering solution.

slices = Import["https://dl.dropbox.com/u/3730003/slicedata.wdx"]
data = Developer`ToPackedArray[Map[ImageData, slices]];

(*3D*)

dims = Dimensions@data;
coords = Table[{x, y, z}, {x, 1, dims[[1]]}, {y, 1, dims[[2]]}, {z, 1,
     dims[[3]]}];
cubes = {Hue @@ #1, EdgeForm[], Cuboid@#2} &
output = MapThread[cubes, {data, coords}, 3];
Graphics3D@output

(*2D*)

dims = Dimensions@data[[45]];
coords = Table[{x, y, 1}, {x, 1, dims[[1]]}, {y, 1, dims[[2]]}];
cubes = {Hue @@ #1, EdgeForm[], Cuboid@#2} &
output = MapThread[cubes, {data[[45]], coords}, 3];
Graphics3D@output
s0rce
  • 9.7k
  • 4
  • 47
  • 79