search for: surface3d

...=TRUE) # Error - Increasing 'x' and 'y' values expected ?? # Next, tried to get the axes with the plot3d command plot3d(mvol2$X1,mvol2$X2,mvol2$value,type="n",xlab="x1",ylab="x2",zlab="value") # then get the surface - why the x,z, y sequence? surface3d(mvol2$X1,vol2,mvol2$X2) I get the following error message for the surface3d command : Error in rgl.surface(x = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, : Bad dimension for cols Additionally, I want points and lines plotted. For example, the points (1,1,200) and (51,2,232). I don't s...

...es: x,y,gridvalue,height. The dataframe has 2.5mio observations (ie grid points), I assign colors through the gridvalue using map_color_gradient thus producing: x,y,gridvalue,height,gridcol as variables of the dataframe. The grid dimensions are 1253 x 2001 (=2507253 data points). My attempts with surface3d fail, mainly because I cannot produce the matrix required for the height input. elev.to.list{CTFS} fails with: "Error in matrix(elevfile$elev, nrow= 1+ydim/gridsize, ncol=1+xdim/gridsize. : attempt to set an attribute on NULL" which I assume means it requires a square grid (=quadrates)....

Hello, Please could someone advise if it's possible to save the graphical output from the surface3d() function? I have tried the dev.copy() function to save as a pdf but an error message says I cannot copy from the null device. Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax: +44 113 343 6716 mail: laura at...

Hi R users, Does anyone know how to change the size of 3d charts? I'm using surface3d in rgl package, opening a new window each time to display the chart. I want it so that the chart fills the whole window, because when I output it to png, I don't want all the white space around the chart (right now, i'm getting this white "border" around the chart because the char...

...of the data. I typically have something like this: X Y Z 0.005 0.023 34.45 0.0035 0.63 28.48 . I've tried looking at the persp and wireframe packages, and the rgl package, but I can't seem to figure out how to use any of them. I tried to take the rgl.surface3d example and use it for myself, and this is what I have so far: setwd(".") data <- read.csv("data.csv",header=T) x <- data$X y <- data$Y z <- data$Z open3d() surface3d(x, y, z) That gives me the following error: Error in rgl.surface(x = c(0...

Hi, I'm using rgl to generate a 3D surface plot and I'm struggling to get the lighting correct. Currently the surface gets plotted, but is very 'shiny'. On rotating the view, I get to see parts of the surface - but overall I don't see much detail because of the spotlight like lighting. I've played around with the specular, ambient and diffuse but I can't

Dear folks, I have created a plot on RGL device :   x = 1:6 y = seq(-12, 5, by=1) z = matrix(0, length(y), length(x)) z[13,3] = 1; z[13,4] = 1.011765 surface3d(x, y, t(z), col=rainbow(1000)) grid3d(c("x-", "y-", "z"))   Now I want to draw 2 lines along x=3 & x=4, over the surface (with different colour). Could somebody help me how to draw that?   Thanks, [[alternative HTML version deleted]]

Dear R-experts, Using persp3d(), I plotted a 3d surface. I would like to colour this surface based on the result values (dependent variable) grouped in, say, three to five different ranges. Later, I would also like to paint the same surface based on the range-values of one of the independent variables. Please help. Thank you. Regards, Suresh -- View this message in context:

...(1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) zlim <- range(y) zlen <- zlim[2] - zlim[1] + 1 colorlut <- terrain.colors(zlen) # height color lookup table col <- colorlut[ z-zlim[1]+1 ] # assign colors to heights for each point #open3d() surface3d(x, y, z, color=col, back="lines", alpha=0.75) @ \end{figure} \end{document}

Hi, I want to make a surface3D plot of a landscape. I have cordinates (x, y, z) recorded with a GPS. The datapoints are not evenly distributed within the rectangular area. To do a fast 3D plot I used following. > library(grid) > library(lattice) > v <- read.table("clipboard") > names(v) <- c("x...

...ns? g.pca <- prcomp(RtG , retx = TRUE) scale <- 0.4 bg3d("black") amp <- 6 zlim <- amp*range(z) zlen <- (zlim[2] - zlim[1])*100 colorlut <- terrain.colors(zlen) # height color lookup table col <- colorlut[ ((z*amp)-zlim[1])*100+1 ] # assign colors to heights surface3d(x,y,z*amp,color=col) for (i in 1:30) { for (i in 1) { pc <- g.pca$rotation[,i] pc.pos <- (pc > 0) if (any(pc.pos)) { spheres3d(xval[pc.pos],yval[pc.pos],6*zsamp[pc.pos], radius=scale*pc[pc.pos],color="red") } if (any(!pc.pos)) { spheres3d(xval[!pc.pos],yval[!pc.p...

Hi All- I would like to create a 3d topographic map using lat/lon and z(height). I have been scouring the R help pages and have not located the package I am looking for. Does anyone have a suggestion of package that will work for this? thanks- sherri

Dear all, I would like to put two perspective plots into one plot. The help page for ?persp shows how one can add points and lines but not another perspective plot. data(volcano) z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) ## Don't draw the grid lines : border = NA persp(x,

Hello, I requested help a couple of weeks ago creating a dipole field in R but receieved no responses. Eventually I opted to create a 3d sinusoidal plot and concatenate this with its inverse as a means for a "next best" situation. It seems that this isn't sufficient for my needs and I'm really after creating a continuous 3d gaussian mesh with a "positive" and

...ame pattern: XXX3d() The functions map to rgl functions internally (with some exceptions). Here's a short list of interesting spots in the interface API: * par3d() can be used to get/set values on the current RGL device. * Primitives functions such as points3d(), lines3d(), surface3d(), ... * First bunch of generic transformation utilities, rendering methods and mesh objects are implemented. The rgl interface using the pattern "rgl.XXX()" is still available. You can decide between more OpenGL/RGL'is low-level functions or more R'ish function...

...ame pattern: XXX3d() The functions map to rgl functions internally (with some exceptions). Here's a short list of interesting spots in the interface API: * par3d() can be used to get/set values on the current RGL device. * Primitives functions such as points3d(), lines3d(), surface3d(), ... * First bunch of generic transformation utilities, rendering methods and mesh objects are implemented. The rgl interface using the pattern "rgl.XXX()" is still available. You can decide between more OpenGL/RGL'is low-level functions or more R'ish function...

Hello all, I use persp3d from the rgl-package to plot a sruface. The typical call is persp3d(x, y, z) With cooridinate-vectros x, y and a function-values matrix z. Now I have different z's, say z_1,...,z_n Question: Is it possible to generate an animation from a sequence of such calls, for different z's? I would like to see how the surface is beeing changed in the time. Thank you

Hi, I would like to plot a matrix which comprises 3 columns. So first column should be the x-axis. Second the y-axis and the third one should be the z-axis. I know that I can plot such data using the scatterplot3d, but I would like to have a surface plot like you can do with wireframe and persp. But both, persp and wireframe, want to have x and y as descending data, at least persp. I am not sure

Dear R-helpers: I, an entry level R user, wonder how make a 3d scatter plot with both error bars and a flexibly fitted surface. Can anyone eligthen me? Many Thanks in advance. -Sean [[alternative HTML version deleted]]

...im[1], to = ylim[2], length.out = resolution) Z <- outer(X, Y, f) ## Create 3d scatter plot and add the regression surface open3d() with(data = data, plot3d(x = vars[2], y = vars[3], z = vars[1], xlim = xlim, ylim = ylim, zlim = zlim)) par3d(ignoreExtent = TRUE) surface3d(X, Y, Z, col = "blue", alpha = .6) par3d(ignoreExtent = FALSE) return(summary(my.model)) } -- Joshua Wiley Ph.D. Student, Health Psychology University of California, Los Angeles http://www.joshuawiley.com/