Interpolação de valor de grade hexadecimal no POSTGIS

Eu tenho uma tabela de hexágonos em mosaico que criei como uma grade de bin hexadecimal para cobrir uma área.

Sobreponha um conjunto de pontos para atribuir um valor a alguns hexágonos, veja o diagrama: -

As formas com o valor "0" não têm valor e eu gostaria de interpolá-las dos valores dos vizinhos, conforme Kriging. Como posso fazer isso no postgres, dada a minha tabela:

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

Não sei de que tipo de interpolação você está falando, mas se a média de todos os vizinhos tiver um bom valor, essa pode ser a solução:

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
  )

Ou após as dúvidas de @MakinFlippyFloppy, se realmente 0 no exemplo significa nulo (sem valor):

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
  )

Ou se os valores nulos ou zero não devem diminuir a média:

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
  )