Perguntas com a marcação «r»

Eu gostaria de realizar uma sobreposição espacial para identificar o polígono no qual um conjunto de pontos se enquadra. No entanto, primeiro quero tamponar e dissolver os polígonos, de modo que os pontos que caem em qualquer lugar dentro dos polígonos mesclados (mas não dentro de quaisquer furos) sejam rotulados …

9 r  sp 

Eu criei dois SpatialLinesobjetos em R: . Esses objetos foram criados desta maneira: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Agora, quero concluir de alguma forma que essa é a mesma linha girada e invertida, e que a diferença entre eles é igual a 0 (ou …

9 r  geometry  similarity 

Estou tentando encontrar alguns centróides de um shapefile para os municípios dinamarqueses e depois encontrar o tempo de condução entre eles. Eu uso Ra readShapePolyfunção de maptoolscombinada com a gCentroidfunção de rgeose tudo funciona. No entanto, recebo pontos espaciais como SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) …

9 shapefile  r  coordinate-system 

Como posso encontrar um valor máximo de célula de uma varredura empilhada. Rmax <- maxValue(RAD1998.all[[1]]) funciona bem, mas Rmax <- maxValue(RAD1998.all[[2]]) dá NA. Claro que não em uma varredura empilhada. Aqui está o meu código: RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names …

9 raster  r 

Eu calculei a área de superfície das distribuições de espécies (mesclando polígonos de shapefiles), mas como essa área pode ser composta por polígonos bastante distanciados, eu gostaria de calcular alguma medida de dispersão. O que fiz até agora é recuperar os centróides de cada polígono, calcular a distância entre eles …

9 polygon  r  distance  centroids 

todo mundo, desculpe pela perturbação, mas estou sendo bastante novo com r enfrentou uma dificuldade crucial: quero criar um mapa animado de Russin com mudanças no desemprego com anos diferentes, como. Na imagem você pode ver dados de um ano require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header …

9 cartography  r  animation 

Estou tentando remover a interseção do Poly2 no Poly1. Abaixo estão as cordas dos 2 polígonos. > coords1 [,1] [,2] [1,] 992.0161 7462.531 [2,] 950.1962 7413.532 [3,] 902.3632 7421.207 [4,] 706.1985 7472.378 [5,] 654.0139 7711.760 [6,] 657.5960 7726.950 [7,] 786.2667 7844.380 [8,] 789.5935 7842.925 [9,] 981.7046 7498.659 [10,] 983.4246 7493.271 …

9 r  polygon  intersection  symmetricdifference 

Estou interessado em criar uma superfície raster que represente a intensidade solar média por unidade de área para uma pequena horta. Eu imagino que as entradas no modelo seriam as seguintes: Área da parcela - por uma questão de simplicidade, vamos assumir que a área do jardim é circular. Altura …

9 arcgis-desktop  python  arcpy  raster  r 

A imagem em anexo mostra uma lacuna na floresta com pinheiro vermelho representado como círculos e pinheiro branco representado como cruzes. Estou interessado em determinar se existe uma associação positiva ou negativa entre as duas espécies de pinheiros (ou seja, se elas estão crescendo ou não nas mesmas áreas). Estou …

9 arcgis-desktop  python  arcpy  r  spatial-statistics 

Na documentação de spvi que o coordsargumento da SpatialPolygonsfunção aceita apenas 2 colunas. Não é possível criar polígonos tridimensionais? Estou fazendo krigagem espaço-temporal usando pacote gstatadicionando tempo como uma terceira dimensão espacial de a SpatialPointsDataFrame. Eu gostaria de estender isso para bloquear a krigagem com polígonos como blocos. Existe alguma …

9 r  interpolation  spatial-statistics  kriging  gstat 

Preciso de uma grade espacial como uma grade principal para diversos mapas temáticos. Como produzo uma Grade espacial a partir de uma varredura descartando todos os pixels de NA?

9 r 

Existem 100 SpatialPointsDataFrames no meu Espaço de Trabalho. Carreguei-os desta maneira: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Eles são chamados assim: XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp etc. Como posso convertê-los em rasters que percorrem o espaço de trabalho? Eu sou um usuário R muito novo e …

9 r  rasterization 

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 …

9 raster  r  georeferencing 

Provavelmente uma pergunta muito simples. Eu tenho uma lista de cerca de mil locais geográficos em potencial (lat lon) e desses preciso selecionar os 200 locais que estão "mais espalhados". Eu acho que seriam os 200 pontos com maior distância média total. Pense em lojas em uma cidade. Existe um …

9 r  spatial-statistics  business 

Esta é uma referência cruzada do stackoverflow . Eu sei muito pouco sobre o software GIS e estou fazendo todo o meu mapeamento no R. Desculpas com antecedência, se essa é uma pergunta muito básica. Digamos que eu tenho dois shapefiles de fontes diferentes, mas com atributos diferentes. Digamos, um …

9 r  overlay