Nội suy giá trị lưới hex trong POSTGIS

Tôi có một bảng các hình lục giác tessellated mà tôi đã tạo như một lưới binning hex để bao phủ một khu vực.

Sau đó, tôi đã phủ lên một điểm được đặt để cung cấp cho một số hình lục giác một giá trị, xem sơ đồ: -

Các hình có giá trị "0" không có giá trị và tôi muốn nội suy chúng từ các giá trị của hàng xóm theo Kriging. Làm thế nào tôi có thể thực hiện điều này trong các postgres khi bảng của tôi trông như sau: -

hex_grid_data
(
 gid SERIAL,
 wkb_geometry (geometry),
 value INTEGER
)

Tôi không biết bạn đang nói về loại nội suy nào nhưng nếu trung bình của tất cả các hàng xóm sẽ có giá trị tốt thì đây có thể là giải pháp:

create table hex_grid_data_av as 
  (
  select gid, wkb_geometry, value,
    case
      when value > 0 then value
      else (select sum(h2.value)/6 from hex_grid_data h2 where ST_Touches(h1.wkb_geometry, h2.wkb_geometry))
    end as int_value 
  from hex_grid_data h1
  )

Hoặc theo dõi @MakinFlippyFloppy nghi ngờ nếu thực sự 0 trong ví dụ có nghĩa là null (không có giá trị):

create table hex_grid_data_av as 
  (
  select gid, wkb_geometry, value,
    case
      when value > 0 then value
      else (select sum(coalesce(h2.value,0))/6 from hex_grid_data h2 where ST_Touches(h1.wkb_geometry, h2.wkb_geometry))
    end as int_value 
  from hex_grid_data h1
  )

Hoặc nếu giá trị null hoặc zero không nên giảm trung bình:

create table hex_grid_data_av as 
  (
  select gid, wkb_geometry, value,
    case
      when value > 0 then value
      when value = 0 and not exists (select 1 from hex_grid_data h2 where ST_Touches(h1.wkb_geometry, h2.wkb_geometry) and h2.value != 0 and h2.value is not null ) then 0
      else (select sum(coalesce(h2.value,0))/(select count(nullif(value,0)) from hex_grid_data h3 where ST_Touches(h1.wkb_geometry, h3.wkb_geometry) ) from hex_grid_data h2 where ST_Touches(h1.wkb_geometry, h2.wkb_geometry))
    end as int_value 
  from hex_grid_data h1
  )