Sistemas de Informação Geográfica

Eu tenho uma grande tabela PostGIS armazenando recursos de linha. Gostaria de mesclar recursos que estejam conectados entre si e que tenham valores de atributo idênticos. Aqui está um exemplo do que eu gostaria de alcançar. Os recursos mais extremos (preto e azul) não são mesclados porque não são contíguos, …

9 postgis  postgis-2.0  merge 

Estou tentando criar o choropleth()mapa correto em R usando o GISToolspacote. Eu uso choro.legend()para exibir minha legenda. No entanto, não consigo criar o posicionamento correto da legenda, E / OU Não entendo o papel do fmtparâmetro na minha choro.legend()função. Parece-me que fmtpode diminuir o espaço entre as cores e a …

9 r  legend  choropleth 

Estou tentando empacotar um script Python que usa GDAL para PyPI. Comecei incluindo uma referência direta no meu setup.py: install_requires=['GDAL==1.11.2'], Dessa forma, o pacote não foi instalado no meu ambiente virtual de teste: extensions/gdal_wrap.cpp:2855:22: fatal error: cpl_port.h: No such file or directory #include "cpl_port.h" ^ compilation terminated. error: Setup script …

9 python  gdal 

Estou trabalhando de maneira a criar projetos hidrológicos usando o QGIS e uma planilha do Excel que possuo. Para fazer isso, quero extrair algumas informações das linhas, incluídas em uma camada vetorial, que representa uma seção de tubo. As informações que eu preciso extrair são: Número de identidade comprimento Coordenadas …

9 qgis  vector  autofields 

Recentemente, observei uma mudança na projeção geográfica dos dados do Sentinel 2. Consulte as capturas de tela abaixo para obter esclarecimentos (somente Granule T33UXP, Banda 8). A primeira imagem foi tirada em 20160314, a segunda (deslocada) em 20160403. Adicionei um polígono como referência. São dados não processados, baixados do hub …

9 coordinate-system  remote-sensing  sentinel-2 

Eu tenho um número de grades hexagonais de 1 km que cobrem vários municípios nos Estados Unidos em um banco de dados postgreSQL / postGIS. Cada grade possui o CRS EPSG: 3857 e a camada de condados possui o EPSG: 3857. Ao visualizar as grades com os condados no QGIS, …

9 qgis  postgis  ogr2ogr  epsg 

Quero suavizar polígonos que compartilham um limite comum. Como posso suavizá-los para que os limites compartilhados permaneçam compartilhados? Estou usando o Arcgis 10.3

9 arcgis-desktop  polygon  arcgis-10.3  smoothing 

Eu tenho um conjunto de dados nacional de pontos de endereço (37 milhões) e um conjunto de dados poligonais de contornos de inundação (2 milhões) do tipo MultiPolygonZ, alguns dos polígonos são muito complexos, o máximo de ST_NPoints é de cerca de 200.000. Estou tentando identificar usando o PostGIS (2.18) …

9 postgis  point-in-polygon 

Uma camada de memória pode ser usada como entrada e / ou saída com o runalg ? Com runandload, eu faço isso funcionar: processing.runandload("qgis:mergevectorlayers","memory:bufferOne", "memory:bufferTwo", r"memory:merged") processing.runandload("qgis:dissolve", r"memory:merged", True, '', r"D:\PythonTesting\dissolved.shp") Os parâmetros 2 e 3 são as camadas de entrada que já estão na memória e o parâmetro 3 …

9 qgis  pyqgis  qgis-processing 

Este pequeno script exclui todos os recursos pequenos em todas as camadas. Isso funciona bem, mas o tamanho dos arquivos .shp não é reduzido. Por favor me diga o porquê... canvas = iface.mapCanvas() expr = QgsExpression("$length<100") for l in canvas.layers(): l.startEditing() it = l.getFeatures(QgsFeatureRequest(expr)) for ft in it: l.deleteFeature(ft.id()) l.commitChanges()

9 qgis  python 

Eu fiz um curso de treinamento do ArcGIS Online há algum tempo e aprendi um pouco mais sobre os pacotes de blocos. Segui as instruções, como criar um pacote de ladrilhos e criei um pacote de ladrilhos para algumas fotografias aéreas. O processo foi concluído e agora tenho um pequeno …

9 arcgis-desktop  arcgis-10.3  mbtiles  tile-package 

Tenho dados de pontos ao longo das ruas, gostaria de transformar esses pontos em simples linhas coloridas. Qualquer ponteiro como esse problema pode ser chamado ou qualquer algoritmo que possa me ajudar a resolver isso? Eu esperava usar PostGISfunções para fazer isso, mas estou aberto a sugestões, esses são dados …

9 postgis  line 

Estou tentando gerar polígonos para faixas orbitais de satélite. Até agora, tenho um método para gerar duas linhas que representam a borda de cada faixa em [lat, long]. Alguns dos trechos cruzam a linha de dados internacional e, por isso, envolvem-se: Consegui resolver isso com ogr2ogr -wrapdateline: ogr2ogr -wrapdateline -f …

9 python  ogr  wgs84  shapely  antimeridian 

Esta é uma questão de geometria do SQL Server aparentemente simples que eu pensei que teria uma solução pronta para uso, mas não estou tendo sorte em encontrar uma. Minha intenção é selecionar todos os registros em uma tabela que tenham polígonos aninhados (contidos) em um polígono maior de outra …

9 sql-server  intersection 

Eu tenho um grande projeto que consiste em cerca de 10x arquivos de projeto QGIS. Cada projeto contém um compositor de impressão com um atlas configurado para gerar até 10 mapas. Qual é a maneira mais fácil de automatizar a abertura de cada projeto e criar um PDF usando o …

9 pyqgis  print-composer  pdf