Se você deseja criar em lote índices em colunas de geometria, tente esta função plpgsql que acabei de chamar:
CREATE OR REPLACE FUNCTION BatchIndex(sn text, tn text, cn text) RETURNS void AS $$
DECLARE i_exists integer;
DECLARE idxname text;
BEGIN
idxname := 'idx_' || tn || '_' || cn;
select into i_exists count(*) from pg_class where relname = idxname;
IF i_exists = 0 THEN
EXECUTE 'CREATE INDEX ' ||idxname || ' ON '
|| sn || '.' || tn
|| ' USING GIST(' || cn || ')';
END IF;
END;
$$ LANGUAGE plpgsql;
Eu não tentei com raiva no meu banco de dados, mas parece fazer o trabalho.
Para usá-lo, basta executar uma SELECT
declaração como esta:
select BatchIndex('public', f_table_name, f_geometry_column) from geometry_columns where f_table_schema = 'public';
Para criar índices em todas as colunas geométricas, você pode usá-lo assim:
select BatchIndex(f_table_schema, f_table_name, f_geometry_column) from geometry_columns;
Depois, execute um VACUUM ANALYZE
para arrumar tudo.