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...
Eu tenho duas camadas no ArcGIS: uma é um ponto vetorial e a outra é uma camada raster. Quero acrescentar o valor da grade de dados raster ao ponto do vetor. O problema é que o ponto do vetor está longe de qualquer grade raster existente com valor; portanto, o comando "Amostra" não retorna nenhum...
Estou interessado em uma maneira de identificar rasters (por exemplo, formato ESRI GRID) ou shapefiles corrompidos que impediriam que eles fossem usados em uma ferramenta de análise. Eu tenho um grande número na árvore de pastas estruturadas e gostaria de poder criar um script para algo que...
Gostaria de importar todos os múltiplos rasters de mosaico (cada mosaico é composto por 7 camadas) da pasta indicada para o meu R. Em seguida, acesse-os como rasters de multicamadas únicas. Então eu fiz: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(),...
Eu tenho um DEM de uma região montanhosa. Gostaria de extrair as áreas do DEM acima das quais estão 60% da área total do DEM. Então, basicamente, eu gostaria de extrair as áreas que representam partes superiores das montanhas. Já ouvi falar dessa linha hipsométrica (linha H60), mas não consegui...
Referindo-me à minha pergunta anterior , como "geocode" / converto - usando R - uma imagem raster gerada para um formato que pode ser mostrado como uma camada em um mapa? Em resumo, pretendo recriar os resultados do software comercial SpatialKey ( captura de tela ). Estou usando uma instância do...
Estou tentando criar uma varredura usando as bibliotecas python gdal e cheguei ao ponto em que os dados estão sendo produzidos, mas os dados de saída são invertidos no eixo x do ponto de origem . Eu sei que devo estar ignorando alguma coisa, mas não consigo descobrir onde estou errado. Alguma...
Eu tenho duas imagens raster, classificadas como booleanas (1 ou 2). Eles se sobrepõem em determinadas áreas. Eu preciso de uma varredura que é apenas as seções sobrepostas. Usando o Arcinfo 10. A resposta pode usar apenas métodos raster, sem conversão em polígonos, pontos,...
Gostaria de passar de uma varredura binária de floresta / não floresta com resolução de 30m para uma varredura com resolução de 240m, com o valor de cada célula sendo a proporção dessa célula florestada - ou seja, os valores de dados para as células de 240m passarão de 0 (onde todas as células na...
Minha formação: programador experiente, historiador amador, iniciante em GIS (estou usando o QGIS). Aprendi a usar a ferramenta de georreferenciador QGIS. A entrada: suponha que eu tenha uma imagem como esta Wikipedia: File: VoyagesOfRabbanBarSauma.jpg . Suponha que seja um arquivo JPG "nu" sem...
Eu tenho uma varredura GeoTIFF no ArcGIS Desktop, que pode ser baixada nos dois links a seguir: http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tfw http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tif Se você carregá-los no ArcGIS Desktop, verá grandes áreas...
Estou tentando projetar uma varredura. Em R existe a projectRaster()função para isso (abaixo de um exemplo totalmente reprodutível): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project...
Desejo carregar arquivos .adf no R. Os dados são desta página: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Eu tentei o seguinte código que encontrei depois de algumas pesquisas na internet. O problema é que, na classe RasterLayer, recebo valores negativos que não deveriam estar...
Gostaria de adicionar atributos a uma varredura armazenada em um banco de dados postgis. Eu uso o seguinte comando para inserir a varredura no banco de dados: raster2pgsql -s 4236 -I -C -M /home/user/mysrc/outlabeled.tif public.test | psql -d template_postgis A varredura outlabeled.tifé de fato...
Trabalhando com o QGIS 2.4.0 Chugiak, tenho uma pergunta sobre o recorte de imagens raster. Tenho uma imagem aérea que quero recortar usando as extensões de um arquivo de forma. Até aí tudo bem, usei a ferramenta raster / extract / clipper para fazer isso. O problema é que o tamanho do arquivo da...
Estou tentando cortar um grande ECW (detalhes abaixo), mas o arquivo raster é muito grande para ser processado completamente. Alguns detalhes da ECW abaixo Driver: Wavelets compactadas ECW / ERDAS (SDK 5.0) tamanho do arquivo: 50gb O tamanho é 450000, 565081 Tamanho do pixel: 0,15 0,15...
Eu tenho uma imagem raster de elevações. Eu adicionei um shapefile de vetor que contém divisões políticas locais. Quero manter as cores da elevação, mas também quero manter as linhas pretas representando as divisões políticas. Reduzi a transparência do arquivo de forma das divisões políticas, mas...
Estou tentando melhorar um processo de vetor / python atualmente extremamente complicado para um modelo de risco natural. No momento, temos um longo script que gera distância / linhas de rumo a partir de um determinado ponto para determinar: o tipo de polígono que cruza (por exemplo: floresta,...
Compilei e instalei o QGIS master atual como de costume e importei alguns dados rasterizados no meu banco de dados do postgres 9.3 com a extensão postgis 2.0.2. Sempre que tento importar os dados rasterizados no qgis usando o gerenciador de banco de dados, recebo esta mensagem: GDAL provider:...
Eu tenho duas camadas .cw sobrepostas da mesma área, embora uma esteja cortada no meio, com uma borda branca. Quando defino o valor de segundo plano como 255.255.255 e o defino como oco, ele apenas oculta esse valor RGB específico; quando, na realidade, existem outros valores "brancos" na minha...