search for: gridtopology

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20080203/5ea1e3ad/attachment.pl

Is there any R software that create an image from Yahoo maps together with points of known UTM coordinates (or lat/long marked? Note that my region of interest is not covered in sufficient detail by Google maps. It actually does not have to be Yahoo maps as long as it has sufficient coverage of my region but that's the one I have found with coverage of my region. The scale I am interested in

...found for interfacing between GRASS 6 and R (library (spgrass6)). I understood the following (I already managed to configured properly the mapset and so on): 1. I have to create a SpatialGridDataFrame object before to write into a Raster file. 2. In order to create a SpatialGridDataFrame I need a GridTopology. 3. After doing it, I am allow to write the file into GRASS using the function writeRAST6, like: writeRAST6 (spterrain, "maps.dem"). How to create the SpatialGridDataFrame and/or the GridTopology? I am doing something wrong. Thanks, Javier Hidalgo. _____________________________________...

...defs")) x2 <- readShapeSpatial("ptna2", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) bb <- bbox(rsa2) cs <- c(1, 1) cc <- bb[, 1] + (cs/2) cd <- ceiling(diff(t(bb))/cs) rsa2_grd <- GridTopology(cellcentre.offset = cc,cellsize = cs, cells.dim = cd) getClass("SpatialGrid") p4s <- CRS(proj4string(rsa2)) x2_SG <- SpatialGrid(rsa2_grd, proj4string = p4s) x2_SP <- SpatialPoints(cbind(x2$X, x2$Y)) v <- variogram(log1p(tsport_ace) ~ 1, x2, cutoff=100, width=9) te<- fit.va...

...@ data :'data.frame': 499105 obs. of 2 variables: .. ..$ PRECIP.pred: num [1:499105] 0.000536 0.000536 0.000536 0.000536 0.000536 ... .. ..$ PRECIP.var : num [1:499105] 1.90e-05 1.90e-05 1.90e-05 1.90e-05 1.90e-05 ... ..@ coords.nrs : num(0) ..@ grid :Formal class 'GridTopology' [package "sp"] with 3 slots .. .. ..@ cellcentre.offset: Named num [1:2] 345 4278 .. .. .. ..- attr(*, "names")= chr [1:2] "x" "y" .. .. ..@ cellsize : Named num [1:2] 0.1 0.1 .. .. .. ..- attr(*, "names")= chr [1:2] "x&quo...

Hi, I have a matrix of data which i can vizualize as an image - for example. I would like to save this image as a geotiff file or at a tiff file with a world file which holds the projection of my data (ultimately the data represent a map of some sort). I know i can save the data as an ESRI grid, but i am not interested in that. I wonder if anybody knows about any code which will help me do

...a) = ~Easting + Northing + SampleElevation class(PAHdata) summary(PAHdata) dimensions(PAHdata) hist(logTotalPAH) lpah.vgm = variogram(logTotalPAH ~1, PAHdata) lpah.vgm plot(lpah.vgm ) lpah.fit = fit.variogram(lpah.vgm, model = vgm( 6, "Sph", 80, 1.5)) lpah.fit plot(lpah.vgm, lpah.fit) x = GridTopology(c(534500,3531400, 20), c(2,2, 0.5), c(176,126, 31)) class(x) coordinatevalues(x) x = SpatialGrid(x) class(x) grd = as(x, "SpatialPixels") lpah.kriged = krige(logTotalPAH~1, PAHdata, x, model = lpah.fit) class(lpah.kriged) lpah.kriged = as(lpah.kriged, "SpatialPixelsDataFrame")...

...aller) do not overlap with the grids I have already created with spplot(). I would like shaded gridboxes, the boxes should be filled with colour. Can anyone help me with this problem? Thank you in advance! This is part of my code: minlon<-4.5 maxlon<-12 minlat<-46.5 maxlat<-52.5 gt=GridTopology(cellcentre.offset=c(5.75,47.5),cellsize=c(2.5,2),cells.dim=c(3,3)) grd=SpatialGrid(gt,proj4string = CRS(as.character(NA))) gridparameters(grd) gd<-as.data.frame(grd,"SpatialGrid") h<-season_djf h<-as.data.frame(h) djf.att=SpatialPointsDataFrame(gd,h) gridded(djf.att)=TRUE spplot...

...ts for my data. I've pasted my code below and I am not sure how to define the grd1 (cellsize and cells.dim). Here is my code: coords <- S at polygons[[1]]@Polygons[[1]]@coords coord <- as.points(coords) cases1.xy <- as.points(mergedis$jit.x.x, mergedis$pre.fsa.shp.Y2) grd1 <- GridTopology(cellcentre.offset=c(50281627.2, 580082), cellsize=c(0.2, 0.2), cells.dim=c(75,100)) k1000 <- spkernel2d(cases1.xy, coord, h0=1000, grd1) k5000 <- spkernel2d(cases1.xy, coord, h0=5000, grd1) if (.sp_lt_0.9()) { df <- AttributeList(list(k1000=k1000, k5000=k5000)) } else { df <- data.f...

Now I have the forest plot data with x, y locations, and I measured the DBH for every indivicuals. The data looks like that: x=runif(100,0,100) y=runif(100,0,100) dbh=runif(100,1,100) rdata=data.frame(x,y,dbh) > rdata[1:5,] x y dbh 1 99.5354145 1.412844 34.10112 2 0.8259361 87.737036 39.12914 3 6.5678613 65.699032 22.55990 4 67.2987881 72.053877 45.83978 5 2.2491372

Dear all, I'd like to change a point pattern to a grid of cells and use one of the variables as the output. e.g. The point pattern is of a window of (500*500) and several features such as pH, SoilType etc. I like to divide it into a grid with cell size 5*5, and use the mean of the point values falling inside the cell as the output. Is there any package in R working with this? Thanks in

R-geographers... I'm trying to solve a problem to implement a line-by-line tiled processing using RGDAL (read 1 line of an image, process the one line, write one line of the image to a binary file). Everything except for the final step I'm able to do using a combination of RGDAL and r-base commands. Below is the basic structure, the input file "elev" can be any image file

Hi, I am stuck in a looping problem. It might be an easy problem for experienced R users but I have been unable to do it. Any kind of help or advice will be great appreciated. I am creating a spatial grid and have a list of spatial points in a folder. I can read the spatial points separately and get a count of the points in the spatial grid in R. The output should be a table showing the counts of

Dear Happy R-users & experts, I am in need of advice, While working with spatial data (x & y coordinates of seed locations) I have come accross the problem that I need to convert my point data into a matrix or grid system. I then need to count how often a point falls into a certain position in the matrix or grid. I have searched all day online, asked collegeas but nothing works. Sadly