Hello @twhiteaker
using http://www.epsg-registry.org/export.htm?wkt=urn:ogc:def:crs:EPSG::6931 as
a start point, we can adjust the datum definition and name as follows
(disclaimer, this is a draft, please check the contents)
```
PROJCRS["NSIDC EASE-Grid 1.0 North",
BASEGEODCRS["GCS_Sphere_International_1924_Authalic",
DATUM["D_Sphere_International_1924_Authalic",
ELLIPSOID["Sphere_International_1924_Authalic",6371228.0,0.0,LENGTHUNIT["metre",1.0]]],
PRIMEM["Greenwich",0.0]],
CONVERSION["US NSIDC EASE-Grid 1.0 North",
METHOD["Lambert Azimuthal Equal Area",ID["EPSG",9820]],
PARAMETER["Latitude of natural
origin",90,ANGLEUNIT["degree",0.01745329252]],
PARAMETER["Longitude of natural
origin",0,ANGLEUNIT["degree",0.01745329252]],
PARAMETER["False easting",0,LENGTHUNIT["metre",1.0]],
PARAMETER["False northing",0,LENGTHUNIT["metre",1.0]]],
CS[cartesian,2],
AXIS["easting
(X)",south,MERIDIAN[90,ANGLEUNIT["degree",0.01745329252]],ORDER[1]],
AXIS["northing
(Y)",south,MERIDIAN[180,ANGLEUNIT["degree",0.01745329252]],ORDER[2]],
LENGTHUNIT["metre",1.0]]
```
note, you can check syntax for errors with a modern version of proj
(e.g.proj-6.2.1)
```
projinfo 'PROJCRS["NSIDC EASE-Grid 1.0 North",
> BASEGEODCRS["GCS_Sphere_International_1924_Authalic",
> DATUM["D_Sphere_International_1924_Authalic",
>
> ELLIPSOID["Sphere_International_1924_Authalic",6371228.0,0.0,LENGTHUNIT["metre",1.0]]],
> PRIMEM["Greenwich",0.0]],
> CONVERSION["US NSIDC EASE-Grid 1.0 North",
> METHOD["Lambert Azimuthal Equal Area",ID["EPSG",9820]],
> PARAMETER["Latitude of natural
> origin",90,ANGLEUNIT["degree",0.01745329252]],
> PARAMETER["Longitude of natural
> origin",0,ANGLEUNIT["degree",0.01745329252]],
> PARAMETER["False easting",0,LENGTHUNIT["metre",1.0]],
> PARAMETER["False northing",0,LENGTHUNIT["metre",1.0]]],
> CS[cartesian,2],
> AXIS["easting
> (X)",south,MERIDIAN[90,ANGLEUNIT["degree",0.01745329252]],ORDER[1]],
> AXIS["northing
> (Y)",south,MERIDIAN[180,ANGLEUNIT["degree",0.01745329252]],ORDER[2]],
> LENGTHUNIT["metre",1.0]]'
PROJ.4 string:
+proj=laea +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +R=6371228 +units=m +no_defs
+type=crs
```
(though also note that projinfo also outputs it's own adapted WKT, which may
have ID entries that you do not want)
I hope this helps
marqh
--
You are receiving this because you are subscribed to this thread.
Reply to this email directly or view it on GitHub:
https://github.com/cf-convention/cf-conventions/issues/278#issuecomment-645302428
This list forwards relevant notifications from Github. It is distinct from
[email protected], although if you do nothing, a subscription to the
UCAR list will result in a subscription to this list.
To unsubscribe from this list only, send a message to
[email protected].
