|
@@ -3,12 +3,10 @@
|
|
|
|
|
|
|
|
<mapper namespace="com.onemap.analyse.mapper.raster.SlopMapper">
|
|
<mapper namespace="com.onemap.analyse.mapper.raster.SlopMapper">
|
|
|
<select id="getAveSlop" resultType="Double">
|
|
<select id="getAveSlop" resultType="Double">
|
|
|
- -- 创建一个多边形
|
|
|
|
|
WITH poly AS (SELECT public.st_geomfromtext(#{geomWkt}, 4326) AS geom)
|
|
WITH poly AS (SELECT public.st_geomfromtext(#{geomWkt}, 4326) AS geom)
|
|
|
- -- 裁剪坡度图层并计算平均坡度
|
|
|
|
|
SELECT AVG((gv).val) AS avg_slope
|
|
SELECT AVG((gv).val) AS avg_slope
|
|
|
FROM (SELECT public.st_clip(rast, 1, geom, false) AS clipped_rast
|
|
FROM (SELECT public.st_clip(rast, 1, geom, false) AS clipped_rast
|
|
|
- FROM sy_slop,
|
|
|
|
|
|
|
+ FROM raster.sy_slop,
|
|
|
poly
|
|
poly
|
|
|
WHERE public.st_intersects(rast, geom)) AS clipped,
|
|
WHERE public.st_intersects(rast, geom)) AS clipped,
|
|
|
LATERAL public.st_pixelaspoints(clipped_rast) AS gv
|
|
LATERAL public.st_pixelaspoints(clipped_rast) AS gv
|
