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 =...
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á 2 anos . Antecedentes Esta é...
Minha pergunta é sobre a extração de valores de rasters no local dos pontos. Com a extração da função, isso é muito fácil, e a função fornece um quadro de dados com os valores de todas as variáveis nos pontos. Eu quero ter nesse quadro de dados as coordenadas de cada ponto. Como posso fazer isso...
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 estava tentando criar uma imagem rasterizada de um banco de dados de pontos com espaçamento irregular. Os dados parecem- > head(s100_ras) x y z 1 267573.9 2633781 213.29545 2 262224.4 2633781 69.78261 3 263742.7 2633781 51.21951 4 259328.4 2633781 301.98413 5 264109.8 2633781 141.72414 6...
Eu quero criar dois polígonos. Um deles é da extensão retangular da minha varredura. Eu sei que posso obter a extensão da varredura usando r <- raster("band5.tif") e <- extent(r) plot(e)me dá o seguinte: Como posso criar UM polígono do limite da varredura, como mostrado abaixo? É assim...
Eu queria converter Raster em Point Vector no QGIS. É possível? O QGIS tem a opção de convertê-lo em polígono, mas não encontrei nenhuma ferramenta para convertê-lo em ponto. Alguém pode me
Eu tenho dois arquivos LiDAR (.las), um é original, digamos, com pontos X. E a outra é a cópia do primeiro arquivo .las, mas com pontos Y, onde Y é menor que X. Agora, quero comparar os DEMs (Digital Elevation Models) desses dois arquivos .las e visualizar como eles são diferentes. Quero obter...
Eu tenho uma grade binária Arc / Info - especificamente, uma varredura de acumulação de fluxo ArcGIS - e gostaria de identificar todas as células com um valor específico (ou em um intervalo de valores). Por fim, eu gostaria de um shapefile de pontos representando essas células. Posso usar o QGIS...
Estou trabalhando em um projeto que envolverá "dados raster" (imagens aéreas e DEMs) e "dados em grade", que eu entendo serem outros tipos de medidas numéricas tomadas em pontos da grade e que podem ser completamente não relacionados à imagem (por exemplo, densidade populacional). Estou tendo...
Eu tenho um conjunto de dados raster no formato ASCII. Usando Python, preciso determinar os valores mine maxdentro do conjunto de dados. Foi-me dito que as informações do cabeçalho são essenciais, que contêm itens como número de linhas / colunas, tamanho da célula e etc. Você não pode simplesmente...
Eu tenho uma variedade de tipos de habitat para uma área específica na Escócia. Preciso criar futuros cenários de habitat com mudanças no habitat para avaliar a viabilidade populacional de uma espécie de ave. Por exemplo, no futuro, pode haver 10% mais florestas na área. Gostaria de alterar o mapa...
Tentando entender o porquê, quando leio em uma varredura de NDVI, o slot @ data @ values não contém os valores reais até que eu os defina manualmente. Por exemplo: NDVI <- raster("./filename.tif", crs="+proj=longlat +datum=WGS84") NDVI@data@values ## returns: logical(0) Isso não aconteceu...
Apenas imaginando, é sempre melhor ter mais de 3 GCPs? Os GCPs devem ser distribuídos uniformemente? Onde posso encontrar alguns recursos gerais sobre a teoria disso? Além de tentar um erro, que simplesmente me ensina o que funciona naquele momento e não o que deve funcionar em determinados casos....
Eu criei uma camada de polígono a partir de rasters e só tenho um único ID de classificação (e acres calculados) para cada polígono. Gostaria de "dissolver" todos os recursos tocantes em polígonos individuais maiores. As 3 áreas circuladas em vermelho (abaixo) contêm polígonos / pixels verdes que...
Existe uma maneira de executar a mesma tarefa que o utilitário gdalbuildvrt usando as ligações GDAL Python? Até agora, não encontrei outra maneira de fazer isso, a não ser criar um vrt de um único conjunto de dados e editar manualmente o xml. Gostaria de criar um vrt a partir de vários rasters...
Eu li o código fonte de vários filtros raster QGis-1.7.4, computando inclinação, aspecto e curvatura. Existe uma fórmula no filtro que calcula a curvatura total que me intriga. O arquivo de origem está na versão atual do QGis, com o seguinte caminho: qgis-1.7.4 / src / analysis / raster /...
Eu tenho que transformá-lo em um UTM, a fim de tornar o buffer funcional. wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <-
Estou tentando converter dados rasterizados em dados poligonais. Anteriormente, eu estava usando o shapefile.py para criar dados de pontos em lat / lon individuais na varredura. Mas isso criará uma enorme quantidade de dados e a consulta será um grande problema. Portanto, seria melhor se eu...
Estou tentando gravar um shapefile em uma varredura usando o RasterizeLayer da GDAL. Eu pré-crio uma área de interesse de varredura a partir de um shapefile diferente, considerando um tamanho de pixel específico. Essa AOI serve como base para todas as seguintes rasterizações (mesmo número de...