jueves, 9 de octubre de 2008

Importar MDT en Arcview 3.x : Primer formato ASCII

Tal y como comentaba en el primer post de la semana, vamos a ver como importamos los distintos formatos de MDT a ArcView 3.x .

El primer formato, y más extraño encontrarlo, es el simple listado de puntos X, Y, Z.
Un ejemplo de ese formato sería:
218338.00 3187217.00 1458.46
218340.00 3187217.00 1457.29
218342.00 3187217.00 1456.11
218344.00 3187217.00 1454.94
..........................................................
Lo primero que habría que fijarse es que elemento se está usando de separador decimal y cuál de separador entre coordenadas, en el ejemplo es el punto decimal y el tabulador y son los que normalmente hay que usar (es un error muy frecuente emplear como separador decimal la coma, y arcview 3.x no lo admite)

Lo segundo sería cuál es el paso de malla, en los datos de ejemplo se ve claramente que es 2 metros (ya que es la diferencia de coordenadas entre un punto y otro). Habría que verificar que en el eje Y también existe el mismo paso de malla.

Si trabajamos con ArcView 3.x, hay que tener además instaladas y cargadas las extensiones de Spatial Analyst y Xtools. El procedimiento de trabajo es el siguiente:
  • Abrir el fichero en un bloc de notas o similar (recomiendo el UltraEdit).
  • Remplazar la "," por "." en el caso de que el separador decimal no sea el punto.
  • Añadir una primera fila con los textos x y z, con el mismo separador que el resto de las columnas (recomiendo que este sea el tabulador).
  • Guardamos el fichero y abrimos ArcView.
  • En un proyeto nuevo, añadimos una tabla en formato TXT, y selecionaremos el fichero creado anteriormente.
  • Verificamos que se han cargado correctamente las tres columnas, y que el formato de las mismas es de número, para ello picamos en XTools/Información de la estructura de la tabla. Nos debe de salir Field_Decimal.
  • En una Vista nueva, cargamos un tema de Evetos(Vista/Añadir tema de eventos). Seleccionando el fichero anterir.
  • Con el tema seleccionado, pulsamos en Tema/Convert to Grid, le damos el nombre de fichero.
  • Seguidamente nos pide la extensión (output grid extent) del GRID, en este caso le damos Same as "Nombre del fichero TXT" y el tamaño de celda (CellSize) que se dará en función de la separación que deben de tener nuestros puntos (en nuestro caso al observar los valores señalados en el ejemplo vemos que es 2).
  • Posteriormente nos pide en que columna tiene los valores que pasamos a la celda (en este caso es la columna Z).
  • Al darle la columna, transforma el fichero de Texto a un Grid.
En el vídeo podemos ver los pasos realizados en Arcview.



En el sucesivos post, veremos el resto de formatos ASCII. Es posible descargar el ejemplo de fichero XYZ desde aquí.

He usado el servidor de archivos de Rapidshare. Para descargarlo se pulsa en FREE, se espera el tiempo señalado y se podrá descargar el fichero. La clave para descomprimir el fichero RAR es "ITACAS" sin las comillas.

2 comentarios:

manolo dijo...

Pido disculpas,la pregunta que realizé para este artículo la puse por error en el artículo"Mapas de Polución en Google Earth y Maps"

Moisés Martín dijo...

En relación a la pregunta que realizas:
una vez que el asigna a las coordenadas xyz el paso de malla, que es este caso es de 2, porque vuelves a seleccionar la z

La respuesta es bien sencilla, el 2 como indicas es el paso de malla, por lo tanto representa la separación que existe entre cada uno de los píxeles, en cambio la selección de la Z es el atributo que será representado en el GRID, en este caso es la altitud.

Ahora bien ese atributo en principio podría ser cualquier variable que cumpla dos condiciones:
1. Que se pueda medir
2. Que sea continua (altura, precipitación, temperatura, salinidad, humedad, etc)

Un saludo