Fundamentos de análisis geográfico con SEXTANTE - La Salle

240 CAPÍTULO 13. INTERPOLACIÓN. GEOESTADÍSTICA
13.2. Interpolación de datos por distancia inversa (IDW)
El más sencillo de los módulos de interpolación es el denominado Distancia Inversa (puedes
encontrar estos módulos en el menú Grid/Raster–Vectorial/Interpolación desde Puntos)
Selecciona la capa de puntos que deseas interpolar y en Campo escoge cuál de sus campos
de atributos contiene la variable a utilizar (sólo puedes interpolar una variable, ya que sólo se
creará una única capa raster). Por supuesto, el campo escogido tiene que ser de tipo numérico.
En el grupo Opciones encontrarás cuatro campos. Uno de ellos es el conocido Dimensiones
del grid resultante, con el que ya hemos trabajado y que no requiere mayor explicación, así como
los parámetros de las sucesivas ventanas de parámetros que se hayan relacionadas con el mismo.
Los tres campos restantes controlan el comportamiento del algoritmo de interpolación. Para
cada celda de la capa resultante, un área de radio fijo es ✭✭rastreada✮✮ alrededor de la misma en
busca de puntos con valores. Puedes ajustar dicho radio (en las unidades propias de la capa),
introduciendo su valor en el campo Radio de búsqueda. Trata de utilizar el radio mínimo que
asegure que alrededor de todas las celdas se encuentra un número suficiente de puntos. Sólo
los puntos dentro del radio de búsqueda se utilizan, por lo que, si éste es insuficiente, puede
darse el caso de que ningún punto caiga dentro y la interpolación no se pueda realizar, dando
como resultado una celda sin datos.
Valores muy elevados del radio de búsqueda hacen que el módulo requiera un mayor tiempo
de ejecución (y en ocasiones pueden no aumentar la precisión del cálculo, sino todo lo contrario),
por que debes encontrar un equilibrio adecuado. Prueba diversos ajustes hasta encontrar
uno correcto.
En el campo Número Máximo de Puntos puedes introducir la cantidad máxima de puntos
que deseas utilizar. si dentro del radio de búsqueda hay más puntos de los especificados (sea
n), entonces únicamente los n más cercanos serán utilizados por el algoritmo de interpolación.
Para explicar el significado del último término es necesario algo de matemática. Si se toman
n puntos alrededor de una celda para proceder a la interpolación, el valor asignado a dicha
celda mediante el método de distancia inversa es el resultante de la siguiente expresión.
z =
n
i=1 zid k i
n
i=1 dk i
(13.1)
donde zi es el valor del punto y di la distancia desde ese punto hasta la celda interpolada
El valor de k se puede ajustar en el campo Exponente. 1 y 2 son valores típicos que sirven
para la mayor parte de los casos.
Interpolando la capa de puntos creada a partir de una tabla de coordenadas (si no lo
recuerdas, vuelve a 11.6 para ver cómo lo hicimos) se obtiene el siguiente resultado (he utilizado
un tamaño de celda de 0.1 y he añadido unas curvas de nivel para que puedas interpretar más
claramente el ✭✭relieve✮✮ en la representación de mapa).