在 R 中的地图上绘制英国邮政编码
Plotting UK postcodes on a map in R
请帮助R新手
我有一份英国邮政编码列表 - 实际上有 300 万条观察结果。使用 R 将它们绘制到地图上的最佳方法是什么?
谢谢
第 1 步 - 下载英国邮政编码地理信息:例如
http://www.doogal.co.uk/UKPostcodes.php
https://data.gov.uk/search?q=postcodes
http://www.freemaptools.com/download-uk-postcode-lat-lng.htm
您可以选择您喜欢的格式。 CSV 易于使用。
将此 table 导入 R.
第 2 步 - 使用列表创建子集
#You have now a data.frame Df_UK containing geoinfo.
#Your initial list is in Df_JVT with variable PostCodes.
list <- as.list(unique(Df_JVT$PostCodes))
#Select your postcodes from Df_UK and choose variable to display on the map
datamap <- subset(Df_UK, Df_UK$POSTNR %in% list, select= c("POSTNR","CITY", "COUNTY", "LAT", "LON"))
row.names(datamap) <- 1:nrow(datamap)
步骤 3 - 创建空间对象并绘制地图
#Transform data.frame in spatial object
require(rgdal)
require(sp)
require(plotGoogleMaps)
datamap_mat<- cbind(datamap$LON,datamap$LAT)
row.names(datamap_mat) <- 1:nrow(datamap_mat)
AACRS <- CRS("+proj=longlat +ellps=WGS84")
UK_Map <- SpatialPointsDataFrame(datamap_mat, datamap, proj4string = AACRS, match.ID = TRUE)
#Map Points on Googlemaps
m <- plotGoogleMaps(UK_Map , filename='MAP_UK.html')
请帮助R新手
我有一份英国邮政编码列表 - 实际上有 300 万条观察结果。使用 R 将它们绘制到地图上的最佳方法是什么?
谢谢
第 1 步 - 下载英国邮政编码地理信息:例如 http://www.doogal.co.uk/UKPostcodes.php
https://data.gov.uk/search?q=postcodes
http://www.freemaptools.com/download-uk-postcode-lat-lng.htm
您可以选择您喜欢的格式。 CSV 易于使用。 将此 table 导入 R.
第 2 步 - 使用列表创建子集
#You have now a data.frame Df_UK containing geoinfo.
#Your initial list is in Df_JVT with variable PostCodes.
list <- as.list(unique(Df_JVT$PostCodes))
#Select your postcodes from Df_UK and choose variable to display on the map
datamap <- subset(Df_UK, Df_UK$POSTNR %in% list, select= c("POSTNR","CITY", "COUNTY", "LAT", "LON"))
row.names(datamap) <- 1:nrow(datamap)
步骤 3 - 创建空间对象并绘制地图
#Transform data.frame in spatial object
require(rgdal)
require(sp)
require(plotGoogleMaps)
datamap_mat<- cbind(datamap$LON,datamap$LAT)
row.names(datamap_mat) <- 1:nrow(datamap_mat)
AACRS <- CRS("+proj=longlat +ellps=WGS84")
UK_Map <- SpatialPointsDataFrame(datamap_mat, datamap, proj4string = AACRS, match.ID = TRUE)
#Map Points on Googlemaps
m <- plotGoogleMaps(UK_Map , filename='MAP_UK.html')