0

I´ve a .gpkg file in Mollweide projection. When I open it in geopandas, it shows me the correct crs information:

df1.crs
>>>>> {'init': 'epsg:54009'}

I´d like to re-project it to EPSG 4326, in a way to cross this layer with a global map in lat/long. If I try:

df1.to_crs(epsg=4326)

I get the error message:

pyproj.exceptions.CRSError: Invalid projection: +init=epsg:54009 +type=crs

How can I proceed to reproject those layers in a way I could intersect them?

3
  • What version of pyproj are you using?
    – joris
    Mar 25, 2019 at 20:59
  • joris, the last one, I just updated it. I think it´s because EPSG 54009 is not directly supported. Apr 1, 2019 at 18:10
  • 1
    Mollweide CRS projection string is esri:54009, not epsg:54009, see. Try setting this string with the set_crs() method.
    – Alex
    Dec 7, 2020 at 10:56

2 Answers 2

2

The Mollweide CRS projection is not "EPSG:54009" but "ESRI:54009". You should try:

df1.to_crs("ESRI:54009") 
0

You can get it from the OGC Well Known Text definition:

PROJCS["Sphere_Mollweide",
    GEOGCS["Unknown datum based upon the Authalic Sphere",
        DATUM["Not_specified_based_on_Authalic_Sphere",
            SPHEROID["Sphere",6371000,0],
            AUTHORITY["EPSG","6035"]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Mollweide"],
    PARAMETER["central_meridian",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["ESRI","53009"]]

Source: https://epsg.io/53009

In .NET, i can then get the projection from the esri-string:

            // Mollweide equal-area projection 
            // DotSpatial.Projections.ProjectionInfo toMollweide = DotSpatial.Projections.ProjectionInfo.FromEpsgCode(53009); // not working 

            // https://epsg.io/53009
            DotSpatial.Projections.ProjectionInfo toMollweide = DotSpatial.Projections.ProjectionInfo.FromEsriString(@"
PROJCS[""Sphere_Mollweide"",
    GEOGCS[""Unknown datum based upon the Authalic Sphere"",
        DATUM[""Not_specified_based_on_Authalic_Sphere"",
            SPHEROID[""Sphere"", 6371000, 0],
            AUTHORITY[""EPSG"", ""6035""]],
        PRIMEM[""Greenwich"", 0],
        UNIT[""Degree"", 0.0174532925199433]],
    PROJECTION[""Mollweide""],
    PARAMETER[""central_meridian"", 0],
    PARAMETER[""false_easting"", 0],
    PARAMETER[""false_northing"", 0],
    UNIT[""metre"", 1,
        AUTHORITY[""EPSG"", ""9001""]],
    AXIS[""Easting"", EAST],
    AXIS[""Northing"", NORTH],
    AUTHORITY[""ESRI"", ""53009""]]
");

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.