<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>
</p>
<br>
<div class="moz-cite-prefix">On 11/08/2017 10:01 PM, Ari Jolma
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:ae0fd2a1-6b6b-c194-5f43-8905c8130054@gmail.com">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<p>Even Rouault kirjoitti 08.11.2017 klo 17:31:
</p>
<blockquote type="cite" cite="mid:19119758.vFUb0KjAUD@even-i700">
<meta name="qrichtext" content="1">
<style type="text/css">
p, li { white-space: pre-wrap; }
</style>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">On mercredi 8 novembre 2017 15:50:48 CET Piero Campalani wrote:</p>
<br>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">> This order is what can be referred to in the GridFunction: +1 --> first</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">> GRID axis, +2 --> second GRID axis, etc.</p>
</blockquote>
<br>
Ok, this tells me that the axisOrder attribute in sequenceRule
determines how to consider the GridEnvelope and the order of the
offsetVectors in the Grid element. <br>
</blockquote>
<br>
the axisLabels attribute determines axis order, and this in turn
refers to the CRS definition behind the URL in srsName. <br>
SequenceRule / GridFunction is just for the internal linearization
(whose description depends, ie: interprets axis order) - and BTW
only applies to GML (and JSON and RDF), not necessarily to binary
encodings which usually determine their own internal structure.<br>
<br>
My team has analysed your report and says:<br>
<br>
<div class="moz-cite-prefix">On 11/08/2017 02:56 PM, Bang Pham Huu
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:456b654e-d758-90c5-ecf7-d07a05339b7c@jacobs-university.de">
<p>I checked the result as I think it is like this
</p>
<p>+ GeoServer is correct with <a class="moz-txt-link-freetext" href="http://www.opengis.net/def/crs/EPSG/0/2393">http://www.opengis.net/def/crs/EPSG/0/2393</a> as this one is X, Y so the offset vector is correct.
</p>
<p>offsetVector1: X (positive), Y (negative)
</p>
<p>offsetVector2: Y (negative), X (positive)
</p>
<p>+ Rasdaman is also correct with <a class="moz-txt-link-freetext" href="http://ows.rasdaman.org/def/crs/EPSG/0/4326">http://ows.rasdaman.org/def/crs/EPSG/0/4326</a> (Lat, Long) so the offset vector is correct.</p>
<p>offsetVector1: Y (negative), X (positive)
</p>
offsetVector2: X (positive), Y (negative) <br>
<br>
and with a CRS with X, Y order like <a
class="moz-txt-link-freetext"
href="http://www.opengis.net/def/crs/EPSG/0/3857">http://www.opengis.net/def/crs/EPSG/0/3857</a>,
Rasdaman returns offset vector with correct order, e.g from a 2D
coverage imported with this CRS:<br>
<br>
<span class="tag" style="box-sizing: border-box; color: rgb(0, 0,
136);"><offsetVector</span><span class="pln"
style="box-sizing: border-box; color: rgb(0, 0, 0);"> </span><span
class="atn" style="box-sizing: border-box; color: rgb(102, 0,
102);">srsName</span><span class="pun" style="box-sizing:
border-box; color: rgb(102, 102, 0);">=</span><span class="atv"
style="box-sizing: border-box; color: rgb(0, 136, 0);"><a
class="moz-txt-link-rfc2396E"
href="http://localhost:8080/def/crs/EPSG/0/3857">"http://localhost:8080/def/crs/EPSG/0/3857"</a></span><span
class="tag" style="box-sizing: border-box; color: rgb(0, 0,
136);">></span><span class="pln" style="box-sizing:
border-box; color: rgb(0, 0, 0);">23613.2868316 0</span><span
class="tag" style="box-sizing: border-box; color: rgb(0, 0,
136);"></offsetVector></span><br>
<span class="pln" style="box-sizing: border-box; color: rgb(0, 0,
0);"></span><span class="tag" style="box-sizing: border-box;
color: rgb(0, 0, 136);"><offsetVector</span><span class="pln"
style="box-sizing: border-box; color: rgb(0, 0, 0);"> </span><span
class="atn" style="box-sizing: border-box; color: rgb(102, 0,
102);">srsName</span><span class="pun" style="box-sizing:
border-box; color: rgb(102, 102, 0);">=</span><span class="atv"
style="box-sizing: border-box; color: rgb(0, 136, 0);"><a
class="moz-txt-link-rfc2396E"
href="http://localhost:8080/def/crs/EPSG/0/3857">"http://localhost:8080/def/crs/EPSG/0/3857"</a></span><span
class="tag" style="box-sizing: border-box; color: rgb(0, 0,
136);">></span><span class="pln" style="box-sizing:
border-box; color: rgb(0, 0, 0);">0 -23613.2868316</span><span
class="tag" style="box-sizing: border-box; color: rgb(0, 0,
136);"></offsetVector></span><br>
<br>
+ Mapserver is not correct for <a class="moz-txt-link-freetext"
href="http://www.opengis.net/def/crs/EPSG/0/4326">http://www.opengis.net/def/crs/EPSG/0/4326</a>
(it should be as same as Rasdaman).<br>
<br>
<br>
</blockquote>
<br>
<br>
-Peter<br>
<br>
<blockquote type="cite"
cite="mid:ae0fd2a1-6b6b-c194-5f43-8905c8130054@gmail.com"> <br>
<blockquote type="cite" cite="mid:19119758.vFUb0KjAUD@even-i700"><br>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; "> And thus MapServer output at
</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><a class="moz-txt-link-freetext" href="http://194.66.252.155/cgi-bin/BGS_EMODnet_bathymetry/ows?SERVICE=WCS&REQUEST=DescribeCoverage&version=2.0.1&coverageid=BGS_EMODNET_CentralMed-MCol" moz-do-not-send="true">http://194.66.252.155/cgi-bin/BGS_EMODnet_bathymetry/ows?SERVICE=WCS&REQUEST=DescribeCoverage&version=2.0.1&coverageid=BGS_EMODNET_CentralMed-MCol</a></p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">could probably be correct, but it would be better for RectifiedGrid.axisLabels to be changed to "long lat" to better</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">reflect what is done.</p>
</blockquote>
<br>
At least it fixes the driver for MapServer in this case (no hack
options needed) since it does not define the axisOrder.<br>
<br>
<blockquote type="cite" cite="mid:19119758.vFUb0KjAUD@even-i700">
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">
</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">And GeoServer output at <a class="moz-txt-link-freetext" href="https://msp.smartsea.fmi.fi/geoserver/wcs?SERVICE=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=smartsea__eusm2016-EPSG2393" moz-do-not-send="true">https://msp.smartsea.fmi.fi/geoserver/wcs?SERVICE=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=smartsea__eusm2016-EPSG2393</a></p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">would either need to remove the GridFunction (ala MapServer) or keep it and invert the order in which its offsetVector appear (ala Rasdaman)</p>
</blockquote>
<br>
Invert the offsetVector order *and* invert the axis in
GridEnvelope (otherwise the size is wrong).<br>
<br>
<blockquote type="cite" cite="mid:19119758.vFUb0KjAUD@even-i700">
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; "> </p>
<p style="-qt-paragraph-type:empty; margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; "> </p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">And for WCS subsetting, when you write something like SUBSET=AXIS_NAME(min,max),</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">where AXIS_NAME should come from ? From the RectifiedGrid.axisLabels I guess ?</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">(to be opposed to CoverageDescription.boundedBy.Envelope@axisLabels)</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;">Then in Ari's test with GeoServer</p>
<p style=" margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; -qt-block-indent:0; text-indent:0px; -qt-user-state:0;"><a class="moz-txt-link-freetext" href="https://msp.smartsea.fmi.fi/geoserver/wcs?SERVICE=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=smartsea__eusm2016-EPSG2393" moz-do-not-send="true">https://msp.smartsea.fmi.fi/geoserver/wcs?SERVICE=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=smartsea__eusm2016-EPSG2393</a>, in theory "i" and "j" should be the axis requested ?</p>
</blockquote>
<br>
That does not work -- it's good that they are different since that
shows it is not ok (invalid axis label error). Maybe the logic is
related to the axis names? It's hard to tell since I don't find
the logic in the GeoServer source code.<br>
<br>
Ari<br>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></pre>
</blockquote>
<br>
<pre class="moz-signature" cols="80">--
Dr. Peter Baumann
- Professor of Computer Science, Jacobs University Bremen
<a class="moz-txt-link-abbreviated" href="http://www.faculty.jacobs-university.de/pbaumann">www.faculty.jacobs-university.de/pbaumann</a>
mail: <a class="moz-txt-link-abbreviated" href="mailto:p.baumann@jacobs-university.de">p.baumann@jacobs-university.de</a>
tel: +49-421-200-3178, fax: +49-421-200-493178
- Executive Director, rasdaman GmbH Bremen (HRB 26793)
<a class="moz-txt-link-abbreviated" href="http://www.rasdaman.com">www.rasdaman.com</a>, mail: <a class="moz-txt-link-abbreviated" href="mailto:baumann@rasdaman.com">baumann@rasdaman.com</a>
tel: 0800-rasdaman, fax: 0800-rasdafax, mobile: +49-173-5837882
"Si forte in alienas manus oberraverit hec peregrina epistola incertis ventis dimissa, sed Deo commendata, precamur ut ei reddatur cui soli destinata, nec preripiat quisquam non sibi parata." (mail disclaimer, AD 1083)
</pre>
</body>
</html>