[R-sig-Geo] spacetime : the challenge of image time series containing holes

2011-04-08 Thread servet cizmeli
Greetings to all, I have monthly mean satellite image time series that come in a 4-column text file each. I have dozens of such files. Every line correspond to a pixel that is projected onto a regular ISIN grid. The first column is the ISIN grid number (same coordinates for a given grid number

Re: [R-sig-Geo] [R-sig-geo] Sum lines when they intersect polygons

2011-04-08 Thread Michael Sumner
If you really just need a grid then pixellate.psp in spatstat can do this, pass in the "length" with the weights argument. On Sat, Apr 9, 2011 at 12:32 AM, Barry Rowlingson wrote: > Slow and dirty method: > > lineLengthIntersects <- function(gridpolys,linegeom){ >  ngrids = length(gridpolys) >  l

Re: [R-sig-Geo] spplot, polygons, and line density fill

2011-04-08 Thread Matthew Landis
In case my request wasn't really clear, let me rephrase: I'd like to overlay some polygons onto a raster, by plotting with a line density fill instead of a solid fill, as shown here http://r-spatial.sourceforge.net/gallery/#fig15.R. Ideally I would be able to modify the angle of the lines. I

[R-sig-Geo] spplot, polygons, and line density fill

2011-04-08 Thread Matthew Landis
Hello again R-sig-geo-ers, Does anyone know how to convince spplot to create a figure like http://r-spatial.sourceforge.net/gallery/#fig15.R I would like to overlay some polygons onto a raster, and I need to do it with spplot for the lattice functionality. I've looked at the help and the fu

Re: [R-sig-Geo] [R-sig-geo] Sum lines when they intersect polygons

2011-04-08 Thread Barry Rowlingson
Slow and dirty method: lineLengthIntersects <- function(gridpolys,linegeom){ ngrids = length(gridpolys) ll = rep(0,ngrids) for(i in 1:ngrids){ geomint = gIntersection(linegeom,gridpolys[i]) if(!is.null(geomint)){ ll[i] = SpatialLinesLengths(geomint,TRUE) } } return(ll)

Re: [R-sig-Geo] [R-sig-geo] Sum lines when they intersect polygons

2011-04-08 Thread Barry Rowlingson
On Fri, Apr 8, 2011 at 9:10 AM, jonathan daeden wrote: > Dear list, > > I am looking for a way to measure lines crossing polygones or a grid. > I have a line shapefile representing a coastline and a regular grid over it. > I would like to measure the coastline in each cells of the grid and keep it

Re: [R-sig-Geo] Support in krige/gstat

2011-04-08 Thread Edzer Pebesma
On 04/08/2011 12:40 PM, piero campa wrote: > Dear list, > I have datasets of different variables with different spatio-temporal > supports, and I'd like to [co]krige them together. > > What I'd like to ask you is: > - I know that with block kriging one could estimate values over a different > su

Re: [R-sig-Geo] Error message in rgeos

2011-04-08 Thread Leonardo Monasterio
Dear Colin, Thank you very much. I've uploaded shapefile to: http://dl.dropbox.com/u/116353/OD87.zip Yours, Leonardo. On Thu, Apr 7, 2011 at 6:39 PM, Colin Rundel wrote: > >> I am getting the following error message when I try: >> >>> lineborder<-gBoundary(map1) >> Erro em TopologyFunc(spgeom,

[R-sig-Geo] Support in krige/gstat

2011-04-08 Thread piero campa
Dear list, I have datasets of different variables with different spatio-temporal supports, and I'd like to [co]krige them together. What I'd like to ask you is: - I know that with block kriging one could estimate values over a different support area wrt the input variables, but is one allowed to h

Re: [R-sig-Geo] [R-sig-geo] Sum lines when they intersect polygons

2011-04-08 Thread Jacob van Etten
Do you want to measure the *length* of the coastline in each grid cell? Or do you want to measure the *presence* of the coastline in each cell? Jacob. --- On Fri, 8/4/11, jonathan daeden wrote: > From: jonathan daeden > Subject: [R-sig-Geo] [R-sig-geo] Sum lines when they intersect polygons >

[R-sig-Geo] [R-sig-geo] Sum lines when they intersect polygons

2011-04-08 Thread jonathan daeden
Dear list, I am looking for a way to measure lines crossing polygones or a grid. I have a line shapefile representing a coastline and a regular grid over it. I would like to measure the coastline in each cells of the grid and keep it in a data.frame. Could not make it with Qgis so i would like to