У меня есть шейп-файл, 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)
Число уникальных sa$id (4500) совпадает с количеством уникальных b$SMALL_AREA (4500 строк).Как (например) и id 22 отображается из sa на правильную небольшую область в b?
есть 56k строк в sa и 4500 строк в b
Любые предложения приветствуются
Я работаю в R