3
$\begingroup$

I have a data set on the form of 3 spatial coordinates x,y,z and an amplitude. Here is a small crop of the data, :

{{0, 0, 0, 1.00448}, {0, 4, 0, 2.98714}, {0, 4, 0.010101, 2.98712}, {0.010101, 4, 0.010101, 2.9871}, {0, 4, 0.020202, 2.98713}, {0.010101, 4, 0.030303, 2.98709}, {0, 4, 0.030303, 2.9871}, {0.010101, 0, 0.050505, 1.00449}, {0, 0, 0.050505, 1.0045}, {0.010101, 4, 0.050505, 2.9871}}

What I would like is a 2D-contour plot on the form x, y, amplitude, where the contour interface is represented by the values of amplitude within, say, 2.0 +/- 0.05. This is all I want, a black-colored contour on a white background.

First, what I do is to extract the data points that satisfy data[[:,3]] within 0 +/- 0.05, which I do using Select: Select[data, -0.05 < #[[3]] < +0.05 &]; Then I extract the points that satisfy that amplitude is within 2.0 +/- 0.05. This is easy, same procedure.

However, when I plot it using ListContourPlot, not only is the domain reduced, but the interface is still multicolored.

Is there a way to obtain a simple contour plot?

Improve this question
$\endgroup$
7
  • 1
    $\begingroup$ Am I reading correctly that you want basically one contour line, whose thickness is such that it would cover the area bounded by contours at 1.95 and 2.05? $\endgroup$ Commented Oct 21, 2015 at 9:55
  • $\begingroup$ @JasonB That is exactly what I want. That will allow be to plot any cross section of the data and include as many data amplitudes as I desire $\endgroup$ Commented Oct 21, 2015 at 9:56
  • $\begingroup$ I just noticed the dimensions of your data, that you have three dimensions and an amplitude but you want to plot only as a function of 2 - what do you want to do with the z coordinate? $\endgroup$ Commented Oct 21, 2015 at 10:24
  • $\begingroup$ @JasonB I extract the data points that satisfy data[[:,3]] within 0 +/- 0.05, I made a note of it in my OP $\endgroup$ Commented Oct 21, 2015 at 11:00
  • 1
    $\begingroup$ Ahhh, so it wasnt' clear what you were looking for. You want a 2D slice of your 3D contour plot I believe. If your data came from a function f[x,y,z], then you want a contour plot of f[x,y,2.0+delta] where delta is summed from plus/minus 0.05. Is that correct? $\endgroup$ Commented Oct 21, 2015 at 11:11

1 Answer 1

Reset to default
4
$\begingroup$

Here I'm just drawing two contour lines at the boundaries you set, and shading the region between them black and the rest of the image white. With more data it may look better and you won't have to use the InterpolationOrder option.

Try this:

ListContourPlot[data, InterpolationOrder -> 2, Contours -> {1.95, 2.05}, ContourShading -> {White, Black}, PlotRange -> All, ContourStyle -> Black]

enter image description here

Edit: You can set it to do multiple contours lines like this

data2 = Flatten[ Table[{x, y, .8 (x + 1)^3 + .4 (y + 1)^2}, {x, -4, 4, .1}, {y, -4,4, .1}], 1]; width = 0.1; ListContourPlot[data2, Contours -> Flatten@Table[n + width {-.5, .5}, {n, -30, 30}], ContourShading -> {White, Black}, PlotRange -> All, ContourStyle -> Black]

enter image description here

Improve this answer
$\endgroup$

Start asking to get answers

Find the answer to your question by asking.

Ask question

Explore related questions

See similar questions with these tags.