> 2008/10/13 Mathew Yeates <[EMAIL PROTECTED]>
>
> > Is there a routine in scipy for telling whether a point is inside a
> > convex 4 sided polygon?
Mathew,
You could use OGR (www.gdal.org)
Example
-------------
import osgeo.ogr as ogr
vert = [(0,0),(0,1),(1,1),(1,0)]
listvert = [" %s %s" % (x,y) for (x,y) in vert]
listvert.append(listvert[0])
geo = ogr.CreateGeometryFromWkt("POLYGON ((%s))" % ','.join(listvert))
querypoint = (0.5, 0.5)
qpt = ogr.CreateGeometryFromWkt("POINT(%s %s)" % querypoint)
assert geo.Contains(qpt)
querypoint = (0.5, 1.5)
qpt = ogr.CreateGeometryFromWkt("POINT(%s %s)" % querypoint)
assert not geo.Contains(qpt)
_______________________________________________
Numpy-discussion mailing list
Numpy-discussion@scipy.org
http://projects.scipy.org/mailman/listinfo/numpy-discussion
