2

Eu tenho uma matriz com duas variáveis numéricas: lat and long. Tipo isto:

> head(pontos_sub)
  id       lat      long
1  0 -22,91223 -43,18810
2  1 -22,91219 -43,18804
3  2 -22,91225 -43,18816
4  3 -22,89973 -43,20855
5  4 -22,89970 -43,20860
6  5 -22,89980 -43,20860

Agora eu faço um round para arrendondar com 3 dígitos decimais:

pontos_sub$long_r <- round(pontos_sub$long, 3)
pontos_sub$lat_r <- round(pontos_sub$lat, 3)

> head(pontos_sub)
  id       lat      long  long_r   lat_r
1  0 -22,91223 -43,18810 -43,188 -22,912 
2  1 -22,91219 -43,18804 -43,188 -22,912
3  2 -22,91225 -43,18816 -43,188 -22,912
4  3 -22,89973 -43,20855 -43,209 -22,900 
5  4 -22,89970 -43,20860 -43,209 -22,900
6  5 -22,89980 -43,20860 -43,209 -22,900

Ahora quero usar o pacote dplyr para encontrar, agrupado por cada long_r lat_r e usando a função distVincentyEllipsoid, a distancia mínima a todos os lat long do grupo correspondente. Algo assim:

> newdata <- pontos_sub %>% 
               group_by(long_r,lat_r) %>% 
               summarise(min_long = special_fun(arg), 
                         min_lat = special_fun(arg))

O que resultaria algo como isto:

> head(newdata)
  long_r   lat_r   min_long  min_lat
1 -43,188 -22,912   xxxxxx   xxxxxxx
4 -43,209 -22,900   xxxxxx   xxxxxxx

Para finalizar, gostaria de saber se esta é a forma mais rápida, pois tenho milhares de linhas... existe alguma outra maneira de fazer isto bem rápido?

2
  • Você quer a distancia dos lat e long sem arredondar para os lat_r e long_r arredondados? Commented 7/04/2017 às 16:43
  • Para cada combinação única de lat_r e long_r, quero a distância de todas as combinações lat e long, correspondentes a tal lat_r e long_r.
    – MAOC
    Commented 7/04/2017 às 16:47

1 Resposta 1

2

Eu já tentei fazer algo parecido há um tempo atras. Fiz da forma que está abaixo, foi o melhor que consegui.

library(dplyr)
library(tibble)
library(tidyr)

pontos_sub <- tribble(
  ~id, ~lat, ~long, ~long_r, ~lat_r,
  0, -22.91223, -43.18810, -43.188, -22.912, 
  1, -22.91219, -43.18804, -43.188, -22.912,
  2, -22.91225, -43.18816, -43.188, -22.912,
  3, -22.89973, -43.20855, -43.209, -22.900, 
  4, -22.89970, -43.20860, -43.209, -22.900,
  5, -22.89980, -43.20860, -43.209, -22.900
)

dist <- pontos_sub %>% 
  dplyr::select(long_r, lat_r) %>% 
  dist() %>% 
  as.matrix()

dist %>% 
  tibble::as_tibble() %>% 
  dplyr::mutate(from = as.numeric(pontos_sub$id)) %>%
  tidyr::gather(to, dist, -from)

# A tibble: 36 × 3
    from    to       dist
   <dbl> <chr>      <dbl>
1      0     1 0.00000000
2      1     1 0.00000000
3      2     1 0.00000000
4      3     1 0.02418677
5      4     1 0.02418677
6      5     1 0.02418677
7      0     2 0.00000000
8      1     2 0.00000000
9      2     2 0.00000000
10     3     2 0.02418677
# ... with 26 more rows

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .