Я пытаюсь найти кратчайшее расстояние между школами и береговой линией. Все школы расположены в восточном и северном формате, береговая линия состоит из точек, также в восточном и северном формате.
Я решил эту проблему, создав цикл, который проходит через каждую школу, и еще один цикл внутри школьного цикла, который сравнивает местоположение школы со всеми точками береговой линии. Это невероятно медленно, так как у меня 40 000 школ и 180 000 точек на карте, и я знаю, что вы никогда не должны использовать циклы в R! Я попытался собрать следующее вместе:
Данные испытаний:
schools <- structure(list(URN = c(100000L, 100008L, 100009L, 100010L, 100011L, 100012L), Easting = c(533498L, 530238L, 524888L, 529912L, 528706L, 528386L), Northing = c(181201L, 182761L, 185067L, 184835L, 186594L, 185209L)), .Names = c("URN", "Easting", "Northing"), row.names = c(NA, 6L), class = "data.frame")
coastline <- structure(list(Easting = c(219588.203816721, 219623.335092579, 219625.861360502, 219661.118975722, 219664.898582579, 219700.155464073 ), Northing = c(607325.869617586, 607324.434359255, 607386.276450707, 607384.83630279, 607477.377010103, 607475.937159766)), .Names = c("Easting", "Northing"), row.names = c(NA, 6L), class = "data.frame")
Код
for (sch in schools$URN){
minimumDistance <- 500000
SEasting <- schools %>% filter(URN == sch) %$% Easting
SNorthing <- schools %>% filter(URN == sch) %$% Northing
mindisance <- coastline %>% mutate(distance =
min(sqrt((SEasting - Easting)^2 +
(SNorthing - Northing)^2))) %$% distance
print(paste(sch, "minDistance = ", mindisance))
}
Но я получаю результат для каждой точки береговой линии:
[1] "100000 minDistance = 529243.315102678" "100000 minDistance = 529243.315102678"
[3] "100000 minDistance = 529243.315102678" "100000 minDistance = 529243.315102678"
[5] "100000 minDistance = 529243.315102678" "100000 minDistance = 529243.315102678"
Я бы хотел
100000 minDistance = 529243.315102678
Любая идея о том, что я делаю неправильно?