[gdal-dev] EPSG code with vertical datum

[Even Rouault]

On jeudi 30 mai 2019 16:27:41 CEST Even Rouault wrote:
> On jeudi 30 mai 2019 10:17:35 CEST Andrew Bell wrote:
> > On Thu, May 30, 2019 at 5:06 AM Even Rouault <even.rouault at spatialys.com>
> > 
> > wrote:
> > > On mercredi 29 mai 2019 20:38:02 CEST Andrew Bell wrote:
> > > > In earlier GDAL versions I could call
> > > > 
> > > > SetFromUserInput("EPSG:4326+4326");
> > > > 
> > > > without error.  I now get the following error:
> > > > 
> > > > (Error) GDAL failure (1) SetCompoundCS() fails, vertical component is
> > > > not
> > > > VERT_CS.
> > > > 
> > > > Was this a failure to detect a problem in older GDAL versions
> > > 
> > > Yes. EPSG:4326+4326 doesn't make any sense from a geodesic point of
> > > view.
> > > If
> > > you need WGS84 as a Geographic 3D CRS, you need to use EPSG:4979
> > 
> > So is the syntax "EPSG:<horizontal crs>+<vertical datum>" supported?
> 
> Yes
> 
> > I've
> > tried, for example, "EPSG:4326+5103" and also gotten a failure.
> 
> Doesn't work either in GDAL 2.4. EPSG:5103 doesn't exist.

More exacly, EPSG:5103 is not a vertical CRS, but a vertical datum.
EPSG:XXXX+YYYY requires XXXX to be a horizontal geographic/projected CRS and 
YYYY to be a vertical CRS (not a vertical datum)

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com