domingo, 15 de enero de 2017

Reproyectar un archivo vectorial mediante GeoPandas

Como se ha señalado con anterioridad, GeoPandas es un módulo python que extiende los 'datatypes' usados por 'pandas' para permitir operaciones sobre geometrías; las cuales son realizadas con el módulo 'shapely'. Además, GeoPandas depende adicionalmente de 'fiona' para el acceso a los archivos y 'descartes' y 'matplotlib' para graficar.



Otra operación que se puede realizar de manera expedita con GeoPandas es la reproyección. El código siguiente permite este proceso en apenas muy pocas líneas:

1
2
3
4
5
6
7
import geopandas as gpd

tmp = gpd.GeoDataFrame.from_file('/home/zeito/pyqgis_data/ward2012/ward.shp')

tmpWGS84 = tmp.to_crs({'proj':'longlat', 'ellps':'WGS84', 'datum':'WGS84'})

tmpWGS84.to_file('/home/zeito/pyqgis_data/ward2012/wardWGS84_2.shp')

El shapefile usado en el ejemplo anterior se puede bajar de aquí:

ward2012.zip

El resultado de la reproyección se encuentra en la imagen siguiente:



No hay comentarios: