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é...
Para obter um único valor de pixel em um ponto de uma varredura usando rasterio, há um exemplo aqui: https://github.com/mapbox/rasterio/pull/275 No entanto, existe uma API direta no rasterio (e não no cli) que pode ser usada para extrair valor em um único ponto em um raster? - EDITAR with...
É possível ler dados brutos em a geopandas GeoDataFrame, a la a pandas DataFrame? Por exemplo, o seguinte funciona: import pandas as pd import requests data =
O script que estou escrevendo reposiciona dois quadros de dados e define suas extensões. Ao fazer isso, ele redesenha todo o Active View 4 vezes, diminuindo substancialmente o script. Pressionar F9 ou clicar no botão 'Pausar desenho' antes de executar o script desativa o desenho e faz com que o...
Eu quero desenvolver uma GUI em Python para o geoprocessamento ArcGIS. Posso usar o PyQT para programação de GUI no ArcGIS? Também considerei usar o Tkinter e o wxPython.
Eu me deparei com um problema estranho ao executar scripts de processamento PyQGIS independentes. A ordem das importações no script afeta sua execução normal. Você pode reproduzir o problema abrindo um console Python e inserindo o seguinte script (eu uso o GNU / Linux, QGIS 2.6.1, plugin de...
Estou servindo blocos de vetor usando o TileStache , tenho tudo configurado como eu quero. Meus dados são armazenados no Postgres e estou usando o provedor VecTiles para servir blocos GeoJSON . Quero armazenar em cache todos os meus blocos para que eles sejam mais rápidos. Estou usando...
parece uma pergunta boba, mas não consigo encontrar uma resposta clara sobre isso: que unidades geopandas / bem torneadas usam o cálculo da distância / área entre os
eu tenho GeoPandas 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Python 2.7.10 instalado no ambiente Anaconda 2-4.1.1-Windows-x86. Eu sou capaz de construir GeoDataFramelendo conjuntos de dados de entrada e operar os dados, mas salvar um conjunto de dados de saída não preserva o sistema de...
Gostaria de chamar a função de plug-in de interpolação QGIS (método TIN) (Raster-> Interpolate) no console python. Não consigo encontrar a função correspondente na API do QGIS ou na lista do algoritmo de processamento. Encontrei o algoritmo de triangulação SAGA, que funciona bem, mas é 5-10...
Existe uma maneira de executar um plug-in personalizado feito para o QGIS que um Shell Python possa aparecer quando o plug-in for executado, como uma tela de depuração? Estou pedindo para saber se existe uma maneira de ver comandos de impressão por meio de um script Python durante a programação de...
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,...
Eu tenho geopandas GeoDataFrame contendo centenas de formas Polygone MultiPolygonformas geométricas. Os polígonos se sobrepõem em muitos lugares. Gostaria de fazer uma nova geometria que contenha quantas delas se sobrepõem. Algo assim: Alguém tem alguma idéia de como abordar isso? Eu nem consigo...
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 =
Fechado . Esta questão precisa ser mais focada . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela se concentre apenas em um problema editando esta postagem . Fechado há 4 anos . Estou procurando um...
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...
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...
Eu sou novo nisso. Estou "construindo" um plugin python para qgis e quero importar módulos do GRASS e SAGA, mas não entendo como fazê-lo. O plugin que quero importar é a interpolação Kriging do SAGA ou v.surf.idw do GRASS ou o plug-in de interpolação do Raster qgis. Eu quero um método de...
Estou tentando alterar um plugin do QGIS 1.8 para 2.0. Na versão 1.8, crio uma tabela com os nomes dos campos assim: layer = vlayer provider = layer.dataProvider() feat = QGSFeatrue() allAttrs = provider.attributeIndexes() provider.select(allAttrs) fieldNames = [] while
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