判断是否在多边形内

判断是否在多边形内
1 生成测验数据
在arcgis中将excel中的点依次相连产生面状范围，选取面内和面外的一些测验点，见ceshi.jpg。
2 在urls.py中添加路径
url(r'^api/gis/inpoly/(?P[+-]?\d.\d)/(?P[+-]?\d.\d)$', InPoly)
3 在world.views中增加函数
def InPoly(request,x,y): path=os.path.abspath(os.path.join(os.path.dirname(file)))+'/data/polygon.xlsx' bk = xlrd.open_workbook(path) table = bk.sheets()[0] nrows = table.nrows ncols = table.ncols plist=[] for i in range(1,nrows): plist.append(GEOSGeometry('POINT(%s %s)' %(table.cell(i,1).value, table.cell(i,2).value))) p = Polygon(plist) print p # pnt = Point(113.885, 22.517) lon=float(x) lat=float(y) pnt = Point(lon, lat) InOrNot=p.contains(pnt) # InOrNot=p.filter(poly__contains=pnt) Bldict={"InOrNot":InOrNot} print json.dumps(Bldict) return HttpResponse(json.dumps(Bldict), content_type="application/json")
4 测验
输入测试的url，例如：http://127.0.0.1:1234/api/gis/inpoly/-113.876/+22.538
返回：{"InOrNot": false}
5 api接口说明
[get] api/gis/inpoly/x/y
参数说明：x为精度，y为纬度。返回说明：返回为json类型，当InOrNot的值为true就在范围内，为false则在范围外。