Тут я демонструю невелику оптимізацію щодо рішення @whuber, і я вказую на "ширину буфера", оскільки це корисно для інтеграції рішення більш загальної проблеми: Чи є зворотна функція st_buffer, яка повертає оцінку ширини?
CREATE FUNCTION buffer_width(
-- rectangular strip mean width estimator
p_len float, -- len of the central line of g
p_geom geometry, -- g
p_btype varchar DEFAULT 'endcap=flat' -- st_buffer() parameter
) RETURNS float AS $f$
DECLARE
w_half float;
w float;
BEGIN
w_half := 0.25*ST_Area(p_geom)/p_len;
w := 0.50*ST_Area( ST_Buffer(p_geom,-w_half,p_btype) )/(p_len-2.0*w_half);
RETURN w_half+w;
END
$f$ LANGUAGE plpgsql IMMUTABLE;
Для цієї проблеми, @celenius питання про ширині вулиці , sw
, рішення
sw = buffer_width(ST_Length(g1), g2)
де sw
"середня ширина", g1
центральна лінія g2
, а вулиця g2
- ПОЛІГОН . Я використовував лише стандартну бібліотеку OGC, протестовану з PostGIS , і вирішував інші серйозні практичні програми з тією ж функцією buffer_width.
ДЕМОНСТРАЦІЯ
A2
- площа g2
, L1
довжина центральної лінії ( g1
) g2
.
Припустимо, що ми можемо генерувати g2
за допомогою g2=ST_Buffer(g1,w)
, а g1
це пряма, так g2
це прямокутник з довжиною L1
і шириною 2*w
, і
A2 = L1*(2*w) --> w = 0.5*A2/L1
Це не та сама формула @whuber, тому що тут w
половина g2
ширини прямокутника ( ). Це хороший оцінювач, але, як ми бачимо тестами (нижче), не є точним, і функція використовує його як підказку, для зменшення g2
площі та як остаточний оцінювач.
Тут ми не оцінюємо буфери з "endcap = square" або "endcap = round", яким потрібна сума до A2
області точкового буфера з однаковою w
.
ДОВІДКИ: на подібному форумі 2005 року В. Хубер пояснює подібні та інші рішення.
ТЕСТИ ТА ПРИЧИНИ
Для прямих ліній результати, як очікувалося, є точними. Але для інших геометрій результати можуть бути невтішними. Основна причина - мабуть, вся модель - для точних прямокутників, або для геометрій, які можна наблизити до «смугового прямокутника». Тут "тестовий комплект" для перевірки меж цього наближення (див. wfactor
У результатах вище).
SELECT *, round(100.0*(w_estim-w)/w,1) as estim_perc_error
FROM (
SELECT btype, round(len,1) AS len, w, round(w/len,3) AS wfactor,
round( buffer_width(len, gbase, btype) ,2) as w_estim ,
round( 0.5*ST_Area(gbase)/len ,2) as w_near
FROM (
SELECT
*, st_length(g) AS len, ST_Buffer(g, w, btype) AS gbase
FROM (
-- SELECT ST_GeomFromText('LINESTRING(50 50,150 150)') AS g, -- straight
SELECT ST_GeomFromText('LINESTRING(50 50,150 150,150 50,250 250)') AS g,
unnest(array[1.0,10.0,20.0,50.0]) AS w
) AS t,
(SELECT unnest(array['endcap=flat','endcap=flat join=bevel']) AS btype
) AS t2
) as t3
) as t4;
РЕЗУЛЬТАТИ:
З ПРАВИЛЬНИМИ (центральна лінія - ПРАВА ЛІНІЯ):
btype | len | w | wfactor | w_estim | w_near | estim_perc_error
------------------------+-------+------+---------+---------+--------+------------------
endcap=flat | 141.4 | 1.0 | 0.007 | 1 | 1 | 0
endcap=flat join=bevel | 141.4 | 1.0 | 0.007 | 1 | 1 | 0
endcap=flat | 141.4 | 10.0 | 0.071 | 10 | 10 | 0
endcap=flat join=bevel | 141.4 | 10.0 | 0.071 | 10 | 10 | 0
endcap=flat | 141.4 | 20.0 | 0.141 | 20 | 20 | 0
endcap=flat join=bevel | 141.4 | 20.0 | 0.141 | 20 | 20 | 0
endcap=flat | 141.4 | 50.0 | 0.354 | 50 | 50 | 0
endcap=flat join=bevel | 141.4 | 50.0 | 0.354 | 50 | 50 | 0
З ІНШИМИ ГЕОМЕТРІЯми (у центрі складена):
btype | len | w | wfactor | w_estim | w_near | estim_perc_error
-----------------------+-----+------+---------+---------+--------+------------------
endcap=flat | 465 | 1.0 | 0.002 | 1 | 1 | 0
endcap=flat join=bevel | 465 | 1.0 | 0.002 | 1 | 0.99 | 0
endcap=flat | 465 | 10.0 | 0.022 | 9.98 | 9.55 | -0.2
endcap=flat join=bevel | 465 | 10.0 | 0.022 | 9.88 | 9.35 | -1.2
endcap=flat | 465 | 20.0 | 0.043 | 19.83 | 18.22 | -0.9
endcap=flat join=bevel | 465 | 20.0 | 0.043 | 19.33 | 17.39 | -3.4
endcap=flat | 465 | 50.0 | 0.108 | 46.29 | 40.47 | -7.4
endcap=flat join=bevel | 465 | 50.0 | 0.108 | 41.76 | 36.65 | -16.5
wfactor= w/len
w_near = 0.5*area/len
w_estim is the proposed estimator, the buffer_width function.
Про btype
дивіться посібник ST_Buffer , з хорошими ілюстрастратинами та LINESTRING, які використовуються тут.
ВИСНОВКИ :
- оцінка
w_estim
завжди краща, ніж w_near
;
- для "близької до прямокутної"
g2
геометрії, її нормально, будь-якаwfactor
- для іншої геометрії (поблизу "прямокутних смуг") використовуйте обмеження
wfactor=~0.01
на 1% помилки на w_estim
. До цього wfactor використовуйте інший оцінювач.
Обережність і профілактика
Чому виникає помилка оцінки? Коли ви використовуєте ST_Buffer(g,w)
, ви «від прямокутної моделі смужок» очікуєте, що нова область, додана буфером ширини, w
становить приблизно w*ST_Length(g)
або w*ST_Perimeter(g)
... Коли ні, зазвичай це накладення (див. Складені лінії) або «стилізація» - це коли оцінка середньої w
несправності . Це головне повідомлення тестів.
Щоб виявити цю проблему в будь-якому короткі буфері , перевірте поведінку генерації буфера:
SELECT btype, w, round(100.0*(a1-len1*2.0*w)/a1)::varchar||'%' AS straight_error,
round(100.0*(a2-len2*2.0*w)/a2)::varchar||'%' AS curve2_error,
round(100.0*(a3-len3*2.0*w)/a3)::varchar||'%' AS curve3_error
FROM (
SELECT
*, st_length(g1) AS len1, ST_Area(ST_Buffer(g1, w, btype)) AS a1,
st_length(g2) AS len2, ST_Area(ST_Buffer(g2, w, btype)) AS a2,
st_length(g3) AS len3, ST_Area(ST_Buffer(g3, w, btype)) AS a3
FROM (
SELECT ST_GeomFromText('LINESTRING(50 50,150 150)') AS g1, -- straight
ST_GeomFromText('LINESTRING(50 50,150 150,150 50)') AS g2,
ST_GeomFromText('LINESTRING(50 50,150 150,150 50,250 250)') AS g3,
unnest(array[1.0,20.0,50.0]) AS w
) AS t,
(SELECT unnest(array['endcap=flat','endcap=flat join=bevel']) AS btype
) AS t2
) as t3;
РЕЗУЛЬТАТИ:
btype | w | straight_error | curve2_error | curve3_error
------------------------+------+----------------+--------------+--------------
endcap=flat | 1.0 | 0% | -0% | -0%
endcap=flat join=bevel | 1.0 | 0% | -0% | -1%
endcap=flat | 20.0 | 0% | -5% | -10%
endcap=flat join=bevel | 20.0 | 0% | -9% | -15%
endcap=flat | 50.0 | 0% | -14% | -24%
endcap=flat join=bevel | 50.0 | 0% | -26% | -36%