FIX - WCPS creates multi dimensional coverage constructor properly

[Bang Pham Huu]

for example, this coverage constructor below should be 2D not 1D

coverage lonLatCoverage
                over $i i(imageCrsDomain($volume_2d, Lon))
                values (
                          coverage latCoverage
                          over $px x(0:1) 
                          value list <4.5846843086; 4.5846848184>
                       )

it comes from this big query on 2D coverages

for $c in WorldDEMNeo_DSM,
    $d in WorldDEM_DSM

let $cutOut := [Long(-62.633:-62.623), Lat(-8.75:-8.7499)],
    $neo := $c[ $cutOut ],
    $dem := crsTransform($d [ $cutOut ], "EPSG:4326", 
                { cubicspline }, 
                { Lat: domain($c, Lat).resolution, Lon: domain($c, Lon).resolution },
                { domain($neo) }
            ),

    $positive_diff := ($neo - $dem) * ( ($neo - $dem) > 0),
    $volume_2d := ($positive_diff * 5),

    $lat_res := coverage lonLatCoverage
                over $i i(imageCrsDomain($volume_2d, Lon))
                values (
                          coverage latCoverage
                          over $px x(0:1) 
                          value list <4.5846843086; 4.5846848184>
                       )
    
return  
  encode(
          add(
              $volume_2d
              *
              $lat_res
             )

    , "json")