ST_HillShade — 使用提供的方位角、高度、亮度及比例輸入,傳回高程網格圖層的假設照明。
raster ST_HillShade(
raster rast, integer band=1, text pixeltype=32BF, double precision azimuth=315, double precision altitude=45, double precision max_bright=255, double precision scale=1.0, boolean interpolate_nodata=FALSE)
;
raster ST_HillShade(
raster rast, integer band, raster customextent, text pixeltype=32BF, double precision azimuth=315, double precision altitude=45, double precision max_bright=255, double precision scale=1.0, boolean interpolate_nodata=FALSE)
;
使用方位角、高度、亮度及比例輸入,傳回高程網格圖層的假設照明。使用地圖代數並將山陰方程套用至相鄰的像素。傳回的像素值介於 0 到 255 之間。
azimuth
是介於 0 到 360 度之間的值,以順時針方向從北方測量。
altitude
是介於 0 到 90 度之間的值,其中 0 度在地平線上,而 90 度在正上方。
max_bright
是介於 0 到 255 之間的值,其中 0 代表無亮度,而 255 代表最大亮度。
scale
是垂直單位與水平單位的比率。對於英尺:經緯度,請使用 scale=370400;對於公尺:經緯度,請使用 scale=111120。
如果 interpolate_nodata
為 TRUE,則在計算山陰照明之前,將使用 ST_InvDistWeight4ma 來內插輸入網格圖層中 NODATA 像素的值。
![]() |
|
有關山陰的更多資訊,請參閱山陰運作方式。 |
可用性:2.0.0
增強功能:2.1.0 使用 ST_MapAlgebra() 並新增了可選的 interpolate_nodata
函式參數
變更:2.1.0 在先前的版本中,方位角和高度是以弧度表示。現在,方位角和高度是以度數表示
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_Hillshade(rast, 1, '32BF')) FROM foo st_dumpvalues ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ ----------------------------------------------------------------------- (1,"{{NULL,NULL,NULL,NULL,NULL},{NULL,251.32763671875,220.749786376953,147.224319458008,NULL},{NULL,220.749786376953,180.312225341797,67.7497863769531,NULL},{NULL,147.224319458008 ,67.7497863769531,43.1210060119629,NULL},{NULL,NULL,NULL,NULL,NULL}}") (1 row)
完整的覆蓋範圍圖塊範例。此查詢僅適用於 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_Hillshade(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;