Mapaverde, cómo lo hice

Hace unos días la gente de Buenos Aires Data (@GCBAdata) publicó el dataset del censo de arbolado de la Ciudad de Buenos Aires del año 2011.

Para crear el mapa de densidad de arbolado tuve que hacer lo siguiente:

  1. Usando Quantum Gis y el plugin MMQGIS cree un grid de polígonos hexagonales de 100 metros en la diagonal sobre el mapa de comunas de la Ciudad (http://data.buenosaires.gob.ar/dataset/comunas).
    detallegrid
  2. Subí los árboles a QGis y usando las herramientas incluidas intente que cuente la cantidad de árboles que había en cada polígono, pero por la cantidad de puntos (árboles) y polígonos (hexágonos) que tenia se hacía eterno.
  3. Después de instalar PostGIS cree una tabla para subir el csv de los árboles y usando shp2sql convertí el grid de hexágonos en otra tabla.
  4. Con una consulta a la base de datos pude ver cuantos árboles tenia cada hexágono y agregar ese dato a la tabla de hexágonos
  5. Usando TileMill cree un mapa con los shapefiles de las comunas y barrios de la Ciudad (http://data.buenosaires.gob.ar/dataset/barrios), a ese mapa le agregué una capa PostGIS que consultaba la tabla de hexágonos
  6. Usando el css de TileMill le di el estilo que quería al mapa (colores, grosor de las lineas, etc.) y compartí la primer imagen del mapa por Twitter.tweet
  7. Para crear la versión interactiva del mapa agregué la leyenda al mapa y lo exporté desde TileMill en formato MBTiles
  8. Usando MBUtil extraje los tiles y los subí al servidor web donde cree una página html, subí el script Wax y usando la versión Leaflet entre los mapas disponibles lo configuré para que encuentre los tiles.

Un tema a tener en cuenta es que todos los datos georeferenciados que provee el GCBA están en el sistema de coordenadas y proyección Gauss Krugger Buenos Aires (GKBA) y ninguno de los programas que use tiene esta opción preconfigurada, para poder transformarlo a WGS84 (la proyección mas conocida) pueden usar el API que ofrece el GCBA y que explican como usar en este post o agregar los datos de la proyección a los distintos programas, en general esto se hace pasando los parámetros que usa proj4, que para el caso de GKBA a WGS84 son:
+proj=tmerc +lat_0=-34.629269 +lon_0=-58.4633 +k=0.9999980000000001 +x_0=100000 +y_0=100000 +ellps=intl +units=m +no_defs

Si todavía no lo vieron, el resultado final lo pueden encontrar en http://unahormiga.com/mapaverde/

6 Comments

Filed under Varios

6 Responses to Mapaverde, cómo lo hice

  1. Te felicito por el trabajo. Te hago una pregunta, el dataset que usaste es el mismo que está actualmente en la página? O puede haber pasado que haya cambiado. Este es el subido actualmente (http://data.buenosaires.gob.ar/dataset/censo-arbolado).
    Te lo comento porque por que observando esta información en ArcMap veo una menor densidad en algunos sectores como por ejemplo en la comuna 10 y 9. Desde ya muchas gracias y nuevamente felicitaciones por el trabajo.

  2. Acá te dejo el enlace a la imagen que estoy obteniendo con este dataset. https://dl.dropboxusercontent.com/u/6811962/CABA.jpg
    Los puntos grandes son simplemente algunas especies que me interesaban remarcar.
    ¿Será el mismo dataset?

    ¡Desde ya muchas gracias!

  3. Martín, el dataset cambió en algún momento. El que yo todavía tengo en disco tiene por encabezado las siguientes campos:
    nombre_barrio,codigo_manzana,alt_ini,alt_fin,calle,coord_x,coord_y,id_especie,nombre_familia,
    nombre_genero,nombre_cientifico,nombre_comun,tipo_follaje,origen

    el que tienen actualmente tiene los campos:
    nombre_barrio;codigo_manzana;alt_ini;alt_fin;calle;coord_x;coord_y;id_especie;nombre_familia;
    nombre_genero;nombre_cientifico;nombre_comun;tipo_follaje;origen;longitude;latitude;geometry

    en cuanto a los arboles que figuran en cada uno por lo visto están en otro orden así que no los puedo comparar rápidamente sin ponerlos dos datasets en un mismo mapa y ver si coinciden. Lo que si te puedo decir es que los dos archivos tienen la misma cantidad de líneas y comparando algunas especies al azar parecen tener la misma información.

    Lo que me llama la atención de tu imagen es que tiene muy marcada una zona de alta densidad y otras con muy pocos arboles. No tengo experiencia con ArcMap, pero puede ser que no te este tomando todos los datos, tal vez algún limite de memoria? Tenes forma de ver cuantos puntos esta mostrando el mapa con ArcMap? Se supone que son 359000 aprox.

    Si te sirve te puedo mandar la copia vieja del dataset.

  4. Juan, muchas gracias por la respuesta. Voy a chequear estas cosas. Es probable que hay un límite de memoria. Voy a intentar haciendo un dataset por comuna y te cuento. ¡Muchas gracias!

  5. Encontré finalmente el problema. Luego de probar muchas cosas, accidentalmente hice zoom out y vi lo siguiente.
    https://dl.dropboxusercontent.com/u/6811962/arbolado.jpg

    Lo que se ve abajo a la izquierda es la capital, el resto son árboles perdidos.
    Parece que hay varios con otra georeferencia. Hasta ahora estoy tomando los campos coord_x, coord_y.
    Tal vez, con otro sistema de coordenadas se puedan usar estos campos nuevos que notaste que agregaron (longitude;latitude). Por ahora no tuve éxito. Te agradecería mucho si me podés pasar la versión que tenés del archivo.

  6. Pude hacerlo andar, usé los nuevos campos (longitude;latitude) y cambié el sistema de coordenadas para esa capa por GCS_WGS_1984 (es el que usan las del IGM, http://www.ign.gob.ar/sig).
    Gracias nuevamente.

Deja un comentario

Tu dirección de correo electrónico no será publicada. Los campos necesarios están marcados *