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.

martes, 17 de enero de 2017

Lista por comprensión en python para eliminar elementos repetidos

Eliminar los elementos repetidos en una lista python es posible mediante el siguiente algoritmo; donde se declara inicialmente una lista que se llena a medida que se cumple la condición establecida en el if del bucle que recorre la lista original.

lunes, 16 de enero de 2017

Reproyectar puntos con módulo pyproj

La reproyección con pyproj es muy fácil y directa. En este caso hay que usar el contructor de la clase 'Proj' y el método 'transform' para la reproyección. El constructor admite varias sintaxis pero, en el ejemplo a continuación, se va a utilizar aquella que hace referencia al código EPSG (European Petroleum Survey Group). Otro aspecto importante a mencionar es que si el datum de origen viene expresado en unidades diferentes a metros, es imprescindible emplear el modificador 'preserve_units = True' en el método 'Proj' para evitar resultados erróneos.