Vondiese Frage Ich habe die großartige Funktion erhaltenareaPolygon() Das gibt mir die Fläche innerhalb eines Polygons von Koordinaten. Wenn ich jedoch mit der Funktion experimentiere, erscheinen die Berechnungen seltsam:

Zunächst erstelle ich einen Satz von zwei Punkten

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

dann überprüfe ich den Abstand zwischen diesen beiden:

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

Erhalten: 0,2827821 Kilometer (das wird die Diagonale des Rechtecks ​​sein)

Ich mache weiter mit dem Erstellen eines Rechtecks

 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,])

Mal sehen, ob das gut aussieht:plot(polygon)

Vierter Schritt: Ich berechne die Fläche innerhalb des Polygons.

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

Ich hätte es jedoch erwartet200*200=40000 m² da sind die seiten meines recangle 200 mal 200 meter. Dies kann über überprüft werden

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

Kommen wir nun zu meiner Frage (endlich), was mache ich falsch? Vielen Dank für Ihre Hilfe!