[PROJ] GOES Geostationary Projection

[Kurt Schwehr]

Thread on geotools-devel at lists.sourceforge.net here
https://sourceforge.net/p/geotools/mailman/message/36704682/

On Thu, Jun 27, 2019 at 10:08 AM Kurt Schwehr <schwehr at gmail.com> wrote:

> 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
>


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