POSTGISでの16進グリッド値の補間

エリアをカバーするために六角ビニンググリッドとして作成したテッセレーションされた六角形のテーブルがあります。

次に、いくつかの六角形に値を与えるためにポイントセットをオーバーレイしました。図を参照してください：

値が「0」のシェイプには値がありません。クリギングに従って、これらを隣接の値から補間したいと思います。私のテーブルが次のようになっているとしたら、どうすればpostgresでこれを達成できますか？

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

どのような補間について話しているのかはわかりませんが、すべての近傍の平均が良い値である場合は、これが解決策になる可能性があります。

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
  )

または、@ MakinFlippyFloppyをフォローして、例で本当に0がnull（値なし）を意味するかどうかを疑います：

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
  )

または、null値またはゼロ値が平均を減少させない場合：

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
  )