[PROJ] GOES Geostationary Projection

[Even Rouault]

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