Perguntas com a marcação «numpy»

Eu sou novo no GIS. Eu tenho algum código que converte imagens infravermelhas de Marte em mapas de inércia térmica, que são armazenados como matrizes numpy 2D. Eu tenho salvado esses mapas como arquivos hdf5, mas eu realmente gostaria de salvá-los como imagens rasterizadas, para poder processá-los...

Eu tenho tentado verificar meus filtros no DEM raster para reconhecimento de padrões e isso sempre resulta em falta das últimas linhas e colunas (como 20) . Eu tentei com a biblioteca PIL, carregamento de imagem. Então com entorpecido. A saída é a mesma. Eu pensei que algo estava errado com meus...

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

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

Comecei a aprender a manipular dados LAS em python e queria ver como outras pessoas lidam com arquivos LAS. Gostaria de ler os pontos (estou usando uma matriz numpy) e filtrar as classes 1 e 2 (não classificadas e fundamentadas) em uma matriz separada. Eu tenho o código a seguir, mas não consigo...

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

Estou usando o ArcGIS 10.1 e quero criar uma nova varredura com base em duas rasters existentes. O RasterToNumPyArray tem um bom exemplo que eu quero adaptar. import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape =...

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

Eu tenho dados rasterizados em uma tabela do postgres que eu quero entrar no python como uma matriz numpy. Estou usando o psycopg2 para conectar ao banco de dados. Consigo baixar os dados, mas eles retornam como uma string (provavelmente um binário serializado). Alguém sabe como pegar essa string...

Estou usando o Python 2.6.5 (32 bits) com o Numpy 1.3 e o Gdal 1.9.1 instalados no Windows 7 de 64 bits. Estou tentando ler um conjunto de dados de varredura Imagine de 800 MB Imagine (.img) em uma matriz Numpy para fazer uma álgebra de varredura, mas assim que executo o código a seguir, o...

Estou tentando extrair valores de quebra naturais de uma varredura usando Python. O processo pretendido vai: Objeto de varredura do Arcpy para NumPy Array (usando RasterToNumPyArray ) até valores de quebra naturais (usando a função PySAL Natural Breaks ). import arcpy, pysal from...

Estou tentando aprender as cordas do processamento de imagens de Sensoriamento Remoto usando ligações Python GDAL e numpy. Como primeira tentativa, estou lendo um arquivo geotiff Landsat8, faça uma manipulação simples e escreva o resultado em um novo arquivo. O código abaixo parece funcionar bem,...

Como posso obter coordenadas projetadas e os valores reais de pixel nessas coordenadas de um arquivo GeoTiff e salvá-los em uma matriz numpy? Eu tenho o arquivo arsenci020l.tif e suas coordenadas estão em metros. Abaixo está a saída abreviada do gdalinfo que eu executei nele. ~$ gdalinfo...

Gostaria de reclassificar um arquivo rasterizado de um raster com 10 classes para um raster com 8 classes usando pyhton, gdal e / ou numpy. As classes são representadas como números inteiros. Eu tentei seguir as etapas deste post Reclassificar rasters usando GDAL e Python , o numpy.equal doc e...

Estou usando o método GDAL ReadAsArray para trabalhar com dados rasterizados usando numpy (especificamente reclassificação). Como minhas rasters são grandes, eu processo as matrizes em blocos, repetindo cada bloco e processando com um método semelhante ao exemplo de GeoExamples . Agora, estou...

alguém conseguiu fazer pandas trabalharem com a instalação Python do QGIS, no Windows? Usei o instalador do OSGeo4W para instalar o QGIS2.2, que vem com o Python 2.7 e o Numpy 1.5.1. Qualquer versão razoavelmente nova do pandas requer o Numpy 1.6.1 ou mais recente. Então esse parece ser o meu...

Hoje, atualizei para o ArcGIS Desktop 10.2.2, bem como a ferramenta de processamento em segundo plano de 64 bits, e notei que tanto o numpy quanto o matplotlib foram instalados com o Python. De fato, na seção " O que há de novo? " Da 10.2.1, afirma: O ArcGIS 10.2.1 foi atualizado para incluir o...

Tenho um problema com a biblioteca GDAL ao abrir um arquivo GeoTiff grande, com tamanho de 32000x32000. Não consigo usar a ReadAsArrayfunção porque o tamanho máximo da matriz numpy em python. Mas estou me perguntando por que o QGIS pode abrir esse arquivo facilmente. Qual é a técnica por...

O problema consiste em modelar a propagação de um sinal (por exemplo, luz ou som etc.) através de uma série de obstáculos, como na figura abaixo. O sinal não pode passar pela superfície inferior (terreno), mas pode atravessar obstáculos. Eu quero contar o número de obstáculos...