Tutorial: Creación de un mapa del nivel de agua subterránea a partir de los datos de pozos y un Modelo Dígital de Elevación (DEM).
6. Calcular la densidad de pozos
Ya tenemos los pozos dentro del área de estudio (o bien dentro de un acuífero) donde calcularemos la densidad de pozos.
1. Abre la caja de procesamiento (Processing Toolbox) desde el menú principal, elige Processing | Toolbox.
2. En la caja de procesamiento (Processing Toolbox) selecciona Vector analysis | Count points in polygon.
3. En la ventana de diálogo Count Points in Polygon selecciona la capa Stampriet_boundary en la sección de Polygons y la capa boreholes_clipped en la sección de Points. Cambia en el campo Count field name a No of boreholes. Guarda la salida al archivo GeoPackage Count Stampriet_Data con el nombre de la capa Stampriet study area.
Usamos ese nombre de archivo de salida, porque la herramienta hace un copia del polígono Stampriet_boundary y agrega el campo "No of boreholes". Aún no tenemos la capa del área de estudio en el archivo GeoPackage, así que lo agregaremos de la siguiente forma.
4. Da click en el botón Run. Cierra (Close ) la ventana después del geoprocesamiento.
Considera que si la capa de entrada tiene múltiples polígonos (e.g. diferentes acuíferos) producirá el conteo por polígono. Puedes seguir los pasos para calcular la densidad de pozos en todos los polígonos.
5. Copia el estilo de la capa Stampriet_boundary a la capa Stampriet study area: en el panel Layers da click derecho sobre la capa Stampriet_boundary y selecciona Styles | Copy Style | All Style Categories. Entonces da un click derecho sobre la capa Stampriet study area y elige Styles | Paste Style | All Style Categories.
6. Borra la capa Stampriet_boundary del panel Layers.
Ahora, revisemos la tabla de estudio de la capa Stampriet study area .
7. En el panel Layers da un click derecho sobre la capa Stampriet study area y selecciona Open Attribute Table.
Observarás que hay 64 pozos en el área de estudio.
Ahora, calcularemos el área de la zona de estudio.
8. En la tabla de atributos da un click en el botón que active la edición.
9. Da un click en el icono que abre la calculadora de campos (Field calculator).
10. En la ventana de la calculadora de campos (Field Calculator) crea un nuevo campo (Create a new field) con el nombre del campo de salida (Output field name) como Area (km2) y el tipo de dato del campo de salida (Output field type) Decimal number (real).
11. Bajo la sección de la Expression (expresión) escoge la función en el panel medio de la parte Geometry | $area da un doble click sobre $area que lo agrega a la sección de la expresión, La variable $area
calcula el área del polígono en unidades del mapa. En nuestro caso en metros cuadrados. Completa la expresión de la siguiente forma:
$area / 1000000
Utiliza el botón que indica la división.
El resultado lucirá como la figura inferior.
El siguiente paso, es el cálculo de la densidad de pozos.
13. Da un click en el icono que abre la calculadora de campos (Field calculator).
14. En la ventana de diálogo de la calculadora de campos (Field Calculator ), verifica el checkbox Create a new field , en el campo Output field name utiliza el nombre Borehole density per km2 y el campo de salida (Output field type) de tipo Decimal number (real). En la sección de la expresión ingresa :
"No of boreholes" / "Area (km2)"
Recuerda que puedes agregar los campos, haciendo doble click en el nombre de los campos bajo la sección Fields and values en el medio del panel. Usa el botón que agrega el símbolo de división.
16. Habilitá la edición, activando el icono y guarda con el botón Save.
Calculemos la densidad de pozos en el área de estudio que es 0.2 per km2.
En la siguiente sección, agregaremos un Modelo Dígital de Elevación (MDE) al proyecto.