将 geocdoed 地址组合到路线图上的困难

问题描述

抱歉，我无法让这段代码重现，但我会尽力而为。

我用这个绘制了多伦多主要干道的图表：

[![好的道路地图][1]][1]

我已经成功地从一个单独的数据集中对 10 个停车违规地址进行了地理编码，其中前 10 个地址如下所示：

"***46474","***46475","***46674","***64245","***80959","***66996","***71445","***66997"),date_of_infraction = c(20200101,20200101,20200101),infraction_code = c(3,29,3,29),infraction_description = c("PARK ON PRIVATE PROPERTY","PARK PROHIBITED TIME NO PERMIT","PARK ON PRIVATE PROPERTY","PARK PROHIBITED TIME NO PERMIT"),set_fine_amount = c(30,30,30),time_of_infraction = c("0001","0002","0004","0005","0006"),location1 = c(NA,"OPP","NR","W/S","AT","NR"),location2 = c("4700 JANE ST","573 CLENDENAN AVE","21 MARIA ST","164 GLENHOLME AVE","SHAW ST","29 ROSEBANK DR","23 CRAYDON AVE","23 CRAYDON AVE"
),location3 = c(NA,NA,"S/O",NA
),location4 = c(NA,"COLLEGE ST",NA),province = c("SK","ON","ON"),address = c("4700 JANE ST,Toronto,Ontario","573 CLENDENAN AVE,"21 MARIA ST,"164 GLENHOLME AVE,"SHAW ST,"29 ROSEBANK DR,"23 CRAYDON AVE,Ontario"
),lon = c(-79.5214055,-79.4734866,-79.47397,-79.4395997,-79.421207,-79.2348412,-79.4976524,-79.4976524),lat = c(43.7692709,43.6668852,43.6661886,43.6806354,43.6573352,43.7911591,43.6906576,43.6906576
)),row.names = c(NA,-10L),class = c("tbl_df","tbl","data.frame"
))

现在，我想我可以用这个添加停车票的积分：

toronto %>% 
#just filter in big roads
  filter(str_detect(FEATURE_00,"Arterial|Expressway")) %>% 
#Draw the roads grph
  ggplot(.)+geom_sf()+
#Add the points,this doesn't work. 
  geom_point(data=geo_test,aes(x=lon,y=lat))

但它实际上最终看起来像这样：

[![坏地图][2]][2]

所以，纬度和经度的比例似乎有些偏差，但我真的不知道如何解决这个问题。 shapefile的原始投影好像是wgs84，但不知道有没有关系。

如果您查看我通过对地址进行地理编码获得的 geo_test 中点的纬度和经度坐标，它们看起来应该适合原始地图上打印的比例尺。再次，很抱歉没有使这完全可重现；我会努力解决这个问题。 [1]：https://i.stack.imgur.com/vwxzX.png [2]：https://i.stack.imgur.com/BnG9U.png

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

ggmap ggplot2 r r