我目前正在尝试检测数据中的疾病簇,并且需要提取每个多边形的质心。我跑步的时候
Shape.dat <- as.data.frame(coordinates(FSAShape)) # extract centroids
但是,它返回不包含多边形名称的形心列表。然后我需要将质心与疾病数据合并,所以我需要知道哪个多边形是哪个多边形,因此需要在输出中保留多边形的名称,它目前将所有多边形名称重新标记为数字,而不是其实际名称名称。我担心质心会与我“独特”运行时的质心顺序不同!
请帮助! 谢谢一百万!
这里是我的多边形的称呼,下面是我得到的输出。我的文件名为“ FSAShape”
R0G R0H R0J R0K R0L R0M R1B R1N R2C R2E R2G R2H R2J R2K
R2L R2M R2N R2P R2R R2V R2W R2X R2Y R3A R3B R3C R3E R3G
V1 V2
0 -98.08717 49.43897
1 -98.51542 50.34608
2 -100.13448 50.58424
3 -99.73617 49.57764
4 -100.29010 52.00918
5 -100.92629 49.71098
class : SpatialPolygonsDataFrame
features : 6
extent : -101.672, -97.09419, 49, 53.01323 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0
variables : 6
names : FSA, PRUID, PRNAME, X, Cases, Pop
min values : R0G, 46, Manitoba, 597, 0, 2
max values : R0M, 46, Manitoba, 598, 0, 5