De Shape a KMl con herramienta Google

Para convertir un archivo Shape a KML puede también emplearse el programa de Google gSpatial Toolbar. Lo puedes encontrar aquí.

La exigencia es que el sistema de coordenadas de origen debe ser WGS84, Lat Long.

Extensión WWF para Hidrología. Arc View

He subido las extensiones de hidrología para Arcview a Media Fire.

Las puedes encontrar aquí .

Si tienes problemas de acceso avisa.

Probando Stitchmaps

Compré el programa stitchmaps, ilusionado con la idea de georreferenciar rápidamente las imágenes de Google Earth... pero mi primera impresión no es buena.

El programa georreferencia directamente en UTM, WGS84... hasta ahí todo era bueno... el problema comenzó cuando hice pruebas a diferentes alturas (escalas). Probé a 2000 m, 3000 m, 5000, 10.0000.... y los resultados son siempre distintos, al punto de haber más de 150 m de diferencia entre la escala más chica y la más grande. Consulté al autor, y según él esa diferencia es normal. Aquí les entrego la respuesta (que no entiendo bien del todo):

You are right, the stitches with eye altitude, higher than 5000 m are no more correct.
Causing is the projection of Google Earth images. They are ever rectangular.
But if you have a big piece of the earth in your image, it is wrong, because it should be rhombohedra. The earth is a roundly, not plain. The north side of the images is shortening the southern side.

Regards
Petr Bezdecka

Supongo que si la georreferenciación fuese en lat/long no debiese haber esa diferencia en el resultado. Le sugerí al Sr. Bezdecka que incorporase esa opcion y .. cri cri cri...

Por lo pronto, los 57 US$ que cuesta el programa son una pérdida de tiempo y dinero, pues ya no me atrevo a usarlo.

Link de Imágenes de Digital Globe

Este es el enlace a las imagenes de Digital Goble

Cómo convertir un archivo KML en SHAPE

El primer paso es asegurarse que el archivo de Google Earth (GE) está efectivamente en formato kml, y no en kmz, que es la otra extensión estándar de GE. Si estuviese en forma de kmz hay que “guardar como” y elegir el formato “kml” dentro del mismo GE. Este nuevo archivo es entonces el punto de partida.

La extensión de arcview que permite la conversión de kml a shape se encuentra aquí . El autor es Antonio Sobral Almeida. Descomprimiendo el zip encontramos el avx. Un vez que tengamos el archivo kml2shp_1.avx situado en el directorio Ext32, se puede activar la extensión.”KML 2 SHP”.

En la vista aparece ahora un botón con una K mayúscula. Pulsándolo entramos a un diálogo en el que podemos elegir el archivo de entrada (kml) y el de salida. Esta operación es similar para puntos, líneas o polígonos. El resultado es un shape 3d en coordenadas geográficas y proyectado en WGS84, el sistema estándar de GE. Misión cumplida, en principio, pues ya tenemos un archivo shape.

En los archivos shape 3d pueden editarse las tablas, pero no así la parte gráfica. He probado incluso a cargar la extensión 3D y aún así el programa no acepta el movimiento de vértices, estando en modo "Start editing".

Plan B: Para convertir un shape 3d a uno 2d se puede emplear un script que viene dentro de los ejemplos (samples) por defecto del arcview ….búscalo en arcview/samples/xyz2xy.ave.

Una vez que corras el script, se convierte el shape 3d en otro 2d y a partir de ahí se puede editar la capa (Star Editing)….. si fuese necesario.

Cómo convertir un archivo SHAPE en KML

Si queremos convertir un archivo shape en kml el script lo facilita Domenico Ciavarella. Puedes bajarlo de aquí. Descomprimir el zip.

Una vez situado el archivo “shape2KML.avx” en el directorio arcview/ext32 ya puedes activarlo.

Recomiendo ver el video demostrativo de esta extensión.

El siguiente paso consiste en proyectar el shape en geográficas (lat/long) y en WGS 84.

Opción 1. La extensión recién cargada tiene una opción para definir el cambio de proyección. Esto se hace con el botón de "preferences". Hay que definir los códigos de ingreso IC ( Input Coordinate System) y de salida OC (Output coordinate system). Por ejemplo, para pasar de UTM 19 S WGS84 a lat/long, se usan los códigos IC : 32719 y OC : 4326. Estos códigos se obtienen de listados que aparecen pulsando en P? y G?.

Opción 2. También se puede proyectar la capa con el Arcview Projection Utility Wizard, que es una extensión habitual de Arc View para este fin. Si para llegar a ese sistema debes, además, hacer un cambio de datum (por ejemplo, de PSAD56 a WGS84) recomiendo comprobar este paso en particular con otro software distinto a Arc view.

Yo hice una conversión a geográficas muy simple, sin cambio de datum y el resultado final fue el mismo por ambas vías.

==== Definido el proceso de conversión ...

En el botón de “preferences” debemos indicar si queremos que los archivos de salida (kml) se desplieguen inmediatamente sobre GE y si lo harán dentro de una carpeta.

El siguiente paso es pulsar el botón azul de la derecha (Shape to kml) y definir los parámetros solicitados en el cuadro de diálogo.

Estando activa una capa shape, esta se convertirá en archivo kml y el despliegue puede ser inmediato. Especialmente interesante es el hecho de que se pueden mantener un campo del shape (nombre, id, etc.) en el archivo kml. A continuación una muestra del resultado con lineas y polígono en kml.


¿Podrá exportarse la base de datos completa en dbf con el kml?.

Suena a próxima meta.

Proyectar un grid de Lat / Long a UTM

En ocasiones es necesario proyectar un grid en otro sistema de proyección. Por ejemplo, si hemos obtenido un mapa de pendientes a partir de un MDT del SRTM, tendremos un grid proyectado originalmente en lat/long. Lo usual en estos casos es que queramos proyectar el mapa de pendientes en UTM. Este es el procedimiento para realizar esta transformación:

1. Cargar la extensón de Arc View Grid and Theme Projector. Esta extensión es shareware y es facilitada por Jennes Enterprise.

2. Activar el Grid que deseamos reproyectar. Pulsar la opción del menú "Create array of points.... " Con ello se generará una matriz de puntos sobre el grid original que tenemos en lat/long

3. Activar la extensión Arcview Projection Utility (ojo que ésta herramienta no nos sirve si vamos a cambiar de datum; da valores erróneos. Ver post a este respecto).

4. Proyectar la capa de puntos, desde la situaciòn original (WGS84 lat/long) a la nueva proyección (por ejemplo WGS84, UTM, Huso 19 S).

5. Situar la nueva capa de puntos en la misma vista ( no importa que no se vea en la imagen).

6. Pulsar Grid Projector 2 (botón con flecha hacia arriba y número 2). Tenemos a la vista tres columnas. Seleccionar sucesivamente (1) el grid de origen, (2) la capa de puntos de origen y (3) la capa de puntos de salida. Para màs información sobre los métodos de interpolación se sugiere leer los archivos de texto que acompañan a esta extensión. Yo suelo dejar la interpolación bilinear y un orden de transformación de 4. El tamaño de celda del grid de salida puede ser el sugerido por el sistema.

7. Completado este procedimiento tendremos un nuevo grid proyectado en el sistema deseado.

Por último insisto que es recomendable leer las instrucciones de esta extensión para comprender cual es su funcionamiento.

Probando formas de hacer metadata

Probando formas rápidas de hacer metadata encuentro dos opciones. He tratado de exportar los resultados de una a otra pero aún no lo he conseguido. Se aceptan sugerencias !!!.

El primero de ellos es CatMDEdit, del Consorcio TEIDE. Este programa español es empleado por diversas instituciones europeas. Para los que trabajamos en sudamèrica las ventajas principales son la presentación en castellano y que ofrece varios Tesauros con términos disponibles en varios idiomas de la UE , para facilitar el proceso de selección de palabras clave. El programa se puede bajar desde aquí.

CatMDEdit permite tener a la vista una tabla con el conjunto de capas que integran un SIG, y facilita el acceso a cada una de ellas. Al seleccionar una de ellas entramos en su metadata. A la izquierda aparece el listado de campos que permiten el cumplimiento del estándar 19.115. Una clave de colores asociada, distingue los campos obligatorios de aquellos que son opcionales.

Ventaja principal: tendrás a la vista todas las capas de tu sistema (selector de registros).

Inconveniente: Para trabajo en grupo, primero hay que confirmar que los formatos de importación y exportación son compatibles entre distintas plataformas, para evitar tener que repetir la tarea de preparar la metadata para cada sistema. Por ejemplo, aún no consigo que CatMDEdit importe con éxito una metadata preparada con Arc Gis.

Para trasladar una metadata desde CatMDEdit a Arc Gis.

Exportar la metadata en CatMDEdit con la opción CSDGM(FGDC) Formato XML (CSDGM). Esta es la tercera opciòn en la lista que aparece a la izquierda del cuadro de diálogo.

Importar la metadata en ArcGis:

Abrir arc catalog, navegar hacia el directorio donde se encuentra el archivo xml. Seleccionarlo
Y està listo, a primera vista lo reconoce sin mayor problema. Es decir o hace falta hacer ninguna importación adicional.

Para que una metadata se asocie con una capa, hay que renombrar el archivo de manera que coincidan, con una extensión, para la metadata, que tiene esta forma xxx.shp.xml

Para trasladar una metadata desde Arc Gis a. CatMDEdit

He probado varias opciones, pero aùn no consigo que CatMDEdit reconozca el archivo de exportación de Arc Gis o archivos estándar tal y como los produce Arc Catalog.

Problema con Arcview Projection Utility

Empleaba Arcview Projection Utility cada vez que la necesitaba.... sin mayor preocupación… hasta que un día un cliente nos avisa de que estamos proyectando mal. Revisamos… aplicamos la prueba de la blancura (proyectar en una dirección, volver a proyectar en sentido inverso y comparar el resultado)..... y…. efectivamente comprobamos que Arcview se equivoca cuando hay un cambio de datum. La diferencia es de varias decenas de metros, es decir, inaceptable.

El error se comete cuando se proyecta desde WGS84 hacia un sistema particular como PSAD 56 o SAD69. Los dos más usuales en Chile.

Arc view NO comete errores, por lo que he visto, en el sentido inverso: desde un sistema local a WGS84.

El error ocurre siempre, es decir, es un problema de software. Busqué en Internet si alguien ha reportado este problema y... sí, existen referencias a este problema.

Solución: no emplear esta herramienta. Elegir Arc Gis o Datum Pro, u otro sistema. Y hacer la prueba de la blancura antes de aceptar el resultado.

Por lo anterior, es recomendable trabajar con puntos de GPS bajados en WGS84, de manera que seamos nosotros quienes revisemos con cuidado las transformaciones de coordenadas que hacemos, antes de que éstas se apliquen de un modo no controlado.

En resumen : No emplear Arcview Projection Utility para trasladar una capa desde el datum WGS84 a un sistema local y como recomendación, hacer la prueba de la blancura (proyectar en una dirección, volver a proyectar en el sentido inverso y comparar los resultados).

COMO GEORREFERENCIAR UNA IMAGEN GOOGLE EARTH.

GE tiene cada día que pasa mejor resolución en las imágenes que aporta. Así que constituye una base excelente para cartografía.

La propuesta esta vez consiste en trasladar la imagen de GE desde su origen hasta un archivo tiff georreferenciado y disponible en Arc View.

Si tienes la licencia de 400 US $ de GE el resultado será mucho mejor, porque el programa guarda las imágenes con mucha mejor calidad (archivos jpg de 1 a 2 Mb). Si no la tienes, búscate un amigo que sì la tenga o confórmate con una imagen estándar (100 kb).

Plan A: Hemos ido a terreno, con la imagen impresa en papel, hemos identificado un conjunto de puntos reconocibles en la imagen (cruce de dos caminos, un árbol grande, un puente sobre un río, la esquina de una casa…), tenemos un GPS decente y hemos tomado las coordenadas de esos puntos en un sistema específico (UTM Huso 19…WGS 84). De vuelta a gabinete, podemos georreferenciar la imagen contando con los puntos obtenidos. Fácil. Pasos E y 1-8, sustituyendo los cruces de grilla por puntos de coordenadas conocidas.


Plan B: No podemos ir a terreno o bien hay que llevar la imagen georreferenciada de manera preliminar, porque de ese modo podemos superponer otros elementos con los que vamos a trabajar (instalaciones, un proyecto, puntos de monitoreo, etc.). Mi experiencia con esta forma de georreferenciaciòn es irregular, es decir, a veces el resultado es bueno y a veces no lo es tanto (errores de más de 20 m), por lo cual siempre es recomendable poder comprobar el resultado con puntos de terreno, para poder reducir el error a la precisión del GPS (5 a 15 m).

Las extensiones que hacen falta en este caso son Arc view projection utility e Image Analysis.

Los pasos a seguir son:

Trabajando en GE

A. Elige el área que te interesa en GE, navegando y ajustando los límites con las opciones del programa. Cuanto más se pueda acercar uno al suelo, mejor será la calidad de la imagen (la resolución, es decir, el tamaño del píxel, es menor) en opciones, selecciona las coordenadas en lat/long en grados en la opción que lleva comas (primera opción). No con la referencia de grados, minutos y segundos.

B. Elimina elementos sobrantes de la imagen, tales como los botones de navegación, escala, etc.
Activa la grilla (ver / cuadrícula).Cuida que la imagen esté bien orientada, es decir, con el Norte en su sitio. Desactiva el terreno, pues así la imagen quedará mejor y la grilla quedará cuadrada. Un ejemplo de qué ocurre si no se desactiva el terreno es la siguiente imagen, donde la grilla está deforme en sectores de montaña. Insisto, hay que trabajar sin relieve.

C. Guarda la imagen (Archivo / Guardar / Guardar imagen) en un archivo jpg. Puedes llamarle Imagen 1.

D. Sin mover la imagen, elimina ahora la grilla (Ver cuadrícula, desactivar).

E. Vuelve a Guardar la imagen (Archivo / Guardar / Guardar imagen) en un archivo jpg. Puedes llamarle Imagen 2. Esta es la imagen que debieses bajar desde un GE con la licencia de 400 $, para que tenga la mejor calidad. Pero no todo el monte es orégano… con esta licencia el resultado a veces es bueno y a veces está lleno de parches que no tienen ningún sentido y hay que repetir el proceso una y otra vez hasta que a GE se le ocurra hacerlo bien… paciencia ..


Trabajando en Arc view


Nuestra imagen 1 viene con las coordenadas en WGS 84 Lat/long, asì que hay que completar el proceso trasladando esta imagen a un sistema local (por ejemplo, WGS 84 UTM Huso 19 S). Para georreferenciar en otro datum recomiendo leer el siguiente post. Es importante !!!!.

Bueno, manos a la obra:

1. Activar la extensión Imagen Analysis

2. Abrir la imagen 1 (Add Theme) cuidando que en el cuadro de diálogo, abajo a la izquierda “Data Soource Type”.. tengamos seleccionada la opciòn “Image Analysis Data Source. Si nos pide construir pirámides, aceptamos.

3. Activar la imagen 1 y pulsar el botòn de Align Tools (cuadrito negro con flechita hacia el noreste. A partir de aquí sólo se puede usar el botón derecho del ratón para acercarse o alejarse de la imagen y el botòn izquierdo para señalar puntos. Prohibido tocar cualquier otro elemento del menú de arriba, otra capa, otro botón,… nada….sólo el botón derecho del ratón para moverse por la imagen y el izquierdo para señalar puntos. Lo que explicamos a continuación hay que hacerlo 4 o 5 veces

4. Identificar un primer cruce de grilla donde sepamos perfectamente cuales son sus coordenadas geográficas (recomendable apuntarlas… el valor de x y el valor de y, tal y como aparecen en la imagen).

5. Acercarse con el botòn derecho del ratón…zoom in., zoom in, zoom in, sobre ese primer punto. Cuando estamos razonablemente bien situados sobre el cruce de grilla del punto seleccionado… sólo entonces señalamos ese punto con el botón izquierdo .... a continuación pulsamos el botón derecho y seleccionamos Enter “to” Coordinates… introducimos las coordenadas….ojo porque es fácil cometer un error. Hay que cuidar que:

- empleemos el símbolo de coma tal y como lo reconoce el sistema, es decir, con punto (no con coma)
- pongamos negativo a la latitud y longitud si estamos trabajando en Sudamérica
- la x es el valor de longitud (eje horizontal) y la y es el valor de latitud (eje vertical).

6. Una vez que tenemos este primer punto señalamos zoom out.. zoom out… o bien zoom to active Theme (siempre con el botón derecho del ratón) y repetimos el proceso para 3 puntos más, ojala formando un rectángulo en cuyo interior quede la zona de trabajo. No muy alejada de esta.

...

7. Justo cuando introduzcamos las coordenadas del cuarto punto el sistema nos entrega información de nuestro error en la parte baja de la pantalla (RMS…) valores de más de 15 o 20 sugieren volver a repetir el proceso. No es tan raro conseguir 5 m de error o menos si la cosa se hace bien y GE fue "bueno" con nosotros al situar su grilla. Ojo es muy fácil equivocarse en un número de manera que todo lo hecho hasta ese momento deja de servir. La prueba es que la imagen se deforma, se invierte o desaparece por completo (aún cuando señalemos Zoom to active Theme). Hay que eliminar la imagen y volver al punto 2. Sorry.

8. Si llegamos con éxito a tener cuatro puntos georreferenciados, una imagen lógica y un valor de RMS decente (menor a 10 m), procedemos a grabar la imagen en formato tif o img, con un nombre que nos permita identificarla más tarde (Save image as).

Hasta ahora hemos trabajado con la imagen 1, es decir, aquella que tenía la grilla encima. A continuación vamos a georreferenciar la imagen sin grilla, que será nuestro producto final.

9. Identificar 6 puntos de la imagen 1 que sea claros de reconocer en la imagen 2, que estén bien distribuidos por la imagen y que formen al menos algún trapecio entre sì. 4 puntos son los que vamos a emplear en la georreferenciaciòn y otros 2 los vamos a emplear como puntos de control. Tomar nota en un papel… punto 1 cruce de calles tal y tal… punto 2 esquina del potrero verde … punto 3… auto sobre carretera…punto 4…..

10. Con Arcview Projection utility u otro sistema externo a Arc View, hay que cambiar la georreferencia, por ejemplo, desde WGS84 lat/long a WGS84 UTM 19 S. El resultado de este paso es una capa de puntos nueva, georreferenciada en la UTM que necesitamos.

11. Generar una nueva vista. Abrir en esta nueva vista la imagen 2 del modo indicado en el punto 2 y la capa de puntos que acabamos de generar.

12. Definir como capa activa la imagen, pulsar entonces el botón de Align tools y recordar las precauciones que hay que tener a partir de ahora: No tocar otros botones del menú… sólo el botón derecho de la imagen para navegar sobre ella.

13. Identificar sobre la imagen el primer punto….cruce de calles tal y tal… zoom in.. zoom in… zoom in… hasta situarse sobre él de forma precisa. Pulsar entonces el botón izquierdo del ratón… pulsar botón derecho.. zoom out zoom out zoom out .. tantas veces como haga falta hasta tener a la vista el punto 1 correspondiente en el shape… acercarse a él… zoom in zoom in zoom in.. hasta tenerlo bien cerca.. pulsar entonces, sobre el punto, por segunda vez, el botón izquierdo. Le acabamos de decir al programa que el punto 1 de la imagen corresponde al punto 1 de la capa shape.

14. Repetimos el proceso del mismo modo, exactamente igual, para los puntos 2, 3 y 4.. ojo que una vez situados los dos primeros (1 y 2), los nùmeros 3 y 4 debiesen quedar prácticamente encima de donde deben estar, es decir, el desplazamiento va a ser muy poco.

15. A partir del 4º punto el programa arroja un valor de error en la parte inferior de la pantalla. RMS. Fíjate en él. Si es del orden de 10 m o menor “estamos bien”. Si fuera superior deberíamos repetir de nuevo todo el proceso y obtener un nuevo RMS mejor que el anterior. El RMS, en suma, define la calidad del resultado y nos permite tomar decisiones en base a esta (hasta donde confiamos en la georreferenciación).

16. Los puntos 5 y 6 nos iban a servir como puntos de control. Este es el momento de verificar si la imagen coincide bien con estos puntos (un cruce, un árbol..). Revisamos.

17. Si aceptamos el trabajo realizado debemos grabar la imagen 2 con un nuevo nombre (Save image as). Debemos decidir si queremos una extensión tif o img. Generalmente me decanto por la primera pero ciertamente no conozco las ventajas de una u otra.

Eso es todo, en adelante tendremos una imagen GE georreferenciada en nuestro sistema de coordenadas.

Cómo delimitar una cuenca hidrográfica.

La delimitación de una cuenca hidrográfica es posible si se dispone de un modelo digital del terreno, por ejemplo el que provee el proyecto SRTM. Para conseguir este MDT puedes ver más abajo. Los programas necesarios son, en orden de aparición: arcview con spatial analysis, basin 1 y wwf.hydro. Esta última extensión no la encuentro en internet, así que la subí a un archivo de mediafire. Hice un zip con todas ellas y las encuentras aquí.

Por cierto, encontré un Manual que explica el mismo procedimiento descrito aqui en la página de Gabriel Ortiz. La única diferencia es que el autor del Manual parte de curvas de nivel, las trasnforma en TIN y de ahí genera un MDT, mientras que yo inicio el texto desde un MDT bajado del SRTM. Pero la explicación para delimitar las cuencas debe ser la misma.

Empecemos.

Recordemos que el SRTM es un MDT con una resolución de 80 m, por lo que el resultado que vamos a obtener es valido a escalas 1:50.000 (o estirando el chicle) o como máximo 1:25.000. En sectores de pendientes escasas o nulas, SRTM lo pasa mal y hay que emplear la información con mucha precaución pues, en base a lo que he visto, puede dar resultados incorrectos. En estos casos es recomendable verificar los resultados con trabajo de terreno (cuando se pueda) y/o revisar otras fuentes (cartografía base oficial).

El procedimiento que explicamos aquí es válido también con un MDT generado a partir de curvas de nivel locales, compradas a un organismo oficial o construidas mediante restitución fotogramétrica a partir de pares estereoscópicos. En este caso la precisión del resultado es mayor.

Es altamente recomendable leer detenidamente la documentación que acompaña a las extensiones que vamos a emplear. De ese modo conoceremos la teoría que hay detrás de este procedimiento.

El punto de partida es el archivo srtm rellenado con el programa “SRTM_fill”. Consideremos que ese MDT ya está proyectado en el sistema que nos interesa (por ejemplo en UTM, WGS 84, Huso 19 S) (proceso explicado más abajo). A continuación se ejecutan estos pasos:

1.- Eliminación de los “sinks”, es decir, las depresiones que pueda haber dentro del MDT. Se pulsa wwf hydro/remove depressions teniendo activo el MDT. Esto lleva varios minutos. Al terminar el proceso el resultado es un archivo temporal, así que hay que crear un grid permanente con “Convert to Grid”.

2.- Generación del grid de dirección del flujo. Para ello pulsar wwf hydro / Flow directions, teniendo activa la capa producida en el punto anterior. Nuevamente el proceso consume tiempo en base a la extensión del archivo y la capacidad del PC (entre 2 y 10 minutos, por lo general).

3.- Generación de la capa de acumulación de flujo. Pulsar wwf hydro / Flow accumulation, teniendo activa la capa producida en el punto anterior. Nuevamente el proceso consume tiempo (entre 2 y 10 minutos, por lo general).

4. Las capas de flow direction y flow accumulation son temporales, por lo que tenemos que grabarlas con otro nombre mediante Convert to grid.

5.Iniciación de la extensión basin 1.Para ello obviamente tenemos que habr cargado la extensión bas1. En primer lugar hay que tener como capa activa el grid que no contiene “sinks”. Pulsar botón “Initiate”. A continuación el programa pregunta el nombre de la capa de “flow direcction” y el de “flow accumulation”. Indicar uno y otro.

5. Si todo se hizo correctamente se activa el botón de la gotita de agua. Eso significa que estamos casi listos.

6. Recomiendo modificar los colores de la capa de “flow accumulation”. Ello nos permitirá tener una idea precisa sobre los cursos de agua de este grid. El proceso consiste en editar la leyenda y asignar colores claros (o mejor ausencia de color) a todas las celdas con valores bajos (pocas celdas aportan agua a ésta), mientras que se asigna una tonalidad oscura a valores altos (muchas celdas le aportan). Este paso es importante porque la definición del punto a partir del cual se crea la cuenca requiere algo de puntería, estos es, señalar un punto del rìo y no de ladera.

7. Pulsando sobre un punto de la red hídrica se generará automáticamente la cuenca aportante de ese punto. Las limitaciones del sistema son:

- generalmente toma unas pocas celdas situadas aguas abajo del punto que señalamos. Por ello recomiendo tomar un punto ligeramente más elevado que el que realmente nos interesa.
- La delimitación de la cuenca emplea en ocasiones trazos rectos, poco “naturales”. Por ello a efectos de presentación a veces hay que retocar el perímetro suavizándolo. Esto lo podemos hacer editando la capa y teniendo como referencia de fondo las curvas de nivel.
- La base de datos asociada hay que revisarla. En ocasiones, los resultados que arroja no son ciertos.

8. El punto 8 puede repetirse tantas veces como uno quiera. El resultado en cada ocasión serà el perímetro de la cuenca aportante del punto seleccionado.


9. Finalmente, si nos interesa tener la red hídrica del MDT podemos obtenerla mediante wwf hydro stream network. Este proceso se alimenta de la capa de dirección de flujo. Hay que responder al sistema con “yes” / “Number of upstream cells” / e indicar el número de celdas mínimo con el que comience a delimitarse la red. Para un resultado muy general se puede indicar 2000 o 3000. Probando con nùmeros màs bajos la red resulta más densa; asì que conviene repetir el proceso hasta que el resultado sea el adecuado para nuestros propósitos. Cuanto más bajo sea el número de celdas aportantes más tiempo precisa el PC para llegar al resultado.

Las limitaciones de esta operación son:

El trazado de un río depende de cuan bueno sea el MDT de partida. Con una resolución de 80 m como la que trae por defecto el SRTM (fuera de USA), no se pueden esperar milagros. Es decir, los ríos que se dibujan solitos van, ... más o menos por ahí.... Conviene contrastar el resultado con una imagen de satélite como la del Google o Landsat y ajustar el trazado de las líneas a la imagen.


Eso es todo … si aún tienen tiempo… pueden jugar con otras opciones contenidas en estas extensiones.

suerte.

Trabajando con un GPS

El objetivo para hoy es hacer que conversen amistosamente un GPS (GARMIN GPSMAP 60CS) y Arcview 3.x. Como de costumbre, hay muchas maneras de conseguir lo mismo, pero después de un dia de prueba y error, esto es lo mejor que he encontrado.

Para empezar cabe recordar qué son los waypoints, las rutas y los tracks:

- Este GPS almacena 1.000 Waypoints, es decir, localizaciones singulares de puntos.
- Las rutas son trayectos que se definen a partir de waypoints. Se permiten 50 rutas diferentes, con 250 puntos cada una.
- Un recorrido particular se registra en una memoria especial, que se llama tracklog y que puede almacenar hasta 10.000 puntos.
- Sin embargo, cuando el tracklog se graba, éste se limita automaticamente a un máximo de 500 puntos. El sistema permite 20 tracks de este tipo. Ah ! ....No tener esto en cuenta conduce a una desagradable sorpresa !!!
- Por el momento no he encontrado cómo subir mapas completos, tales como shapes de ríos, caminos, etc. que puedan completar el mapa que trae de fábrica el aparato.


Las operaciones que se quieren realizar son:

1.- Bajar los Waypoints (puntos) desde el GPS hacia el Arcview.
2.- Idem los Tracks y el Tracklog.
3.- Subir puntos desde el Arcview hacia el GPS.
4.- Idem los Tracks (que se asemejen a limites, caminos, ríos, etc.)


Los programas que se van a emplear son:

1.- MapSource 6.9, actualización del programa facilitado por el fabricante. Ojo que no funciona si se instala de forma independiente (sin comentarios...).
2.- Vector Transformation 1.0, una utilidad que permite la generalización de polilineas. Se encuentra en una página llena de extensiones interesantes. El archivo avx se baja pulsando con el botón derecho del ratón.
3.- DNR Garmin. Del Departamento de Recursos Naturales (DNR) de Minnesota. La página original dónde estaba contenida este link es muy útil, está aquí. Las extensiones del DNR no tienen desperdicio, hay varias muy útiles. Felicitaciones y Muchas Gracias !!

Bueno, con estos tres programas - y un poco de paciencia - se consiguen bastantes cosas. Vayamos de lo fácil a lo difícil. Animo !!



1.- Bajar los Waypoints (puntos) desde el GPS hacia el Arcview.

A. Esta operación es bastante trivial si tenemos muy pocos puntos, pues se puede bajar la información en Mapsource, grabarla como dxf y abrirla seguidamente con Arcview; teniendo activa la extensión de Cad Reader. Aquí se pierde la base de datos que traigan los puntos.

B. Sin embargo, si tenemos más puntos, es muy probable que necesitemos conservar esa base de datos asociada, so pena de dedicar el tiempo que no tenemos a reconstruirla. En este caso, se puede grabar el archivo desde Mapsource como txt, abrirlo cuidadosamente en Excel, definiéndole una estructura adecuada y grabarla acto seguido como dbf para abrirla como tabla en Arcview. Un recorrido rápido y sencillo si se tienen ciertas precauciones.

C. Una tercera vía para realizar la operación es contar con la extensión del DNR. Para ello se abre el Arcview, se activa la extensión (ojo que hay que situar el archivo avx en el directorio Ext32) y se bajan los datos a la tabla. Desde ahí se puede generar una cobertura de Arcview con solo archivar la tabla en ese formato.

2.- Idem bajar los Tracks y el Tracklog desde el GPS hacia el Arcview.

En general, esta operación es simple, pues los tracks suelen ser escasos y se puede aplicar la operación descrita en 1.A. ya en Arcview se genera un shape con el clásico "Convert to Shapefile" y se trabajan las polilineas según convenga con las herramientas propias de Arcview. Se aconseja copiar aquellos elementos de interés, pegarlos en una nueva sesión de Mapsource y grabar cada track por separado, tal como deciamos, como dxf.

La alternativa es trabajar con la extensión del DNR

Nota: En terreno, pienso que hay que procurar trabajar siempre con tracklog, para evitar tener que grabar esta memoría, más generosa, en los pobres bloques de 500 puntos; de otro modo nuestro trabajo puede perder precisión.

3.- Subir puntos desde el Arcview hacia el GPS.

Esta es la situación inversa, es decir, teniendo un conjunto de puntos definidos en un shape, se pretende subir esta información al GPS (pueden ser puntos de muestreo por visitar, por poner un ejemplo). Con Mapsource no encontré cómo hacerlo !! , (...) , así que apliqué la extensión DNR. El objetivo es trasladar hasta el GPS un identificador de cada punto, el que provenga de la base de datos original del shape. DNR lo permite, pues en el proceso de "subir" los puntos,... el sistema solicita el nombre del campo que actuará como tal identificador ¡...vamos bien..!. Y con la conexión USB el proceso de generar los Waypoints se completa rápidamente.

4.- Idem los Tracks (que se asemejen a limites, caminos, ríos, etc.)

A falta de un buen mapa base cargado en el GPS (grrrr....), ésta es la única opción válida. Lamentablemente, en este modelo sólo se premiten 20 tracks, con un máximo de 500 puntos cada uno. Así que es necesario seleccionar aquellos elementos (caminos, límites prediales, etc.) que realmente necesitamos transportar en el GPS y grabarlos como tracks GRABADOS; Dejando libre memoria para almacenar el tracklog (los 10.000 puntos de ensueño) y disponible algún que otro track que tenga que ser grabado en terreno.

El proceso consiste, entonces, en:

- Elegir una linea en arcview (camino, límite, río, etc.)
- Grabarlo como un shape independiente
- Aplicarle una generalización con el programa Vector Transformation 1.0, de manera que no sobrepase los 500 puntos, pero tampoco pierda más precisión de la necesaria. Una serie de pruebas con 10-20 m de distancia nos permiten llegar a una solución adecuada.
- Trasladar el shape de lineas hasta la tabla de datos del DNR.
- Subir los puntos del track al GPS.
- Grabar el tracklog como TRACK GRABADO. Ponerle un nombre y un color.
- Borrar el contenido del tracklog, dejándolo nuevamente vacio
- Volver a repetir la misma operación con otra línea.

......

PROXIMAMENTE... USANDO EL GPS Y EL ARCVIEW EN TIEMPO REAL...

saludos.

Mapas de pendientes ....a partir de SRTM

Este es el proceso para generar mapas de pendientes (y/o exposiciones) a partir de los archivos del SRTM. El primer paso es abrir el archivo del SRTM (archivo.hgt) mediante la extensión Gridmachine. Para ello hay que:


- activar la extensión “Gridmachine” en el menú de extensiones de Arc view

- seleccionar “Import Machine”, en el menú de Gridmachine

- indicar “SRTM data” como formato de entrada, “Input: Sort of data”

- permitir que la búsqueda se añada a la vista (opción “Load theme into active view”)

- indicar el nombre del archivo de seleccionado, por ejemplo archivo.hgt

- indicar el nombre del directorio de trabajo

Cumplidos estos pasos un tema con aspecto de grid se habrá incorporado a la vista.



Este grid estará georreferenciado en coordenadas geográficas. Puede ser convertido en un grid de pendientes aplicando la función “Derive Slope” cuyo formato de salida es nuevamente un archivo TEMPORAL, que (OJO) debe ser convertido en un grid permanente, nuevamente mediante “Create Grid”.



En la base de datos encontraremos la pendiente en grados. El proceso para generar mapas de exposición (aspect, en inglés) es muy similar al aquí descrito. La diferencia estriba en que debe seleccionarse la opción “Derive Aspect” en vez de aquella de “Derive Slope”


Si necesitamos el mapa de pendientes en porcentaje hay que usar la extensión DEMAT.

Para transformar esta capa en coordenadas geográficas a coordenadas UTM, mira en este blog el post que lo explica.

Suerte .

Creando curvas de nivel ....a partir de SRTM

Algunas veces el problema es que no disponemos de curvas de nivel en el área de trabajo o éstas tienen un alto costo económico. Por si les fuera de utilidad, aquí presento el detalle de un proceso para generar curvas de nivel en cualquier parte del mundo, a partir de Modelos Digitales del Terreno (MDT) bajados del SRTM, para utilizarlas, por ejemplo, en ArcView.

1. Bajar los archivos desde el servidor de SRTM, por ejemplo desde aquí:

El nombre de cada archivo está relacionado con la latitud y longitud. En nuestro caso, interesados en la Isla de Tierra del Fuego, bajamos los siguientes archivos:

S54W069.hgt.zip
S54W070.hgt.zip
....
....

Información relevante sobre SRTM puede encontrarse aqui.

2. Los archivos pueden ser descomprimidos con WinZip.

De esta manera quedan como sigue:

S54W069.hgt
...
...

3. Existe una utilidad que permite corregir los vacios de información que se encuentren en los archivos originales. Esta utilidad se llama SRTMfill y fue creada por 3dNature. El programa se descarga aqui.

Una vez procesado el archivo, este aparece con las siglas"_fill" en el nombre.

4. Otra utilidad, pero que NO es necesaria para nuestro proceso en ArcView, es la que crea una cabecera para los usuarios de ArcInfo. Puede encontrarse en esta dirección.

5. Ya tenemos los archivos del MDT. Estos pueden importarse en ArcView empleando el programa "shareware" Gridmachine, de esta página

En la parte izquierda de la pantalla señalar "Desktopgis" y buscar en el listado el nombre de "Gridmachine 6.63". bajar el programa, que tiene una duración de 90 días sin necesidad de comprarlo (hay un truco para extender su funcionaliad: cerrando el cuadro de diálogo que exige la compra).

6. Instalar el programa en el directorio Arcview / Ext32 (dónde se sitúan las extensiones de Arcview)

7. Abrir el Arcview y activar las extensiones Projection, Gridmachine, Image Analyst y Spatial Analyst

8. En el menú superior de Arcview, aparecerá GridMachine, hay que seleccionar Import Machine.

9. En el área de "Input: Sort of Data", indicar SRTM; más abajo, en browse seleccionar el archivo de entrada.

10. En el último browse (esquina inferior izquierda) se indica el nombre del archivo de salida y su localización.

11. Indicar que el tema aparezca en la vista

12. Si todo ha ido bien, una imagen con aspecto de grid aparecerá en la vista

13. Es el momento de crear las curvas de nivel . Ello se hace

- activando el tema en la vista
- Menu Surface / Create Contour
- Indicar el intervalo entre líneas (50 m, por ejemplo).

Nota: Creo que podría bajarse a 20-25 m, pero no estoy muy seguro de ello.

14. A continuación pueden proyectarse las curvas en UTM. En nuestro caso en el Huso 19S. Para ello abrir la ventana de dialogo de

"Arcview Projection Utility"

*** En este momento las curvas se encuentran en Geographic / WGS 1984 / Units Degree

*** Hay que cambiarlas a Projected / WGS 1984_UTM_Zone 19S / Units Meters

Indicar el destino de los archivos de salida y salir de esta utilidad


15. Es necesario limpiar los archivos producidos de "impurezas", especialmente en las áreas de mar y lagos.



=====

Otras opciones posibles son la generación de mapas de exposición y pendiente, creación de TIN o Grids. Sobre ello estoy aún avanzando.


Este es un ejemplo de cómo se comportan las curvas de nivel creadas comparadas con las adquiridas en el IGM, teniendo de fondo una imagen satelital Mr. Sid . Sector sur de Lago Lynch, en la isla de Tierra del Fuego.

El color rojo se corresponde con la cartografía 3D comprada al IGM, mientras que las curvas en amarillo proceden del SRTM.

Suerte !