判斷是否在多邊形內(nèi)
1 生成測(cè)驗(yàn)數(shù)據(jù)
在arcgis中將excel中的點(diǎn)依次相連產(chǎn)生面狀范圍，選取面內(nèi)和面外的一些測(cè)驗(yàn)點(diǎn)，見(jiàn)ceshi.jpg。
2 在urls.py中添加路徑
url(r'^api/gis/inpoly/(?P<x>[+-]?\d.\d)/(?P<y>[+-]?\d.\d)$', InPoly)
3 在world.views中增加函數(shù)
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 測(cè)驗(yàn)
輸入測(cè)試的url，例如：http://127.0.0.1:1234/api/gis/inpoly/-113.876/+22.538
返回：{"InOrNot": false}
5 api接口說(shuō)明
[get] api/gis/inpoly/x/y
參數(shù)說(shuō)明：x為精度，y為緯度。返回說(shuō)明：返回為json類(lèi)型，當(dāng)InOrNot的值為true就在范圍內(nèi)，為false則在范圍外。