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.

• Como 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.