Estou tentando alterar o valor do campo de um shapefile. No entanto, isso não altera o valor. O que há de errado com o meu código? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' dataSource = driver.Open(fn, 0) layer = dataSource.GetLayer() feature =...
No Python, usando GDAL, extraí a projeção de uma varredura como uma sequência WKT da seguinte maneira: wkt = dataset.GetProjection() # wkt is
Eu tenho um shapefile com polígonos. E eu tenho um arquivo raster global. Quero sobrepor os polígonos do shapefile na grade de varredura e calcular o valor médio da varredura para cada polígono. Como posso fazer isso usando GDAL, escrevendo os resultados no
Estou tentando processar alguns arquivos raster bioclimáticos, como os que podem ser baixados em http://www.worldclim.org/current (conjunto bioclim). Eles parecem ter valores nodata definidos de -3.4e+38acordo com o QGIS (olhando para a saída do gdalinfo, é -3.39999999999999996e+38). Parece que as...
Vejo muitas vezes no código Python GDAL que as pessoas fecham os conjuntos de dados no final de seu script. Por que faz sentido fechar um conjunto de dados no Python GDAL? Existem consequências se eu não fizer isso? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds =...
Estou preparando dados para inclusão em um nó da geonetwork, onde os mapas são exibidos via geoserver. O problema é que as células que "não têm dados" estão sendo exibidas em preto no geotiff e eu gostaria que a cor fosse definida como branca ou transparente. Estou usando o ArcView e não tenho...
Estou reprojetando rasters em python usando GDAL. Preciso projetar vários tiffs, desde as coordenadas geográficas WGS 84 até o WGS 1984 Web Mercator (Auxiliary Sphere), para usá-los posteriormente em Openlayers, juntamente com o OpenStreetMap e talvez com o Google maps. Estou usando Python 2.7.5 e...
Eu tenho uma imagem .tiff que estou tentando sobrepor no google earth. Para isso, eu preciso georreferenciá-lo primeiro. A projeção da imagem é UTM com NAD83 como referência. Conheço o lat / long dos quatro cantos da imagem. Eu estou fazendo da seguinte maneira: gdalinfo para obter as quatro...
Atualmente, estou trabalhando para importar dados climáticos do CANGRID (fornecidos como Surfer Grid ascii, arquivos ".grd") para o ArcGIS. A grade é de 95 linhas por 125 colunas de tamanho. Os metadados fornecem lat / lon de origem (canto inferior esquerdo), tamanho da célula (50 km), bem como...
Eu tenho várias fontes de bloco geradas usando gdal2tiles que gostaria de apresentar no mesmo mapa. Portanto, ao apresentar um bloco, preciso determinar de qual fonte servir, verificando os limites. Alguém sabe como calcular os limites posteriores para um único bloco baseado puramente em zoom, x e...
Acabamos de receber um grande conjunto de DEMs no trabalho e gostaria de gerar contornos a partir deles. Os DEMs têm uma resolução de 1m e um tamanho de 1kmx1km. Saída do gdalinfo: Driver: AAIGrid/Arc/Info ASCII Grid Files: 380000_6888000_1k_1m_DEM_ESRI.asc Size is 1000, 1000 Coordinate System is...
Alguém sabe qual é a maneira mais fácil de instalar o GDAL 2.0 em um MacOSX? A distribuição binária kingchaos suspeita usual e homebrew parecem estar presos no GDAL 1.11 Preciso de suporte para as propriedades do recurso TopoJSON que parecem ter sido incluídas no 2.0 EDIT 2017/01/07: a...
Desculpas se a pergunta a seguir é um pouco estúpida, mas eu sou muito MUITO novo nessa coisa toda de GIS. Estou tentando converter algumas imagens geoTiff projetadas para WGS84 usando gdal em python. Eu encontrei um post que descreve o processo para transformar pontos nos GeoTiffs projetados...
Como exemplo, vou pegar o seguinte bloco http://a.tile.openstreetmap.org/3/4/2.png e salvá-lo como "4_2.png". As coordenadas WGS84 desta peça pode ser calculado ou ler lá clicando na telha correspondente: 0 66.51326044311185 45 40.97989806962013 (West North East South) Como georreferenciar o...
Desejo desenvolver um portal da web que mostre as imagens produzidas dinamicamente pelo QGIS, mas o QGIS forneceu imagens tiff muito grandes e não é possível mostrá-las no site usando visualizadores de imagens. É por isso que decidi converter essas imagens tiff em jpg / jpeg para reduzir o tamanho...
Eu uso o GDAL em Python para importar o shapefile. Quero saber os nomes dos campos do arquivo, minha maneira atual é: fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) Mas, dessa maneira,...
é possível dividir um shapefile por recurso em python? (o melhor seria uma solução em que eu possa salvar temporariamente os objetos vetoriais resultantes na memória e no disco). O motivo: quero usar a função gdal rasterizeLayer com vários subconjuntos diferentes de shapefile. A função requer um...
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...
Ok, estou cansado de diferentes soluções de varredura no banco de dados que tentei nas últimas 2 semanas que não funcionam atualmente (varredura PostGIS) ou que não conseguem lidar com meu volume de dados (Rasterlite), então Estou lendo sobre o formato VRT. Fique comigo, gostaria das suas opiniões...
Breve: Quero alterar a resolução de uma varredura e suavizar a rampa cinza, como mostrado nas imagens abaixo. A preferência é usar GDAL, PIL ou Numpy. Descrição: estou transformando pontos em rasters com uma resolução de saída de 20 metros na Biblioteca Geoestatística de Alto Desempenho . Não...