Skip to content

Possible bug in coord_sf reprojection of vector data (Natural Earth coastline to Robinson projection) #3949

@AndrewCunliffe

Description

@AndrewCunliffe

Reprojection of the small scale Natural Earth coastline dataset to a Robinson projection returns an erroneous line feature on the plot. (Apologies if this is the wrong place to report, but it appears to be an issue with coord_sf rather than the Natural Earth vector dataset).

Medium-scale coastlines reproject okay to a World Robinson projection (EPSG:54030)
coastlines <- ne_coastline(scale = "medium", returnclass = "sf")
ggplot(data = coastlines) + geom_sf(colour = "black") + coord_sf(crs = st_crs(54030)) + theme_bw()
Global map from R 1

Small-scale coastlines display okay at the default projection
coastlines <- ne_coastline(scale = "small", returnclass = "sf")
ggplot(data = coastlines) + geom_sf(colour = "black") + theme_bw()
Global map from R 3

Small-scale coastlines reprojected to the World Robinson projection (EPSG:54030) include the erroneous line feature at ca. 70 degrees north
coastlines <- ne_coastline(scale = "small", returnclass = "sf")
ggplot(data = coastlines) + geom_sf(colour = "black") + coord_sf(crs = st_crs(54030)) + theme_bw()
Global map from R 2

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions