Prev: How to plot shapefile onto a regular x-y plot
Next: programmatically find nested vs. subfunctions
From: James Ramm on 23 Jun 2010 09:36
"Adrian McCallum" <abm35(a)cam.ac.uk> wrote in message <hov2so$b54$1(a)fred.mathworks.com>...
> "Mark Shore" <mshore(a)magmageosciences.ca> wrote in message <hotcfb$r9l$1(a)fred.mathworks.com>...
> > "Adrian McCallum" <abm35(a)cam.ac.uk> wrote in message <hot16s$h8k$1(a)fred.mathworks.com>...
> > > Apologies, I am new to MATLAB.
> > > I am trying to create a 3D image of borehole data.
> > > I have say 10 boreholes at different x,y locations, each with a data value at ~1000 different depth (z) values.
> > > How can I plot this in 3D and interpolate between the boreholes to define layering etc?
> > > I have been advised that I can use interp3, however I am unsure how I can easily read-in so many z data points and create multidemensional arrays.
> > > Any advice is appreciated, thank you.
> >
> > I strongly suspect that isosurface won't be that useful to you, given the strongly asymmetrical nature of resistivity in most geological environments which typically amount to an assemblage of different log-normal distributions. Even if your particular variation is unusually smooth, you'll have to consider down-sampling the borehole data to a grid size sufficiently small for a reasonable XYZ array size. And MATLAB probably uses an isometric interpolation method in 3D, which will not respect any real cross-hole structures.
> >
> > It may be worth a try, but there is a reason specialized (and horrendously expensive) tools such as Gocad exist.
>
> Hi Mark, thanks very much for that. I'll cx out Gocad else I might smooth/simplify the data a bit more before playing further. Thanks very much for your help. Cheers, Adrian.

Hi, if you are still having trouble with this, this may help:
For interpolation, I reccomend interpolating in 2D (1 for each layer) and using gstat or some other kriging package to do that. I have a rudimentary matlab interface to gstat which im still working on, but can email it to you if you like.
After gridding your layered data, then I reccomend you looking at
http://www.mathworks.com/matlabcentral/fileexchange/27965-surf3d
Which is a function I've created for visualizing layered data in 3D. Get in touch for help on things like thresholding.
Good luck!
 | 
Pages: 1
Prev: How to plot shapefile onto a regular x-y plot
Next: programmatically find nested vs. subfunctions