Sistemas de Información Geográfica

¿Es posible crear una capa virtual a través de un script de Python? Por ejemplo, tengo una capa "road" y me gustaría realizar sql "SELECT * FROM road WHERE type = 'Expressway'" ¿Será esto posible? ¿Hay algún ejemplo al que me pueda referir?

10 pyqgis  virtual-layer 

Me gustaría anular el estilo predeterminado de la información sobre herramientas de Leaflet 1.0.0, especialmente la burbuja / marco. No veo ninguna opción o método para esto. ¿Cómo me engancho al CSS? Quiero cambiar las capas individualmente, así que necesito seleccionar la información sobre herramientas de cada capa individualmente con …

10 leaflet  css 

Para obtener las coordenadas exteriores, necesito convertir un MultiPolygon bien formado en un Polígono. Lo hago así: if poly.geometry.type == 'Polygon': x, y = poly.geometry.exterior.xy elif poly.geometry.type == 'MultiPolygon': allparts = [p.buffer(0) for p in poly.geometry] poly.geometry = shapely.ops.cascaded_union(allparts) x, y = poly.geometry.exterior.xy # here happens the error Esto tiene …

10 polygon  shapely  multipart  singlepart 

¿Hay alguna herramienta SWIPE en QGIS similar a la que tiene en ArcGIS? Lo que estoy buscando es una herramienta para intercambiar rápidamente entre dos capas. Estoy comparando el mosaico del sonar de barrido lateral con la batimetría para problemas de posicionamiento.

10 qgis  tools  swipe 

Estoy usando el folleto para representar un mapa. Creé un mapa con marcadores y no sé cómo implementar el escucha de eventos 'onClick' en cada marcador. Mi código var stops = JSON.parse(json); var map = new L.Map('map', { zoom: 12, minZoom: 12, center: L.latLng(41.11714, 16.87187) }); map.addLayer(L.tileLayer('https://api.tiles.mapbox.com/v4/{id}/{z}/{x}/{y}.png?access_token={accessToken}', { attribution: 'Map …

10 javascript  leaflet 

Leí una publicación sobre mapas interactivos con R usando el leafletpaquete. En este artículo, el autor crea un mapa de calor como este: X=cbind(lng,lat) kde2d <- bkde2D(X, bandwidth=c(bw.ucv(X[,1]),bw.ucv(X[,2]))) x=kde2d$x1 y=kde2d$x2 z=kde2d$fhat CL=contourLines(x , y , z) m = leaflet() %>% addTiles() m %>% addPolygons(CL[[5]]$x,CL[[5]]$y,fillColor = "red", stroke = FALSE) No …

10 r  leaflet  heat-map 

Cuando se usa vrt para mosaico sobre la marcha grandes colecciones de imágenes, ¿es mejor agregar vistas generales (pirámides) a los rásteres individuales o al .vrt final, o ambos? (Donde "mejor" significa encontrar el punto de equilibrio entre no usar más almacenamiento del necesario al mismo tiempo que proporciona una …

10 vrt  image-pyramids  raster-catalog 

Necesito rasterizar un archivo de formas realmente simple un poco como este http://tinyurl.com/odfbanu . ¿Qué es solo un archivo de forma que contiene condados en los EE. UU. He visto esta respuesta anterior: GDAL RasterizeLayer no quema todos los polígonos en Raster? pero me preguntaba si hay una manera de …

10 python-2.7  fiona  geopandas 

Estoy escribiendo un script de Python usando la API QGIS. Tengo que obtener resultados de la tabla de base de datos postgres. La tabla tiene un campo de geometría y otros campos que no son de geometría. Quiero realizar consultas SQL en campos no geométricos y obtener resultados en mi …

10 qgis  postgresql  pyqgis 

Normalmente construyo mapas topográficos estáticos, que en su mayoría se imprimen. A menudo uso OSM y rasterfiles (principalmente DEM diferentes). Las áreas que estoy mapeando rara vez exceden los 20,000 km². Estoy usando exclusivamente software de código abierto como QGIS, SAGA y GRASS. Conozco los problemas con los archivos de …

10 postgis  shapefile  spatialite 

Mis disculpas si esta es una pregunta repetida, busqué aquí y no pude encontrar una respuesta. Tengo datos JSON que exporté de Openpaths.cc que contienen valores de lat y lon. Sin embargo, los datos no están en formato GeoJSON, por lo que OGR2OGR no puede leerlos. ¿Alguien podría señalarme en …

10 leaflet  geojson  ogr2ogr  json 

Soy nuevo en R. ¿Cómo puedo reemplazar los NA con 0 para mis datos ráster que tienen información espacial?

10 raster  r  missing-data  replace 

Estoy editando un archivo de forma en QGIS y quiero dibujar en un rectángulo con dimensiones específicas (por ejemplo, digamos 101 pies por 305 pies). ¿Hay alguna manera en QGIS de hacer que dibujar un rectángulo sin tener que usar el mouse para colocar manualmente cada uno de los puntos …

10 qgis 

Estoy teniendo muchos problemas con GDAL. Además de la falta de documentación, parece haber poco soporte en Python. De todos modos, esto es solo un script de reproyección por lotes que va así: source_file = gdal.Open(filepath.encode("ascii")) source_wkt = source_file.GetProjectionRef() source_srs = osr.SpatialReference() source_srs.ImportFromWkt(source_wkt) reproj_file = gdal.AutoCreateWarpedVRT(source_file, source_wkt, dest_wkt) gdal.ReprojectImage(source_file, reproj_file, …

10 python  gdal  coordinate-system 

Muchos artículos científicos informan la ubicación del sitio en unidades de grado / min / seg, pero para el análisis con software que no es SIG, es más fácil usar grados decimales. ¿Hay alguna razón para que los grados / min / seg sean la convención más utilizada? ¿Cuáles son …

10 coordinate-system