GIS

A continuación se muestran una serie de peticiones que se pueden hacer a una API STAC para realizar determinadas acciones comunes. Éste documente pretende ser una referencia rápida para realizar peticiones de manera inmediata. Las especificaciones completas de STAC API están disponibles en GitHub. Ver Colecciones de datos curl -X GET "https://my-stac-api.com/collections" \ -H "Authorization: Bearer <API_TOKEN> Crear nueva Colección de datos curl -X POST "https://my-stac-api.com/collections" \ -H "Authorization: Bearer <API_TOKEN>" \ -d '{"id": "<NOMBRE_COLECCIÓN>", \ "type": "Collection", \ "license": "proprietary", \ "description": "<DESCRIPCIÓN_COLECCIÓN>", \ "stac_version": "1.0.0", \ "stac_extensions": [], \ "extent": {"spatial": {"bbox": [[-180.0, -90.0, 180.0, 90.0]]}, "temporal": {"interval": [["1970-01-01T00:00:00Z", null]]}}}'

Suponiendo que tengamos un fichero, o simplemente un string, con GeoJSON, como por ejemplo: { "type": "MultiPoint", "coordinates": [ [-1.6402482999999999, 37.411843399999995], [-1.645345, 37.413398199999996], [-1.6453499999999999, 37.4129673], [-1.6453347, 37.41289], [-1.6455814, 37.4127428], [-1.6455834, 37.4136486], [-1.6402482999999999, 37.411843399999995] ] } y una base de datos PostGIS donde tengamos una tabla llamada field con un campo geometry, podríamos hacer una consulta para obtener todos los fields que intersecten con la geometría definida en el GeoJSON: SELECT field.id, field.name FROM field WHERE ST_Intersects(field_zone.geometry, ST_GeomFromGeoJSON(' { "type": "MultiPoint", "coordinates": [ [-1.6402482999999999, 37.411843399999995], [-1.645345, 37.413398199999996], [-1.6453499999999999, 37.4129673], [-1.6453347, 37.41289], [-1.6455814, 37.4127428], [-1.6455834, 37.4136486], [-1.6402482999999999, 37.411843399999995] ] }')); Los métodos de PostGIS utilizados son ST_Intersects, que comprueba intersección entre 2 objetos geometry o geography, y ST_GeomFromGeoJSON, que construye un objecto geometry a partir de su representación en GeoJSON. ...

A continuación se muestra una serie de ejemplos útiles de uso de la librería Python GeoPandas. Activar soporte para ficheros KML GeoPandas utiliza la librería Fiona para cargar datos desde ficheros. Por defecto, Fiona tiene el soporte para ficheros KML desactivado. Por lo tanto si usamos la función read_file() de GeoPandas con un fichero KML nos devuelve un error. Para activar el soporte de estos ficheros KML, podemos ejecutar la siguiente sentencia en nuestro programa Python: ...

A continuación se muestra una pequeña referencia sobre la instalación y algunos ejemplos de uso de la librería GDAL. Instalación Necesitaremos la librería base del sistema, escrita en C++, la cual es usada por GDAL-Python durante su compilación. La podemos instalar mediante apt en sistemas Debian: sudo apt-get install libgdal1-dev Una vez instalada la librería del sistema ya podemos instalar GDAL-Python. Es mucho mejor utilizar pygdal, que es una versión actualizada más amigable con el sistema de instalación de pip: ...