Преобразование области в пределах многоугольника координат в метры
Заданный набор координат
lat <- c(47.2325618, 47.2328269, 47.2330041, 47.2330481, 47.2330914,
47.2331172, 47.2331291, 47.2331499)
lon <- c(11.3707441, 11.3707791, 11.3708087, 11.3708031, 11.3707818,
11.3707337, 11.3706588, 11.370284)
coords <- cbind(lon,lat)
Я хочу рассчитать площадь многоугольника. Я использую функцию areapl()
из пакета splancs:
library(splancs)
areapl(coords)
# [1] 1.4768e-07
это оставляет меня с измерением в градусах в квадрате (?). поэтому мой вопрос: как я могу преобразовать это в метры / километры?
Заранее большое спасибо.
1 ответ
Решение
Либо преобразовать в декартову сеточную систему (например, зону UTM), используя spTransform
в sp
пакет или попробуйте areaPolygon
в geosphere
пакет.
> areaPolygon(coords)
[1] 7688.568