0

我在 Postgres 中有一个 PostGIS 点数据库，我想将几​​个地理上不同区域的点提取到 CSV 文件中，每个区域一个文件。

我已经建立了一个area带有区域多边形和区域标题的表格，我想有效地循环通过该表格，使用 Postgis 之类的东西st_intersects()来选择要进入每个 CSV 文件的数据，并从 CSV 文件中获取文件名面积表中的标题。

我对执行交集代码和设置 CSV 输出的细节感到满意——我不知道如何为每个区域执行此操作。是否可以通过某种加入来做这样的事情？或者我是否需要使用存储过程来完成，并在 plpgsql 中使用循环构造？

4

2 回答 2

0

您可以使用 plpgsql 函数或内联 do （如果您只需要执行一次并且不想存储函数。）

do $body$
  DECLARE i int;
  BEGIN FOR i IN SELECT DISTINCT city FROM table
    LOOP RAISE
      NOTICE 'foo';
      EXECUTE format($$COPY (SELECT * FROM foo WHERE x='%s') TO /tmp/%s$$, i, i); 
    END LOOP;
    RETURN;
  END;
$body$ LANGUAGE plpgsql;
于 2017-04-08T20:38:51.190 回答
0

您可以area在 plpgsql 中遍历表中的行。但请注意正确引用标识符和值：

假设这个设置：

CREATE TABLE area (
  title text PRIMARY KEY
, area_polygon geometry
); 

CREATE TABLE points(
  point_id serial PRIMARY KEY
, the_geom geometry);

你可以使用这个 plpgsql 块：

DO
$do$
DECLARE
   _title text;
BEGIN
   FOR _title IN
      SELECT title FROM area
   LOOP
      EXECUTE format('COPY (SELECT p.*
                            FROM   area   a
                            JOIN   points p ON ST_INTERSECTS(p.the_geom, a.area_polygon)
                            WHERE  a.title = %L) TO %L (FORMAT csv)'
                   , _title
                   , '/path/to/' || _title || '.csv');
   END LOOP;
END
$do$;

使用 format with %L(for string l iteral) 来获得正确引用的字符串，以避免语法错误和可能的SQL 注入。您仍然需要在area.title该工作中使用字符串作为文件名。）

还要小心地引用整个文件名，而不仅仅是它的标题部分。

您必须将整个命令连接为string。“实用命令”COPY不允许变量替换。这只能通过核心 DML 命令SELECT、INSERT、UPDATE和DELETE. 看：

所以不要area.area_polygon在循环中读出。必须将text其转换为将其连接到查询字符串中，其中文本表示将被转换回geometry（或您实际未公开的数据类型是什么）。这很容易出错。

相反，我只阅读area.title以唯一标识行并在内部处理查询中的其余部分。

于 2017-04-09T20:02:15.123 回答