Análisis Vectorial

gdal vector convert

gdal vector convert se puede utilizar para convertir datos entre diferentes formatos.

Ejemplos

Ejemplo 1: Convertir el archivo poly.shp a un GeoPackage

gdal vector convert poly.shp output.gpkg

Ejemplo 2: Añadir una nueva capa del archivo line.shp a un GeoPackage existente y renombrarla "lines"

gdal vector convert --update --output-layer=lines line.shp output.gpkg

Ejemplo 3: Añadir entidades del archivo poly2.shp a una capa poly existente de un GeoPackage, con barra de progreso

gdal vector convert --append --output-layer=poly --progress poly2.shp output.gpkg

gdal vector concat

gdal vector concat concatena varios conjuntos de datos de origen.

Ejemplos

Ejemplo 1: Creación de un GeoPackage que apila todos los shapefiles de entrada en capas separadas.

gdal vector concat --stack *.shp out.gpkg

Ejemplo 2: Adición de un campo para indicar la capa de origen y reproyección a un único CRS

Concatenar el contenido de france.shp y germany.shp en merged.shp, reproyectándolos a ETRS89, y añadir un campo "country" a cada entidad cuyo valor sea "france" o "germany", según su origen:

gdal vector concat --single --source-layer-field-name=country --dst-crs=EPSG:4258 france.shp germany.shp merged.shp

gdal vector filter

El filtro vectorial gdal se puede utilizar para filtrar un conjunto de datos vectoriales a partir de su extensión espacial o una cláusula WHERE de SQL.

Ejemplos

Ejemplo 1: Seleccionar entidades de un archivo GeoPackage que intersecan el cuadro delimitador desde la longitud 2, latitud 49, hasta la longitud 3, latitud 50 en WGS 84

gdal vector filter --bbox=2,49,3,50 in.gpkg out.gpkg --overwrite

gdal vector geom

Ejercicios

1. Utilizando el conjunto de capas del archivo Geo_CR.zip (disponible en el repositorio de datos del curso), extraiga todos los poblados de la provincia de Cartago. Exporte el resultado a un archivo GeoJSON.

2. Utilizando el mismo conjunto de capas, ejecute una unión de los distritos de cada cantón. Exporte el resultado a un archivo GeoJSON y comprímalo en ZIP.

3. Genere un buffer (círculo) de 5km (aproximadamente) alrededor de los poblados que son del tipo CIUDAD. Exporte el resulta a un archivo tipo GeoJSON.