search for: cellcentr

Displaying 9 results from an estimated 9 matches for "cellcentr".

Did you mean: cellcentre
2011 Jan 05
1
Prediction error for Ordinary Kriging
...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.variogram(v,...
2010 Oct 19
1
Problems with a specific calculate.
.... ..$ 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" "y" .. .. ..@ cells.dim : Named int [1:2] 2...
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
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
2007 Jul 06
0
import DTM with readRAST6()
...coords.x1 1611600 1729700 coords.x2 5058500 5158000 Is projected: TRUE proj4string : [+proj=tmerc +lat_0=0.0000000000 +lon_0=9.0000000000 +k_0=0.9996000000 +x_0=1500000.0000000000 +y_0=0.0000000000 +a=6378388 +rf=297 +no_defs +towgs84=-225.000,-65.000,9.000] Number of points: 2 Grid attributes: cellcentre.offset cellsize cells.dim 1 1611650 100 1181 2 5058550 100 995 Data attributes: DTM_trentino_100m Min. :1.192e-319 1st Qu.:1.065e-314 Median :1.066e-314 Mean :1.285e-312 3rd Qu.:2.722e-312 Max. :4.774e-312 As you can see, the import has been fully do...
2010 Nov 02
0
spatial plots maps-ssplot
...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(djf.att,...
2009 Aug 18
0
kernel density estimation for univariate data using splancs
...a. 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.frame(k1000...
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