Poligonos de thiessen en ArcGis

jul. 11, 2011 3 min de lectura

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:

Construido el shape de puntos, tendremos la vista que se muestra en la siguiente figura.

F1. Datos  Thiessen

Nos vamos a ArcToolbox seguido de Analysis Tools luego Proximity y por ultimo Create Thiessen Polygons.

F2. Crear poligonos de Thiessen

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.
F3. Ventana crear poligonos de thiessen

Damos clic en Ok para obtener el siguiente resultado.

F4. Poligonos creados

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….

F5. Spatial analyst - convert

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.
F6. thiessen convertido a raster

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.

F7. thiessen raster

Pedro Villegas

Ingeniero Civil, Master en Ingeniería con énfasis en hidrogeología