Interpolación de valor de cuadrícula hexadecimal en POSTGIS

Tengo una tabla de hexágonos teselados que creé como una cuadrícula de agrupamiento hexadecimal para cubrir un área.

Luego superpuse un conjunto de puntos para dar un valor a algunos de los hexágonos, vea el diagrama:

Las formas con el valor "0" no tienen valor y me gustaría interpolarlas a partir de los valores de los vecinos según Kriging. ¿Cómo puedo lograr esto en postgres dado que mi tabla se ve así: -

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

No sé de qué tipo de interpolación estás hablando, pero si el promedio de todos los vecinos será un buen valor, esta podría ser la solución:

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
  )

O siguiendo las dudas de @MakinFlippyFloppy si realmente 0, por ejemplo, significa nulo (sin 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
  )

O si los valores nulos o cero no deberían disminuir el promedio:

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
  )