Recientemente, me vi en la necesidad de instalar un nuevo sistema GNU/Debian estable (Debian 8) porque la actualización de su predecesor transcurrió con un corte de energía eléctrica y los daños fueron irrecuperables. Posteriormente, cuando en consola (como superusuario) me dispuse a instalar fiona y geopandas con easy_install, el cual se encuentra disponible en el paquete python-setuptools de Debian, se produjo un error inesperado. Es de hacer notar que el único paquete SIG que había instalado hasta el momento era QGIS a través de su repositorio respectivo.
Investigando la posible causa encontré que fiona requiere compilar código C++ de GDAL\OGR para instalarse (al igual que geopandas), por lo que resultaba evidente la posibilidad de dependencias insatisfechas, es decir, la necesidad de un compilador g++. Una investigación algo más exhaustiva me llevó a que la instalación del compilador y sus dependencias estaba garantizada con:lunes, 15 de mayo de 2017
miércoles, 1 de febrero de 2017
Intersecciones con método 'overlay' de GeoPandas
El método 'overlay' de GeoPandas permite realizar intersecciones (y otras operaciones de geoproceso) en la misma capa (auto intersección) o en capas diferentes; aunque en la misma capa puede ser difícil de interpretar. Por tanto, para comenzar, se va a realizar la operación con los shapefiles de la imagen siguiente:
domingo, 29 de enero de 2017
Intersecciones en la misma capa con GeoPandas
sábado, 21 de enero de 2017
Cómo crear un geodataframe con GeoPandas empleando rasgos aislados
Una de las principales facetas del módulo python 'GeoPandas' es la facilidad con la que se pueden realizar ciertas tareas. Sin embargo, en el manual a veces no aparecen las cosas que se nos ocurren y toca averiguarlas en otra parte. Por ejemplo, con PyQGIS, si queremos visualizar las capas producto de nuestros geoprocesos se puede recurrir a una memory layer. No obstante, esto puede considerarse como un procedimiento algo laborioso que amerita conocer varios métodos de diferentes clases. Con GeoPandas prácticamente se realiza en una sola línea.
viernes, 20 de enero de 2017
Cómo crear nuevos campos en una geodataframe con GeoPandas y modificar sus valores
GeoPandas es un módulo python que facilita la manipulación de datos espaciales a través de instrucciones muy concisas. En este caso tocará introducir un nuevo campo al vectorial de la imagen siguiente (con un sólo campo en su tabla de atributos) cuando lo carguemos en memoria como geodataframe.
jueves, 19 de enero de 2017
Cómo seleccionar polígonos que contienen al menos un punto con 'sjoin' de GeoPandas
Si se tiene instalado el módulo GeoPandas de Python es muy sencillo realizar este proceso mediante el método 'sjoin' y su parámetro (opción) 'intersects'. Con el ejemplo de la imagen siguiente:
miércoles, 18 de enero de 2017
Union espacial con GeoPandas (método 'sjoin')
Las uniones espaciales en GeoPandas son posibles mediante el método 'sjoin'. Sin embargo, esto requiere además el concurso del módulo Python 'Rtree' (instalado con easy_install) y la librería libspatialindex tomada de los repositorios de mi Debian.
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.