Petascope_WCS_Check subsettingCrs functionality

[Bang Pham Huu]

From the CRS Extension OGC

subsettingCrs: CRS Identifier indicating the CRS in which the request subsetting coordinates are expressed.

However, it seems it does not work when requested subsetting is not native CRS (e.g: request in 3857 and native coverage is 4326)

http://localhost:8080/rasdaman/ows?service=WCS&version=2.0.1&
request=GetCoverage&coverageId=test_eobstest1&
subset=Lat(2999080.943,3503549.844)&subset=Long(2894306.761,3339584.724)&
subset=t("1950-01-01")
&subsettingCrs=http://www.opengis.net/def/crs/EPSG/0/3857&format=image/tiff

and return error

Lat axis: subset (2999080.943:3503549.844) is out of bounds.

This works with WCPS

http://localhost:8080/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=
for c in (test_eobstest1) return encode(c[t:"CRS:1"(0),
 Long:"http://www.opengis.net/def/crs/EPSG/0/3857"(2894306.761:3339584.724), 
Lat:"http://www.opengis.net/def/crs/EPSG/0/3857"(2999080.943:3503549.844)] , 
"tiff")

The WCS query with subsetting in 3857 is equivalent with this subsetting 4326 query

http://localhost:8080/rasdaman/ows?service=WCS&version=2.0.1&request=GetCoverage&coverageId=test_eobstest1&
subset=Lat(26,30)&subset=Long(26,30)
&subset=t:"CRS:1"(0)
&outputCrs=http://www.opengis.net/def/crs/EPSG/0/3857&format=image/tiff

[Bang Pham Huu]

Merged with Alex's example

 See this request
 {{{
 http://138.201.18.85/tb12/rasdaman/ows?request=GetCoverage&service=WCS&version=2.0.1&coverageId=SFTest&subset=E,http://www.opengis.net/def/crs/EPSG/0/3857(-1.3637472939075228E7,-1.3636585328807762E7)&subset=N,http://www.opengis.net/def/crs/EPSG/0/3857(4540876.150377409,4541515.136893244)&format=image/tiff&outputCrs=http://www.opengis.net/def/crs/EPSG/0/3857
 }}
 on the testbed server that fails because it tries to use the subsets given
 in EPSG:3857 as if they were in EPSG:3210
}}}