我不熟悉使用R并使用DMTI Canmap数据集中的形状文件:http://www.dmtispatial.com/canmap/#top
我使用readOGR成功将形状文件导入为'SpatialLineDataFrame'。我在映射数据方面没有成功。这是我到目前为止所尝试的。注意:我已将数据框命名为“trans”。
library(maps)
library(rgdal)
1绘制经度和纬度坐标
plot(trans$LONGITUDE,trans$LATITUDE)
此命令将坐标绘制为点而不是线。我想绘制线条。
2使用ggmap绘制
library(ggmap)
can <- get_map(location="canada",zoom=3)
ggmap(can) + geom_point(data=trans, aes(x="LONGITUDE",y="LATITUDE"))
这将返回错误,指出“错误:提供给连续比例的离散值”
ggmap(can) + geom_path(data=trans, aes(x="LONGITUDE",y="LATITUDE"))
这也会返回错误“错误:提供给连续比例的离散值”
作为单挑,当我询问存储在'trans'中的数据的名称时,我得到以下名称:“UID”“NAME”“TYPE”“LOCATION”“NUM_LINES”“MUNICIPAL”“PROV”“ LONGITUDE“”LATITUDE“”BEGELEV“”ENDELEV“”ACCURACY“”ACQ_TECH“”VALDATE“”PROVIDER“”SHAPE_Leng“
您是否建议尝试使用R来映射存储在此.shp文件中的行?