Индексирование запроса PostgreSQL с использованием геометрического поля (gist), но, похоже, оно не работает

Я вспомнил, чтобы удалить и воссоздать индекс, чтобы исключить вероятность того, что что-то изменилось

gidrop index if exists incidents;
create index idx_geo on incidents using gist(geo);

Тогда это:

vacuum analyze incidents;

Вот запрос, который все еще O(n^2) и невероятно медленный. Обратите внимание, что я запустил его с и без enable_seqscan. Нет разницы.

SET enable_seqscan TO off;

drop table if exists BSC;

create table BSC(
  most_recent_id int not null,
  incident_id int not null
);

insert into BSC(most_recent_id, incident_id)
select *
from (
  select
    (select max(id)
     from incidents i2
     where i2.geo_mesh && i.geo_mesh
       and ST_DWithin(i2.geo_mesh, i.geo_mesh, 0)
       and i2.id in (select most_recent_id from temp_unique)
    ) as most_recent_id,
    id as incident_id
  from incidents i
  where i.id in (select most_recent_id from temp_unique)
) t
where t.most_recent_id <> t.incident_id;

SET enable_seqscan TO on