R: Найти/сопоставить код NUTS в шейп-файле

Question

Где я могу найти код NUTS ЕС в шейп-файле в R?

NUTS = Номенклатура территориальных единиц для статистики

Каждый NUTS регион имеет NUTS код:

.

Вы также можете найти подробные коды NUTS в this Excel file.

Когда я загружаю шейп-файл NUTS, я не могу найти этот код NUTS.

Я скачал шейп NUTS_2013_20M_SH.zip от:

ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts#nuts13

После загрузки, я загрузить шейп в R, используя код:

library(sp) 
library(rgdal) 

dir <- setwd(getwd()) 

NUTS_shape = readOGR(dsn = (dsn = "~/NUTS_2013_20M_SH/data", layer = "NUTS_BN_20M_2013") 

Примечание:

• dsn = - путь к папке, в которой вы сохранили шейп-файл.
• layer = это имя файла без расширения файла (например, .shp)

Если я использую summary(NUTS_shape), я только найти «NUTS_BN_ID», но нет NUTS кода. Также, если я использую str(NUTS_shape[4461,]), чтобы увидеть структуру списка из примера-объекта 4461, похоже, что существует код NUTS.

Кто-нибудь знает, содержит ли файл формы NUTS код NUTS или если NUTS_BN_ID каким-то образом связан с этим кодом?

Answer 1

Я нашел ответ:

В случае, если кто также работает с орехом данными, то вы найдете NUTS код не в слое «NUTS_BN_20M_2013», но вы должны использовать еще один слой, то есть: «NUTS_RG_20M_2013 ».

Следовательно, используйте код:

NUTS_shape = readOGR(dsn = (dsn = "~/NUTS_2013_20M_SH/data", layer = "NUTS_RG_20M_2013") 

Затем, если вы пишете NUTS_shape@data вы получаете:

NUTS_ID STAT_LEVL_ SHAPE_AREA SHAPE_LEN 
0  AT   0 10.04269653 22.922441 
1  AT1   1 2.84477225 10.876468 
2 AT11   2 0.47903755 5.591853 
3 AT111   3 0.08480488 1.178272 
4 AT112   3 0.21836213 2.649698 
5 AT113   3 0.17587054 2.276286 

NUTS_ID равен код.