Tengo un shapefile con polígonos. Y tengo un archivo ráster global. Quiero superponer los polígonos del archivo de forma en la cuadrícula de trama y calcular el valor medio de la trama para cada polígono. ¿Cómo puedo hacer esto usando GDAL, escribiendo los resultados en el...
Estoy reproyectando rásteres en python usando GDAL. Necesito proyectar varios tiffs desde las coordenadas geográficas WGS 84 a WGS 1984 Web Mercator (esfera auxiliar), para poder usarlos más tarde en Openlayers junto con OpenStreetMap y tal vez los mapas de Google. Estoy usando Python 2.7.5 y GDAL...
Disculpas si la siguiente pregunta es algo estúpida, pero solo soy MUY nuevo en todo este asunto SIG. Estoy tratando de convertir algunas imágenes geoTiff proyectadas a WGS84 usando gdal en python. He encontrado una publicación que describe el proceso para transformar puntos dentro de GeoTiffs...
Soy realmente nuevo en Python y me gustaría saber si hay un método rápido para obtener valores de celda de un ráster píxel por píxel y las coordenadas (mapa XY del centro de cada píxel) usando Python en ArcGIS 10. Para describir esto más, necesito obtener el mapa X, el mapa Y y el valor de celda...
Dado que geopandas GeoDataFramecontiene una serie de polígonos, me gustaría obtener el área en km² de cada entidad en mi lista. Este es un problema bastante común, y la solución habitual sugerida en el pasado ha sido usar shapelyy pyprojdirectamente (por ejemplo, aquí y aquí ). ¿Hay alguna manera...
Quiero crear mediante programación un archivo GeoJSON usando polígonos de un archivo de forma pero agregando atributos desde mi propia aplicación. Esto se hace fácilmente para un shapefile: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist...
Estoy trabajando con un complemento de ArcMap en C #. Desde el código C #, he ejecutado algunos scripts de Python. Ahora, para ejecutar ese script, tengo una ruta de python codificada. Pero esto no es portátil. Entonces, quiero obtener la ruta del ejecutable de Python desde el código y...
Actualmente estoy trabajando para importar datos climáticos CANGRID (proporcionados como Surfer Grid ascii, archivos ".grd") en ArcGIS. La cuadrícula tiene 95 filas por 125 columnas de tamaño. Los metadatos proporcionan lat / lon de origen (esquina inferior izquierda), tamaño de celda (50 km), así...
Tengo un software (no Arc) que se ejecuta todas las noches en una máquina que no es de producción y que actualiza las características de los archivos de forma desde una base de datos externa. Periódicamente me gustaría eliminar todas las características en el archivo de forma (no el archivo en sí,...
El paquete de geopy de Python presenta dos técnicas de medición de distancia: las fórmulas Great Circle y Vincenty . >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) -
Cerrada . Esta pregunta necesita estar más centrada . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que se centre en un problema solo editando esta publicación . Cerrado hace 9 meses . Excluyendo...
Estoy interesado en aprender a utilizar matrices NumPy para optimizar el geoprocesamiento. Gran parte de mi trabajo implica "big data", donde el geoprocesamiento a menudo lleva días para realizar ciertas tareas. No hace falta decir que estoy muy interesado en optimizar estas rutinas. ArcGIS 10.1...
Tengo una lista de marcos de datos, todos de la misma estructura y todos son polígonos. ¿Cómo puedo agregarlos todos juntos usando
Puntos de datos dados con longitud, latitud y un tercer valor de propiedad de este punto. ¿Cómo puedo agrupar puntos en grupos (subregiones geográficas) en función del valor de la propiedad? Busqué en Google y descubrí que este problema parece llamarse "agrupamiento restringido espacial" o...
¿Alguien sabe de una biblioteca pura de Python (sin dependencias de C o Java) para realizar operaciones de geometría como búfer, unión, intersecciones, etc. Estoy buscando algo que sea compatible con Google
Intentando importar gdal para Python 2.7.3 en Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File...
He estado tratando de crear una función que básicamente haga lo mismo que la función de "disolución" de QGIS. Pensé que sería súper fácil pero al parecer no. Entonces, por lo que he reunido, el uso de fiona con bien proporcionado debería ser la mejor opción aquí. Empecé a jugar con archivos...
¿Es posible dividir un archivo de forma por función en Python? (Lo mejor sería una solución donde pueda guardar temporalmente los objetos vectoriales resultantes en la memoria en lugar de en el disco). La razón: quiero usar la función gdal rasterizeLayer con varios subconjuntos diferentes de...
Estoy usando los datos SIG de mi condado para actualizar las carreteras en OpenStreetMap. Tengo una copia que publicaron el año pasado y una que acaba de publicarse, y me gustaría encontrar todas las LineStrings que hayan cambiado sus atributos o su geometría. Esto me ayudará a asegurarme de que...
Estoy tratando de capturar todos los polígonos no superpuestos que se indican a continuación usando Shapely (dados los polígonos A, B y C). Además, espero hacerlo sin iteraciones, pruebas de intersección, etc. La respuesta aceptada a esta pregunta expresa el método PostGIS, pero parece que 'unión'...