search for: contourlin

I want to use contourLines() to get contour line coordinate vectors, but I don't want to make a plot. However contourLines() insists on opening a graphics device. Is there a way tell it not to do this? > version platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status Patch...

...ur plots (of a spatially-smooted density produced by kde2d()), with very satisfactory results. I now want access to the coordinates of the points on the contours, and it would seem that contour() does not return a value, so there is nothing from which these could be extracted. However, apparently contourLines() does, and it seems to be invoked in essentially the same way (so far as the x,y,z and levels parameters are concerned). However, I am not sure whether the contour lines generated by contourLines() are the same numerically as those generated (at any rate plotted) by contour(). By over-plotting...

Hello, I have searched for documentation on the function contourLines's algorithm but cannot find a thing. I am about to submit a paper to a journal but cannot yet do so because I need to provide some reference for this function. Does anyone know what algorithm is used for this function? Thanks, Derek Lacoursiere -- View this message in context: http://www...

A follow-up to my previous query. Say one of the results returned by contourLines() is C.W <- contourLines(....) Then C.W is a list of (in this case 28) lists, each of which is a list with components $level (a single number), $x (the x coords of the points on a contour at that level) and $y (the y coords). I can of course get these individually with, for the 5th one fo...

Hello everybody, I search a function which returns the contour of map with levels like contourLines, but I would like this function return the border of the map too, because the function contourLines cannot consider the corner of the map and it is not adapted to fill polygon after that. Thanks in advance Pierre Bruyer

Full_Name: Dominik Heide Version: 2.3.1 OS: Linux (Suse 9.3) Submission from: (NULL) (134.76.220.200) The contour functions has sommewhere (I think in .../main/plot3d.c) a maximum number of line segments that is set using a define. The contour line in my data too long for this and therefore not correct analysed (I get the message 'contour(): circular/long seglist -- bug.report()!').

Hello, Is it possible to draw level lines on the surface of a bivariate function? In the following example, to draw surface and levels lines for a multivariate normal law, I use persp, trans3d, contourLines and lines, but if the lines are correctly drawn, some parts of them are, of course, visible even if they are drawn on a non visible "face". Any suggestion to avoid this problem ? Thank you Etienne Example : trans3d <- function(x,y,z, pmat) { tr <- cbind(x,y,z,1) %*% pmat...

I have a list of polygons generated by the contourLines() command (each object of the list is a list in itself with two objects: a vector of x values, and a vector of y values for each vertex). I wish to convert that list into a gpc.poly object of multiple contours. How do I do this? gpclib apparently has no method of coercing lists into the gpc.po...

...from list of colors according to its z-value. (The range should be from blue (z=0) to red (z=1).) The grid should then be omitted. [I have seen "terrain.colors", but don't know how to use it for this purpose.] (2) Add two lines to the surface for, say z=0.8 and z=0.3. [Can contour or contourLines be used?] --- x <- seq(0, 1, length = 50) y <- x f <- function(x,y) { sin ((1-x)*y) } z <- outer(x,y,f) persp(x, y, z, theta = 30, phi = 30, shade = 0.3, col = "red" ) --- Finally, I would also produce a flattened 2D map of the same function, i.e. a map in which ea...

Hi everybody, I want to add a single contourline to a levelplot. While everything works fine if the value at which the line should be drawn is positive, there is an error if the value is negative: library(lattice) my.panel <- function(..., at, contour=FALSE, labels=NULL) { panel.levelplot(..., at=at, contour=contour, labels=labels)...

..., and I don't want to evaluate it at a dense grid of points (i.e. I don't want a huge array of values). Say I have a vector of x-values such as x <- 0.1*c(1:10), and the same for y <- 0.1*(0:10). I then evaluate a 10*10 matrix z of values of z = f(x,y). so I can then do CL <- contourLines(x,y,z) and get a list of contour-lines, each list containing the x and y coordinates. When I look at the plot from contour(x,y,z), I see that the lines are (as expected) made up of straight-line segments. But the break-points between the segments seem to either have x-coords on the x-values in...

...NULL,ylab= NULL,col.regions=heat.colors(100)[100:1]),split=c(1,1,1,1),more=T) and then to add the contour plot (I only want a contour at level 5): print(levelplot(d~x+y,data=t,cuts=20,scales=list(draw=F),xlab=NULL,ylab= NULL,contour=T,at=c(0,5),reg=F,labels=F),split=c(1,1,1,1)) this shifts the contourline to the right of the levelplot. If I use region=T then it looks better but it is still shifted slighlthy, plus have an extrakey overlayed to the original key with different numbers on it. If I draw the first plot, but only put one contourline on it, I loose all the other colors...Any ideas...

I'm trying to make a set of contour plots of bivariate kernel density estimates, showing three such plots overlaid, similar to this plot http://euclid.psych.yorku.ca/SCS/Private/Test/ridge-boot2.pdf except that I would like to have the contours *filled* (using transparent colors). To make this reproducible, I've saved the results of KernSmooth::bkde2D() in the following file:

Hi all, I'm trying to draw a contour plot with rounded (smoothed) contour lines instead of sharp angles; when the grid consists of only a few points in x- and y- axis, the resulting contour is in facts rather ugly since very sharp angles may appear. I did not find any way to do it, by using either "contour" or "contourplot" (from the lattice package), I wonder if there

...e image matrix but I can't see how to over come this other than loading in a second image that has been rotated 90 degrees and plotting the contours on the non-rotated one (obviously a bad solution) 2) I want to be able to select a single contour after they have been drawn. I know there is contourLines but I can't see how to relate this to the contours plotted. In matlab there is a function to label contours so that they can be called by their contour label in order to get the coordinates of the contour. Can this be done in R. Sorry if these are easy / stupid questions. I'm new to a...

...ce) cloud(z~x+y|g,data=d,layout=c(2,1), type='h', lwd=3, par.box=list(lty=0), scales=list(z=list(arrows=F,tck=0)), panel.3d.cloud = function(x, y, z,rot.mat, distance, zlim.scaled, nlevels=20,...){ add.line <- trellis.par.get("add.line") clines <- contourLines(surf[[packet.number()]],nlevels = nlevels) for (ll in clines) { m <- ltransform3dto3d(rbind(ll$x-.5, ll$y-.5, zlim.scaled[1]), rot.mat, distance) panel.lines(m[1,], m[2,], col = add.line$col, lty = add.line$lty,...

...up the conour. Maybe, there is something I miss in here. What went wrong there? Also, is there any diffeence between contourplot(lattice) and counterLine(base)? These 2 functions seem to give difference contour from the same data set. contourplot(lattice) seem to give better looking contour than contourLine() or contour(). ps. I krige 1700+ simulated observations into 4000+ regular spaced data to get contour. Thanks in advance, simon [[alternative HTML version deleted]]

Hi all, I used contour() to add contour lines to a plot. Now I?m wondering if there is a way to get an output of the calculated x- and y- coordinates of the contour lines? Tonja

Hello, I'm trying to assess the similarity of two lines that are represented as points (output of differential equation solvers). Is there a function or a package that deals with things like this? Thanks, John. ==================================================== Johannus (John) Janmaat Assistant Professor of Economics Barber School of Arts and Sciences University of British Columbia -

...39;CobbGRAF.tex', standAlone = TRUE, width=5, height=5) # La gr?fica de la funci?n de utilidad f = function(x, y) ((y)^1*(x)^1) x = seq(0,5,len=40) y = seq(0,5,len=40) z = outer(x, y, f) showsurface = function(x, y, z) persp3d(x,y,z, col="blue", alpha=0.3, axes= F)+{ contours = contourLines(x,y,z) for (i in 1:length(contours)) { with(contours[[i]], lines3d(x, y, level, col="darkred")) } } open3d() showsurface(x,y,z) # Cierro el device dev.off() # Compilo el archivo tex tools::texi2dvi('CobbGRAF.tex',pdf=F) -- View this message in context: http://r.78969...