Desdeesta pregunta Obtuve la gran funciónareaPolygon() Lo que me da el área dentro de un polígono de coordenadas. sin embargo, cuando experimento con la función, los cálculos parecen impares:

Empiezo creando un conjunto de dos puntos.

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

Luego verifico la distancia entre estos dos:

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

obteniendo: 0.2827821 kilómetros (que será la diagonal del rectángulo)

Sigo creando un rectángulo.

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

ver si esto se ve bien:plot(polygon)

Cuarto paso: calculo el área dentro del polígono.

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

sin embargo, habría esperado200*200=40000 m² ya que los lados de mi recangle son 200 por 200 metros. Esto se puede comprobar a través de

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

así que ahora que voy a mi pregunta (finalmente), ¿qué estoy haciendo mal? ¡muchas gracias por su ayuda!