Deessa questão Eu obtive a grande funçãoareaPolygon() que me dá a área dentro de um polígono de coordenadas. no entanto, quando experimento com a função, os cálculos parecem estranhos:

Eu começo criando um conjunto de dois pontos

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

então eu verifico a distância entre estes dois:

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

obtenção: 0,2827821 quilômetros (que será a diagonal do retângulo)

Eu continuo criando um retâ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,])

veja se isso parece bom:plot(polygon)

Quarta etapa: eu calculo a área dentro do polígono.

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

no entanto, eu teria esperado200*200=40000 m² desde os lados da minha lembrança são 200 por 200 metros. isso pode ser verificado via

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

Então, chegando agora à minha pergunta (finalmente) o que estou fazendo errado? muito obrigado pela sua ajuda!