Guest User

Untitled

a guest
Oct 17th, 2018
107
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.71 KB | None | 0 0
  1. canvas = qgis.utils.iface.mapCanvas()
  2. csql='ST_Intersects("p1".geom, "p2".geom)'
  3. uri = QgsDataSourceURI()
  4. uri.setConnection("localhost", "5432", "ZBE", "postgres", "postgres")
  5. uri.setDataSource("public","Sample", "geom",csql,"dn")
  6. vlayer = QgsVectorLayer(uri.uri(), "test", "postgres")
  7. if not vlayer.isValid():
  8. print "not valid"
  9. QgsMapLayerRegistry.instance().addMapLayer(vlayer)
  10. canvas.refresh()
  11.  
  12. csql='SELECT field1,field2,etc
  13. FROM "Trainingset" t, "Polygon" p
  14. WHERE ST_Intersects(t.geom, p.geom)'
  15.  
  16. CREATE VIEW boo AS SELECT * FROM p1, p2 WHERE ST_Intersects(p1.geom, p2.geom)
  17.  
  18. csql = "p1.id IN (SELECT p1.id FROM p1, p2 WHERE ST_Intersects(p1.geom, p2.geom)"
Add Comment
Please, Sign In to add comment