Prev: getting subimages
Next: PDF as image
From: James Ramm on 10 May 2010 09:11
I'm a little stuck with a problem;
I have x,y,z data with associated data values R. R is arranged as layers of thickness T, for each xyz point. For example, a simplified data file may look like this:

X Y Z R1 R2 R3 T1 T2 T3
794843.2 9915493 100 55 70 500 5 10 14
794826.1 9915463 90 60 80 545 5 10 14
794798.5 9915441 120 76 74 600 5 10 14
794765.3 9915430 75 45 100 560 5 10 14

Where Rn is the data value of layer n at the point described by x,y,z, down to the bottom of the layer at this point. Layer n has an associated thickness, Tn.

This means that, e.g R2 will describe the data value of layer 2 (which has a thickness T2), for the range Z-T1 to Z-(T1+T2).

I would like to express this data as a volumetric grid, which can then be exported for visualisation/volume rendering (in a program such as paraview).

I know meshgrid can be used to set up a structured grid using the co-ordinate limits,
But how can I grid R across this, taking into account the layer thicknesses?

For example, if I define a 3D grid which has a Z spacing of 1, and the first layer has a thickness N, not only am I doing a 'normal' gridding of R (e.g with griddata3), I also need every point in the z-direction down to the thickness of the first layer (in this example that would be 5 points) to have the data value for R1 at this 5 points, then R2 at the points underneath (again, for the whole thickness of the layer) and so on and so on...

Thanks for your help
 | 
Pages: 1
Prev: getting subimages
Next: PDF as image