Export contourlines to serve them through WMS

Question

I am trying to do automatic mapping. I want to have a uncertainty map , which uncertainties are displayed by isolines. I want to export these isolines and the values associated with them and then display them through web map service. Is it possible? by following code you can produce the contour lines

require(gstat)
require(lattice)
data(meuse)
data(meuse.grid)
coordinates(meuse) = ~x+y
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m)
l=as.data.frame(x)
contourplot(var1.var ~ x+y,data=l,cuts=10)

Why convert contour LINES to polygons? Why not serve them as lines through your web map service? — Spacedman, Jan 21 '13 at 11:47

score 1 · Answer 1 · answered Jan 21 '13 at 10:35

1

This is not a direct answer for your question, but maybe it can point you in the right direction:

you can generate an object from your contourplot

CP<- contourplot(var1.var ~ x+y,data=l,cuts=10)

And try to look at its structure

fix(CP)

The data about isolines is stored here:

x <- CP$panel.args.common[1]$x
y <- CP$panel.args.common[2]$y
z <- CP$panel.args.common[3]$z

But I cannot figure out now how to extract it. Maybe this can help you a bit.

answered Jan 21 '13 at 10:35

Andrey Dyachenko

551
5
8

1

Dear Andrey, thanks for your answer But these values are not contour values. These are the out put of object x. – Saba Eshraghi Jan 21 '13 at 11:46

Export contourlines to serve them through WMS

1 Answers1