Perguntas com a marcação «gdalwarp»

Eu tenho uma imagem raster geotiff que possui um sistema de coordenadas com longitudes de 0 a 360. O centro horizontal da imagem é 180 longitude. Veja a imagem abaixo: Quero transformá-lo em EPSG: 4326 SRS com -180 180 de longitude. E eu quero que o centro da imagem esteja no meridiano de...

Estou tentando recortar meu arquivo Raster de acordo com uma geometria do shapefile. Estou usando o seguinte código gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Mas está resultando uma cor preta fora da extensão da geometria do shapefile. Estou dando o exemplo aqui. Na primeira imagem que eu...

Estou tentando reprojetar / reamostrar com as ligações python GDAL, mas estou obtendo resultados ligeiramente diferentes em comparação aos do utilitário de linha de comando gdalwarp. Veja a atualização abaixo para obter um exemplo mais curto Este script ilustra a abordagem do Python: from osgeo...

Depois de usar gdalwarppara projetar e alinhar à grade (via -tap) um número de rasters, notei que os rasters de saída eram significativamente maiores que os rasters originais. Uma pesquisa na Web bastante completa mostrou esse problema do Trac : Frank Warmerdam explicou o motivo: "Em uma...

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 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...

No QGIS 2.XX, as ferramentas da barra de menus suspensa tinham uma caixa de chamada editável na parte inferior, onde os comandos não disponíveis na GUI podiam ser inseridos manualmente (clicando em um botão de lápis). No 3.0, as ferramentas da barra de menus são as mesmas que as ferramentas GDAL...

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...

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 usando o gdalwarp para manipular blocos SRTM perto da linha de dados (ou seja, 180 °, também conhecido como antimeridiano). Os blocos SRTM apresentam uma sobreposição muito pequena (1/2 pixel) com o meridiano. Você pode ver isso usando gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT...

Estou tentando remapear de geoestacionário para lambert conforme gdalwarp. Meus dados de entrada estão em netcdf e em coordenadas geográficas (graus) e eu gostaria de enviar os dados remapeados para o netcdf. Eu criei um arquivo vrt correspondente para os dados netcdf de entrada. Gdalwarp produzirá...

Quero georreferenciar uma varredura usando pythone GDAL. Minha abordagem atual é chamar gdal_translatee gdalwarpusar os.systemuma lista feia de pontos de controle de solo. Eu realmente gostaria de uma maneira de fazer isso nativamente por dentro python. Este é o processo atual que estou...

Eu tenho me esforçado um pouco para gerar blocos para uma imagem de alta resolução que temos. A imagem atual que temos é uma imagem muito grande (+ 20 GB), salva como um arquivo GeoTiff. Gostaria de gerar os blocos usando o utilitário de linha de comando gdal2tiles e depois abri-los e...

Estou usando python e QGIS 2.0. Estou tentando cortar rasters em uma pasta por um recurso de polígono. É a primeira vez que eu uso (digamos) "PyQGIS", eu já estava acostumado ao arcpy antes. De qualquer forma, não consigo que meu script simples funcione, qualquer sugestão seria muito...

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...

Estou tentando agregar alguns dados rasterizados (todos .tif) até uma resolução mais grossa (de 0,05 a 0,25 graus) usando gdalwarp em python, mas o comando não está funcionando. Em vez de obter uma saída com uma ampla faixa de valores, todos os valores da saída são 0. A resolução e a profundidade /...