I'm using the MCNP transport code to create the data. I'm using mathematica to analyze it. In my case I used MCNP to track 100 keV electrons through a material and score the energy deposited in each cell. My primary goal was to find the isodose curves at any slice of the data set, which I was able to do. I was hoping to find some interesting results with image3d and was greeted with a red blob.

As far as what the data is like? I'm not sure if I understand your question. Could you be specific? The data is pretty much equivalent to any study on electron beams, so I can only refer you to those. Here's one that shows electron tracks. https://www.researchgate.net/publication/269112923_YW-PolymerPhysics-2013/figures?lo=1

My exact mathematica code I used to create the plots are:

Manipulate[ ListContourPlot[totaldose[[1 ;; 100, 1 ;; 100, m]], PlotRange -> All, Contours -> 40, ColorFunction -> "TemperatureMap"], {m, 1, 100, 1}]

logged = Table[Table[Table[ If[ totaldose[[i, j, k]] == 0, 0, Log10[totaldose[[i, j, k]]] ] , {i, 1, 100}], {j, 1, 100}], {k, 1, 100}]

Image3D[logged]

"totaldose" is my data set. I created the "logged" term because I wasn't sure if the issue was the huge difference in minimum and maximum values (from 10 to 10^7). It didn't seem to matter because the results were essentially identical as if I had just used "totaldose" with image3d. I tried using the slider options on the image to adjust my output. I'm not familiar with image3d.

Dear Henrik I have a one .txt file with 4 columns (not .mx file) i want to plot my data same as your plotting. with Image3D and ListSliceCountorPlot3D. can you help me more. the files size is more than 20 MB and i could not attached it.( my mesh cells are 100 (x) *100 (y) * 200 (z) ). The data with 4 columns and 2000,000 rows.

Best regards

Saeid

Dear Saeid,

to give an unspecific answer to your unspecific question: To begin with I would try to create a 3D array from your data. To be able to give more help, one needs at least a small but representative sample of your data.

Regards -- Henrik

Dear Saeid,

the first three rows of your data appear to be the coordinates and are Integers - excellent! (This is why I wished to see a sample!) Try this:

ClearAll["Global`*"]

textData = Import[" < your file > ", "Table"];
(* dimenstions taken from the last row: *)
dims = Most[Last@textData] + 1
(* generate empty array: *)
data3D = ConstantArray[0, dims];
(* values are filled in: *)    
Set[data3D[[#1 + 1, #2 + 1, #3 + 1]], #4] & @@@ textData;

Attributes[logScale] = Listable;
logScale[data_] := Piecewise[{{0, data < 0.001}}, Log[1. + data]];

logData3D = logScale[data3D];

Now you can do e.g.:

Regards -- Henrik