сопоставить правильный идентификатор уникальному идентификатору в координатах файла shp - PullRequest
0 голосов
/ 24 мая 2019

У меня есть шейп-файл, http://census.cso.ie/censusasp/saps/boundaries/Census2011_Small_Areas_generalised20m.zip

и я хочу извлечь long / lat, но я не уверен, как сопоставить правильную координату с правильной маленькой областью.

mycodeis:

require(ggplot2)
require(proj4)
require(rgdal)

a=readOGR(....shp)
dublin = a[a$NUS3NAME=='Dublin',]
dublin=spTransform(dublin,CRS('=proj=longlat +ellps=WGS84 +datum=WGS84'))
b=data.frame(dublin)

sa=fortify(dublin,SA='SMALL_AREA')
pj=project(sa[,1:2],proj4string(dublin),inverse=TRUE)
latlon=data.frame(latdeg=pj$y,londeg=pj$x)
sa=data.frame(cbind(latlon,sa)

enter image description here

Число уникальных sa$id (4500) совпадает с количеством уникальных b$SMALL_AREA (4500 строк).Как (например) и id 22 отображается из sa на правильную небольшую область в b?

есть 56k строк в sa и 4500 строк в b

Любые предложения приветствуются

Я работаю в R

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...