我的学习区域中有一个带有等压线的shapefile,我很容易将其作为SpatialLinesDataFrame导入到R中,>
sdf.isobath <- readOGR("./gis","isobath")
sdf.isobath <- spTransform(sdf.isobath, CRS("+proj=longlat +datum=WGS84"))
proj4string(sdf.isobath)
plot(sdf.isobath)
尽管如此,我想在ggplot映射中使用。为此,我“加强”了它。
# add id column with row names sdf.isobath$id <- rownames(sdf.isobath@data) # create dataframe dat.isobath_spc <- fortify(sdf.isobath,region="id") head(dat.isobath_spc) # add map data dat.isobath_spc <- merge(dat.isobath_spc,sdf.isobath@data,by="id") head(dat.isobath_spc) summary(dat.isobath_spc)
但是当我用ggplot等值线对其进行绘制时感到困惑
ggplot() + geom_line(data = dat.isobath_spc, aes(x=long, y=lat, group=id), color="black",size=0.1)
我尝试了使用等值线的多个shapefile,并且遇到了相同的问题。我要去哪里错了?
谢谢。
我的学习区域中有一个带有等压线的shapefile,我很容易将其作为带有sdf.isobath的SpatialLinesDataFrame导入到R中]
我得到了答案。我必须使用geom_path