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
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 =...
O pacote geopy do Python apresenta duas técnicas de medição de distância: as fórmulas de Great Circle e Vincenty . >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) -
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...
Desejo criar programaticamente um arquivo GeoJSON usando polígonos de um shapefile, mas adicionando atributos do meu próprio aplicativo. Isso é feito facilmente para um shapefile: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist is the...
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...
Sou realmente novo em Python e gostaria de saber se existe um método rápido para obter valores de célula de um pixel raster por pixel e as coordenadas (mapear coordenadas XY do centro de cada pixel) usando Python no ArcGIS 10? Para descrever isso ainda mais, preciso obter o mapa X, o mapa Y e o...
Eu tenho um software (não o Arc) que é executado todas as noites em uma máquina de não produção que atualiza os recursos de shapefile de um banco de dados externo. Periodicamente, eu gostaria de excluir todos os recursos do shapefile (não o arquivo em si, que deve permanecer) e deixar o software...
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á 8 meses . Excluindo o ArcPy,...
Estou interessado em aprender a utilizar matrizes NumPy para otimizar o geoprocessamento. Grande parte do meu trabalho envolve "big data", onde o geoprocessamento geralmente leva dias para realizar determinadas tarefas. Escusado será dizer que estou muito interessado em otimizar essas rotinas. O...
Dado que geopandas GeoDataFramecontém uma série de polígonos, gostaria de obter a área em km² de cada elemento da minha lista. Este é um problema bastante comum, e a solução sugerida no passado foi usar shapelye pyprojdiretamente (por exemplo, aqui e aqui ). Existe uma maneira de fazer isso de...
Estou trabalhando com um suplemento do ArcMap em c #. No código C #, eu executei alguns scripts Python. Agora, para executar esses scripts, tenho um caminho python codificado. Mas isso não é portátil. Então, eu quero pegar o caminho do executável Python a partir do código e usá-lo. Questão: Como...
Estou tentando capturar todos os polígonos não sobrepostos indicados abaixo usando Shapely (dados polígonos A, B e C). Além disso, espero fazê-lo sem iteração, teste de interseção etc. A resposta aceita a essa pergunta expressa o método PostGIS, mas parece que 'união' significa coisas diferentes...
Dados dados pontos com longitude, latitude e um terceiro valor de propriedade desse ponto. Como posso agrupar pontos em grupos (sub-regiões geográficas) com base no valor da propriedade? Eu procurei pelo google e descobri que esse problema parece ser chamado de "agrupamento espacial restrito" ou...
Estou usando os dados GIS do meu município para atualizar as estradas no OpenStreetMap. Tenho uma cópia que eles publicaram no ano passado e uma que acabou de ser publicada e gostaria de encontrar todas as LineStrings que tiveram seus atributos ou sua geometria alterada. Isso me ajudará a garantir...
Tentando importar o gdal para Python 2.7.3 no Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File...
Alguém conhece uma biblioteca pura de Python (sem dependências de C ou Java) para realizar operações de geometria como buffer, união, interseção etc.? Estou procurando por algo que seja compatível com o Google
Eu tenho uma lista de quadros de dados, todos da mesma estrutura e todos são polígonos. Como posso adicioná-los todos juntos usando
é 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...
Eu tenho tentado criar uma função que faz basicamente a mesma coisa que o QGIS "dissolve" a função. Eu pensei que seria super fácil, mas aparentemente não. Então, pelo que entendi, o uso de fiona com bem torneado deve ser a melhor opção aqui. Eu comecei a mexer com arquivos vetoriais, então esse...