1

与えられた座標のセット

 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

これにより、平方度（？）の寸法が残ります。だから私の質問は：これをメートル/キロメートルに変換するにはどうすればよいですか？

よろしくお願いします。

4

1 に答える 1

5

パッケージで使用してデカルトグリッドシステム（UTMゾーンなど）に変換するかspTransform、spパッケージで試しareaPolygonてくださいgeosphere。

> areaPolygon(coords)
[1] 7688.568
于 2012-09-27T15:03:04.407 に答える