[PROJ] GOES Geostationary Projection

[Kurt Schwehr]

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.

On Thu, Jun 27, 2019 at 9:32 AM Even Rouault <even.rouault at spatialys.com>
wrote:

> On jeudi 27 juin 2019 08:43:06 CEST Kurt Schwehr wrote:
> > Thanks Even!
> >
> > Since the end target isn't geotiff and it does have WKT
>
> Which WKT are you talking about ?
> I consider the WKT1 with PROJ strings you'll find reported by GDAL on the
> netCDF or in GeoTIFF to be hack, but I don't know a standardized form of
> WKT
> for that.
>
> With GDAL 3 / PROJ 6, the WKT2 is probably a bit better (see just below).
> But
> the method name and parameter names are still my own invention (hopefully
> reasonable enough)
>
> $ gdalsrsinfo "+proj=geos +sweep=x +h=35786023"
>
> PROJ.4 : +proj=geos +sweep=x +lon_0=0 +h=35786023 +x_0=0 +y_0=0
> +datum=WGS84
> +units=m +no_defs
>
> OGC WKT2:2018 :
> PROJCRS["unknown",
>     BASEGEOGCRS["unknown",
>         DATUM["World Geodetic System 1984",
>             ELLIPSOID["WGS 84",6378137,298.257223563,
>                 LENGTHUNIT["metre",1]],
>             ID["EPSG",6326]],
>         PRIMEM["Greenwich",0,
>             ANGLEUNIT["degree",0.0174532925199433],
>             ID["EPSG",8901]]],
>     CONVERSION["unknown",
>         METHOD["Geostationary Satellite (Sweep X)"],
>         PARAMETER["Longitude of natural origin",0,
>             ANGLEUNIT["degree",0.0174532925199433],
>             ID["EPSG",8802]],
>         PARAMETER["Satellite Height",35786023,
>             LENGTHUNIT["metre",1,
>                 ID["EPSG",9001]]],
>         PARAMETER["False easting",0,
>             LENGTHUNIT["metre",1],
>             ID["EPSG",8806]],
>         PARAMETER["False northing",0,
>             LENGTHUNIT["metre",1],
>             ID["EPSG",8807]]],
>     CS[Cartesian,2],
>         AXIS["(E)",east,
>             ORDER[1],
>             LENGTHUNIT["metre",1,
>                 ID["EPSG",9001]]],
>         AXIS["(N)",north,
>             ORDER[2],
>             LENGTHUNIT["metre",1,
>                 ID["EPSG",9001]]]]
>
>
> --
> Spatialys - Geospatial professional services
> http://www.spatialys.com
>


-- 
--
http://schwehr.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20190627/67fbf578/attachment.html>