Понимание вывода `areaPolygon ()` из пакета геосфера
Отэтот вопрос Я получил отличную функциюareaPolygon()
что дает мне площадь в пределах многоугольника координат. однако, когда я экспериментирую с функцией, вычисления кажутся странными:
I start by creating a set of two points
require(fields)
coords <- c(11.3697193956209, 47.233380520521, 11.3723606043791,
47.235179479479)
coords <- matrix(coords, nrow=2, ncol=2, byrow=TRUE)
then I check the distance between these two:
rdist.earth(coords,coords,miles=FALSE)[1,2]
получение: 0,2827821 километров (что будет диагональю прямоугольника)
I go on with creating a rectangle
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,])
посмотрим, хорошо ли это выглядит:plot(polygon)
Fourth step: I calculate the area within the polygon.
geosphere::areaPolygon(polygon)
[1] 31.99288 #from the help file I know this ought to be square metres.
однако я бы ожидал200*200=40000 m²
так как стороны моего переплета 200 на 200 метров. это можно проверить через
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
так что теперь перейду к моему вопросу (наконец), что я делаю не так? Большое спасибо за Вашу помощь!