Lo que estoy tratando de hacer: bucle a través de un shapefile de puntos y seleccionar cada punto que cae en un polígono. El siguiente código está inspirado en un ejemplo de consulta espacial que encontré en un libro: mitte_path = r"D:\PythonTesting\SelectByLocation\mitte.shp" punkte_path =...
Encontré este artículo que describe la "forma de QGIS GUI" para calcular las longitudes de línea. http://www.qgistutorials.com/en/docs/calculating_line_lengths.html ¿Cómo puedo calcular la longitud de una línea (característica) en una capa vectorial con código Python? Como ejemplo, quiero saber...
Estoy tratando de entender los códigos EPSG. He estado en el sitio spatialreference.org, y en su enlace "acerca de", dice: "Actualmente hay tres catálogos distintos de sistemas de referencia espacial en este sitio web: EPSG, ESRI y definido por el usuario (llamado sr-org)". ¿Cuál es la diferencia...
Usando arcpy, ¿cómo verifica si dos clases de entidad tienen la misma referencia espacial? Solo verificar si los dos son iguales no funciona: >>> import arcpy >>> fc1 = r"C:\Users\e1b8\Desktop\E1B8\GIS_Stackexchange\data.gdb\test" >>> sr1 = arcpy.Describe...
Me gustaría usar algunas proyecciones en QGIS que no puedo ver soportadas directamente. Un ejemplo sería algo así como esta proyección de cubo conforme. http://www.csiss.org/map-projections/Polyhedral_Globes/Cube_Gnomonic.pdf ¿ Alguna pista sobre cómo encontrar detalles y utilizar esta o...
Quiero importar GeoJSON disponible en el siguiente volcado en PostGIS pero no puedo importarlo. Ya importé los datos geousa disponibles en este basurero público a MongoDB sin problemas. Avíseme, ¿hay alguna manera de importar estos datos directamente o a través de MongoDB? Como estos datos son...
Estoy usando el método GDAL ReadAsArray para trabajar con datos ráster usando numpy (específicamente reclasificación). Como mis rásteres son grandes, proceso las matrices en bloques, iterando en cada bloque y procesando con un método similar al ejemplo de GeoExamples . Ahora estoy buscando la...
Los barrios reina y torre son dos formas comunes de calcular estadísticas para una celda focal. También se les conoce como vecindarios Moore y (von) Neumann. Creo que también existe una "versión obispo", es decir, desde la esquina inferior izquierda a la superior derecha. Pero, ¿por qué estos dos...
Estoy tratando de mapear un CSV con más de un millón de puntos lat / lon en los Estados Unidos sin metadatos. Nunca he intentado asignar tantos datos y quiero saber qué pasos debo seguir. Sé que son demasiados datos para QGIS, así que intenté dividirlos en CSV de 100,000 líneas, pero luego obtuve...
Necesito encontrar un algoritmo o método que pueda detectar latitude longitude puntos atípicos en una trayectoria durante el procesamiento posterior , que luego se puede arreglar (volver a la trayectoria de la trayectoria en función de sus vecinos). Como ejemplo del tipo de puntos atípicos que me...
Estoy desarrollando un complemento QGIS y me gustaría usar paquetes de flecha o solicitudes , ya que estoy acostumbrado a manipular fechas y API REST. QGIS no viene con estos paquetes (con menos flecha), así que cuando intento acceder a mi complemento en QGIS tengo: "ImportError: ningún módulo...
¿Cómo cambiar la expresión para definir el número de decimales? layer.dataProvider().addAttributes([QgsField("fieldname", QVariant.Double)]) Como leí en la documentación de la API de QGIS: Parameters name Field name type Field variant type, currently supported: String / Int / Double typeName...
Intento desarrollar un complemento para agregar capas postgis a qgis desde un complemento como parte del aprendizaje de Python. El código para agregar una capa postgis se toma del libro de cocina y funciona bien si lo ejecuto desde la consola de Python dentro de qgis (si se presiona el botón OK en...
Quiero usar este archivo json, todavía no es un archivo GeoJSON, pero noté que contiene varias características y no una que me parece confusa. Quería preguntarle si conoce alguna herramienta en la que pueda fusionar todas las características / colecciones de características en un archivo GeoJSON...
Deseo recortar un conjunto de polilíneas (líneas negras en la imagen a continuación) en el límite exterior de un polígono. Cualquier vacío dentro del polígono debe ser ignorado. Mi salida ideal son las líneas amarillas discontinuas. Las líneas iniciales pueden o no ser rectas. La imagen es un...
Tengo una copia maestra de algunos datos en los que estoy trabajando en un servidor remoto. A veces hago una copia en mi disco local para acelerar los tiempos de procesamiento, pero las fuentes de datos en mis MXD todavía apuntan al servidor remoto. Me pregunto si hay un truco rápido y fácil para...
Tengo un gran proyecto que consta de alrededor de 10x archivos de proyecto QGIS. Cada proyecto contiene un compositor de impresión con un atlas configurado para generar hasta 10 mapas. ¿Cuál es la forma más fácil de automatizar la apertura de cada proyecto y crear un PDF utilizando el atlas...
Estoy ejecutando un servidor Windows Server 2012. PostgreSQL 9.3.5, 64 bits POSTGIS 2.1.3 QGIS Server 2.6.1-2 QGIS Desktop 2.8.3 y 2.12 Estoy usando una tableta Microsoft Surface Pro 4 con QGIS Desktop 2.12 para editar algunas capas almacenadas en la base de datos Postgresql. Las capas en la...
Estoy desarrollando un complemento PyQGIS usando python. Me enfrento a un problema durante la implementación de subprocesos en mi complemento. He escrito un script de Python que contiene una clase que envuelve todos los métodos de geoprocesamiento en uso. El problema que estoy enfrentando es...
Esta es una pregunta de geometría de SQL Server aparentemente simple que pensé que tendría una solución lista para usar, pero no tengo suerte para encontrarla. Mi intención es seleccionar todos los registros dentro de una tabla que tengan polígonos anidados (contenidos) dentro de un polígono más...