Cómo reproyectar un shapefile en QGIS 1.8?

Ayer, en un comentario de un post del Blog, me señalan la imposibilidad de abrir un kml en ArcGis (no me especifican la versión) y me agradecen una solución alternativa en QGIS que incluyera también la reproyección de WGS 84 long/lat a WGS 84 UTM zona 13. Aunque tenía más de un año que no usaba ArcGis (solía emplear las versiones 9.2 ó 9.3), como recientemente me instalé la versión 10.1 (también sin uso), me entró la curiosidad por averiguar que tan difícil podría ser abrir un kml con este software; algo sencillo de hacer en QGIS.

Point shapefile, con librería ogr de Python, a partir de un archivo de texto plano

En un artículo anterior, se consideró un script de Python que crea un shapefile de puntos y los añade de manera interactiva a través de la cónsola. Sin embargo, si se dispone, por ejemplo, de gran cantidad de puntos es preferible que éstos sean añadidos a partir de un archivo de texto plano. El considerado (datos1) para probar el script fue el siguiente:

Script de Python para generar aleatoriamente palabras “colgadero” para “Supermemoria”

Últimamente, le he dedicado bastante atención a los procedimientos necesarios para recordar grandes volúmenes de información. Para ello me he basado en los sugeridos por Harry Lorayne en su libro “COMO ADQUIRIR UNA SUPERMEMORIA”. En éste se plantea el uso de palabras “colgadero”, que se asocian a un número y una imagen, que al aprenderlas de memoria podrán ser usadas para asociarlas en cadena para recordar otros objetos o números de cierta extensión.

Dividir shape circular en 4 cuadrantes

En el foro de Gabriel Ortiz una usaria tenía la necesidad de una extensión para ArcGIS o ArcView que pudiese, de manera automática, dividir un polígono circular en cuatro sectores iguales NW, NE, SW y SE. Parece obvio pensar que usando grass-script de Python la tarea es fácilmente realizable. Sin embargo, me entró curiosidad por saber que tan tedioso podía ser usando sólo QGIS.

Los números de 2012 en WordPress.com

Los duendes de las estadísticas de WordPress.com prepararon un informe sobre el año 2012 de este blog.

Aquí hay un extracto:

Unos 55,000 turistas visitan a Liechtenstein cada año. Este blog ha sido visto cerca de 250.000 veces en 2012. Si fuera Liechtenstein, se necesitarían alrededor de 5 años para que todos lo vean. Tu blog tuvo más visitas que un pequeño país en Europa!

Haz click para ver el reporte completo.

Conversión de coordenadas UTM a geográficas (datum WGS84) usando los elipsoides Hayford o WGS84 con módulo pyproj

En un artículo precedente se elaboró una aplicación escrita en C/C++ destinada a convertir coordenadas UTM en geográficas (datum WGS84) usando los elipsoides Hayford o WGS 84 y tomando como base para la codificación las ecuaciones de Coticchia-Surace. Codificar este mismo tipo de cálculos en Python es mucho más sencillo porque se cuenta con el módulo pyproj; tal como se presenta a continuación:

Determinación de los valores promedios de un ráster usando el módulo GDAL de Python

En un artículo precedente se calculó la media aritmética por columnas en un ráster, en formato ascii, usando un script de python. En el presente, se va a usar el código de la página 51 de Python Geospatial Development, ligeramente modificado para adaptarlo a mis condiciones locales, para determinar el promedio total del mismo ráster. Este código usa, primero, el módulo gdal de Python para obtener y leer la banda simple del ráster a través de una línea de barrido y luego el módulo struct para leer los valores individuales fuera de la línea referida. La ventaja, con relación al primer código, es que permite usar directamente el ráster en su formato *.tif original; sin necesidad de convertirlo previamente a ascii. Además, tiene métodos que permiten extraer el número de filas y columnas del ráster y su tiempo de ejecución es considerablemente menor.