Primeiro de tudo, eu sou um novato em GIS. Eu tenho uma nuvem de pontos lidar com dados de GPS e telemetria associados. Em seguida, transformo as coordenadas de todos os pontos em um sistema de coordenadas ECEF e, finalmente, transformo as coordenadas angulares em um sistema de coordenadas...
Desejo carregar arquivos .adf no R. Os dados são desta página: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Eu tentei o seguinte código que encontrei depois de algumas pesquisas na internet. O problema é que, na classe RasterLayer, recebo valores negativos que não deveriam estar...
Estou tentando projetar uma varredura. Em R existe a projectRaster()função para isso (abaixo de um exemplo totalmente reprodutível): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project...
Então, eu estou trabalhando em um projeto que está em desenvolvimento há mais de um ano, mas só recentemente me envolvi. Minha tarefa envolve converter aproximadamente 2.500 arquivos KML em outro formato, importar esses arquivos para um Geoserver e copiar as informações de estilo associadas a cada...
Fechado . Esta pergunta precisa de detalhes ou clareza . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Adicione detalhes e esclareça o problema editando esta postagem . Fechado há 5 meses . Estou tentando criar um buffer para...
Eu criei um data.frame no R. Quero obter SpatialLinesDataFrame como resultado, para que eu entenda que preciso converter o data.frame em linhas, as linhas em SpatialLines e as SpatialLines em SpatialLinesDataFrame. eu fiz isso filedata.frame=data.frame(matrix(file),ncol=14,
Trabalhando com o QGIS 2.4.0 Chugiak, tenho uma pergunta sobre o recorte de imagens raster. Tenho uma imagem aérea que quero recortar usando as extensões de um arquivo de forma. Até aí tudo bem, usei a ferramenta raster / extract / clipper para fazer isso. O problema é que o tamanho do arquivo da...
Quero criar um conjunto de dados de vetores de pontos de 10000 pontos (ou maiores) em uma varredura binária, em que os pontos devem ser restritos a áreas onde o valor da varredura é 1. Eu tentei os seguintes passos. Polygonize raster QGIS: Vetor -> Ferramentas de pesquisa -> Pontos...
Estou usando o método GDAL ReadAsArray para trabalhar com dados rasterizados usando numpy (especificamente reclassificação). Como minhas rasters são grandes, eu processo as matrizes em blocos, repetindo cada bloco e processando com um método semelhante ao exemplo de GeoExamples . Agora, estou...
Bairros da rainha e da torre são duas maneiras comuns de calcular estatísticas para uma célula focal. Também são conhecidos os bairros Moore e (von) Neumann. Eu acho que também existe uma "versão bispo", ou seja, do canto inferior esquerdo para o canto superior direito. Mas por que esses dois...
Preciso encontrar um algoritmo ou método que possa detectar latitude longitude pontos extremos em uma trajetória durante o pós-processamento , que poderá ser corrigido (trazido de volta ao caminho da trajetória com base em seus vizinhos). Como exemplo do tipo de pontos extremos que gostaria de...
Eu tento desenvolver um plugin para adicionar camadas postgis ao qgis a partir de um plugin como parte do aprendizado de python. O código para adicionar uma camada postgis é retirado do livro de receitas e funciona bem se eu executá-la no console python no qgis (se o botão OK for pressionado no...
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 atlas do...
Estou executando um servidor Windows Server 2012. PostgreSQL 9.3.5, 64 bits POSTGIS 2.1.3 QGIS Server 2.6.1-2 QGIS Desktop 2.8.3 e 2.12 Estou usando um tablet Microsoft Surface Pro 4 usando o QGIS Desktop 2.12 para editar algumas camadas armazenadas no banco de dados Postgresql. As camadas no...
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...
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...
Estou usando a biblioteca maptools do R para analisar o shapefile em uma lista de polígonos. A função readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") me fornece um objeto SpatailPolygonDataFrame . No meu exemplo, meu SpatialPolygonsDataFrame possui as seguintes...
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, tudo parece ótimo. Mas ... para...
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...
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...