TransWikia.com

Where do I find CRS axis names?

Geographic Information Systems Asked on May 1, 2021

When using the SUBSETTINGCRS argument of a WCS query, one needs to provide the axis labels in the subsetting CRS.

This answer to GeoServer does not work WCS parameters "subsettingcrs" and "outputcrs" states that

The axis names can be checked from http://www.opengis.net/def/crs/EPSG/0/3857 and they appear to be

Cartesian 2D CS. Axes: easting, northing (X,Y). Orientations: east, north. UoM: m.

However I don’t see this mention or any mention related to axis on the provided link, for this particular CRS or others.

How do I check what are the axis names for a particular CRS ?

One Answer

It seems that the format of the XML report has been changed and it is now more complicated for beginners.

If you just want to see the axis order with your eyes it is easiest to make a text search with "3857" at https://epsg.org. The result is one hit that leads to

https://epsg.org/crs_3857/WGS-84-Pseudo-Mercator.html?

The Coordinate System Details section gives the desired information:

Cartesian 2D CS. Axes: easting, northing (X,Y). Orientations: east, north. UoM: m.

When using the original link or the URL of the new REST API of epsg.org user must know how the system is working.

  1. EPSG:3857 is a Projected Coordinate System
      <gml:ProjectedCRS xmlns:xlink="http://www.w3.org/1999/xlink" 
      xmlns:epsg="urn:x-ogp:spec:schema-xsd:EPSG:2.2:dataset" 
      xmlns:gml="http://www.opengis.net/gml/3.2" gml:id="epsg-crs-3857">
  1. From the bottom of the document you can find that 3857 is based on a Cartesian Coordinate System with EPSG code 4449
<gml:cartesianCS xlink:href="https://apps.epsg.org/api/v1/CoordSystem/4499/export?format=gml"/>
  1. The provided link https://apps.epsg.org/api/v1/CoordSystem/4499/export?format=gml gives details of this Cartesian system, includina axis names, axis abbreviations, and axis order.
<gml:CartesianCS xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:epsg="urn:x-ogp:spec:schema-xsd:EPSG:2.2:dataset" xmlns:gml="http://www.opengis.net/gml/3.2" gml:id="epsg-cs-4499">
<gml:metaDataProperty>
<epsg:CommonMetaData>
<epsg:type>cartesian</epsg:type>
<epsg:informationSource>OGP</epsg:informationSource>
<epsg:revisionDate>2001-04-29</epsg:revisionDate>
<epsg:changes/>
<epsg:isDeprecated>false</epsg:isDeprecated>
</epsg:CommonMetaData>
</gml:metaDataProperty>
<gml:identifier codeSpace="EPSG">4499</gml:identifier>
<gml:name>Cartesian 2D CS. Axes: easting, northing (X,Y). Orientations: east, north. UoM: m.</gml:name>
<gml:remarks>Used in projected and engineering coordinate reference systems.</gml:remarks>
<gml:axis>
<gml:CoordinateSystemAxis gml:id="epsg-axis-41" uom="9001">
<gml:descriptionReference xlink:href="https://apps.epsg.org/api/v1/CoordinateAxisName/9906/export?format=gml"/>
<gml:identifier codeSpace="EPSG">41</gml:identifier>
<gml:axisAbbrev>X</gml:axisAbbrev>
<gml:axisDirection codeSpace="EPSG">east</gml:axisDirection>
</gml:CoordinateSystemAxis>
</gml:axis>
<gml:axis>
<gml:CoordinateSystemAxis gml:id="epsg-axis-42" uom="9001">
<gml:descriptionReference xlink:href="https://apps.epsg.org/api/v1/CoordinateAxisName/9907/export?format=gml"/>
<gml:identifier codeSpace="EPSG">42</gml:identifier>
<gml:axisAbbrev>Y</gml:axisAbbrev>
<gml:axisDirection codeSpace="EPSG">north</gml:axisDirection>
</gml:CoordinateSystemAxis>
</gml:axis>
</gml:CartesianCS>

It is also possible to see the axis names with one request into the EPSG API https://apps.epsg.org/api/swagger/ui/index by using WKT as outputformat.

https://apps.epsg.org/api/v1/CoordRefSystem/3857/export/?format=WKT

The response contains these strings

...AXIS["Easting (X)",east,ORDER[1]],AXIS["Northing (Y)",north,ORDER[2]]...

If you have Proj and its utility programs available on your computer you can also use projinfo for finding the axis info. It should be exactly same WKT than the one from the EPSG service but more conveniently formatted. GDAL utility gdalsrsinfo is also doing about the same thing.

projinfo epsg:3857 -o WKT2:2019
WKT2:2019 string:
PROJCRS["WGS 84 / Pseudo-Mercator",
    BASEGEOGCRS["WGS 84",
        DATUM["World Geodetic System 1984",
            ELLIPSOID["WGS 84",6378137,298.257223563,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4326]],
    CONVERSION["Popular Visualisation Pseudo-Mercator",
        METHOD["Popular Visualisation Pseudo Mercator",
            ID["EPSG",1024]],
        PARAMETER["Latitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["unknown"],
        AREA["World - 85┬░S to 85┬░N"],
        BBOX[-85.06,-180,85.06,180]],
    ID["EPSG",3857]]

Correct answer by user30184 on May 1, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP