Parece que isso deve ser direto, mas não consigo descobrir pela vida. Dado um GeoTiff existente contendo uma banda única, como você pode adicionar outra banda ao GeoTiff? A única maneira de fazer isso é usar GDALCreate () para criar um novo arquivo de saída com uma banda adicional e copiar a banda...
Estou tentando instalar o GDAL 2.X no CentOS7 sem construir a partir do código-fonte, porque precisarei configurá-lo em outros servidores e instalar pelo código-fonte será uma tarefa muito cansativa e demorada. Encontrei alguns repositórios como o opengeo e o ELGIS, mas eles são para o CentOS7. No...
Estou tentando otimizar meu código python usando pypy . Embora eu consiga importar a maioria das bibliotecas e execute com êxito, não consigo importar as bibliotecas GDAL. As mesmas bibliotecas quando importadas no Python 2.7 funcionam bem. Eu tentei adicionar o local GDAL ao PYTHONPATH e até...
Eu tenho uma varredura inteira para a qual gostaria de criar uma tabela de atributos de varredura usando Python e GDAL. Posso criar uma tabela de atributos de varredura GDAL em Python da seguinte maneira: >>> rat = gdal.RasterAttributeTable() Isso funciona bem, como podemos...
Eu tenho um arquivo de imagem tif e um arquivo mundial .tfw (criado a partir de dados do OSM com o Mapertive). Como posso (usando a linha de comando linux) converter isso em um arquivo GeoTIFF? (ou seja, eu quero ter um arquivo, não 2). Tenho certeza que é possível com gdal ou algo...
Eu sou extremamente novo no GIS. Estou usando gdalpara ler em um mapa de uso e cobertura de terra e preciso escolher a lat / lng de certos tipos de cobertura da terra para indexar em um conjunto de dados diferente, expresso apenas em lat / lng. Unfortuantely, eu não entendo a forma das...
Estou escrevendo um utilitário simples para cortar lotes de arquivos raster geotiff de várias bandas na mesma área (menor). Usando gdalwarp, posso cortar facilmente um arquivo usando um shapefile de recorte de polígono único: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif No...
Eu queria mesclar cerca de 20 dem tiffs cada um com cerca de 100 MB e usei o gdalwarp (eu sei sobre a opção de arquivo virtual): gdalwarp $(list_of_tiffs) merged.tiff No entanto, esse comando levou muito tempo e, após 15 minutos, estava processando a 5ª imagem do conjunto. Eu termino. Descobri...
Estou mesclando alguns shapefiles e tive alguns problemas ao fazê-lo no QGIS, por isso estou usando ogr2ogr diretamente. Estou fazendo isso (em um lote): ogr2ogr -overwrite %destination% %n1% ogr2ogr -update -append %destination% %n2% -nln all_new ogr2ogr -update -append %destination% %n3% -nln...
Eu quero criar dois polígonos. Um deles é da extensão retangular da minha varredura. Eu sei que posso obter a extensão da varredura usando r <- raster("band5.tif") e <- extent(r) plot(e)me dá o seguinte: Como posso criar UM polígono do limite da varredura, como mostrado abaixo? É assim...
Instalei o GDAL e verifiquei na linha de comando, ele está funcionando e quando executo o gdal2tiles.py na linha de comando, ele funciona, mas quando executo o mesmo arquivo gdal2tiles.py com python IDLE, recebo "Nenhum módulo chamado 'osgeo '"erro. Você pode
Acabo de exportar uma tabela do Postgis para o shp usando essas dicas, mas não consigo importar um shp para o Postgis usando a mesma biblioteca (ogr). Qualquer ideia? Muito obrigado
Eu sou um cara de vetores principalmente, mas para um projeto atual eu preciso trabalhar com um monte de velhos mapas de papel digitalizados (mapas de dano de bomba ww2 para Londres, se você estiver interessado!) Digitalizamos e georreferenciamos os mapas e agora queremos criar uma camada composta...
Existe uma maneira de definir todos os pixels de valor 0 em uma varredura para nodata usando GDAL? Gostaria de usar fillnodata para preencher alguns furos, mas os furos estão definidos como 0, então preciso excluí-los
Estou tentando processar alguns dados rasterizados usando ogr / gdal e não consigo obter a utilização completa de todos os núcleos da minha máquina. Quando apenas executo o processo em um único núcleo, obtenho 100% de utilização desse núcleo. Quando tento dividir em multicore (no exemplo abaixo,...
Veja este link para mais detalhes. O problema: Quero fazer um loop através de uma varredura contínua (que não possui tabela de atributos), célula por célula, e obter o valor da célula. Eu quero pegar esses valores e executar condicionais neles, emulando as etapas de álgebra do mapa detalhadas...
Na documentação do gdal_calc , é apresentada a calculadora raster de linha de comando com sintaxe numpy . Mais tarde, existem poucos exemplos em que um deles: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - significa definir valores de zero e abaixo...
Eu regularmente crio meus próprios rasters GeoTIFF com GDAL em Python, por exemplo: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band =
Eu tenho dados rasterizados em uma tabela do postgres que eu quero entrar no python como uma matriz numpy. Estou usando o psycopg2 para conectar ao banco de dados. Consigo baixar os dados, mas eles retornam como uma string (provavelmente um binário serializado). Alguém sabe como pegar essa string...
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...