WMS - WCPS fragment to compute difference between two slices in WMS style

[Bang Pham Huu]

From rasdaman user groupd with 6D coverage demo web page: ​http://mapvue-test.s3-website-us-west-2.amazonaws.com/, the WMS style of this coverage iem_ar5_cruts31_temp_precip_wms works incorrectly for X and Y bounding box when WCPS fragment has the slicing coordinates.

# style 3

$c[period(1), scenario(0), season(0), model(0)].tas - $c[period(0), scenario(0), season(0), model(0)].tas

which can be reproduced by importing 2 tiff files to build a 3D coverage in the attachment and create a WMS style as below

$c[ansi("2020-12-30T23:54:58.500Z")] -  $c[ansi("2021-01-04T00:15:04.500Z")] + 30

http://localhost:8080/rasdaman/ows/?bbox=70586.06568358302,969097.8829550613,594874.0656835827,1493385.8829550543&crs=EPSG:3338&format=image/png&height=256&layers=test_wms&request=GetMap&service=WMS&srs=EPSG:3338&styles=subset_temperature_new1&transparent=true&version=1.3.0&width=256&time="2020-12-30T23:54:58.500Z"

The generated rasql query is

RasUtil@137: Executing rasql query: SELECT ENCODE( (  SCALE ( ( ( c0[0,0:88,0:70] - c0[1,54:66,31:49] + 30 ) ), [0:255, 0:255] ) ) , "png", "{\"nodata\":[0]}") FROM  test_wms AS c0
rasj[0] RasNetImplementation.getResponse: query failed, errNo=343, lineNo=0, colNo=0, token=Execution error: Spatial domains of the binary induce operands are incompatible.