Re: [R-sig-Geo] point.in.polygon() on massive datasets

2007-12-14 Thread Galkowski, Jan
Hi Markus, I have had success with a somewhat smaller point-in-polygon problem using a relational database (PostgreSQL) with suitable indexes to determine containment. The database serves as a large lookup table for this purpose, using RODBC. I use R to build a series of SQL INSERT statements to

Re: [R-sig-Geo] overlay(SpatialPolygonsDataFrame, SpatialPointsDataFrame)

2007-12-14 Thread Roger Bivand
On Fri, 14 Dec 2007, Ingo Holz wrote: > Hi, > > I have a SpatialPolygonsDataFrame (SPolyDF) and a > SpatialPointsDataFrame(SPoiDF). The SPolyDF is a (very) large set of > polygons (some 30.000) and the SPoiDF is actually a grid of some 60.000 > points. > > I would like to use overlay(SPolyDF, SPoi