search for: gridtopology

Displaying 14 results from an estimated 14 matches for "gridtopology".

2008 Feb 04
1
how to get points from SpatialPolygonsDataFrame
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
2007 Mar 03
1
R software to place points on Yahoo maps
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
2009 Jul 01
0
The step before interfacing to GRASS
...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. _____________________________________...
2011 Jan 05
1
Prediction error for Ordinary Kriging
...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...
2010 Oct 19
1
Problems with a specific calculate.
...@ 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...
2007 Sep 05
1
geotiff or tiff files with world files
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
2008 Aug 07
0
3d kriging et al
...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")...
2010 Nov 02
0
spatial plots maps-ssplot
...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...
2009 Aug 18
0
kernel density estimation for univariate data using splancs
...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...
2007 Dec 12
1
How to plot the grid figure using R?
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
2005 Nov 17
8
Point pattern to grid
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
2008 Nov 05
1
Using RGDAL to "copy" header information...
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
2012 Dec 17
2
looping through spatial points and getting counts of spatial points in spatial grid in R
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
2007 Apr 18
3
Matrix or grid conversion of spatial data
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