El Blog colaborativo sobre Geotecnologías, GIS y Mapas para la comunidad GEO
  • Política de cookies
  • Más información
  • Contacto
geomapik
  • Proyecto
  • Contacto
Sin resultados
Ver todos los resultados
geomapik

Interpolación espacial en QGIS: métodos, procesos y evaluación

Raúl Estévez por Raúl Estévez
14 mayo, 2019
en Análisis GIS, Tutoriales
Home Análisis GIS
12k
Vistas
CompartirTwitter

La interpolación espacial es un proceso crucial en múltiples campos de trabajo e investigación. Al trabajar con GIS y tener que evaluar variables con componente espacial, el hecho realizar cualquier tipo de interpolación espacial correctamente es crítico para la calidad del estudio.

Este artículo se centra en definir qué es la interpolación espacial, qué métodos existen, cómo se clasifican y cuáles son los más adecuados. Asimismo, se proponen distintos métodos disponibles de interpolación espacial en QGIS para comparar los resultados, usando como ejemplo un conjunto de datos climáticos.

¿Qué es la interpolación espacial?

Usualmente se considera a la interpolación espacial como parte integral del campo de la geoestadística. La interpolación 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.

qgis interpolacion espacial geoestadística

Los procesos de interpolación espacial son válidos para cualquier variable ambiental continua sobre un territorio concreto. Como ejemplos podemos encontrar variables de temperatura, humedad, precipitación, presión, contaminación acústica, contaminación lumínica…

Las variables climáticas son las que más frecuentemente requieren de esta clase de análisis. En estos casos, la interpolación espacial es algo crucial ya que no existen infinitas estaciones meteorológicas que permitan conocer el valor exacto de temperatura, humedad o precipitación sobre todos y cada uno de los puntos de un territorio.

Para estimar los valores intermedios entre puntos de medición se llevan a cabo procesos de interpolación 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.

El resultado obtenido mediante el proceso de interpolación se denomina habitualmente como superficie estadística. Se trata, pues, de una superficie continua con valores interpolados a partir de otros previamente conocidos.

Clasificación de métodos de interpolación espacial

A pesar de que todos los métodos de interpolación espacial tienen el mismo principio básico, no todos los métodos 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:

  • Globales o locales, según si utilizan todos los valores del área evaluada o sólo una parte de ella (subconjunto).
  • Graduales o abruptas, 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.
  • Determinísticos o estocásticos, según si incorporan o no variaciones aleatorias (incertidumbre) en la superficie interpolada. Los métodos determinísticos son aplicables cuando hay mediciones suficientes para describirla matemáticamente, mientras que los estocásticos incorporan el concepto de aleatoriedad por una insuficiencia de ellas.

Métodos para la interpolación espacial en QGIS

Las posibilidades en cuanto a métodos de interpolación en QGIS son bastante amplias, aunque mayormente pertenecientes al grupo de los métodos determinísticos definido anteriormente.

Asimismo, en QGIS podemos adecuar cada uno de los métodos y ajustar sus parámetros según las necesidades al realizar el geoproceso para generar las capas de interpolación en base a nuestros puntos con valores de medición.

Los principales métodos de interpolación posibles reunidos en las distintas librerías de las que se nutre este GIS de escritorio son los siguientes.

Vecino más cercano (Nearest Neighbor)

La interpolación mediante el método de vecino más cercano 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. El método se basa únicamente en la distancia euclidiana, obviando cualquier tipo de valor asignado a los puntos de muestreo.

interpolacion espacial nearest neighbor QGIS
Resultado obtenido mediante interpolación con método Nearest Neighbor en QGIS.

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 más cercano. El perímetro de cada una de las regiones generadas es equidistante a todos los puntos vecinos

Finalmente, el método asigna a cada polígono el valor del punto que contiene y a partir del cual se ha generado. Al tratarse de un método basado únicamente en la distancia, las variables interpoladas pueden ser tanto cualitativas como cuantitativas.

Este método resulta útil, por ejemplo, en estudios de geomarketing para conocer posibles áreas de influencia, o para realizar regionalizaciones o divisiones territoriales proporcionales.

TIN (Triangulated Irregular Network)

Éste método de interpolación 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.

QGIS interpolacion espacial TIN
Resultado de la interpolación mediante el método TIN en QGIS.

El resultado obtenido, la superfície TIN, es una malla o red de triángulos interconectados, donde cada uno de ellos representa una zona homogénea en lo que a la variable estudiada se refiere. 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.

Es muy habitual su uso sobretodo para modelos del terreno en base a mediciones de elevación conocida, aunque puede aplicarse a otras mediciones cuantitativas de distintas variables ambientales.

Interpolación IDW (Inverse Distance Weighting)

Mediante el método de interpolación 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.

QGIS metodo interpolacion idw
Resultado de la interpolación utilizando el método IDW en QGIS. Variando el coeficiente P se obtendrán superficies de interpolación ligeramente distintas.

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 énfasis en el conjunto de la muestra de valores, dando como resultado superficies más suavizadas.

Generalmente se utiliza en procesos de interpolación donde el conjunto de datos disponible para la interpolación es abundante, se reparte homogéneamente por el espacio y no existen grandes distancias entre sus localizaciones.

Interpolación mediante Spline

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

interpolacion espacial spline QGIS
Resulttado de la interpolación mediante Spline en QGIS. Simbología ráster categorizada según diversas clases.

Como decíamos antes en el apartado de clasificación, este método,  juntamente con Kriging, es uno de los métodos exactos de interpolación existentes que no admite aproximaciones o suavizados de los valores de entrada.

Matemáticamente, 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.

Cómo realizar una interpolación en QGIS a partir de una capa de puntos

Todos los métodos de interpolación disponibles en QGIS pertenecen a las diversas librerías que trae instaladas por defecto: GDAL, GRASS y SAGA.

Los geoprocesos de interpolación utilizados en este artículo para cada uno de los métodos, juntamente con los parámetros especificados en las imágenes, son los siguientes:

  • Método Nearest Neighbor: Cuadrícula (vecino más próximo) de la librería GDAL
  • Método IDW: Interpolación IDW de la librería de QGIS.
  • Método TIN: Interpolación TIN de la librería QGIS.
  • Método Spline: v.surf.bspline de la librería GRASS

Para ejecutar cualquiera de estos geoprocesos de interpolación podemos buscarlos directamente en la caja de herramientas de geoprocesos.

Configuración de parámetros para el método Spline en QGIS.

Cada uno de los métodos de interpolación de QGIS tiene su propio formulario único para especificar las capas de entrada y salida, así como toda una serie de parámetros a configurar para adecuar el modelo.

Parámetros dados para la interpolación con método IDW en QGIS.

Puedes descargar aquí los datos climáticos utilizados en este artículo para realizar pruebas de interpolación con las distintas variables almacenadas. En el caso de este artículo, hemos trabajado la interpolación con la variable de horas de luz solar en Reino Unido.

Se trata de un conjunto de valores medios de precipitación, temperaturas y horas de Sol en un mes concreto registradas en las principales estaciones meteorológicas del Reino Unido. Los datos han sido obtenidos de la Agencia Estatal de Meteorología MetOffice.

¿Cómo elegir el mejor método de interpolación en QGIS?

La elección del método de interpolación más apropiado dependerá de la propia naturaleza del conjunto de datos de muestreo que deseemos interpolar.

Así pues, deberíamos conocer previamente:

  • Tipo de variable a interpolar: cuantitativa o cualitativa y la lógica o necesidad de aplicación de cada método para cada tipo de variable.
  • Características estadísticas de la muestra: valores máximos, mínimos, media y mediana, desviación estándar… Puedes consultar el siguiente artículo para saber cómo realizar un análisis estadístico exploratorio con QGIS.
  • Distribución espacial de la variables: homogeneidad espacial de los muestreos, distancia media entre puntos de muestreo…
  • Existencia de valores anómalos, clústeres y hotspots o coldspots que puedan inferir o alterar la superficie de interpolación.

Evaluación del método y estimación del error de interpolación

Posteriormente será de gran utilidad realizar una evaluación del método y el error obtenido en la interpolación.

Podría realizarse mediante, por ejemplo, la comparación de determinados valores reales sustraídos de la muestra con los valores estimados obtenidos mediante la interpolación en esos mismos puntos.

Posteriormente, debería reajustarse la configuración de los parámetros del método para acercar el modelo lo máximo posible a los valores de las mediciones reales.

Etiquetas: Análisisaprender GISClimageoestadísticaGeomarketingInterpolaciónQGISSistemas de Información Geográfica
Compartir383Tweet240
AVISO
Post anterior

Cómo obtener valores de una capa ráster en QGIS

Post siguiente

Análisis de cuenca visual en QGIS

Raúl Estévez

Raúl Estévez

Geógrafo y Geoinformático. 8 años de experiencia. Apasionado de los GIS y la geotecnología, el análisis espacial y la cartografía digital. Especializado en gestión, análisis, integración y visualización de geoinformación. No entiendo el mundo sin los mapas.

Post siguiente
calcular cuenca visual qgis

Análisis de cuenca visual en QGIS

Mapas con Python utilizando geopandas y matplotlib

Comentarios 4

  1. Rosana says:
    2 años atrás

    Buenos días,
    Quisiera saber cómo coincidir la interpolación con una zona geográfica georeferenciada, tal como muestra en la imagen «Resultado de la interpolación utilizando el método IDW en QGIS. Variando el coeficiente P se obtendrán superficies de interpolación ligeramente distintas», ubicada más arriba en esta web. Muchas gracias

    Responder
    • Raúl Estévez says:
      2 años atrás

      Hola Rosana,

      gracias por tu comentario y tu interés. Entiendo por tu pregunta que lo que requieres es «recortar» el área para que se ajuste a unos límites geográficos o administrativos concretos.

      Si ese es el caso, deberías utilizar un geoproceso de QGIS denominado «Cortar ráster por capa de máscara» en caso de disponer de una capa vectorial que defina dichos límites, o bien «Cortar ráster por extensión». Puedes acceder a ellos desde el menú ráster > extracción.

      Saludos y espero que resuelva tu duda.

      Responder
  2. Jessica Cohen says:
    2 años atrás

    Hola buenas tardes. Me gustaría saber cual herramienta se puede usar para crear en qgis un análisis estadístico multivariado de tal modo que pueda llevar un registro mes a mes de las variables involucradas, para ver evolutivamente los cambios. Espero haberme explicado y gracias de antemano.

    Responder
  3. Matias Castro says:
    2 meses atrás

    Buenas tardes, que recomendarías para un estudio de ruido ambiental de un municipio? Suponiendo unos 70 puntos de muestreo a 400m o 500m de distancia entre cada punto.

    Responder

Deja una respuesta Cancelar la respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Posts Recomendados

GBIF datos en QGIS descargar

El plugin de GBIF en QGIS para descargar datos de distribución de especies

27 marzo, 2019
1.1k
vistas-mapas-3d-con-QGIS-portada

Cómo generar y visualizar mapas 3D con QGIS

3 enero, 2019
10.8k
Geotecnologías GIS Mapas Cartografía Digital

Geotecnologías, GIS y Mapas

Información, recursos de aprendizaje y novedades del sector GEO para los usuarios especializados en el manejo de la información geográfica, su análisis y su representación.

Posts Nuevos

python visualizacion datos

Librerías Python para visualización de datos

21 febrero, 2023
116
geopandas leer archivos GIS

Cómo leer archivos GIS con Geopandas

20 febrero, 2023
105
Join espacial PostGIS

Cómo realizar un join espacial en PostGIS

19 febrero, 2023
85
Overture Maps. ¿Rival para Google Maps?

Overture Maps ¿una nueva plataforma de mapas con datos abiertos?

19 diciembre, 2022
271
El futuro de los GIS según OpenAI

El futuro de los GIS y el sector geoespacial según OpenAI

9 diciembre, 2022
255
Geotecnologías GIS Mapas Cartografía Digital

El blog colaborativo sobre Geotecnologías, GIS y Mapas para la comunidad GEO

Información, recursos de aprendizaje y novedades del sector GEO para los usuarios especializados en el manejo de la información geográfica, su análisis y su representación.

Todo reunido en un blog colaborativo y abierto. ¿Te animas a colaborar con geomapik?

Categorías

  • Análisis GIS
  • Desarrollo GIS
  • Novedades GIS
  • Recursos GIS
  • Spatial Data Science
  • Tutoriales
  • Uncategorized
  • Webmapping
Sin resultados
Ver todos los resultados
Spatial Data Science

Librerías Python para visualización de datos

21 febrero, 2023
116
Spatial Data Science

Cómo leer archivos GIS con Geopandas

20 febrero, 2023
105
Análisis GIS

Cómo realizar un join espacial en PostGIS

19 febrero, 2023
85
Novedades GIS

Overture Maps ¿una nueva plataforma de mapas con datos abiertos?

19 diciembre, 2022
271

Etiquetas

Análisis (16) Análisis de datos (6) Análisis geoespacial (3) aprender GIS (5) ArcGIS (5) Bases de Datos Espaciales (5) Cartografía (29) Data Analysis (6) Data Science (6) Datos (9) Desarrollo GIS (16) Diseño cartográfico (6) ESRI (3) geoestadística (2) Geoinformática (8) GeoJSON (2) Geopackage (2) Geopandas (4) Geoproceso (2) Geotecnología (5) GIS (44) Javascript (4) Jupyter Notebook (3) Mapas (20) Matplotlib (2) MDT (3) Medio Ambiente (2) OpenStreetMap (3) pgAdmin (2) plugins (7) PostGIS (8) PostgreSQL (7) Programación (17) Python (10) QGIS (25) R (5) Shapefile (2) Sistemas de Información Geográfica (34) Spatial Data Science (4) SQL (7) Tutorial (14) visor de mapas (3) visualización (3) Webmapping (11) WKT (2)

© 2018 geomapik - Geotecnologías, GIS y Mapas

Sin resultados
Ver todos los resultados
  • Proyecto
  • Contacto

© 2018 geomapik - Geotecnologías, GIS y Mapas

Este sitio web utiliza cookies para que usted tenga la mejor experiencia de usuario. Si continúa navegando está dando su consentimiento para la aceptación de las mencionadas cookies y la aceptación de nuestra política de cookies, pinche el enlace para mayor información.plugin cookies

ACEPTAR
Aviso de cookies