Estoy trabajando con los archivos TRMM 3B31 sobre precipitación proporcionados por la NASA. Cuando los cargo en QGIS seleccionando la proyección WGS84 junto con otros archivos de formas, los mapas se muestran verticalmente, es decir, QGIS interpreta como latitud lo que en cambio es longitud y...
Quiero crear dos polígonos. Uno es de las extensiones rectangulares de mi trama. Sé que puedo obtener la extensión de la trama usando r <- raster("band5.tif") e <- extent(r) plot(e)me da esto: ¿Cómo puedo crear UN polígono del límite del ráster como se muestra a continuación? Así es...
Tengo un dataset ráster en caché / formato mixto en ArcGIS. Necesito exportar esto a un tiff georreferenciado o cualquier otro formato de imagen ráster para usarlo como mapa base en algún otro software SIG basado en escritorio como QGIS. Hasta ahora, solo he encontrado una herramienta en ArcGIS...
La calculadora de campo tiene una gran función en "Registros" llamada "$ rownum" que devuelve (asigna) el número de la fila actual. ¿Hay alguna manera de que funcione junto con la clasificación? Digamos que quiero ordenar la tabla por "longitud" de columna y obtener el orden (secuencia) de las...
He intentado varios ejemplos de código utilizando bibliotecas como shapefile, fiona y ogr para intentar verificar si un punto (x, y) se encuentra dentro de los límites de un multipolígono creado con ArcMap (y, por lo tanto, en formato shapefile). Sin embargo, ninguno de los ejemplos funciona bien...
Tengo una tabla con 8 columnas y ~ 16.7 millones de registros. Necesito ejecutar un conjunto de ecuaciones if-else en las columnas. He escrito un script usando el módulo UpdateCursor, pero después de algunos millones de registros se queda sin memoria. Me preguntaba si hay una mejor manera de...
Tengo una capa vectorial con millones de polígonos que hacen una cobertura continua. Necesito clasificarlos según su forma. Ya estoy usando varios índices de formas de la ecología del paisaje como la compacidad ( 4piA / P ^ 2 ), el ancho medio ( 2A / P ), el número de forma ( P / sqrt (A) ),...
Parece que las opciones de Procesamiento en QGIS 3.0 ya no incluyen Scripts R (no están listados como disponibles junto con SAGA, GRASS y GDAL). ¿Es este un cambio en esta versión de QGIS o he tenido un problema con la
usando Gdal / OGR, ¿hay alguna forma de convertir una capa de línea shp, kml o PostGIS en una red de nodos y segmentos para que pueda usarse en paquetes como
¿Hay una buena manera de filtrar conjuntos de datos con aproximadamente 50,000 características en Geoserver? Tendría que seleccionar entre 100 y 200 características basadas en identificadores. Estos identificadores son diferentes de los FeatureIds. No se pudo encontrar un operador "IN" similar a...
Estoy buscando una función similar a GRASS r.resample. Tengo una carpeta llena de archivos TIF y quiero reducir su resolución ("reducirlos") para reducir el tamaño del
Tengo dos archivos LiDAR (.las), uno es original, digamos con X puntos. Y el otro es una copia del primer archivo .las pero con puntos Y, donde Y es menor que X. Ahora, quiero comparar los Modelos Digitales de Elevación (DEM) de estos dos archivos .las y visualizar cuán diferentes son. Quiero...
Me pregunto si hay una fuente "unificada" para las convenciones de color del mapa, que los creadores de mapas SIG deben seguir. EDITAR Para aclarar un poco mi pregunta. Estoy más interesado en el color con respecto al objeto que representa. Por ejemplo, cuando mira un mapa y ve "azul", la...
Estoy tratando de calcular el valor promedio de lluvia a partir de varios puntos usando la herramienta Distancia inversa ponderada en ArcGIS 9.3. Mi problema es que: cada punto tiene su propia serie de tiempo, por lo tanto, el proceso de interpolación debería poder llevarse a cabo durante todos...
Tengo un shapefile de subdivisiones de país y me gustaría extraer una serie de límites de longitud y longitud para cada división ... ¿es posible
¿Cómo puedo hacer selecciones espaciales o filtros en R? Puedo hacer esto en PostGIS usando (por ejemplo): SELECT * FROM mypoint, mypoly WHERE ST_Within(mypoint.geom, mypoly.geom); y de manera similar es pan comido en la mayoría de los demás software SIG. Así que esperaría que algo en R se vea...
Estoy tratando de configurar un Python virtualenv aislado para trabajar en proyectos SIG en mi Mac OS X. Parece que necesitaré instalar GEOS, PROJ, GDAL / OGR desde el sitio de Kyngchaos aquí http://www.kyngchaos.com/ software / frameworks y no fácil_instalarlo en mi virtualenv. ¿Es eso...
Estoy tratando de encontrar una manera de subdividir un mapa basado en vectores en bloques más pequeños. Puedo colocar el mapa en mosaico en función de diferentes estrategias, como q-tree. Sin embargo, es inevitable cortar algunas polilíneas y polígonos que cruzan los límites de los bloques. Al...
Estoy tratando con algunos conjuntos de datos grandes y me gustaría minimizar la cantidad de datos enviados a través de la red. Estoy usando la estrategia BBox de OpenLayers para transferir solo los datos que se necesitan, pero me pregunto si puedo simplificar la salida para poder mostrar niveles...
Tengo una nube de puntos en mi base de datos (espacial de SQL Server 2008). Eso es alrededor de 6 millones de registros. Hay 3 columnas: id, valor, geom. ¿Cuál es la forma óptima de obtener el 'valor' en la entrada lat largo? Soy nuevo en las consultas espaciales en SQL Server 2008. ¿Alguien puede...