使用sp :: Polygon @ labpt进行意外的质心计算

Question

我最近计算了一个地图项目的质心，我注意到了一些可疑的位置。例如，日本的labpt输出远远超过其纬度范围。我的第一个问题是这里发生了什么？我的第二个问题是质心计算是否有一个很好的替代方案可以正确地获得质心。

library(sp)
library(maps)

japan <- map_data("world",region="japan")

> max(japan$lat)
[1] 45.50952
> min(japan$lat)
[1] 24.26606

max(japan$long)
[1] 145.833
min(japan$long)
[1] 123.6789

# note 63 is calculated Latitude for Japan
sp::Polygon(japan[c("long","lat")])@labpt
[1] 130.32951  63.62904

Answer 1

谷歌搜索了一下，我想出了以下内容：

library(rgeos)

rgeos::gCentroid(SpatialPoints(japan[c('long','lat')]))

SpatialPoints:
     x        y
1 135.2924 35.47393