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...
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...
Este é um processo que leva apenas alguns segundos no software GIS. Minha tentativa de fazê-lo no R usa uma grande quantidade de memória e falha. Existe algo errado no meu código ou isso é algo que R não pode fazer? Eu li que R pode trabalhar dentro do Grass, posso usar uma função Grass de dentro...
Qual a importância de corrigir o Landsat para atmosfera, neblina e converter DN em refletância? Que impacto isso tem nos resultados, digamos, para muitos dos índices de
Eu tenho um arquivo raster em que todas as células têm valores entre -3 e 7, nenhum valor de dados é -9999. Como posso calcular o número de células com um valor específico, por exemplo, 6? Posso usar a
Estou criando alguns rasters de ponto flutuante no GDAL com o driver GeoTIFF. Quando carrego as imagens resultantes no QGIS ou no Arc, a simbologia padrão coloca o mínimo em -3.40282e38 e o máximo em 3.40282e38, para que a varredura pareça exatamente cinza. Existe uma maneira de escrever o...
Eu tenho um conjunto de rasters (.vrt) com dados diários de umidade do solo. Quero somar os pixels de todos os rasters para ter uma medida por mês. No entanto, o mundo inteiro não é coberto todos os dias, o que resulta em valores nodata em locais onde as medidas existem nos outros dias. O que eu...
Eu importei uma imagem raster tif no PostgreSQL com este comando: raster2pgsql -s 32643 -I -M filepath.tif -F -t 100x100 public.databassename > filepath.sql E importou o arquivo SQL de saída dentro do banco de dados PostgreSQL executando o seguinte: psql -U postgres -d databasename -f...
Estou escrevendo um método Python que adiciona dois rasters e gera uma única saída raster. Por razões fora do meu controle, as extensões dos rasters de entrada são diferentes, mas elas se sobrepõem. É possível operar exclusivamente nos pixels de varredura de entrada sobrepostos nas 2 regiões...
Para um estudo sobre a expansão das redes de estradas na floresta, estou tentando extrair estradas das imagens do Landsat. Já temos composições nítidas e sem nuvens, nas quais as estradas são claramente visíveis a olho nu, mas a extração delas nas características da linha está se mostrando difícil,...
Eu estou tentando fazer a correspondência de histograma usando Python para melhorar o processo de mosaico de várias rasters sobrepostas. Estou baseando meu código no encontrado em: http://www.idlcoyote.com/ip_tips/histomatch.html Até o momento, consegui recortar a área sobreposta de duas rasters...
Fiz o download de vários blocos de varredura SRTM DEM em https://dwtkns.com/srtm30m/ e criei versões de colinas a partir dos blocos. No entanto, após uma inspeção mais detalhada, os blocos de colina agora têm linhas de borda que eu não sei como ocultar ou remover (como na imagem abaixo, que mostra...
Eu tenho tentado converter arquivos de velocidade do vento em rasters sem sucesso. Foi discutido em um fórum da Arcgis que eu segui. Isso envolveu a substituição de ponto-e-vírgula por espaços e a adição de um cabeçalho. Esperava incluir um dos arquivos .asc originais junto com uma versão .csv que...
Eu tenho que fazer uma análise de vários critérios para responder à pergunta: "qual é o melhor lote a ser desenvolvido". Alguns dos critérios são: distância do ponto de ônibus mais próximo (camada de pontos com pontos de ônibus) distância da loja mais próxima (camada pontual com lojas) qual é o...
Eu tenho um mosaico de um DEM derivado de 1m LiDar. Preciso enviar um subconjunto dos dados como um DEM de 10m. Atualmente, estou usando a ferramenta agregada no ARCGIS 10 para produzir um valor médio para cada novo pixel de 10 m. Algum conselho sobre se essa é a melhor técnica para tarefas como...
São 4 seções de imagens (NE, SE, NW, SW). Eles têm alguma sobreposição entre eles. Trabalhando com ortoimagens de 1m e 3 bandas. Às vezes, partes diferentes são transportadas em diferentes épocas / fontes, daí a diferença nas épocas do ano. As imagens estão no formato .tif. As propriedades de...
No meu trabalho diário, sou constantemente solicitado a calcular áreas de conjuntos de dados raster globais em projeção geográfica com resolução de 30 segundos de arco. Esses conjuntos de dados normalmente são o resultado de uma operação Combinar (um exemplo típico é uma classe de vegetação...
... em vez de manter os rasters como GEOTIFFs armazenados localmente? Eu tenho um projeto grande que exigirá disponibilizar aproximadamente 50.000 GEOTIFFs de 10.000 MB como WMS para um mapa da Web baseado em OpenLayers. Somente 1 WMS será exibido a qualquer momento, portanto, o desempenho não é...