Poligonos de thiessen en ArcGis
El método de los polígonos de Thiessen consiste en delimitar áreas de influencia (unidades discretas) a partir de un conjunto de puntos. El tamaño y la configuración de los polígonos depende de la distribución de los puntos originales. Una limitante que tiene el método es que no se puede estimar el error asociado, pues el valor para cada polígono se obtiene a partir de un solo punto.
Para ilustrar la aplicación de los polígonos de Thiessen utilizaremos la precipitación correspondiente al año 1985 registradas en estaciones del IDEAM en el golfo de Urabá, tal como se muestra en la siguiente Tabla.
| ESTACION | NOMBRE | LATITUD (N) | LONGITUD (W) | PP 1985 |
| 1201007 | EL CASCO | 7,88 | -76,6 | 1679,6 |
| 1201502 | UNIBAN | 7,82 | -76,7 | 2301,6 |
| 1201006 | LA LORENA | 7,85 | -76,7 | 2167 |
| 1201507 | APTO LOS CEDROS | 7,82 | -76,7 | 2993,5 |
| 1201506 | TULENAPA | 7,77 | -76,7 | 2326,5 |
| 1201009 | TRIGANA | 7,73 | -76,7 | 2197 |
| 1201005 | TOSCANA | 7,8 | -76,7 | 2353,5 |
| 1201509 | LA PALMERA | 7,53 | -76,6 | 4031,9 |
| 1201003 | BARRANQUILLITA | 7,57 | -76,7 | 3292 |
| 1202504 | TURBO | 8,08 | -76,7 | 1919,7 |
| 1202001 | PUEBLO BELLO | 8,2 | -76,5 | 1293 |
| 1201017 | NUEVO ORIENTE | 7,45 | -76,7 | 3684,4 |
| 1201012 | SANTA MARTA | 7,88 | -76,6 | 2314 |
| 1201011 | PRADOMAR | 7,98 | -76,6 | 2260,1 |
| 1201001 | EL TORMENTO | 7,75 | -76,9 | 2989 |
En primer lugar se debe realizar un shape de puntos con los datos registrados en las estaciones. Sobra ilustrar la metodología para crear un shape de puntos con ArcGis, los interesados pueden revisar los siguientes artículos:
- Shape de puntos en ArcGIS a partir de coordenadas geográficas.
- Cómo crear un shape de puntos a partir de coordenadas geográficas.
- Shape de puntos a partir coordenadas XYZ en ArcGIS
Construido el shape de puntos, tendremos la vista que se muestra en la siguiente figura.
Nos vamos a ArcToolbox seguido de Analysis Tools luego Proximity y por ultimo Create Thiessen Polygons.
Aparecerá la siguiente ventana donde.
- Input Features: introducimos el Feature que contiene los datos de precipitación
- Output Feature Class: Introducimos el nombre del archivo de salida
- Output Fiels (optional): Indicamos All.
Damos clic en Ok para obtener el siguiente resultado.
Para obtener los valores de precipitación media de estos polígonos, los tenemos que convertir un raster, para ello seleccionamos Spatial Analyst seguido de Convert finalmente features to Raster….
En la ventana que aparece, especificamos la siguiente información
- Input Features: Elegimos PoligonosThiessen
- Field: elegimos el campo que contiene los datos de precipitación, en este caso PP1985.
- Ouput cell size. Especificamos el tamaño de la celda del raster de salida. Colocamos 100.
- Output raster: colocamos el nombre del archivo de salida. Colocamos thiessen.
El resultado debe ser una red que da en cada celda el valor de precipitación anual de los correspondientes polígonos de Thiessen, tal como se muestra en la siguiente figura.