Currently the scale_grid/array() function used by our WCS 2 service endpoint assumes the underlying data uses the WGS84 CRS (by default, it may take the CRS as a parameter). WGS84 presumes a certain axis order (latitude, longitude). This corresponds to (y, x). How should this be addressed by our WCS 2? Should the scale function transpose data? or should we add a transpose function for just this purpose?
From the above link:
... there are no EPSG codes that specify geographic coordinates in (lon,lat) order.
The OGC, however, did define a code for geographic coordinates in (lon,lat) order, specified in OGC document http://portal.opengeospatial.org/files/?artifact_id=30575
The full URN for the (lon,lat) ordering of WGS84 coordinates is:
Of course, whatever software you are using would have to support this. PyCSW will supports this in version 1.10.1 and higher.