[postgis-users] Help with ST_CLIP

Simon SPDBA Greener simon at spdba.com.au
Thu Mar 11 22:37:21 PST 2021


"POSTGIS=""2.4.3 r16312"" PGSQL=""100"" GEOS=""3.6.2-CAPI-1.10.2 
4d2925d6"" PROJ=""Rel. 4.9.3, 15 August 2016"" GDAL=""GDAL 2.2.3, 
released 2017/11/20"" LIBXML=""2.9.4"" LIBJSON=""0.12.1"" 
LIBPROTOBUF=""1.2.1"" RASTER"

I have a polygon that defines the boundary of the Great Barrier Reef.

A) That boundary is stored as a vector geography column (4326) in a table.
B) I have loaded a raster dataset of depths (4326) and checked its 
statistics.
     The data loaded correctly and I can see valid depth values.

Now I want to clip the data in B with the polygon in A.

(I can't use a raster function with a geography object so I cast it to 
geometry.)

create table gbr100x100clip
as
select ST_Clip(g.rast, 1, m.geog::geometry, false) as rast
   from gbr_marine_park_boundary as m,
        gbr100x100 as g;

This works but the value (depth) of the whole (slipped) rasters is NULL.

I expected to see depth values but this returns null values.

SELECT (stats).*
   FROM (SELECT ST_SummaryStats(rast, 1) As stats FROM gbr100x100clip) 
As foo
  order by max desc;

Can anyone give me an idea of what I am doing wrong?

regards
Simon



More information about the postgis-users mailing list