Билинейная интерполяция

[Svility]

Добрый день.
Есть ли уже реализованный алгоритм билинейной интерполяции на Postgre/Postgis?
Собственно необходим алгоритм, который вернет высотную отметку с растра (загруженного в PG) по входным данным - X,Y.
Алгоритм: https://ru.wikipedia.org/wiki/Билинейная_интерполяция

Спойлер

ST_Value(rastr, geom) - отбирает высоты по методу "ближайший сосед" - быстро, но не точно (целые значения).

[oldbay]

Посмотрите в сторону SAGA GIS http://www.saga-gis.org
В качестве источника там можно использовать postgis - а в инструментах посмотрите на "Shapes" -> "Grid Tools"

[Svility]

Модуль rasterinterpolation выполняет билинейную и бикубическую интерполяцию... А вот как бы функцию описать на pg?

[antonv]

Svility, для билинейной интерполяции нужны значения четырёх пикселей, верно? Функция ST_Value умеет не только по координатам точки возвращать значение, но и по номеру столбца и строки ячейки в растре. Верните себе четыре значения вокруг точки и проинтерполируйте сами формулой из Википедии:)

[MBG]

См. https://postgis.net/docs/RT_ST_Resample.html - интерполирует целый растр, в том числе билинейная интерполяция поддерживается. Если определить новый растр из одной точки, получите требуемое.

[Svility]

Svility, для билинейной интерполяции нужны значения четырёх пикселей, верно? Функция ST_Value умеет не только по координатам точки возвращать значение, но и по номеру столбца и строки ячейки в растре. Верните себе четыре значения вокруг точки и проинтерполируйте сами формулой из Википедии:)

А нет справки никакой или примера. Теория стала понятной, но с реализацией беда.

[Svility]

Ну вот определяю я границы пикселя, вокруг искомой точки, но что то не пойму что дальше нужно сделать.. как то пиксель изменить или с границами пикселя какую то операцию провести нужно...
Кто подскажет в каком направлении дальше двигаться?

Код: Выделить всё

SELECT rid, ST_AsText(ST_Envelope(rast)) As geom, ST_Value(rast, ST_SetSRID(ST_POINT(40.074284, 44.702956),4326)) As z, '219.061719108487' z_RI /*эталон*/
FROM dem b
where ST_Intersects(ST_SetSRID(ST_POINT(40.074284, 44.702956),4326), b.rast)

Снимок экрана (60).png (12.93 КБ) 10241 просмотр