الاسم
ST_Slope — تُعيد الانحدار (بالدرجات افتراضيًا) لنطاق رستر ارتفاعات. مفيدة لتحليل التضاريس.
المختصر
raster ST_Slope(raster rast, integer nband=1, text pixeltype=32BF, text units=DEGREES, double precision scale=1.0, boolean interpolate_nodata=FALSE);
raster ST_Slope(raster rast, integer nband, raster customextent, text pixeltype=32BF, text units=DEGREES, double precision scale=1.0, boolean interpolate_nodata=FALSE);
الوصف
تُعيد الانحدار (بالدرجات افتراضيًا) لنطاق رستر ارتفاعات. تستخدم الجبر الخرائطي وتطبّق معادلة الانحدار على البكسلات المجاورة.
يحدد units وحدات الانحدار. القيم الممكنة: RADIANS وDEGREES (الافتراضية) وPERCENT.
scale هي نسبة الوحدات الرأسية إلى الأفقية. للنسبة Feet:LatLon استخدم scale=370400، وللنسبة Meters:LatLon استخدم scale=111120.
إذا كانت interpolate_nodata تساوي TRUE، تُستكمَل قيم بكسلات NODATA من الرستر المُدخل باستخدام ST_InvDistWeight4ma قبل حساب انحدار السطح.
|
|
|
لمزيد من المعلومات عن الانحدار والاتجاه وتظليل التضاريس، راجع ESRI - كيف يعمل تظليل التضاريس ودليل ERDAS الميداني - صور الانحدار. |
متاحة منذ: 2.0.0
حُسّنت في: 2.1.0 تستخدم ST_MapAlgebra() وأُضيفت معاملات الدالة الاختيارية units وscale وinterpolate_nodata
تغيّرت في: 2.1.0 في الإصدارات السابقة كانت القيم المُعادة بالراديان. أصبحت القيم المُعادة الآن بالدرجات افتراضيًا
أمثلة: الصيغة 1
WITH foo AS (
SELECT ST_SetValues(
ST_AddBand(ST_MakeEmptyRaster(5, 5, 0, 0, 1, -1, 0, 0, 0), 1, '32BF', 0, -9999),
1, 1, 1, ARRAY[
[1, 1, 1, 1, 1],
[1, 2, 2, 2, 1],
[1, 2, 3, 2, 1],
[1, 2, 2, 2, 1],
[1, 1, 1, 1, 1]
]::double precision[][]
) AS rast
)
SELECT
ST_DumpValues(ST_Slope(rast, 1, '32BF'))
FROM foo
st_dumpvalues
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
---------------------------------------------------------------------
(1,"{{10.0249881744385,21.5681285858154,26.5650520324707,21.5681285858154,10.0249881744385},{21.5681285858154,35.2643890380859,36.8698959350586,35.2643890380859,21.5681285858154},
{26.5650520324707,36.8698959350586,0,36.8698959350586,26.5650520324707},{21.5681285858154,35.2643890380859,36.8698959350586,35.2643890380859,21.5681285858154},{10.0249881744385,21.
5681285858154,26.5650520324707,21.5681285858154,10.0249881744385}}")
(1 row)
أمثلة: الصيغة 2
مثال كامل لبلاطات تغطية. يعمل هذا الاستعلام فقط مع PostgreSQL 9.1 أو أحدث.
WITH foo AS (
SELECT ST_Tile(
ST_SetValues(
ST_AddBand(
ST_MakeEmptyRaster(6, 6, 0, 0, 1, -1, 0, 0, 0),
1, '32BF', 0, -9999
),
1, 1, 1, ARRAY[
[1, 1, 1, 1, 1, 1],
[1, 1, 1, 1, 2, 1],
[1, 2, 2, 3, 3, 1],
[1, 1, 3, 2, 1, 1],
[1, 2, 2, 1, 2, 1],
[1, 1, 1, 1, 1, 1]
]::double precision[]
),
2, 2
) AS rast
)
SELECT
t1.rast,
ST_Slope(ST_Union(t2.rast), 1, t1.rast)
FROM foo t1
CROSS JOIN foo t2
WHERE ST_Intersects(t1.rast, t2.rast)
GROUP BY t1.rast;