Skip to content

CRS for rotated pole projection works with rgdal, but not with sf #651

@adrfantini

Description

@adrfantini

The EURO-CORDEX climate models use a rotated pole projection described in this document (page 10, table 1)

I've been trying to find the correct CRS to use with SF, but I cannot find it. I tried for example: +proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25 +lon_0=180 +ellps=WGS84, which is accepted by rgdal::CRS and raster::projection, but sf::st_crs fails:

st_crs("+proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25 +lon_0=180 +ellps=WGS84")
Coordinate Reference System: NA
Warning message:
In CPL_crs_from_proj4string(x) :
  Cannot import crs from PROJ.4 string `+proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25 +lon_0=180', missing crs returned

rgdal::checkCRSArgs("+proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25 +lon_0=180  +ellps=WGS84")
[[1]]
[1] TRUE

[[2]]
[1] "+proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25 +lon_0=180 +ellps=WGS84"

Why is rgdal ok with this, but not sf?

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