Re: ListContourPlot for irregular sampling
- To: mathgroup at christensen.Cybernetics.NET
- Subject: [mg251] Re: [mg218] ListContourPlot for irregular sampling
- From: "Roberto Abraham at the Institute of Astronomy, Cambridge" <abraham at mail.ast.cam.ac.uk>
- Date: Tue, 29 Nov 1994 18:43:16 +0000 (GMT)
Hi Stefan,
I've encountered situations similar to yours on a number of occasions
(it's a common problem in astronomy), and I've found a number of
solutions. The one that seems to work best for me is to write a small
function has as its arguments the x,y position you want to sample and the
list of data. The function simply takes the list of data and returns a
value at the x,y position that's the weighted average of the nearest n
data points, where each point is weighted according to the square of its
distance from to the point at which you're evaluating the function. I've
found that using the nearest 3-5 neighbors seems to work well, depending
on how sparse your data is, but you should experiment with different
values. The final function can be plotted up using ContourPlot or Plot3D
or whatever. If you need to evaluate the function really quickly then
it's easy to resample the grid evenly and define function values using
the Interpolation function.
If you can wait a couple of days I can send you a package with some code
(I've just moved to Cambridge and all my Mathematica stuff is still on my
Mac at home, rather than on my machine at work). However, it's very easy
to code up the algorithm.
Cheers,
Bob
--------------------------------------------------------------------
Dr. R. G. Abraham abraham at mail.ast.cam.ak.uk
Institute of Astronomy
Cambridge University
Madingley Rd., Cambridge
CB3 OHA
United Kingdom
-------------------------------------------------------------------
On Wed, 23 Nov 1994, Dr. Stefan P. Mueller wrote:
> I need to generate contour plots of data which are not sampled on an evenly
> spaced Cartesian grid. ListContourPlot ony accepts data in the form of an
> nxm matrix of heights z which are evenly sampled, i.e.
>
> data = {z11,z12,z13, ...,z1n},
> {z21,z22,z23, ...,z2n},
> ...
> {zm1,zm2,zm3, ...,zmn}};
>
> whereas my data have the following structure
>
> data = {{x1,y1,z1},{x2,y2,z2}, ..., {xm,ym,zm}};
>
> i.e. each height is also assciated with its own x and y coordinate.
> Does anybody know about a package to plot data with this structure?
>
> Thank you in advance!
>
> Stefan
> ____________________________________________________________________________
> Dr. Stefan P. Mueller, Nuklearmedizin, Universitaetsklinikum, D-45122 Essen,
> FR Germany, e-mail: onm010 at aixrs1.hrz.uni-essen.de, FAX: <49> (201) 723-5964
>
>
>