Sistemas de Informação Geográfica

Estou manipulando e processando rasters globais com resolução de 30m. O tamanho total da varredura é geralmente [1.440.000 560.000]. Como tenho acesso a um supercomputador, escrevi um código que me permite dividir rasters globais em blocos gerenciáveis, executar alguns cálculos em paralelo e gravá-los no disco rapidamente. Eu bati em …

10 qgis  gdal  performance  display  gdaladdo 

Sou novo no fascinante mundo do SIG e tenho pesquisado o assunto nos últimos dias. Meu objetivo é fornecer mapas em todo o mundo com dados OSM em meus próprios servidores. Minha pilha é: Postgres + Postgis - Para armazenar, transformar e processar dados GIS. Mapnik - Biblioteca para renderização …

10 coordinate-system  openstreetmap  mapnik 

Estou tentando escrever um script de processamento usando o r.neighborsalgoritmo de grama para rasters. Ao usar o Gui para esse algoritmo, você pode deixar o campo de extensão da camada GRASS em branco para usar a extensão mínima. Pelo contrário, quando usado no console python, não consigo encontrar como inserir …

10 qgis  pyqgis  qgis-processing  extents  qgis-python-console 

Eu vi que o gdalwarp tem a -multiopção de usar vários threads. Mas não vejo essa opção para gdal_translate. O gdal_translate suporta multiencadeamentos?

10 gdal-translate 

Estou usando o QGIS e estou procurando um script ou um plug-in que possa criar um grande número de linhas perpendiculares de uma camada de ponto a linhas em uma camada de linha separada. Até agora, tentei usar a função Distância do Hub no MMQGIS (transformar as linhas em pontos …

10 qgis  pyqgis  line  perpendicular 

Quero saber se a versão mais recente do Openlayers (v4.0.1) suporta o Google Maps como uma camada de mosaico. Não consigo encontrar nenhum documento sobre isso. Se o Openlayers não suportar o Google Maps, alguém poderia me dizer se existe alguma maneira de fazer isso?

10 javascript  openlayers  google-maps  google-maps-api 

É possível criar uma camada virtual através de um script python. Por exemplo, tenho uma camada "road" e gostaria de executar o sql "SELECT * FROM road WHERE type = 'Expressway'" Isso será possível? Existe algum exemplo ao qual eu possa me referir?

10 pyqgis  virtual-layer 

Gostaria de substituir o estilo padrão das dicas de ferramentas do Leaflet 1.0.0, especialmente a bolha / quadro. Não vejo opções ou métodos para isso. Como eu me conecto ao CSS? Quero alterar as camadas individualmente, por isso preciso selecionar as dicas de cada camada individualmente com CSS.

10 leaflet  css 

Para obter as coordenadas externas, preciso converter um MultiPolygon bem torneado em um polígono. Eu faço assim: 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 Isso é …

10 polygon  shapely  multipart  singlepart 

Existe alguma ferramenta SWIPE no QGIS semelhante à que você possui no ArcGIS. O que estou procurando é uma ferramenta para trocar rapidamente entre duas camadas. Estou comparando o mosaico do sonar sidescan com a batimetria para problemas de posicionamento.

10 qgis  tools  swipe 

Estou usando o Leaflet para renderizar um mapa. Criei um mapa com marcadores e não sei como implementar o Event Listener 'onClick' em cada Marcador. Meu 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 data', …

10 javascript  leaflet 

Eu li um post sobre mapas interativos com R usando o leafletpacote. Neste artigo, o autor cria um 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) Eu …

10 r  leaflet  heat-map 

Ao usar o vrt para criar mosaicos dinâmicos de grandes coleções de imagens, é melhor adicionar visões gerais (pirâmides) às rasters individuais ou ao final .vrt, ou a ambos? (Onde "melhor" significa encontrar o ponto de equilíbrio entre não usar mais armazenamento do que o necessário e, ao mesmo tempo, …

10 vrt  image-pyramids  raster-catalog 

Eu preciso rasterizar um shapefile realmente simples um pouco como este http://tinyurl.com/odfbanu . O GDAL RasterizeLayer não queima todos os polígonos em Raster? mas eu queria saber se existe uma maneira de fazer isso usando Geopandas ou fiona e talvez rasterio para a parte da escrita tiff. Portanto, meu objetivo …

10 python-2.7  fiona  geopandas 

Estou escrevendo um script python usando a API QGIS. Eu tenho que obter resultados da tabela de banco de dados postgres. A tabela possui campo de geometria e outros campos que não são de geometria. Eu quero executar consultas SQL em campos não-geométricos e obter resultados no meu código. Como …

10 qgis  postgresql  pyqgis