<div dir="ltr">The GDAL in my tool chain is circa Oct-2018, so the hack.  But since I'm writing a custom converter, I should be able to force the WKT to be what ever GeoTools wants for GEOS 16 & 17.  I'll start up a discussion on the GeoTools mailing list and link back to here.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 27, 2019 at 9:32 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On jeudi 27 juin 2019 08:43:06 CEST Kurt Schwehr wrote:<br>
> Thanks Even!<br>
> <br>
> Since the end target isn't geotiff and it does have WKT<br>
<br>
Which WKT are you talking about ?<br>
I consider the WKT1 with PROJ strings you'll find reported by GDAL on the <br>
netCDF or in GeoTIFF to be hack, but I don't know a standardized form of WKT <br>
for that. <br>
<br>
With GDAL 3 / PROJ 6, the WKT2 is probably a bit better (see just below). But <br>
the method name and parameter names are still my own invention (hopefully <br>
reasonable enough)<br>
<br>
$ gdalsrsinfo "+proj=geos +sweep=x +h=35786023"<br>
<br>
PROJ.4 : +proj=geos +sweep=x +lon_0=0 +h=35786023 +x_0=0 +y_0=0 +datum=WGS84 <br>
+units=m +no_defs<br>
<br>
OGC WKT2:2018 :<br>
PROJCRS["unknown",<br>
    BASEGEOGCRS["unknown",<br>
        DATUM["World Geodetic System 1984",<br>
            ELLIPSOID["WGS 84",6378137,298.257223563,<br>
                LENGTHUNIT["metre",1]],<br>
            ID["EPSG",6326]],<br>
        PRIMEM["Greenwich",0,<br>
            ANGLEUNIT["degree",0.0174532925199433],<br>
            ID["EPSG",8901]]],<br>
    CONVERSION["unknown",<br>
        METHOD["Geostationary Satellite (Sweep X)"],<br>
        PARAMETER["Longitude of natural origin",0,<br>
            ANGLEUNIT["degree",0.0174532925199433],<br>
            ID["EPSG",8802]],<br>
        PARAMETER["Satellite Height",35786023,<br>
            LENGTHUNIT["metre",1,<br>
                ID["EPSG",9001]]],<br>
        PARAMETER["False easting",0,<br>
            LENGTHUNIT["metre",1],<br>
            ID["EPSG",8806]],<br>
        PARAMETER["False northing",0,<br>
            LENGTHUNIT["metre",1],<br>
            ID["EPSG",8807]]],<br>
    CS[Cartesian,2],<br>
        AXIS["(E)",east,<br>
            ORDER[1],<br>
            LENGTHUNIT["metre",1,<br>
                ID["EPSG",9001]]],<br>
        AXIS["(N)",north,<br>
            ORDER[2],<br>
            LENGTHUNIT["metre",1,<br>
                ID["EPSG",9001]]]]<br>
<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">--<div><a href="http://schwehr.org" target="_blank">http://schwehr.org</a></div></div>