使用tmap制作英国的犯罪数量地图

Question

我正在通过使用区域代码组合数据在R中按tmap制作地图。并非所有犯罪数据都可用。将区域数据和犯罪数据合并后，我将无法绘制地图

uk_la1 <- readOGR(dsn = "./infuse_dist_lyr_2011", layer = "infuse_dist_lyr_2011")

Totalcrime <- read.csv('Total_crime_in_each_area_full.csv', header = TRUE)

# Calculate the nnumber of crime per 1000 people
Totalcrime <- transform(Totalcrime, Crime_per_1000_people = Total / Population * 1000)

Totalcrimeno<-Totalcrime %>% select(geo_code, Crime_per_1000_people)

uk_la1@data<-left_join(uk_la1@data, Totalcrimeno,
                               by=c('geo_code'))

qtm(uk_la, fill="Crime_per_1000_people")

但是我得到一个错误：

  

$<-.data.frame（*tmp*中的错误，“几何”，值= list（list（list（：     替换有404行，数据有405

Answer 1

您需要使用WGS84坐标系转换经度和纬度：

library(rgdal)
library(spdplyr)
library(geojsonio)
library(rmapshaper)

uk_la <- readOGR(dsn = "./infuse_dist_lyr_2011", layer = "infuse_dist_lyr_2011")

wgs84 <- "+proj=longlat +datum=WGS84"
uk_la_trans <- spTransform(uk_la, CRS(wgs84))

#Convert from Spatial Dataframe to GeoJson 
uk_la_trans_json <- geojson_json(uk_la_trans)
uk_la_trans_sim <= ms_simplify(uk_la_trans_json)

有关WGS84 Coordinate System的更多信息