Realizar Consultas Espaciales (QGIS3)

Las consultas espaciales son centrales a varios tipos de análisis SIG. Las consultas espaciales le permiten seleccionar objetos espaciales en una capa por su relación espacial (intersecta, contiene, toca, etc.) con objetos espaciales de otra capa. en QGIS, esta funcionalidad está disponible mediante las herramientas de Procesamiento Seleccionar por Ubicación y Extraer por Ubicación.

Vista general de la tarea

Estaremos trabajando con 2 capas de datos para la ciudad de Melbourne, Australia. Dadas las capas de datos para las tabernas y bares en la ciudad y las ubicaciones de todas las estaciones de metro, vamos a descubrir todos los bares y tabernas dentro de los 500 metros de una estación de metro.

Otras habilidades que aprenderá

  • Escoger una proyección apropiada y reproyectar datos vector.

  • Crear zonas de influencia.

  • Trabajar con el formato de datos geopackage (.gpkg).

Obtener los datos

La Plataforma de Datos Abiertos de la Ciudad de Melbourne provee varios conjuntos de datos listos para el SIG para la ciudad.

Download the Metro Train Stations with Accessibility Information dataset by Metro Trains Melbourne. Export the data in the Original format.

../../_images/data117.png

Download the Bars and pubs, with patron capacity dataset by City of Melbourne’s Census of Land Use and Employment (CLUE). Export the data as a CSV.

../../_images/data27.png

Para su comodidad, puede descargar directamente una copia de los conjuntos de datos del enlace abajo:

metro_stations_accessbility.zip

Bars_and_pubs__with_patron_capacity

Fuente de Datos: [CITYOFMELBOURNE]

Procedimiento

  1. Localice el archivo metro_stations_accessbility.zip en el Explorador QGIS y expándalo. Seleccione el archivo metro_stations_accessbility.shp y arrástrelo a la pantalla. Una nueva capa metro_stations_accessbility será cargada en el panel Layers.

../../_images/1137.png
  1. La capa de datos para bares y tabernas está en el formato CSV. Para cargarla en QGIS, vaya a Capa ‣ Agregar Capa ‣ Agregar Capa de Texto Deliminatdo…. ( Ver Importar Hojas de Cálculo o archivos CSV (QGIS3) para más detalles en la importación de archivos CSV)

../../_images/2108.png
  1. In the Data Source Manager | Delimited Text dialog, browse and select the downloaded Bars_and_pubs__with_patron_capacity.csv file as File name. The X field and Y field columns should be auto selected to x coordinate and y coordinate respectively. Click Add.

../../_images/356.png

Nota

QGIS saves the last-selected options for this dialog, so if you previously configured this dialog with other options, they may still be selected. Specifically, make sure the Record and Fields Options section has Number of header lines to discard option set to 0.

  1. Verá una nueva capa Bars_and_pubs__with_patron_capacity añadida al panel Capas`. Ambas capas de entrada están en el sistema de referencia de coordenadas geográficas (SRC) EPSG:4326 WGS84. Para realizar análisis espaciales, se recomienda utilizar un Sistema de Referencia de Coordenadas Proyectado (SRC). Así que ahora vamos a reproyectar ambas capas a un SRC regional apropiado que minimice las distorsiones y nos permita trabajar en unidades de distancia como metros en lugar de grados. Vaya a :menuselección:`Procesos --> Caja de herramientas`.

../../_images/428.png
  1. Busque y localice la herramienta Vector general ‣ Reproyectar capa. Haga doble-clic para iniciarla.

../../_images/526.png
  1. Seleccione Bars_and_pubs__with_patron_capacity como la Capa de entrada. Clic el botón Seleccionar SRC``junto al :guilabel:`SRC Destino.

../../_images/626.png
  1. Cuando seleccione un sistema de referencia de coordenadas para su análisis, el primer lugar para mirar es un SRC regional para el área de interés. Para Australia, la Cuadrícula de Mapa de Australia (MGA) 2020 es un sistema de cuadrícula basado en UTM que es usado para mapeo local y regional. Melbourne cae en la Zona 55 UTM, por lo que podemos seleccionar el SRC GDA 2020 / MGA zone 55 EPSG:7855.

../../_images/725.png

Nota

Si no está seguro de un SRC local para la región en la que está trabajando, seleccione un SRC para la zona UTM basado en el datum WGS84 que es una elección segura. Puede descubrir el numero de zona UTM de su región usando las Zonas de Cuadrícula UTM del Mundo.

  1. A continuación, haga clic en el botón junto a Reprojected y seleccione Guardar en GeoPackage. Geopackage es el formato recomendado para datos espaciales abiertos y el formato de intercambio de datos predeterminado para QGIS3. Un solo archivo GeoPackage .gpkg puede contener múltiples capas vectoriales y ráster.

../../_images/825.png
  1. Nombre el geopackage como spatialquery y clic Guardar.

../../_images/925.png
  1. Cuando se le pida un nombre de capa, ingrese bars_and_pubs y clic Aceptar. Clic Ejecutar para reproyectar la capa.

../../_images/1032.png
  1. La ventana se cambiará a la pestaña Bitácora y verá que el algoritmo se ejecuta y crea la nueva capa de salida bars_and_pubs.

../../_images/1138.png
  1. Ahora vamos a reproyectar la capa metro_stations_accessbility. Cámbiese de vuelta a la pestaña Parámetros en la ventana Reproyectar capa. Seleccione metro_stations_accessbility como la Capa de entrada. Mantenga el mismo SRC Destino. A continuación, clic en el botón junto a Reprojected y seleccione Guardar a Geopackage. Seleccione el mismo archivo de salida spatialquery (Recuerde que un solo archivo geopackage puede contener múltiples capas, por lo que guardaremos la nueva capa en el mismo archivo geopackage). Ingrese metro_stations como el Nombre de capa. Clic en Ejecutar.

../../_images/1233.png
  1. De vuelta en la ventana principal QGIS, veremos 2 nuevas capas agregadas en el panel Capas: bars_and_pubs y metro_stations. Puede desactivar la visibilidad de las capas originales. Ahora, estamos listos para hacer la consulta espacial. Como estamos interesados en seleccionar bares y tabernas a 500 m de las estaciones de metro, el primer paso es crear un área de influencia alrededor de las estaciones de metro que represente nuestra área de búsqueda. Busque y localice la herramienta Geometría vectorial ‣ Buffer en la Caja de herramientas de Procesos y haga doble-clic para iniciarla.

../../_images/1331.png
  1. En el cuadro de diálogo Buffer, seleccione metro_stations como Capa de entrada. Establezca 500 metros como Distancia. Guarde el resultado en el mismo geopackage spatialquery e introduzca metro_stations_buffers como Nombre de capa`. Haga clic en Ejecutar.

../../_images/1429.png
  1. Verá una nueva capa metro_stations_buffers cargada en el panel Capas. Ahora podemos averiguar qué puntos de la capa bars_and_pubs caen en los polígonos de la capa metro_station_buffers. Localice la herramienta Selección vectorial ‣ Extraer por ubicación desde la Caja de herramientas de Procesos y haga doble-clic para iniciarla.

../../_images/1524.png

Nota

Extraer por Ubicación creará una nueva capa con los objetos espaciales que coincidan con la consulta espacial. Si sólo quiere seleccionar los objetos espaciales, use la herramienta Seleccionar por Localización.

  1. En el diálogo Extraer por ubicación, seleccione bars_and_pubs como Extraer objetos de. Marque Intersectan como el Predicado geométrico. Defina metro_station_buffers como Comparando los objetos de. Guarde la salida al geopackage spatialquery como la capa selected. Clic en Ejecutar.

../../_images/1623.png
  1. Una vez que termine el procesamiento, verá la capa selected agregada al panel Capas. Note que esta capa sólo contiene puntos de bars_and_pubs que caen dentro de los polígonos del área de influencia.

../../_images/1724.png
  1. Nuestro análisis está completo. Puede que note que los polígonos del área de influencia lucen una forma oval. Esto es debido a que nuestro SRC de Proyecto aún está definido como EPSG:4326 WGS84. Para visualizar mejor los resultados, puede ir a Proyecto ‣ Propiedades ‣ SRC y seleccionar GDA 2020 / MGA zone 55 EPSG:7855 que usamos para el análisis. Una vez que se defina este SRC, el área de influencia aparecerá con la forma correcta.

../../_images/1821.png

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)