The second half of this blog post has a code example that will help you code the test for the points being in the polygons
On Wed, Nov 20, 2013 at 2:54 PM, Abhishek <cheerful...@gmail.com> wrote: > Hello, > > First, I'm very new to doing gis programming in python (and GIS in > general) so apologies if this has been answered before. > > I was wondering if there was a straightforward way for me to do the > following: > 1. Read in a shapefile (all counties in the UK) > 2. Read in a csv with a list of lat / lon (this I know how to do) > 3. Spit out a csv with the same list of points and the polygon id from > the shapefile, indicating that the point lies within the polygon (i.e. > in the county) > > Ive got this to work by manually constructing polygons from a csv > export of the shapefile, but I'm afraid that Im not going to deal very > well with holes in counties, or edge cases. A standard out-of-the-box > solution would be great. > > Thanks! > > Abhishek > ------------ > http://abhishek.mit.edu > _______________________________________________ > Community mailing list > Community@lists.gispython.org > http://lists.gispython.org/mailman/listinfo/community >
_______________________________________________ Community mailing list Community@lists.gispython.org http://lists.gispython.org/mailman/listinfo/community
