Database Reference
In-Depth Information
) AS gv
FROM prism
JOIN sfpoly sf
ON ST_Intersects(prism.rast,
ST_Transform(sf.geom,
4322))
WHERE prism.month_year =
'2012-01-01'::date
)
SELECT
(gv).val,
ST_AsText((gv).geom) AS geom
FROM geoms;
The output is as follows:
val |
geom------------------+--------------------------------7.3899998664856
| POLYGON(...)7.3600001335144 |
POLYGON(...)...7.30000019073486 |
POLYGON(...)(7 rows)
Forthesakeofbrevity,theprecedingqueryresultshavebeentrimmed.Whatismost
important about the results is that thereare seven rows.
Now,
replace
the
function
with
ST_DumpAsPolygons()
ST_PixelsAsPolyons()
.
WITH geoms AS (
SELECT (
ST_PixelAsPolygons(
ST_Union(
ST_Clip(prism.rast, 2,
ST_Transform(sf.geom, 4322), TRUE)
