[gdal-dev] WCS GetCoverage with AxisOrder swap

jratike80 jukka.rahkonen at maanmittauslaitos.fi
Wed Nov 8 07:57:09 PST 2017


Even Rouault-2 wrote
> So, it would seem that there are 2 more or less equivalent way of
> achieving the same end result ?

My understanding is that these two are equivalent and define a grid in
EPSG:4326 with pixel size of 1 degree (if it is an image). First way is to
turn the offset vectors, this is also used in an example in the OGC GML in
JPEG 2000 (GMLJP2) standard.

offsetVector srsName="http://ows.rasdaman.org/def/crs/EPSG/0/4326">0
1</offsetVector
offsetVector srsName="http://ows.rasdaman.org/def/crs/EPSG/0/4326">-1
0</offsetVector

The second way is to turn the axes order.

<GridFunction>
  <sequenceRule axisOrder="+2 +1">Linear</sequenceRule>
</GridFunction>
offsetVector srsName="http://ows.rasdaman.org/def/crs/EPSG/0/4326">-1
0</offsetVector
offsetVector srsName="http://ows.rasdaman.org/def/crs/EPSG/0/4326">0
1</offsetVector

It is also possible to write the first variant as complete so it is not
necessary to know what are the default values of sequenceRule

<GridFunction>
  <sequenceRule axisOrder="+1 +2">Linear</sequenceRule>
</GridFunction>
offsetVector srsName="http://ows.rasdaman.org/def/crs/EPSG/0/4326">0
1</offsetVector
offsetVector srsName="http://ows.rasdaman.org/def/crs/EPSG/0/4326">-1
0</offsetVector


-Jukka Rahkonen-




--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html


More information about the gdal-dev mailing list