Eu tenho o NE lat lnge SW lat lng. Meu objetivo é simples, estou usando o RStudio e quero criar uma caixa delimitadora a partir dos dois lats e longos acima e, em seguida, gerar Polygon e quero adicionar pontos aleatórios no polígono quadrado. Eu tenho NELat / Lng e SW-Lat / Lng em um arquivo csv...
Eu tenho 12 imagens que representam as condições climáticas do dia 1 ao dia 12. Estou usando o gdal2tiles para produzir blocos de mapa de cada uma dessas imagens. O gdal2tiles também produz um arquivo KML que pode ser usado no Google Earth para exibir os blocos. No entanto, quero criar uma camada...
Eu tenho um conjunto de rasters (.vrt) com dados diários de umidade do solo. Quero somar os pixels de todos os rasters para ter uma medida por mês. No entanto, o mundo inteiro não é coberto todos os dias, o que resulta em valores nodata em locais onde as medidas existem nos outros dias. O que eu...
O script básico para iterar recursivamente sobre subpastas e mesclar todos os shapefiles em um único é: #!/bin/bash consolidated_file="./consolidated.shp" for i in $(find . -name '*.shp'); do if [ ! -f "$consolidated_file" ]; then # first file - create the consolidated output file ogr2ogr -f...
Estou procurando escalar uma imagem tiff de 16 bits não assinada para uma imagem de 8 bits usando GDAL (OSGeo4W). No entanto, quero restringir os valores de pixel para que eles fiquem dentro do intervalo das imagens de 8 bits. O que quero dizer é que quero garantir que a proporção visual na imagem...
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...
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...
Tenho uma grande quantidade de imagens do RapidEye (+300 Gb) que quero usar como mapa base para um aplicativo online. Até agora, eu consegui servi-los como um ImagePyramid no geoserver, seguindo as etapas mostradas em uma apresentação chamada "Geoserver on steroids". O problema é que não consegui...
Eu tenho o arquivo GeoTIFF de 1 banda. A interpretação das cores da banda é cinza. Quero uma tabela de cores para transformar meu Col.tiffem Col.jpegou Col.png. Eu tentei colocar uma tabela de cores neste tiff através de um .vrt: gdal_translate -of VRT Col.tiff Col.vrt adicionando no...
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",
Carreguei com sucesso um monte de shapefiles em um banco de dados PostGIS. Eu gostaria de poder recuperar geometria e camadas específicas como arquivos dxf com base em consultas SQL. Estou usando PostGIS, GDAL, PostgreSQL e Python da suíte OSGeo4W no Windows. Estou imaginando que eu deveria estar...
Eu tenho usado GDAL com python para manipular shapefiles. Ao abrir um arquivo de origem, primeiro defino o driver como “ESRI Shapefile” antes de chamar o método "Open" do driver com o shapefile necessário. Eu não tive problemas com isso; no entanto, estou curioso para saber qual é a diferença...
Estou tendo muitos problemas com a GDAL. Além de, às vezes, falta documentação, parece haver pouco suporte no Python. De qualquer forma, este é apenas um script de reprojeto em lote que segue assim: source_file = gdal.Open(filepath.encode("ascii")) source_wkt =
Eu segui vários guias de ajuda on-line sobre como fazer o ecw para o QGIS 1.8 funcionar no Ubuntu 12.04. Eu tinha tudo trabalhando no Ubuntu 10.04, mas desde a migração para o 12.04 não tive essa sorte. Existe um guia que funcione para o Ubuntu 12.04? No momento, parece mais fácil instalar a...
Estou tentando criar um mapa base de blocos. Eu tenho uma imagem da qual quero criar o bloco e estou usando o GDAL. Consegui criar os blocos, mas quando visualizo o resultado, parece que o mapa lado a lado é: É invertido As proporções não estão certas Esta é a imagem de sobreposição do...
Quero fazer uma imagem multiespectral de cero para fazer alguns testes nela. Algo realmente simples, como 5 bandas completamente uniformes com ruído de sal e pimenta ou um quadrado de valores diferentes no centro. Claramente, isso seria apenas uma pilha de matrizes, uma matriz multidimensional, que...
Como posso tirar uma foto e ortorretificá-la automaticamente com software de código aberto? As informações associadas a cada imagem são: especificações da câmera latitude, longitude, elevação da câmera orientação da câmera: de frente para um UAV, a inclinação do rolo e a guinada são dadas Eu...
É possível capturar mensagens de erro de saída do PostgreSQL no OGR após uma falha na consulta? Por exemplo conn.ExecuteSQL('SELECT * FROM non_existing_table;') No Postgres, isso retorna uma mensagem de erro: "Error: relation 'non_existing_table' does not exist." No ogr, obtemos o...
Preciso criar os contornos de polígonos de várias imagens raster de camada única, não a caixa de extensão / delimitação, mas a área sem os valores nodata, como mostrado aqui: Criando shapefile mostrando pegadas de Rasters? . Na resposta à pergunta acima, o plugin Image Boundary foi mencionado, no...
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...