Estou tentando mesclar 60 arquivos geotiff em um arquivo enorme para poder produzir blocos desse enorme arquivo geotiff resultante usando gdal2tiles. O processo de mesclagem funcionou, mas, por algum motivo, produz um arquivo com um tamanho de 78 GB em vez dos 19 GB esperados. Não tenho certeza do...
Estou usando gSimplify(pacote rgeos) para simplificar as geometrias de um shapefile. A função funciona bem, mas agora não consigo gravar a saída em um novo shapefile. Eu tentei algumas maneiras: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') eu recebo obj deve ser um...
Estou tentando configurar um virtualenv python isolado para trabalhar em projetos GIS no meu Mac OS X. Parece que vou precisar instalar o GEOS, PROJ, GDAL / OGR no site kyngchaos aqui http://www.kyngchaos.com/ software / frameworks e não é fácil_instalá-lo no meu virtualenv. Isso é preciso? Os...
Eu tentei usar o código do QGis Save Raster como imagem renderizada . ié uma camada raster: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs =...
Estou fazendo mosaico de algumas imagens com gdal e gostaria de melhorar o resultado final usando uma camada alfa gradual / desbotada na borda de cada imagem para remover as bordas afiadas no meio do mosaico. O problema que estou tendo é que a parte de cada imagem individual com a camada alfa...
Estou trabalhando em uma projeção de seleção de locais em que não poderei realizar nenhuma observação in situ. Um componente-chave do processo de seleção é a robustez do terreno subjacente. Planejo usar o GDALdem para criar os mapas de robustez derivados. Eu li tanto o Wilson, et al. Documento de...
Houve muitas perguntas sobre isso, mas nenhuma delas resolveu meu problema. Eu tenho uma máquina com: Windows 7 x64 Python 3.4.3 x64 gdal instalado: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; gdal foi instalado a partir de uma roda GDAL-1.11.4-cp34-none-win_amd64.whlde aqui O...
O ArcGIS Desktop possui um conjunto de dados de mosaico de varredura , o QGIS e seus amigos têm tabelas de varredura virtual . Um pouco da documentação revela muita semelhança. Ambos criam ou habilitam uma representação alternativa de dados que fica em outro lugar sem ter que criar uma nova...
Estou tentando converter dados rasterizados em dados poligonais. Anteriormente, eu estava usando o shapefile.py para criar dados de pontos em lat / lon individuais na varredura. Mas isso criará uma enorme quantidade de dados e a consulta será um grande problema. Portanto, seria melhor se eu...
Estou trabalhando com o GDAL em Python e queria que o GDAL levantasse exceções para que eu pudesse fazer algo útil com elas. Por exemplo, se um arquivo não puder ser aberto com gdal.Open (), eu gostaria de lançar uma exceção em vez de gdal retornar None. Gostaria apenas de uma maneira de capturar...
Estou tentando gravar um shapefile em uma varredura usando o RasterizeLayer da GDAL. Eu pré-crio uma área de interesse de varredura a partir de um shapefile diferente, considerando um tamanho de pixel específico. Essa AOI serve como base para todas as seguintes rasterizações (mesmo número de...
Estou tendo problemas para usar a função Polygonize em python. O exemplo do livro de receitas para isso pode ser encontrado aqui . A parte relevante do meu código é: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile =...
Ao tentar instalar o Shapely no meu computador com Windows de 64 bits, não consigo fazer com que a biblioteca GEOS funcione. Até agora, eu executei o instalador OSGeo4W a partir do qual instalei o GDAL (acredito que a biblioteca geográfica esteja incluída nesse pacote). Depois disso, verifiquei e...
Eu tenho uma grade binária Arc / Info - especificamente, uma varredura de acumulação de fluxo ArcGIS - e gostaria de identificar todas as células com um valor específico (ou em um intervalo de valores). Por fim, eu gostaria de um shapefile de pontos representando essas células. Posso usar o QGIS...
Existe uma maneira de executar a mesma tarefa que o utilitário gdalbuildvrt usando as ligações GDAL Python? Até agora, não encontrei outra maneira de fazer isso, a não ser criar um vrt de um único conjunto de dados e editar manualmente o xml. Gostaria de criar um vrt a partir de vários rasters...
Estou usando o Python 2.6.5 (32 bits) com o Numpy 1.3 e o Gdal 1.9.1 instalados no Windows 7 de 64 bits. Estou tentando ler um conjunto de dados de varredura Imagine de 800 MB Imagine (.img) em uma matriz Numpy para fazer uma álgebra de varredura, mas assim que executo o código a seguir, o...
Movendo-se da lista de discussão gdal-dev: Em segunda-feira, 2 de setembro de 2013 às 19:09, David Shean escreveu: Olá, estou tentando empacotar séries temporais de rasters GTiff com projeção / extensão / resolução idênticas como um único arquivo NetCDF para distribuição. Passei a última hora...
Às vezes, escrevo breves scripts / instruções de trabalho para colegas que não são muito habilidosos em GDAL / programação em geral. Para evitar problemas com os scripts, ao mesclar / converter geotiffs com o GDAL, tento usar a opção BIGTIFF = YES, independentemente do tamanho do arquivo a ser...
Estou tentando usar o ogrinfo para obter alguns detalhes de um shapefile que baixei. Atualmente, a única maneira de saber como fazer isso é carregá-lo no QGIS e clicar manualmente em volta para encontrar qualquer informação, como abrir a tabela de atributos. Eu só quero ver qualquer metadado...
Eu sei que essa pergunta é bastante vaga, mas por favor, tenha paciência comigo. Estou tentando ter uma idéia de que tipo de desempenho do produto - especificamente o tempo - as pessoas viram por várias metodologias usadas para criar blocos de mapas do Google / Bing. Existem vários métodos para...