Estou tentando gravar várias camadas no mesmo GeoPackage no R, mas recebo um erro Creation of output file failed,. Eu tentei procurar documentação de leitura e gravação em arquivos .gpkg com RGDAL, por exemplo, para descobrir se writeOGR()realmente suporta várias camadas, com pouco sucesso. Isso é...
Já vi essa pergunta antes, mas as respostas que estou tentando não estão funcionando. Eu estou olhando para consultar todas as maneiras OSM que estão dentro de uma caixa delimitadora. Os dados do OSM foram importados com o mercator esférico padrão. Estou consultando usando LAT / LON, portanto, a...
Alguém tem uma maneira eloquente de empilhar vários arquivos .tif em uma pilha de várias bandas usando o Rasterio e / ou o GDAL? Estou procurando uma maneira de evitar o uso de um comando de subprocesso como gdal_merge.py e, em vez disso, tê-lo como parte do meu script python. Eu sei que tanto o...
I têm uma enorme quantidade de dados na forma de CSV contendo coordenadas UTM como Xe Ye um valor de elevação como Zinformação. Preciso converter esses dados em um DEM como GeoTIFF para análise posterior. Nesse caso, uma quantidade enorme significa 16 m. linhas, com um ponto em X, Ye Zpor linha. Os...
Estou criando uma ferramenta que baixa imagens do "Sentinel-2 na AWS" Estou com problemas para traduzir coordenadas Long / Lat para um bloco Sentinel específico. Eu sei que o Sentinel 2 usa o MGRS e tentei simplesmente converter Long / Lat em MGRS e, com baixa precisão, obtenho o ID do bloco S2A....
Preciso saber se é cientificamente válido fundir / afiar os produtos de refletância de superfície Landsat 8 com uma banda de pan dessa respectiva banda? Os detalhes do produto de refletância Landsat podem ser encontrados aqui . É necessário mencionar que é necessário solicitar o produto de...
Estou usando o QGIS 2.18. Eu tenho shapefile que contém várias linhas. Após o processamento de certos dados, preciso desenhar linhas paralelas a diferentes distâncias das linhas existentes. Por exemplo, eu tenho 32 linhas no shapefile existente. Quero desenhar linhas paralelas a uma distância de...
Que tipo de visualizações / estatísticas posso obter no índice GiST do PostGIS? Se olharmos para a entrada da árvore R da Wikipedia , existem diferentes variações e visualizações da árvore R nas caixas delimitadoras de cada uma delas usando o ELKI. O índice GiST é semelhante ao índice R-Tree....
Eu tenho um shapefile na projeção da British National Grid: Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936",
Já faz algum tempo que estou tentando usar o r.mapcalc na QGIS Processing Toolbox. No entanto, o algoritmo não me fornece nenhuma saída. Eu tentei várias camadas de varredura e expressões diferentes, mas o algoritmo nunca gera nada para o "Diretório de Resultados". Abaixo há uma imagem que mostra...
Digamos, por exemplo, se houver dados pontuais com árvores antigas, quando você usar a ferramenta de identificação, poderá ver a fotografia da árvore Ou ainda mais um nível, os rótulos que exibem essas imagens imediatamente, sem nenhum clique necessário Eu acredito que o termo é streaming de...
Gostaria de alinhar vários pontos a uma linha, ou linhas dentro de uma camada, usando uma tolerância especificada ou buffer em torno dos objetos de linha. Consulte o exemplo de desenho em anexo. Para fins de exemplo, os pontos mais próximos da linha na imagem ANTES estão dentro de 5 unidades de...
Adicionei dois elementos de escala à minha composição de mapa. Para harmonização de layout, quero os carrapatos e o texto do rótulo abaixo da barra de escala. Padrão é o seguinte: O que eu quero é isso (alterado manualmente na imagem): Existe uma opção no QGIS para isso?...
Quando abro o QGIS, o painel Camadas é padronizado com uma largura que não mostra o botão "Remover camada / grupo". Então eu tenho que redimensionar o painel para mostrar o botão toda vez que eu abrir o QGIS. É possível definir uma largura padrão para este painel para que eu não precise...
Gostaria de acelerar o processo correspondente a este comando: gdalwarp -t_srs EPSG:4326 -overwrite input.ntf output.tif Existe uma maneira de usar o processamento paralelo no GDAL? Semi-offtopic: Se não, você recomenda uma solução não GDAL para acelerar um processo semelhante ao gdalwarp? Aqui...
Quero gerar IDs automaticamente enquanto crio novos recursos. Portanto, toda vez que eu criar um novo ponto, seu ID aumentará em 1 como o seguinte: Feature 1: id 1 Feature 2: id 2 Feature 3: id 3 Como posso atacar
Toda vez que crio um script para manipular vetor (ponto, linha e polígono) em python, encontro essa dúvida. Qual é a diferença entre caixa delimitadora (bbox), envelope, extensão,
Eu tenho um projeto QGIS que contém um compositor. É claro que se eu entrar na GUI, posso exportá-lo como um pdf, etc. Para o meu script make, eu quero que isso seja feito através de um script python. Da leitura online aqui está a minha tentativa: #!/usr/bin/env python import sys from qgis.gui...
Às vezes, vejo srsName="EPSG:4326"nos documentos GML (por exemplo , https://postgis.net/docs/ST_GeomFromGML.html , embora minha fonte específica não seja pública). A intenção (exceto, talvez, a ordem dos eixos que não é importante aqui) é clara e provavelmente interoperável em vários sistemas, mas...
Estou trabalhando em um projeto que requer um buffer de 30m em torno de um recurso linear. O recurso representa a costa e preciso que o buffer mostre 30m para o interior da linha que representa a costa. Existe uma maneira de fazer isso usando o