Question

在从这个网站收到很多帮助后，我终于创建了一组功能，完成了我想要的工作。谢谢你，但似乎有一个最后的问题，即高效的内存使用。这是问题：
在postgreSQL 9.3中，使用postGIS 2.1和pgRouting 2.0，我创建了一个函数，使用pgrouting函数pgr_trsp计算2个点之间的路径，并返回几何（Linestring）值。这是代码：

CREATE OR REPLACE FUNCTION fm_pgr2geom(edge1 integer, pos1 double precision, edge2 integer, pos2 double precision)
  RETURNS geometry AS
$BODY$
--We have to do a routing query. And declare a cursor for it
DECLARE resc CURSOR FOR
SELECT * FROM pgr_trsp (
    'SELECT * FROM th_2po_4pgr',
    $1, $2, $3, $4, false, true);
doline geometry[];
temp_point geometry;
geom geometry;
temp_rec RECORD;
n integer;
BEGIN

--Append all the edges
FOR temp_rec IN SELECT * FROM pgr_trsp (
    'SELECT * FROM th_2po_4pgr',
    $1, $2, $3, $4, false, true) LOOP
        doline := array_append(
        doline, (SELECT map.geom_way FROM th_2po_4pgr map WHERE map.id = temp_rec.id2));
END LOOP;
--Remove 1st and last edge
n := array_length (doline, 1);
doline := doline [2:n-1];
--Find startpoint and append to doline
doline := array_prepend(
    ST_LineInterpolatePoint((SELECT map.geom_way FROM th_2po_4pgr map WHERE map.id = $1),$2),doline);
--Append the endpoint
doline := array_append(
    doline,ST_LineInterpolatePoint((SELECT map.geom_way FROM th_2po_4pgr map WHERE map.id = $3),$4));
geom := ST_MakeLine(doline);
RETURN geom;
EXCEPTION
WHEN SQLSTATE 'XX000' THEN RETURN NULL;
WHEN SQLSTATE '38001' THEN RETURN NULL;
END;
$BODY$
  LANGUAGE plpgsql VOLATILE
  COST 100;
ALTER FUNCTION fm_pgr2geom(integer, double precision, integer, double precision)
  OWNER TO postgres;

此函数在另一个函数中用于批量更新具有几何结果的大表（800k +）。这里仅供参考：

CREATE OR REPLACE FUNCTION fm_seqrouting()
  RETURNS integer AS
$BODY$
--Declarations
DECLARE
    r record;
    i integer;
BEGIN
--CODE to calculate routes and update table
    i := 0;
FOR r IN 
    SELECT  veh_id
        ,dt
        ,map_edge_id                AS map_id1
        ,map_edge_pos               AS map_pos1
        ,lead(map_edge_id)      OVER w  AS map_id2
        ,lead(map_edge_pos)     OVER w  AS map_pos2
    FROM taxilocs
    WINDOW w AS (ORDER BY veh_id, dt)
    LOOP

        UPDATE taxilocs
            SET geom_route = fm_pgr2geom (r.map_id1,r.map_pos1,r.map_id2,r.map_pos2)
            WHERE r.veh_id = taxilocs.veh_id AND r.dt=taxilocs.dt;
        i := i + 1;

    END LOOP;
RETURN i;
END;
$BODY$
  LANGUAGE plpgsql;

异常是绝对必要的，因为它们处理某些情况，其中上述表中缺少某些数据，或者找不到路由路径。但是，似乎它们导致了问题，即更新查询在几分钟后崩溃。我执行几分钟后收到的消息是：

错误：内存不足
   SQL状态：53200
   细节：根据要求失败   大小640000。

所以问题是：我如何有效地使用或重新编码此函数来更新我想要的表？任何想法？
提前谢谢！

Answer 1

尝试使用以下命令从命令行运行sql psql -h＆lt;主机名＆GT; -d＆lt; DBNAME＆GT; -U＆lt;用户名＆gt; -p＆lt;端口名＆GT; -f＆lt;文件名＆GT;

PostgreSQL函数和内存问题

1 个答案: