2014-01-13 81 views
2

我想优化下面的查询。postgresql hashaggregate查询优化

select cellid2 as cellid, max(endeks) as turkcell 
from (select a.cellid2 as cellid2, b.endeks 
    from (select geom, cellid as cellid2 from grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000) a join (select endeks, st_transform(geom, 2320) as geom_tmp from turkcell_data) b on st_intersects(a.geom, b.geom_tmp)) x 
group by cellid2 limit 5 

和解释分析回报

"Limit (cost=81808.31..81808.36 rows=5 width=12) (actual time=271376.201..271376.204 rows=5 loops=1)" 
" -> HashAggregate (cost=81808.31..81879.63 rows=7132 width=12) (actual time=271376.200..271376.203 rows=5 loops=1)" 
"  -> Nested Loop (cost=0.00..81772.65 rows=7132 width=12) (actual time=5.128..269753.647 rows=1237707 loops=1)" 
"    Join Filter: _st_intersects(grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000.geom, st_transform(turkcell_data.geom, 2320))" 
"    -> Seq Scan on turkcell_data (cost=0.00..809.40 rows=3040 width=3045) (actual time=0.031..7.426 rows=3040 loops=1)" 
"    -> Index Scan using grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_geom_gist on grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000 (cost=0.00..24.76 rows=7 width=124) (actual time=0.012..0.799 rows=647 loops=3040)" 
"     Index Cond: (geom && st_transform(turkcell_data.geom, 2320))" 
"Total runtime: 271387.499 ms" 

存在着几何列和小区id列的索引。我读过,而不是使用max,desc命令并限制1更好。然而,由于我有分组的条款,我认为这是行不通的。有没有办法做到这一点或其他方式来提高性能?

表定义:

CREATE TABLE grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000 
(
    regionid numeric, 
    geom geometry(Geometry,2320), 
    cellid integer, 
    turkcell double precision 
) 
WITH (
    OIDS=FALSE 
); 
ALTER TABLE grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000 
    OWNER TO postgres; 

-- Index: grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_cellid 

-- DROP INDEX grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_cellid; 

CREATE INDEX grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_cellid 
    ON grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000 
    USING btree 
    (cellid); 

-- Index: grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_geom_gist 

-- DROP INDEX grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_geom_gist; 

CREATE INDEX grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000_geom_gist 
    ON grd_90098780_7c48_11e3_8876_f0bf97e0dd001000000000 
    USING gist 
    (geom); 

CREATE TABLE turkcell_data 
(
    gid serial NOT NULL, 
    objectid_1 integer, 
    objectid integer, 
    neighbourh numeric, 
    endeks numeric, 
    coorx numeric, 
    coory numeric, 
    shape_leng numeric, 
    shape_le_1 numeric, 
    shape_area numeric, 
    geom geometry(MultiPolygon,4326), 
    CONSTRAINT turkcell_data_pkey PRIMARY KEY (gid) 
) 
WITH (
    OIDS=FALSE 
); 
ALTER TABLE turkcell_data 
    OWNER TO postgres; 

-- Index: turkcell_data_geom_gist 

-- DROP INDEX turkcell_data_geom_gist; 

CREATE INDEX turkcell_data_geom_gist 
    ON turkcell_data 
    USING gist 
    (geom); 
+1

如果您希望我们帮助优化查询,您需要向我们显示表格和索引定义**以及每个表格的行数。也许你的表格定义不好。也许索引没有正确创建。也许你没有一个你认为你做过的那个专栏的索引。没有看到表和索引定义,我们不能说。我们还需要行计数,因为这会大大影响查询优化。 –

+0

我已经添加了必要的定义。 – adaminasabi

+0

您的嵌套循环具有N = 7和N = 3040的子查询,并导致N = 1237707行。这比汽车产品还要糟糕! – joop

回答

2

无论是存储你的数据重新投影到2320,索引列,并在您的加入使用它,或者在几何在turkcell_data变换投影创建索引。我通常更喜欢后者:

CREATE INDEX turkcell_data_geom_gist2320 
    ON turkcell_data 
    USING gist 
    (st_transform(geom, 2320)); 

另一个问题可能是,如果你的几何形状是非常复杂的 - 如果你的任何多边形具有相对大量的你可能会卡住捣鼓远的交点。不过,请先尝试索引。

+0

我把这个索引放进去了,但并没有改变那么多。 – adaminasabi

+0

首先你可以“设置enable_seqscan = false”,然后运行并发布解释分析(以验证索引是否正确创建)?然后“设置enable_seqscan = true”,并再次运行并发布解释分析。你还可以描述你的图层吗?它看起来像一层大约600个特征,另一个大约3000个。是否有任何几何体具有过多的点数?几乎每一个第一层的几何都与第二层几乎每一个几何相交? – yieldsfalsehood