codepad
[
create a new paste
]
login
|
about
Language:
C
C++
D
Haskell
Lua
OCaml
PHP
Perl
Plain Text
Python
Ruby
Scheme
Tcl
vlPolygons = QgsVectorLayer("C:\\Users\Windows\\Documents\\python_code\\test_contains\\us.shp", "poligono", "ogr") featPoly = QgsFeature() vlPoints = QgsVectorLayer("C:\\Users\\Windows\\Documents\python_code\\test_contains\\quote.shp", "punto", "ogr") featPoints = QgsFeature() if not vlPolygons.isValid(): print "Layer failed to load!" if not vlPoints.isValid(): print "Layer failed to load!" iterPoly = vlPolygons.getFeatures() iterPoints = vlPoints.getFeatures() vlPoints.select( [] ) list_poly = [] for i in iterPoly: list_poly.append(i) results_true = [] results_false = [] for feature in list_poly: iterPoints = vlPoints.getFeatures() geom = feature.geometry() bbox = feature.geometry().boundingBox() pointBBox = vlPoints.select( feature.geometry().boundingBox(), True) #print "geom is %d" % feature.id() for points in iterPoints: #print "geom is %d" % feature.id() if geom.contains( points.geometry()) == True: results_true.append('US%d contains Quota %d' % (feature.id(),points.id())) print results_true
Private
[
?
]
Run code
Submit