En esta lección, serás guiado a través de un análisis SIG completo en QGIS.
Nota
Lección desarrollada por Linfiniti y S Motala (Universidad Tecnológica de Península del Cabo)
Eres el encargado de encontrar áreas en Península del Cabo que sean adecuadas para el hábitat de especies de plantas fynbos raras. La extensión de tu área de investigación en Península del Cabo es: Sur de Melkbosstrand, Oeste de Strand. Botánicos te han proporcionado las preferencias exhibidas por las especies en cuestion:
Crecen en pendientes orientadas al este.
Crecen en pendientes con gradiente entre 15 y 60%.
Crecen en áreas que tienen un total de precipitación anual de > 1200 mm.
Solo se encuentran a al menos 250 m de distancia de los asentamientos humanos.
El área de vegetación en la que ocurre debería ser de al menos 6000m2.
Como voluntario de Naturaleza del Cabo, has acordado buscar la planta en el terreno adecuado más cercano a tu casa. Utiliza tus habilidades SIG para determinar dónde deberías ir a ver.
In order to solve this problem, you will have to download data from https://www.dropbox.com/s/q5evvkizuunrcs0/more_analysis.zip?dl=0 and place it under exercise_data/more_analysis folder. You will use it to find the candidate area that is closest to your house. If you don’t live in Cape Town (where this problem is based) you can choose any house in the Cape Town region. The solution will involve:
análisis del MDE para encontrar las pendientes orientadas al este y las pendientes con el gradiente correcto;
análisis del ráster de precipitación para encontrar las áreas con cantidad de precipitación correcta;
análisis de las zonas de capa vectorial para encontrar áreas que tengan la distancia a enclaves humanos y sean del tamaño correctos.
Haz clic en el botón “estado SRC” en la esquina inferior derecha de la pantalla. En la pestaña SRC de la pestaña que aparece, verás la casilla Sistema de referencia de coordenadas del mundo.
En esta casilla, navega hasta Sistemas de coordenadas proyectadas ‣ Universal Transverse Mercator (UTM).
Guarda el mapa clicando en el botón de la barra de herramientas Guardar proyecto como, o utiliza el elemento del menú Proyecto ‣ Guardar como....
Para procesar los datos, necesitarás cargar las capas necesarias (nombres de calles, zonas, precipitación, MDE) en la vista del mapa.
Haz clic en Abrir.
The dialog closes and shows the original dialog, with the file path specified in the text field next to the Browse button. This allows you to ensure that the correct file is selected. It is also possible to enter the file path in this field manually, should you wish to do so.
Navega hasta el archivo apropiado, selecciónalo y haz clic en Abrir.
Guarda el mapa.
Para ver correctamente qué está pasando, se necesita cambiar la simbología para las capas.
Selecciona Propiedades del menú que aparece.
Haga clic sobre el botón etiquetado Cambiar, con un cuadro mostrando el color actual de la capa Streets.
Selecciona un color nuevo en el cuadro de diálogo que aparece.
Haz clic en Aceptar.
Vuelve a hacer clic en Aceptar en el cuadro de diálogo Propiedades de la capa. Esto cambiará el color de la capa Streets.
Sigue un proceso similar para la capa Zoning y elige un color adecuado para ella.
La simbología de capas ráster es algo diferente.
Abre el cuadro de diálogo Propiedades para el ráster Precipitación.
Asegúrate de que el botón Media +/- desviación estándar está seleccionado.
En el título Mejora de contraste, cambia el valor del menú desplegable a Estirar a MinMax.
Haz clic en Aceptar. El ráster “Precipitación”, si está visible, debería cambiar de colores, permitiendote ver diferentes valores de brillo para cada píxel.
Ahora que todos los datos están cargados y son adecuadamente visibles, el análisis puede comenzar. Es mejor si las operaciones de recorte se hacen primero. Así no se malgasta potencia procesando áreas que no serían utilizadas de todas formas.
En el menú que aparece, selecciona el elemento del menú Filtrar.... El cuadro de diálogo Constructor de Consultas aparecerá.
Ahora construirás una consulta para seleccionar solo la siguiente lista de distritos:
En la lista Campos, haz doble clic en el campo NAME_2. Aparece en el campo de texto inferior Expresión de filtrado específica del proveedor.
Haz clic en el botón Todos bajo la lista Values (actualmente vacío). Después de una corta espera, se rellenará la lista Values con los valores del campo seleccionado (NAME_2).
Haz doble clic en el valor Bellville en la lista Values. Como antes, será añadido a la consulta SQL.
In order to select more than one district, you’ll need to use the OR boolean operator.
Haz clic en el botón O y será añadido a la consulta SQL.
Utilizando un proceso similar al anterior, añade lo siguiente a la consulta existente SQL:
"NAME_2" = 'Cape'
Add another OR operator, then work your way through the list of districts above in a similar fashion.
La consulta final debería ser
"NAME_2" = 'Bellville' OR "NAME_2" = 'Cape' OR "NAME_2" = 'Goodwood' OR
"NAME_2" = 'Kuils River' OR "NAME_2" = 'Mitchells Plain' OR "NAME_2" =
'Simons Town' OR "NAME_2" = 'Wynberg'
Haz clic en Aceptar. Los distritos mostrados en tu mapa están limitados a los de la lista anterior.
Ahora que tienes un área de interés, puedes recortar los ráster a esa área.
Asegúrate de que las únicas capas visibles son las de MDE, Rainfall y Districts.
Los Districts deben estar en primer lugar para ser visibles.
Abre el cuadro de diálogo para recortar seleccionando el elemento del menú Ráster ‣ Extracción ‣ Clipper.
En la lista del menú desplegable Archivo de entrada (ráster), selecciona la capa MDE.
Especifica una localización de salida en el campo de texto Archivo de salida clicando en el botón Seleccionar....
Introduce un nombre de archivo.
Guarda el archivo. Deja la casilla de verificación Valor de sin datos vacía.
Utiliza el modo de corte Extensión asegurándote que el botón radial correcto está seleccionado.
Haz clic y arrastra un área en el mapa, de modo que el área que incluye los distritos esté seleccionada.
Haz clic en Aceptar.
Cuando la operación de corte se complete, NO CIERRES el cuadro de diálogo Clipper. (Cerrarlo te haría perder el área a cortar que ya habías definido)
Selecciona el ráster Precipitación en el menú desplegable Archivo de entrada (ráster) y elige un nombre de archivo de salida diferente.
Después de que la segunda operación de corte haya terminado, puedes cerrar el cuadro de diálogo Clipper.
Guarda el mapa.
Haz clic derecho en esas capas y selecciona Eliminar.
Esto no borrará los datos de tu dispositivo de almacenamiento, solamente lo quitará de tu mapa.
Desactiva las etiquetas de la capa Streets.
Haz clic en el botón Etiquetado.
Quita la marca de la casilla Etiquetar esta capa con.
Haz clic en Aceptar.
Muestra todas las Streets de nuevo:
Selecciona Filtrar.
En el cuadro de diálogo que aparece, haz clic en el botón Limpiar, luego haz clic en Aceptar.
Espera mientras se cargan los datos. Todas las calles serán visibles ahora.
Guarda el mapa.
In order to create the hillshade, you will need to use a plugin that was written for this purpose.
This plugin is included by default in QGIS 1.8. However, it may not be immediately visible. To check if it is accessible on your system:
Haz clic en Aceptar.
You will now have access to this plugin via the Raster ‣ Terrain analysis menu item.
Remember that plugins may sometimes depend on certain Python modules being installed on your system. Should a plugin refuse to work while complaining of missing dependencies, please ask your tutor or lecturer for assistance.
Especifica una localización apropiada para la capa de salida y nómbrala hillshade.
Haz clic en Aceptar.
Espera a que termine de procesar.
La nueva capa hillshade a aparecido en tu Lista de capas.
Haz clic derecho en la capa hillshade en tu Lista de capas y abre el cuadro de diálogo Propiedades.
Observa el efecto cuando el sombreado de relieve transparente se sobre impone sobre el MDE recortado.
Haz clic en el elemento del menú Ráster ‣ Análisis de terreno.
Selecciona el tipo de análisis Pendiente, con el MDE acortado como la capa de entrada.
Especifica un nombre de archivo y localización apropiados para propósitos de salida.
Haz clic en Aceptar.
La imagen pendiente ha sido calculada y añadida al mapa. Sin embargo, como siempre, solo es un rectángulo gris. Para ver adecuadamente qué está pasando, cambia la simbología como sigue.
Abre el cuadro de diálogo Propiedades de la capa (como siempre, con el menú de clic derecho en la capa).
Donde dice Unibanda gris (en el menú desplegable Tipo de renderizador), cámbialo a Unibanda pseudocolor.
Asegúrate de que el botón radial Media +/- desviación estándard está seleccionado.
Utiliza el mismo enfoque que para calcular la pendiente, pero selecciona Orientación en el cuadro de diálogo inicial.
Recuerda guardar el mapa periódicamente.
Haz clic en el elemento del menú Ráster ‣ Calculadora ráster....
In the Raster bands list on the left, you will see all the raster layers in your Layers list. If your Slope layer is called slope, it will be listed as slope@1.
The slope needs to be between 15 and 60 degrees. Everything less than 15 or greater than 60 must therefore be excluded.
Utiliando la lista de elementos y botones en la interfaz, construye la siguiente expresión:
((slope@1 < 15) OR (slope@1 > 60)) = 0
Ajusta el campo Capa de salida a un nombre y localización adecuados.
Haz clic en Aceptar.
Now find the correct aspect (east-facing: between 45 and 135 degrees) using the same approach.
Construye la siguiente expresión:
((aspect@1 < 45) OR (aspect@1 > 135)) = 0
Find the correct rainfall (greater than 1200mm) the same way. Build the following expression:
(rainfall@1 < 1200) = 0
Having reclassified all the rasters, you will now see them displayed as gray rectangles in your map (assuming that they have been added to the map correctly). To properly display raster data with only two classes (1 and 0, meaning true or false), you will need to change their symbology.
Bajo el título Cargar valores min / max, selecciona el botón radial Real (más lento).
Haz clic en el botón Cargar.
The Custom min / max values fields should now populate with 0 and 1, respectively. (If they do not, then there was a mistake with your reclassification of the data, and you will need to go over that part again.)
En el título Mejora de contraste, ajusta la lista del menú desplegable Sin realce a Estirar a MinMax.
Haz clic en Aceptar.
Hazlo para los tres ráster reclasificados, ¡y recuerda guardar tu trabajo!
The only criterion that remains is that the area must be 250m away from urban areas. We will satisfy this requirement by ensuring that the areas we compute are 250m or more from the edge of a rural area. Hence, we need to find all rural areas first.
Hide all layers in the Layers list.
Haz visible la capa vectorial Zoning.
Haz clic derecho en ella y abre el cuadro de diálogo Filtrar.
Construye la siguiente consulta:
"Gen_Zoning" = 'Rural'
Mira antes las instrucciones para construir la consulta ara la capa Streets si te quedas atascado.
Cuando hayas terminado, cierra el cuadro de diálogo Constructor de consultas.
Deberías ver una colección de polígonos de la capa Zoning. Necesitarás guardarlos en un nuevo archivo de capa.
En el menú de clic derecho en Zoning, selecciona Guardar como....
Guarda tu capa en el directorio Zoning.
Haz clic en Aceptar.
Añade la capa a tu mapa.
Haz clic en el elemento del menú Vectorial ‣ Herramientas de geoproceso ‣ Disolver.
Selecciona la capa rural como tu capa vectorial de entrada, mientras dejas la casilla Usar sólo objetos espaciales seleccionados sin marcar.
Guarda tu capa en el directorio Zoning.
Cierra el cuadro de diálogo Disolver.
Borra las capas rural y Zoning.
Guarda el mapa.
Now you need to exclude the areas that are within 250m from the edge of the rural areas. Do this by creating a negative buffer, as explained below.
Haz clic en el elemento del menú Vectorial ‣ herramientas de geoproceso ‣ Buffer(s).
En el cuadro de diálogo que aparece, selecciona la capa rural_dissolve como tu capa vectorial de entrada (Usar sólo objetos espaciales seleccionados no debería estar seleccionado).
Selecciona la casilla Disolver resultados de buffer.
Ajusta el archivo de salida al mismo directorio que los otros archivos vectoriales rurales.
Haz clic en Guardar.
Haz clic en Aceptar y espera a que se complete el proceso.
Selecciona Sí en el cuadro de diálogo que aparece.
Cierra el cuadro de diálogo Buffer.
Borra la capa rural_dissolve.
Guarda el mapa.
Para incorporar las zonas rurales al mismo análisis con los tres ráster existentes, será necesario rasterizarlo. Pero para que los ráster sean compatibles para el análisis, necesitarán tener el mismo tamaño. Además, antes de rasterizar, necesitarás recortar el vector a la misma área que los ráster. Un vector solo puede ser recortado por otro vector, así que primero necesitaras crear un polígono de encuadre del mismo tamaño que los ráster.
En el título Tipo, selecciona el botón Poligono.
Haz clic en Aceptar.
Haz clic en Aceptar en el cuadro de diálogo Nueva capa vectorial.
Guarda el vector en el directorio Zoning.
Oculta todas las capas excepto la nueva capa bbox y uno de los ráster reclasificados.
Navega hasta el elemento del menú Ver ‣ Barras de herramientas y asegúrate de que Digitalización está seleccionado. Deberías ver en una barra de herramientas el icono de un lápiz o un koki. Es el botón Conmutar edición.
Haz clic en el botón Conmutar edición para entrar al modo edición. Esto te permite editar una capa vectorial.
Haz clic en el botón Añadir objeto espacial, que debería estar cerca del botón Conmutar edición. Puede que esté oculto tras un botón de flechas dobles; si lo está, haz clic en las flechas dobles para mostrar los botones ocultos de la barra de herramientas Digitalización.
Con la herramienta Añadir objeto espacial activada, haz clic izquierdo en las esquinas del ráster. Puede que necesites ampliar el zum con la rueda del ratón para asegurarte que es exacto. Para desplazarte sobre el mapa en este modo, haz clic y arrastra el mapa con el boton central o la rueda del ratón.
Para el cuarto y último punto, haz clic derecho para finalizar la forma.
Introduce un número arbitrario para la forma ID.
Haz clic en Aceptar.
Haz clic en el botón Guardar cambios.
Haz clic en el botón Conmutar edición para terminar tu sesión de edición.
Guarda el mapa.
Ahora que tienes un cuadro delimitador, puedes utilizarlo para recortar la capa rural buffer.
Asegúrate de que solo son visibles las capas bbox y rural_buffer, con el último por encima.
Haz clic en el elemento del menú Vectorial ‣ Herramientas de geoproceso > Cortar.
En el cuadro de diálogo que aparece, ajusta la capa vectorial de entrada a rural_buffer y la capa para cortar a bbox, con ambas casillas usar sólo objetos espaciales seleccionados sin marcar.
Pon el archivo de salida en el directorio Zoning.
Haz clic en Aceptar.
Cuando se requiera añadir la capa al TDC, haz clic en Sí.
Cierra el cuadro de diálogo.
Compara los tres vectores y ve los resultados por ti mismo.
Borra las capas bbox y rural_buffer, luego guarda tu mapa.
Ahora ya está lista para ser rasterizada.
Necesitarás especificar un tamaño de píxel para un nuevo ráster que crees, así que primero necesitarás saber el tamaño de uno de los ráster existentes.
Abre el cuadro de diálogo Propiedades de cualquiera de los tres ráster existentes.
Cambia a la pestaña Metadatos.
Cierra el cuadro de diálogo Propiedades.
Haz clic en el elemento del menú Ráster ‣ Conversión ‣ Rasterizar. Puede que recibas una alerta sobre conjunto de datos no soportado. Haz clic de todos modos he ignoralo.
Selecciona rural_clipped como tu capa de entrada.
Ajusta la localización de archivo de salida dentro del directorio Zoning.
Marca la casilla Cargar en la vista del mapa cuando se termine.
Haz clic en Aceptar.
El nuevo ráster se debería mostrar en tu mapa una vez se ha computado.
El nuevo ráster se verá como un rectángulo gris - puedes cambiar el estilo de la presentación como hiciste para reclasificar ráster.
Guarda tu mapa.
Now that you have all four criteria each in a separate raster, you need to combine them to see which areas satisfy all the criteria. To do so, the rasters will be multiplied with each other. When this happens, all overlapping pixels with a value of 1 will retain the value of 1, but if a pixel has the value of 0 in any of the four rasters, then it will be 0 in the result. In this way, the result will contain only the overlapping areas.
Haz clic en el elemento del menú Ráster ‣ Calculadora ráster....
Construye la siguiente expresión (con los nombres apropiados para tus capas, dependiendo de cómo las llamaste)
[Rural raster] * [Reclassified aspect] * [Reclassified slope] *
[Reclassified rainfall]
Set the output location to the Rasterprac directory.
Name the output raster cross_product.tif.
Ensure that the Add result to project box is checked.
Haz clic en Aceptar.
Cambia la simbología del nuevo ráster del mismo modo que ajustaste el estilo para los otros rásters reclasificados. El nuevo ráster ahora representa adecuadamente las áreas donde todos los criterios se satisfacen.
To get the final result, you need to select the areas that are greater than 6000m^2. However, computing these areas accurately is only possible for a vector layer, so you will need to vectorize the raster.
Haz clic en el elemento del menú Ráster ‣ Conversión ‣ Poligonizar.
Haz clic en Aceptar.
Cierra el cuadro de diálogo cuando el proceso esté completo.
All areas of the raster have been vectorized, so you need to select only the areas that have a value of 1.
Abre el cuadro de diálogo Query para el vector nuevo.
Construye esta consulta:
"DN" = 1
Haz clic en Aceptar.
Create a new vector file from the results by saving the candidate_areas vector after the query is complete (and only the areas with a value of 1 are visible). Use the Save as... function in the layer’s right-click menu for this.
Save the file in the Rasterprac directory.
Nombra al archivo candidate_areas_only.shp.
Guarda tu mapa.
Abre el menú de la nueva capa vectorial con clic derecho.
Seleccióna Abrir tabla de atributos.
Haz clic en el botón Conmutar el modo edición en la parte inferior de la tabla, o pulsa Ctrl+E.
Haz clic en el botón Abrir calculadora de campos en la parte inferior de la tabla, o pulsa Ctrl+M.
Under the New field heading in the dialog that appears, enter the field name area. The output field type should be an integer, and the field width should be 10.
En tipo Expresión de la calculadora de campos:
$area
Esto significa que la calculadora de campo calculará el área de cada polígono en la capa vectorial y creará una columna nueva (llamada area) con el valor calculado.
Haz clic en Aceptar.
Hacer lo mismo para otro nuevo campo llamado id. En tipo Expresión de la calculadora de campos:
$id
Eso asegura que cada polígono tiene una ID única para su identificación.
Haz clic de nuevo en Conmutar el modo edición, y guarda tus ediciones si es requerido.
Ahora que las áreas son conocidas:
Build a query (as usual) to select only the polygons larger than 6000m^2. The query is:
"area" > 6000
Save the selection as a new vector layer called solution.shp.
Ahora tienes tus áreas para tu solución, de las que tendrás que elegir la más proxima a tu casa.
Crea una nueva capa vectorial como antes, pero esta vez, selecciona el valor Tipo como Punto.
¡Asegúrate de que es el SRC correcto!
Acaba de crear la nueva capa.
Entra en el modo de edición (mientras la capa nueva está seleccionada).
Haz clic en el punto donde tu casa o otra residencia actual esté, utilizando las calles como guía. Puede que tengas que abrir otras capas para ayudarte a encontrar tu casa. Si no vives cerca, haz clic en algún lugar donde una casa podría ser concebible.
Introduce un número arbitrario para la forma ID.
Haz clic en Aceptar.
Guarda tu edición y sal del modo de edición.
Guarda el mapa.
Necesitarás encontrar los centroides (“centros de masa”) para los polígonos solución para decidir cual está más cerca de tu casa.
Especifica la capa de entrada como solution.shp.
Arrastra la nueva capa al principio de la lista para poder verla.
El tipo de matriz de salida debería ser lineal.
Ajusta una localización y nombre adecuadas para la salida.
Haz clic en Aceptar.
Abre el archivo en un editor de textos (o importalo a nuna hoja de cálculo). Observa qué ID está asociada con la Distance más corta. Puede que haya más de una con la misma distancia.
Esta es la respuesta final a la pregunta investigada.
Para tu presentación, incluye la capa semi transparente del relieve sombreado encima del ráster de tu elección (como MDE o slope raster, por ejemplo). Incluye también el polígono del área(s) solución más cercano, así como tu casa. Sigue las mejores prácticas de cartografía para crear tu mapa de salida.