Custom Query (2765 matches)

For example this WCPS query:

for $co in (air_quality_monitoring_data_co)

let $sub := [ansi(153403), X(147522.21879115:351672.21879115004083), Y(2422004.7731021:2799154.7731021)],
    $bands := ( ($co[ $sub ] - min($co[ $sub ])) / (max($co[ $sub ]) - min($co[ $sub ])) ) * 255 
return encode(
         scale(
           (unsigned char) ($bands),
           { 
             X:"CRS:1"(0:600), Y:"CRS:1"(0:700)
           } 
        ),   
       "png"
)

return wrong grid domains 0:3,0:600,0:700

SELECT encode((char)  ( SCALE(c_0, [0:3,0:600,0:700]) ) , "png" , "{\"geoReference\":{\"crs\":\"EPSG:3826\",\"bbox\":{\"xmin\":147522.21879115,\"ymin\":2422004.7731021,\"xmax\":351672.21879115004083,\"ymax\":2799154.7731021}},\"nodata\":[-3.4028234663852886E38]}")  FROM air_quality_monitoring_data_co_4 AS c_0

But without LET clause, this works fine, e.g.

for $co in (air_quality_monitoring_data_co)

return encode(
         scale( $co[ansi(153403), X(147522.21879115:351672.21879115004083), Y(2422004.7731021:2799154.7731021)]
               - min($co[ansi(153403), X(147522.21879115:351672.21879115004083), Y(2422004.7731021:2799154.7731021)])
                , { 
             X:"CRS:1"(0:20), Y:"CRS:1"(0:30)
           } ),  
       "png"
)

returns rasql query:

SELECT encode(SCALE(c_0[0,0:510,0:942] -  min_cells(c_0[0,0:510,0:942]) , [0:20,0:30]), "png" , "{\"geoReference\":{\"crs\":\"EPSG:3826\",\"bbox\":{\"xmin\":147522.21879115,\"ymin\":2422004.7731020999733,\"xmax\":351672.2187911500274,\"ymax\":2799154.7731021}},\"nodata\":[-3.4028234663852886E38]}")  FROM air_quality_monitoring_data_co_8 AS c_0

In WSClient, it should show generated GET requests for the most used WCS / WMS requests.

Because CRS:31467 has XY axes order but it has NorthEast orientation, hence, this query below returns wrong result:

TM coverage has metadata:

Axis name	Type	Geo Extent	         Grid Extent	Resolution	UoM
ansi	Regular Axis	["1989-12-31T12:00:00.000Z", "2023-05-02T12:00:00.000Z"]	[-10592, 1582]	1	d
E	Regular Axis	[3280414, 3934414]	[0, 653]	1000	metre
N	Regular Axis	[5237501, 6103501]	[0, 865]	-1000	metre

for $c IN (TM)
	let
		$p := POLYGON((927819.663673534 6862280.24807332,927599.907217214 6858209.97631714,940087.806713305 6857493.37917696,939285.21791631 6866895.13365604,927533.024817465 6867506.62988232,927819.663673534 6862280.24807332))
	return encode (
		clip(
			$c[ansi("2022-08-01")], $p, "http://crs.rasdaman.com/def/crs/EPSG/0/3857"),
        "tiff"
)