Hello everybody, I have an Ascii-grid (100x100m) with population density over Switzerland (local coordinations system, CH LV03). Now I would like to have a nice graphic to present the data. I would like to make a plot that will show the density as the z-variable and raise for example a stacking pillar showing how high the density is in that hectare. I attached a map that's similar to what I would like, to give you an idea incase I am not making myself clear. <http://r.789695.n4.nabble.com/file/n4658517/Bildschirmfoto_2013-02-14_um_10.08.51.png> Can you tell me how I can do this? I tried persp and persp3d, but there the x and y-vectors have to be increasing values, where as the coordinates in a grid don't necessarily increase for both x and y, as they follow the line horizontally... Your help would be greatly appreciated :), Thank you jacqueline -- View this message in context: http://r.789695.n4.nabble.com/3D-plots-of-2D-grids-tp4658517.html Sent from the R help mailing list archive at Nabble.com.
On Feb 14, 2013, at 1:10 AM, jas wrote:> > Hello everybody, > > I have an Ascii-grid (100x100m) with population density over > Switzerland > (local coordinations system, CH LV03). Now I would like to have a nice > graphic to present the data. I would like to make a plot that will > show the > density as the z-variable and raise for example a stacking pillar > showing > how high the density is in that hectare. I attached a map that's > similar to > what I would like, to give you an idea incase I am not making myself > clear. > <http://r.789695.n4.nabble.com/file/n4658517/Bildschirmfoto_2013-02-14_um_10.08.51.png > > > > Can you tell me how I can do this? I tried persp and persp3d, but > there the > x and y-vectors have to be increasing values, where as the > coordinates in a > grid don't necessarily increase for both x and y, as they follow the > line > horizontally... >I thought you said you had a 'grid' presumably with values at the intersections of the grid? That is the sort of input that persp will take. You would supply x=seq( ..., length=100), y= seq( ... , length=100) and your 'grid' (which I assume is stored as a matrix) as z values. In fact if you just supply a matrix as the z value to persp you will get an automatic display, just not with the labels you might have hoped for. -- David Winsemius, MD Alameda, CA, USA
calc_ext3.asc <http://r.789695.n4.nabble.com/file/n4658564/calc_ext3.asc> I supplied a small extent of the grid, that I am trying to visualize. It s a ascii-grid with 800 cells. My goal is to display the Grid-Values as the z-variable. The point that I dont understand, is what is my input as x and y-Variable in the persp-function? Do I have to make 8 vectors with 100 entries and plot them individually? Your help is greatly appreciated, jas pop.grid pop.mat <- as.matrix(pop.grid) # order is not right here already z <- pop.mat x <- 1:nrow(z) y <- 1:nrow(z) persp(x,y,z, phi=30, border=NA) # Error n persp.default(x, y, z, phi = 30, border = NA) : invalid 'z' Argument -- View this message in context: http://r.789695.n4.nabble.com/3D-plots-of-2D-grids-tp4658517p4658564.html Sent from the R help mailing list archive at Nabble.com.