-->

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:

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

15 comentarios :

  1. Sensacional blog!!!

    Muchas gracias por tu dedicación!!

    ResponderEliminar
  2. Hola, una pregunta, si no tengo disponible la opción de crear polígonos es cuestión de mi programa (tengo la versión 9.3) o mejor dicho, que arcgis tienes y si puedes darme una solución a este detallito te lo agradeceré mucho.

    ResponderEliminar
  3. Muy buen artículo guía, muy ilustrativo y claro,
    sigan adelante

    ResponderEliminar
  4. Muchas gracias por tu comentario!

    ResponderEliminar
  5. definitivamente esta muy bueno

    ResponderEliminar
  6. Hola... hice exactamente como lo haces... solo que me me sale el error 999998 unexpected error. failed to execute (createThiessenPolygons).... podrías decirme que esta afectando o que podría hacer.

    ResponderEliminar
    Respuestas
    1. Hola Alex, no tengo información sobre este error, chequea bien los datos, verifican que no existen datos en blanco, verifica que todos sean numéricos...

      Eliminar
  7. Saludos!
    Como puedo ver una tabla donde se muestren las áreas de cada poligono

    ResponderEliminar
    Respuestas
    1. Lo primero que debes hacer es convertir los poligos de thiessen de raster a shape. luego revisas la tabla de atributos del shape y alli debe aparecer el área

      Eliminar
  8. Como calculo la precipitacion media de la cuenca teniendo el raster ??

    ResponderEliminar
    Respuestas
    1. Revisa estos dos articulos quizá te puedan ayudar
      http://www.aguaysig.com/2011/07/calculando-caudales-medios-en-arcgis.html
      http://www.aguaysig.com/2011/07/estimacion-de-caudales-medios-con.html

      Eliminar
    2. Peter tendras informacion de metodo de isoyetas y calcular la precipitacion promedio

      Eliminar
  9. Exelente trabajo, gracias por compartir tu conocimiento. un abrazo desde Bogota.

    ResponderEliminar
  10. Hola, hago todo lo anterior y no se me crean los polígonos, me dice error 000537, error in converting points to voronoi polygons. Failed to execute.

    ResponderEliminar
  11. tendras informacion de metodo de isoyetas.gracias

    ResponderEliminar

Tu comentario es muy importante...¿Cómo te pareció el tema?