9.4. Lesson: Ejercicio Suplementario

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)

9.4.1. Planteamiento del Problema

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.

9.4.2. Esquema de la Solución

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.

9.4.3. Preparando el Mapa

  • 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).

  • Select the entry WGS 84 / UTM zone 33S (with the EPSG code 32733).
  • Click OK. The map is now in the UTM33S coordinate reference system.
  • 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....

  • Save the map in a directory called Rasterprac that you should create somewhere on your computer. You will save whatever layers you create in this directory as well.

9.4.4. Cargar Datos dentro del Mapa

Para procesar los datos, necesitarás cargar las capas necesarias (nombres de calles, zonas, precipitación, MDE) en la vista del mapa.

9.4.4.1. For vectors ...

  • Click on the Add Vector Layer button, or use the Layer ‣ Add Vector Layer... menu item.
  • In the dialog that appears, ensure that the File radio button is selected.
  • Click on the Browse button.
  • In the dialog that appears, open the exercise_data/more_analysis/streets directory.
  • Select the file Street_Names_UTM33S.shp.
  • 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.

  • Click Open. The vector layer will load in your map. Its color is automatically assigned. It will be changed later.
  • Rename the layer to Streets.
  • Right-click on it in the Layers list (by default, the pane along the left-hand side of the screen).
  • Click Rename in the dialog that appears and rename it, pressing the Enter key when done.
  • Repeat the vector adding process, but this time select the Generalised_Zoning_Dissolve_UTM33S.shp file in the Zoning directory.
  • Rename it to Zoning.

9.4.4.2. For rasters ...

  • Click on the Add Raster Layer button, or use the Layer ‣ Add Raster Layer... menu item.
  • Navega hasta el archivo apropiado, selecciónalo y haz clic en Abrir.

  • Do this for each of the two raster files. The files you want are DEM/reproject/DEM and Rainfall/reprojected/rainfall.tif.
  • Rename the rainfall raster to Rainfall (with an initial capital). Initially when you load them, the images will be gray rectangles. Don’t worry, this will be changed later.
  • Guarda el mapa.

Para ver correctamente qué está pasando, se necesita cambiar la simbología para las capas.

9.4.5. Cambio de simbología de capas vectoriales

  • In the Layers list, right-click on the Streets layer.
  • Selecciona Propiedades del menú que aparece.

  • Switch to the Style tab in the dialog that appears.
  • 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.

9.4.6. Cambio de simbología de capas ráster

La simbología de capas ráster es algo diferente.

  • Abre el cuadro de diálogo Propiedades para el ráster Precipitación.

  • Switch to the Style tab. You’ll notice that this style dialog is very different from the version used for vector layers.
  • Asegúrate de que el botón Media +/- desviación estándar está seleccionado.

  • Change the value in the associated box to 2.00 (it should be set to 0.00 by default).
  • 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.

  • Repeat this process for the DEM, but set the standard deviations used for stretching to 4.00.

9.4.7. Cambio de orden de capas

  • In the Layers list, click and drag layers up and down to change the order they appear in on the map.
  • Newer versions of QGIS may have a Control rendering order checkbox beneath the Layers list. Ensure that it is checked.

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.

9.4.8. Encuentra los Distritos Correctos

  • Load the vector layer admin_boundaries/Western_Cape_UTM33S.shp into your map.
  • Rename it to Districts.
  • Right-click on the Districts layer in the Layers list.
  • 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:

  • Bellville,
  • Cape,
  • Goodwood,
  • Kuils River,
  • Mitchells Plain,
  • Simons Town, and
  • Wynberg.
  • 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.

  • Click the = button; an = sign is added to the SQL query.
  • 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.

9.4.9. Recorta los Ráster

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....

  • Navigate to your Rasterprac directory.
  • 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.

  • Check the Load into canvas when finished box.
  • 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.

  • Do not change any other options. Do not alter the existing clipping area which you drew previously. Leave everything the same and click OK.
  • Después de que la segunda operación de corte haya terminado, puedes cerrar el cuadro de diálogo Clipper.

  • Guarda el mapa.

9.4.10. Limpia el mapa

  • Remove the original Rainfall and DEM layers from the Layers list:
  • 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:

    • Right-click on the layer in the Layers list.
    • 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.

  • Change the raster symbology as before (see Changing the symbology of raster layers).
  • Guarda el mapa.

  • You can now hide the vector layers by unchecking the box next to them in the Layers list. This will make the map render faster and will save you some time.

In order to create the hillshade, you will need to use a plugin that was written for this purpose.

9.4.11. Activating the Raster Terrain Analysis plugin

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:

  • Click on the menu item Plugins ‣ Manage Plugins....
  • Ensure that the box next to Raster Terrain Analysis plugin is selected.
  • 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.

9.4.12. Crear el sombreado del relieve

  • In the Layers list, ensure that the DEM is the active layer (i.e., it is highlighted by having been clicked on).
  • Click on the Raster ‣ Terrain analysis ‣ Hillshade menu item to open the Hillshade dialog.
  • Especifica una localización apropiada para la capa de salida y nómbrala hillshade.

  • Check the Add result to project box.
  • 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.

  • Click on the Transparency tab and set the transparency slider to 80%.
  • Click OK on the dialog.
  • Observa el efecto cuando el sombreado de relieve transparente se sobre impone sobre el MDE recortado.

9.4.13. Pendiente

  • 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.

  • Check the Add result to project box.
  • 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).

  • Click on the Style tab.
  • 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.

9.4.14. Orientación

  • 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.

9.4.15. Reclasificar rásters

  • Haz clic en el elemento del menú Ráster ‣ Calculadora ráster....

  • Specify your Rasterprac directory as the location for the output layer.
  • Ensure that the Add result to project box is selected.

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.

9.4.16. Ajuste del estilo para capas reclasificadas

  • Open the Style tab in the layer’s Properties dialog as usual.
  • 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.

9.4.17. Encontrar áreas rurales

  • 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.

  • Name the output file rural.shp.
  • 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.

  • Under Dissolve field, select — Dissolve all —.
  • Guarda tu capa en el directorio Zoning.

  • Click OK. A dialog will appear asking whether you want to add the new layer to the TOC (“Table of Contents”, referring to the Layers list).
  • Click Yes.
  • 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.

9.4.18. Crear un buffer negativo

  • 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).

  • Select the Buffer distance button and enter the value -250 into the associated field; the negative value means that the buffer must be an internal buffer.
  • Selecciona la casilla Disolver resultados de buffer.

  • Ajusta el archivo de salida al mismo directorio que los otros archivos vectoriales rurales.

  • Name the output file rural_buffer.shp.
  • Haz clic en Guardar.

  • Haz clic en Aceptar y espera a que se complete el proceso.

  • Selecciona 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.

9.4.19. Crear un encuadre vectorial

  • Click on the menu item Layer ‣ New ‣ New Shapefile Layer....
  • En el título Tipo, selecciona el botón Poligono.

  • Click Specify CRS and set the coordinate reference system WGS 84 / UTM zone 33S : EPSG:32733.
  • Haz clic en Aceptar.

  • Haz clic en Aceptar en el cuadro de diálogo Nueva capa vectorial.

  • Guarda el vector en el directorio Zoning.

  • Name the output file bbox.shp.
  • Oculta todas las capas excepto la nueva capa bbox y uno de los ráster reclasificados.

  • Ensure that the bbox layer is highlighted in the Layers list.
  • 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.

9.4.20. Recortar una capa vectorial

  • 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.

  • Name the output file rural_clipped.
  • Haz clic en Aceptar.

  • Cuando se requiera añadir la capa al TDC, haz clic en .

  • 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.

9.4.21. Rasterizar una capa vectorial

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.

  • Make a note of the X and Y values under the heading Dimensions in the Metadata table.
  • 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.

  • Name the output file rural_raster.tif.
  • Check the New size box and enter the X and Y values you made a note of earlier.
  • Marca la casilla Cargar en la vista del mapa cuando se termine.

  • Click the pencil icon next to the text field which shows the command that will be run. At the end of the existing text, add a space and then the text -burn 1. This tells the Rasterize function to “burn” the existing vector into the new raster and give the areas covered by the vector the new value of 1 (as opposed to the rest of the image, which will automatically be 0).
  • 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.

9.4.22. Combinación de rásters

  • 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.

9.4.23. Vectorizar el ráster

  • Haz clic en el elemento del menú Ráster ‣ Conversión ‣ Poligonizar.

  • Select the cross_product raster.
  • Set the output location to Rasterprac.
  • Name the file candidate_areas.shp.
  • Ensure that Load into canvas when finished is checked.
  • 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.

9.4.24. Cálgulo del área para cada polígono

  • 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.

9.4.25. Selección de áreas para un tamaño dado

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.

9.4.26. Digitaliza 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!

  • Name the new layer house.shp.
  • 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.

9.4.27. Calcula los centroides de los polígonos

  • Click on the Vector ‣ Geometry Tools ‣ Polygon centroids menu item.
  • Especifica la capa de entrada como solution.shp.

  • Provide the output location as Rasterprac.
  • Call the destination file solution_centroids.shp.
  • Click OK and add the result to the TOC (Layers list), then close the dialog.
  • Arrastra la nueva capa al principio de la lista para poder verla.

9.4.28. Calcula qué centroide está más cerca de tu casa

  • Click on the menu item Vector ‣ Analysis Tools ‣ Distance matrix.
  • The input layer should be your house, and the target layer solution_centroids. Both of these should use the id field as their unique ID field.
  • 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.

  • Build a query in QGIS to select only the solution areas closest to your house (selecting it using the id field).

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.