Zto pytanie Otrzymałem świetną funkcjęareaPolygon() co daje mi obszar w obrębie wielokąta współrzędnych. jednak gdy eksperymentuję z funkcją, obliczenia wydają się dziwne:

Zaczynam od stworzenia zestawu dwóch punktów

 require(fields)
 coords <- c(11.3697193956209, 47.233380520521, 11.3723606043791, 
 47.235179479479)
 coords <- matrix(coords, nrow=2, ncol=2, byrow=TRUE)

potem sprawdzam odległość między tymi dwoma:

 rdist.earth(coords,coords,miles=FALSE)[1,2]

uzyskanie: 0,2827821 kilometrów (co będzie przekątną prostokąta)

Kontynuuję tworzenie prostokąta

 polygon <- matrix(coords, nrow=2, ncol=2)
 polygon <- rbind(polygon, polygon)
 polygon[4,2] <- polygon[1,2]
 polygon[4,1] <- polygon[2,1]
 polygon[3,2] <- polygon[2,2]
 polygon[3,1] <- polygon[1,1]
 polygon <- rbind(polygon, polygon[1,])

sprawdź, czy to wygląda dobrze:plot(polygon)

Czwarty krok: obliczam pole w obrębie wielokąta.

geosphere::areaPolygon(polygon)
[1] 31.99288 #from the help file I know this ought to be square metres.

jednak oczekiwałbym200*200=40000 m² ponieważ boki mojej pętli są 200 na 200 metrów. można to sprawdzić za pośrednictwem

rdist.earth(polygon,coords,miles=FALSE)

           [,1]         [,2]
 [1,] 0.0000000 2.827821e-01
 [2,] 0.2827821 9.504539e-05
 [3,] 0.2002671 1.996434e-01
 [4,] 0.1996501 2.002671e-01

więc teraz przychodzę do mojego pytania (w końcu) co robię źle? bardzo ci dziękuje za pomoc!