Skip to Content

Nous adressons toutes nos pensées à la famille de notre ami Jérôme !

http://www.forumsig.org/showthread.php/43488-Disparition-de-Phoenix

collection géospatiale à partir d'un shapefile point - modification du script de Paolo Corti

  1. from pymongo import Connection, GEO2D
  2. from osgeo import ogr
  3.  
  4. # lecture du shapefile
  5. driver = ogr.GetDriverByName('ESRI Shapefile')
  6. shape_path="points.shp"
  7. ds = driver.Open(shape_path, 0)
  8. lyr = ds.GetLayer()
  9. feat = lyr.GetNextFeature()
  10.  
  11. # connexion et création de la collection géospatiale
  12. db = Connection().geo_example
  13. db.places.create_index([("loc", GEO2D)])
  14.  
  15. # insertion du premier élément dans MongoDB
  16. geom = feat.GetGeometryRef()
  17. db.places.insert({"loc": [ geom.GetX(), geom.GetY()]})
  18.  
  19. # etc. avec la boucle du script de Paolo Corti