使用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