Displaying 9 results from an estimated 9 matches for "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,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] 26...
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 don...
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, c...
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