Estou tentando realizar uma união em um campo comum após mesclar dois shapefiles adjacentes. Os arquivos de forma acabam com pelo menos um pedaço fino de espaço entre eles. Quando tento uma união, recebo o seguinte erro de orfão órfão: Erro no createPolygonsComment (p): rgeos_PolyCreateComment:...
Quero publicar um serviço de mapas no servidor ArcGIS, que será sobreposto no Google Maps ou Bing Maps. Para isso, o Serviço de Mapa precisa estar no EPSG: 3857/900913. A projeção do Google Maps deve assumir uma projeção esférica da terra. Então eu fui para Sistemas de coordenadas \ Sistemas de...
Eu tenho uma classe Featureclass, com um intervalo de 0,2 metros. Desejo selecionar apenas os contornos, que estão no intervalo de 1 m? Tentei usar a seguinte sintaxe na caixa de diálogo 'Selecionar por atributo' (no ArcMap), mas ela seleciona todos os contornos: Mod("ELEVATION", 1)=0 Que...
Olhando para esta página de informações um conjunto de dados hospedado no esri.com: http://fema-services2.esri.com/arcgis/rest/services/2012_Sandy/ImageCat_NLT/MapServer/layers Tentando descobrir como posso obter acesso aos dados brutos (coordenadas lat / lng brutas ou SHP). Parece que os dados...
Eu tenho experimentado técnicas de mapeamento de coropletas no R, tendo colhido os benefícios do ggplot2 por criar belos gráficos em um poderoso pacote de análise de dados. Quando se trata de mapeamento, não consegui produzir resultados comparáveis aos que eu sempre recebo do QGIS. Usando dados...
Eu tenho uma lista de algumas centenas de cidades com latitude / longitude. Dada outra localização (também em lat / long), preciso encontrar a cidade mais próxima. Como eu não uso nenhum SIG, o algoritmo óbvio agora é fazer um loop para todas as cidades, calculando a distância entre os...
Atualmente, estou usando o script a seguir para adicionar alguns dados de atributo de uma tabela a vários shapefiles individuais: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in 1:dim(specieslist)[1])...
Estou usando o ArcGIS 10.2 for Desktop, e eu tenho procurado neste site, e parece que ainda não consigo descobrir minha resposta. Muito novo para VBA e Python, mas muitos anos com ArcGIS. Eu sei que posso fazer isso da maneira mais lenta com Selecionar por atributos, mas é demorado. Estou tentando...
É possível criar visualizações do PostgreSQL sobre uma tabela que contém dados rasterizados? Consegui criar a exibição, mas a entrada raster_columnsnão está correta (todas as informações da varredura estão ausentes). Além disso, não consegui abrir a visualização SQL usando o QGIS (provavelmente...
Estou tentando obter um número contando quantos pontos existem em um shapefile. E isso funciona, exceto que estou tendo problemas para usar esse número em outro lugar. Eventualmente, usarei essa contagem em algumas contas (calculadora de campo), mas, durante a depuração, encontro um erro que...
Eu tenho um DEM de terra nua feito a partir de dados do Lidar (A). Funciona bem e passou em todos os testes de precisão. É usado em modelos de sobretensão. O problema que estamos enfrentando está nas passagens inferiores. Preciso permitir que a água flua através das passagens inferiores quando uma...
Por costume, não quero dizer uma estrutura diferente do arquivo KML, mas na <description></description>seção tenho dados com tags HTML desse tipo: <?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <Document> <Placemark>...
Estou tendo um pequeno problema com a API OGR do Python. O que estou tentando fazer é obter todas as coordenadas de cada vértice do anel externo de um polígono. Isto é o que eu tenho até agora: import osgeo.ogr import glob path = "/home/woo/maps/" out = path + 'output.txt' file =...
Esta postagem é um Wiki da Comunidade . Edite as respostas existentes para melhorar esta postagem. No momento, não está aceitando novas respostas. Estou desenvolvendo algum código em Python e quero usar as ligações python GDAL / OGR. Qual é o mínimo que preciso instalar...
Quero converter um GeoPDF que possua entidade de vetores nele para um formato vetorial e preservar as coordenadas reais desses vetores. O GeoPDF foi exportado do ArcMap e possui camadas. Eu tentei o GDAL, mas ele rasterizou o PDF antes de exportá-lo. O que não é bom para mim. Também não quero uma...
Portanto, a GDAL adicionou recentemente um novo recurso que permite a leitura aleatória dos arquivos de bucket do S3. Estou procurando cortar imagens GDAL de vários blocos de uma imagem sem ter que baixar o arquivo inteiro. Eu só vi documentação muito esparsa sobre como configurar e acessar um...
Estou tentando alterar o valor do campo de um shapefile. No entanto, isso não altera o valor. O que há de errado com o meu código? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' dataSource = driver.Open(fn, 0) layer = dataSource.GetLayer() feature =...
Com muita ajuda de gis.stackexchange.com, criei linhas de líderes para símbolos de pontos. As linhas são baseadas em uma simbologia de gerador de geometria: Encurtando a geometria da linha usando o diálogo de expressão no QGIS? Agora, quero implementar linhas de liderança para etiquetas...
Eu gostaria de generalizar a SpatialPolygonsDataFrame em R para tornar a plotagem mais eficiente. Os polígonos que resultam em enormes gráficos em PDF são lentos para carregar e difíceis de integrar em um documento. Tentei usar a gSimplifyfunção do rgeospacote, mas infelizmente ela não conserva a...
Estou tentando descobrir como carregar uma varredura em um banco de dados PostGIS2.0 (eu fiz perguntas anteriores sobre este tópico aqui e aqui ). Estou tentando usar o raster2pgsql.exeprograma fornecido com o PostGIS2.0. Depois de descobrir que o prompt de comando no Windows precisa ser executado...