Autores
El objetivo del presente trabajo de investigación es presentar cálculos de valores desconocidos de una variable espacial a partir de otros valores conocidos, aplicados a variables ambientales continuas sobre un territorio concreto, este proceso es conocido como la interpolaci´on espacial de datos. La interpolación espacial es una parte de la geoestadística,útil en análisis de datos espaciales y modelos predictivos. Existen varios métodos de interpolación espacial los cuales permiten generar superficies continuas a partir de medidas en localizaciones puntuales. A demás en Interpolacion Espacial mediante Sistemas de Información Geográfica existen métodos como:
Referente a las estructuras espaciales tenemos dos análisis que evidencian la estructura espacial de una área de estudio a partir de disitintas distancias entre localizaciones las cuales son los análisis de vecindad y análisis por cuadrantes, en las estructuras espaciales se presenta el tema de autocorrelación espacial, la cual se basa en la afirmación que considera que en el espacio geográfico todo se encuentra relacionado con todo, pero los espacios más cercanos están más relacionados entre sí que con los más alejados.
INTERPOLACIÓN ESPACIAL:
Usualmente se considera a la interpolación espacial como parte integral del campo de la geoestadística. La interpolacion espacial se basa en el cálculo o la estimación de valores desconocidos de una variable espacial a partir de otros valores cuyo valor es conocido. Los procesos de interpolación espacial son válidos para cualquier variable ambiental continua sobre un territorio concreto. Como ejemplos podemos encontrar va- riables de temperatura, humedad, precipitaci´on, presi´on, contaminación acútica, contaminación lumínica. Para estimar los valores intermedios entre puntos de medición se llevan a cabo procesos de interpolaci´on espacial. Se utilizan los valores de mediciones de la red de estaciones y se realiza un cálculo del valor para el resto del territorio.
Clasificación de Métodos de Interpolación Espacial No todos los m´etodos trabajan igual ni aceptan las mismas premisas, pudiendo devolver resultados muy dispares según sus propias características. La clasificación de los métodos de interpolación puede realizarse en base a múltiples criterios. Así, podríamos hablar de métodos de interpolación: Método de Interpolación Determinista:Dichos métodos generan superficies continuas mediante el grado de similitud o suavizado.
Método de Interpolación Geoestadísticos:Dichos métodos generan superficies continuas a partir de las propiedades estadísticas de los datos de partida.
Globales o locales: Se utilizan todos los valores del área evaluada o solo una parte de ella (subconjunto).
Graduales o abruptas: Se dan según la continuidad y suavidad de la superficie resultante.
Exactos o aproximados: Según si respetan los valores de mediciones exactas de entrada para la interpolación o si, por el contrario, pueden ser alterados o suavizados para ajustarlos al modelo del conjunto.
Univariantes o multivariantes: según si admiten o no valores de múltiples variables de entrada para generar el modelo y la superficie de interpolación. En GIS, generalmente la distancia es la variable admitida para métodos de interpolación univarientes.
Métodos para la interpolación espacial en QGIS
Vecino más cercano (Nearest Neighbor): Se basa en la generación de polígonos de Voronoi. Los polígonos de Voronoi constituyen el método más básico y simple de interpolación vectorial.Este método de interpolación permite dividir el espacio en áreas equivalentes de dominio o influencia para cada uno de los puntos de medición de entrada. Los polígonos de Voronoi o de Thiessen quedan definidos mediante líneas que delimitan la región perteneciente al punto mas cercano. El perímetro de cada una de las regiones generadas es equidistante a todos los puntos vecinos.
TIN (Triangulated Irregular Network): Devuelve una superficie de triángulos formada a partir de la localización de una serie de vértices cuyos valores son conocidos. Los vértices se conectan mediante aristas para generar dicha red triangular interconectados.El método TIN tratará por tanto, de generar un conjunto de triángulos sobre el espacio que maximicen la relación área/perímetro.
IDW: Los puntos de muestreo se ponderan durante la interpolación. De esta manera, la influencia de un punto en relación con otros se reduce o disminuye a medida que aumenta la distancia entre ellos. En el método de interpolación IDW puede establecerse un valor de potencia, denominado coeficiente P de distancia que por defecto es 2. A mayor valor de P, mayor énfasis o peso asignado a los puntos cercanos a evaluar, resultando en una superficie estadística más abrupta. A menor valor de P, mayor ´enfasis en el conjunto de la muestra de valores, dando como resultado superficies más suavizadas.
Interpolación mediante Spline: Spline utiliza un método de interpolación que estima valores usando una función polinómica que minimiza la curvatura general de la superficie, lo que resulta en una superficie suave que pasa exactamente por los puntos de entrada. La herramienta Splines emplea funciones polinómicas distintas más acordes para cada tramo, adecuándose así a una superficie más suave, menos abrupta y uniforme («Interpolación espacial en QGIS», 2019).
Dentro del programa QGIS, localizamos y arrastramos las capas a utilizar.
Figura 1
Desde la Caja de herramientas Procesos, busque y localice la herramienta Interpolación ‣ Interpolación TIN. Haga doble clic para iniciarla. En la caja de diálogo Interpolación TIN, seleccione Arlington_Soundings_2007_stpl83 como la Capa vectorial,Elevation como el Atributo de interpolación. Luego clic en el icono + y ejecutar. Ahora, seleccione Islands_2004_550_stpl83 como la Capa vectorial, Elevation como el Atributo de interpolación. Luego clic en el icono Agregar. Ahora cambie el Tipo de la capa a Lineas de ruptura. Luego en Extensión clic en el … y seleccion el Boundary2004_550_stpl83. En Tamaño ráster salida, define el Tamaño píxel X y Tamaño píxel Y como 5. Luego clic en … bajo Interpolado para guardar la capa como elevation_tin.tif. Clic Ejecutar. Ahora se agregará una nueva capa elevation_tin al lienzo.
Figura 2: Capa Elevación
De la Caja de herramientas Procesos, busque y localice la herramienta GDAL ‣ Extracción ráster ‣ Cortar ráster por capa de máscara. Doble clic para iniciarla. En la caja de diálogo Cortar ráster por capa de máscara, seleccione elevation_tin como la Capa de entrada, Boundary2004_550_stpl83 como la Capa máscara. Luego clic en … bajo Clipped (mask) para guardar la capa como elevation_tin_cortado de mascara.tif. Clic en Ejecutar. Ahora se agregará una nueva capa elevation_tin_cortado de mascara al lienzo. Clic en el icono Abrir el panel de estilo de Capa.
Figura 3: Cortado Máscara
Ahora vamos a la Simbología como Pseudocolor banda única. Clic en la flecha en Rampa de color y seleccione Invertir rampa de color. Ingrese 0 en Precisión de etiqueta. Clic Clasificar. Luego vamos a Raster- elegimos Extracción y curvas de nivel
En la caja de diálogo seleccione la elevation_tin_cortado de máscara como la Capa de entrada, coloque 5 en el Intervalo entre línea de contorno. Luego clic en … bajo Contornos para guardar la capa como contorno.gpkg. Clic Ejecutar.
Ahora se agrega una nueva capa Contorno al lienzo. Clic en Propiedades- y colocar Etiquetas.
Figura 4:Interpolación General