Petascope return wrong domainSet

[barboni@…]

In my opinion, with a wcs2 request, petascope return an invalid domainSet for the temporal axis

---

The collection sdom:
[-159:160,-80:80,1200919:1202376]

---

A wcs2 query with subset T larger than raster sdom:
&subsetT=t,CRS:1(1199944,1203840)

---

The <boundedBy> Tag in the wcs2 response(this is CORRECT !!!):

<boundedBy>

<Envelope srsName="​http://www.opengis.net/def/crs/EPSG/0/4326" axisLabels="x y t" uomLabels="" srsDimension="3">

<lowerCorner>53.553381 32.217016 1200919.0</lowerCorner>
<upperCorner>53.553381 32.217016 1202376.0</upperCorner>

</Envelope>

</boundedBy>

---

The <domainSet> Tag in the wcs2 response(for me this is WRONG):

<domainSet>

<RectifiedGrid dimension="3" gml:id="MACC_OMAOD550_4326_1125-grid">

<limits>

<GridEnvelope>

<low>47 -29 1199944</low>
<high>47 -29 1203840</high>

</GridEnvelope>

</limits>
<axisLabels>x y t</axisLabels>

<gml:origin>

<gml:Point gml:id="MACC_OMAOD550_4326_1125-origin" srsName="http://www.opengis.net/def/crs/EPSG/0/4326">

<gml:pos>53.553381 32.217016 1200919.0</gml:pos>

</gml:Point>

</gml:origin>
<gml:offsetVector srsName="http://www.opengis.net/def/crs/EPSG/0/4326">1.125 0 0</gml:offsetVector>
<gml:offsetVector srsName="http://www.opengis.net/def/crs/EPSG/0/4326">0 1.125 0</gml:offsetVector>
<gml:offsetVector srsName="http://www.opengis.net/def/crs/EPSG/0/4326">0 0 1.0</gml:offsetVector>
</RectifiedGrid>

</domainSet>

---

In my opinion, in <GridEnvelope>, the expected output for the t axis is [1200919:1202376] but petascope return [1199944:1203840] that is exactly the requested subset.

My assumption is based on the old versions of petascope that filled the <GridEnvelope> tag in this way.

Ciao

Damiano

# rpm -qa | grep rasdaman
rasdaman-petascope-8.4.0-1.el6.noarch
rasdaman-8.4.0-1.el6.x86_64

[abeccati]

Thanks for reporting Damiano, indeed it looks like a wrong domainset bounds.

[Piero Campalani]

Hi there,
a patch is submitted for this error:

Name	Piero Campalani <cmppri@unife.it>
Subject	[PATCH] Fix wrong grid envelope when axis is subset via grid coordinates (ticket #341).
Date	Mon, 8 Apr 2013 18:54:06 +0200

I would leave this ticket open anyway, because there is still an issue (which I believe is there from year 0) when subsetting via CRS:1 (grid) coordinates: the grid coordinates are not back-transformed to domain coordinates so the the domainSet is ok while the boundedBy gets broken.

Es.:

http://localhost:8080/petascope/wcs2?
   version=2.0.0&
   service=WCS&
   request=GetCoverage&
   coverageid=eobstest&
   subset=x(0,3111)&
   subset=y,CRS:1(10,15)&

Returns:

[...]
  <boundedBy>
    <Envelope (...)>
      <lowerCorner>0.0 25.0 **10.0**</lowerCorner> WRONG Y BOUND
      <upperCorner>5.0 75.0 **15.0**</upperCorner> WRONG Y BOUND
    </Envelope>
  </boundedBy>

  <domainSet>
    <RectifiedGrid (...)>
      <limits>
        <GridEnvelope>
          <low>0 0 10</low>
          <high>5 100 15</high>
        </GridEnvelope>
      </limits>
    [...]
    </domainSet>
[...]

[abeccati]

I think if this other aspect does not get fixed in the patch its probably better to factor it out to a separated ticket, unless we can have also the second patch in short time. So we can release the patch version 8.4.1 solving the specific user need (and closing this ticket accordingly), while deferring the other issue.

Alan

[abeccati]

Fixed by Piero, giving now the expected output for the t axis of [1200919:1202376]

[Piero Campalani]

Fixed in changeset:eff4fa8.