<div>Dear all,<br></div><div><br></div><div>my apologies in advance if I am missing something and bothering you for naught. I noticed some strange results with the Homolosine projection using the WKIDs from ESRI:<br></div><div><br></div><div>$ cat wgs.md<br></div><div>-15 15<br></div><div>175 0<br></div><div>170 35<br></div><div>160 55<br></div><div>-140 50<br></div><div>$ cs2cs +init=epsg:4326 +to ESRI:54052 wgs.md<br></div><div>-1612895.57 1669792.36 0.00<br></div><div>19480910.89 0.00 0.00<br></div><div>15501890.04 3896182.18 0.00<br></div><div>11323880.61 6049814.33 0.00<br></div><div>-10648864.61    5536706.32 0.00<br></div><div><br></div><div>These results do not match those obtained with an explicit formulation of this CRS:<br></div><div><br></div><div>$ cs2cs +init=epsg:4326 +to +proj=igh +lat_0=0 +lon_0=0 +datum=WGS84 +units=m +no_defs wgs.md<br></div><div>-1499101.98 1669792.36 0.00<br></div><div>19480910.89 0.00 0.00<br></div><div>16105847.11 3896182.18 0.00<br></div><div>12540237.72 6049791.84 0.00<br></div><div>-14174481.82    5536683.82 0.00<br></div><div><br></div><div>The differences are all with easting coordinates, making me suspect that something was amiss with lobe definitions. Therefore I went checking in the database:<br></div><div><br></div><div>spatialite> SELECT text_definition FROM projected_crs WHERE code=54052;<br></div><div>PROJCS["World_Goode_Homolosine_Land",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Goode_Homolosine"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",    0.0],PARAMETER["Option",1.0],UNIT["Meter",1.0]]<br></div><div><br></div><div>The "Goode_Homolosine" projection name refers to this class:<br></div><div><a href="https://github.com/OSGeo/PROJ/blob/46c47e9adf6376ae06afabe5d24a0016a05ced82/src/projections/goode.cpp">https://github.com/OSGeo/PROJ/blob/46c47e9adf6376ae06afabe5d24a0016a05ced82/src/projections/goode.cpp</a><br></div><div><br></div><div>It defines no lobes and does not correspond to any of the projections proposed by Goode. I believe this CRS WKT should instead refer to the actual Homolosine for land masses, identifiable with the long name "Interrupted Goode Homolosine" and the short id "igh":<br></div><div><a href="https://github.com/OSGeo/PROJ/blob/f3561e5749fefd957f67cc631378d1c39161a310/docs/source/operations/projections/igh.rst">https://github.com/OSGeo/PROJ/blob/f3561e5749fefd957f67cc631378d1c39161a310/docs/source/operations/projections/igh.rst</a><br></div><div><br></div><div>Also inquiring the database, one can note the Oceanic version also referring to the rogue lobe-less projection:<br></div><div><br></div><div>spatialite> SELECT text_definition FROM projected_crs WHERE code=54053;<br></div><div>PROJCS["World_Goode_Homolosine_Ocean",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",637137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Goode_Homolosine"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian"    ,-160.0],PARAMETER["Option",2.0],UNIT["Meter",1.0]]<br></div><div><br></div><div>In this case the right projection should possibly be the one named "Oceanic view":<br></div><div><a href="https://github.com/OSGeo/PROJ/blob/2e5470387df8c713af18e601c0e6a4b352294556/docs/source/operations/projections/igh_o.rst">https://github.com/OSGeo/PROJ/blob/2e5470387df8c713af18e601c0e6a4b352294556/docs/source/operations/projections/igh_o.rst</a><br></div><div><br></div><div>I hope this makes sense. Thank you.<br></div><div><br></div><div>--<br></div><div>Luís<br></div>