¿Cómo cargar datos de puntos en formato JSON desde una URL en Qgis utilizando la consola de python?

Question

A través de una URL: http://www.odaa.dk/api/action/datastore_search?resource_id=c3097987-c394-4092-ad1d-ad86a81dbf37 Tengo lo que parece ser un archivo en formato JSON que contiene puntos que quiero en QGIS directamente desde la URL, preferiblemente con la consola de Python. Dentro del archivo JSON columnas para la longitud y latitud están presentes.

¿Alguien sabe cómo se puede hacer esto?

Answer 1

Estos parecen (partes de) carreteras, siempre hay 2 coordenadas de puntos: POINT_1_LAT , POINT_1_LNG , POINT_2_LAT , POINT_2_LNG

Escribí un pequeño script (usando el ejemplo del Libro de cocina de QGIS ) que creará una capa y cargará los datos:

import urllib, json
from PyQt4.QtCore import *
import qgis
from qgis.core import *

# create memory layer
vl = QgsVectorLayer("LineString", "point_data", "memory")
pr = vl.dataProvider()

# add fields
pr.addAttributes([QgsField("_id", QVariant.Int),
                    QgsField("ROAD_TYPE",  QVariant.String),
                    QgsField("ORGANISATION", QVariant.String)])
vl.updateFields() # tell the vector layer to fetch changes from the provider

# download data and parse it
url = 'http://www.odaa.dk/api/action/datastore_search?resource_id=c3097987-c394-4092-ad1d-ad86a81dbf37'
response = urllib.urlopen(url)
data = json.loads(response.read())

# iterate points and add to layer
for point in data['result']['records']:
    linestring_feature = QgsFeature()
    point1 = QgsPoint(float(point['POINT_1_LNG']), float(point['POINT_1_LAT']))
    point2 = QgsPoint(float(point['POINT_2_LNG']), float(point['POINT_2_LAT']))
    linestring = QgsGeometry.fromPolyline([point1, point2]);
    linestring_feature.setGeometry(linestring)
    linestring_feature.setAttributes([point["_id"], point["ROAD_TYPE"], point["ORGANISATION"]])
    pr.addFeatures([linestring_feature])

# update layer's extent when new features have been added
vl.updateExtents()

# sdd layer to layers panel
QgsMapLayerRegistry.instance().addMapLayer(vl)

Puedes utilizarlo a través de la caja de herramientas de Processing - Scripts, o simplemente pegarlo en la consola de Python.

Seleccione WGS 84 (EPSG:4326) como sistema de referencia cuando aparezca la ventana.