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:
Publicar un comentario