В настоящее время я пытаюсь обнаружить кластеры болезней в моих данных, и мне нужно извлечь центроид каждого полигона. Когда я бегу;
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
