La construction d'un GeoJSON avec Python

Je veux générer dynamiquement un geoJSON avec un nombre variable de polygones. Exemple pour 2 polygones:

{
    "type": "FeatureCollection", 
    "features": [
      {"geometry": {
          "type": "GeometryCollection", 
          "geometries": [
              {
                  "type": "Polygon", 
                  "coordinates": 
                      [[11.0878902207, 45.1602390564],
                       [0.8251953125, 41.0986328125], 
                       [7.63671875, 48.96484375], 
                       [15.01953125, 48.1298828125]]
              }, 
              {
                  "type": "Polygon", 
                  "coordinates": 
                      [[11.0878902207, 45.1602390564], 
                       [14.931640625, 40.9228515625], 
                       [11.0878902207, 45.1602390564]]
              }
          ]
      }, 
      "type": "Feature", 
      "properties": {}}
    ]
}

J'ai une fonction qui me donne la liste des coordonnées de chaque polygone, donc je peux créer une liste de polygones, donc je suis en mesure de construire le geoJSON itération avec une boucle for.

Le problème est que je ne vois pas comment le faire facilement (je pense par exemple au retour de la liste comme une chaîne de caractères, mais la construction de la geoJSON comme une chaîne se présente comme une mauvaise idée).

J'ai été suggéré cette très pythonic idée:

geo_json = [ {"type": "Feature",,
              "geometry": {
                  "type": "Point",
                  "coordinates": [lon, lat] }}
              for lon, lat in zip(ListOfLong,ListOfLat) ]

Mais depuis que je suis l'ajout d'une variable de nombre de Polygones au lieu d'une liste de points, ces solutions ne semble pas adapté. Ou au moins, je ne sais pas comment l'adapter.

J'ai pu construire comme une chaîne de caractères, mais j'aimerais le faire de façon plus intelligente. Une idée?

OriginalL'auteur Roman Rdgz | 2013-06-04

6

Si vous pouvez obtenir les librairies installées, django a quelques bons outils pour faire face à la géométrie des objets, et ces objets ont une geojson attribut, vous donnant accès à la GeoJSON représentation de l'objet:

https://docs.djangoproject.com/en/2.0/ref/contrib/gis/install/
```
>>> from django.contrib.gis.geos import Polygon, Point, MultiPoint, GeometryCollection
>>>
>>> poly = Polygon( ((0, 0), (0, 1), (1, 1), (0, 0)) )
>>> gc = GeometryCollection(Point(0, 0), MultiPoint(Point(0, 0), Point(1, 1)), poly)
>>> gc.geojson
u'{ "type": "GeometryCollection", "geometries": [ { "type": "Point", "coordinates": [ 0.0, 0.0 ] }, { "type": "MultiPoint", "coordinates": [ [ 0.0, 0.0 ], [ 1.0, 1.0 ] ] }, { "type": "Polygon", "coordinates": [ [ [ 0.0, 0.0 ], [ 0.0, 1.0 ], [ 1.0, 1.0 ], [ 0.0, 0.0 ] ] ] } ] }'
```
GeometryCollection peut également accepter une liste d'objets de géométrie:
```
>>> polys = []
>>> for i in range(5):
...     poly = Polygon( ((0, 0), (0, 1), (1, 1), (0, 0)) )
...     polys.append(poly)
...
>>> gc = GeometryCollection(polys)
```
Il y a un problème avec cette approche: "je ne sais pas combien de polygones sont en train d'être créée pour chaque appel. GeometryCollection ne peut pas être créé vide comme ce que j'ai fait des recherches. Que dois-je faire?
juste un dump de votre polygones dans une liste et de le nourrir à GeometryCollection.
J'ai essayé de nourrir les GeometryCollection avec une liste de Polygones, mais il ne fonctionne pas. Il fonctionne si je crée le GeometryCollection avec le premier polygone, puis ajouter le reste d'entre eux, mais semble laid. Je suis de l'acceptation de votre réponse, mais s'il vous plaît fournir un exemple de la nourrir avec une liste si vous le pouvez
voir la mise à jour, il illustre la création d'un GeometryCollection() à l'aide de la liste, polys.
Le lien est cassé

OriginalL'auteur monkut

Il est le python-bibliothèque geojson (https://github.com/frewsxcv/python-geojson), ce qui semble rendre cette tâche beaucoup plus facile. Exemple à partir de la page de la bibliothèque:

>>> from geojson import Polygon

>>> Polygon([[(2.38, 57.322), (23.194, -20.28), (-120.43, 19.15), (2.38,   57.322)]])  
{"coordinates": [[[2.3..., 57.32...], [23.19..., -20.2...], [-120.4..., 19.1...]]], "type": "Polygon"}

OriginalL'auteur crisscross

Puisque vous avez déjà savoir comment construire une point, c'est assez similaire à la construction d'un polygon objet.
vous pouvez utiliser json.décharges pour convertir un objet python à la chaîne

Quelque chose comme:

geos = []
for longs,lats in LongLatList
    poly = {
        'type': 'Polygon',
        'coordinates': [[lon,lat] for lon,lat in zip(longs,lats) ]
    }
    geos.append(poly)

geometries = {
    'type': 'FeatureCollection',
    'features': geos,
}

geo_str = json.dumps(geometries) // import json

OriginalL'auteur gongzhitaao

Vous devez vous connecter pour publier un commentaire.