AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • Início
  • system&network
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • Início
  • system&network
    • Recentes
    • Highest score
    • tags
  • Ubuntu
    • Recentes
    • Highest score
    • tags
  • Unix
    • Recentes
    • tags
  • DBA
    • Recentes
    • tags
  • Computer
    • Recentes
    • tags
  • Coding
    • Recentes
    • tags
Início / dba / Perguntas / 334021
Accepted
idrizza
idrizza
Asked: 2023-12-13 22:12:37 +0800 CST2023-12-13 22:12:37 +0800 CST 2023-12-13 22:12:37 +0800 CST

Dividir linhas em polígonos em um gatilho PostGIS

  • 772

Eu tenho uma tabela de linhas nomeada segmentarmazenada no referenceesquema e uma tabela de polígonos nomeada communesem um donnees_externesesquema. Tento criar um gatilho que divida os recursos dos segmentrecursos communesna atualização do geomcampo ou na inserção de um novo arquivo segment. Duas situações:

  • inserção: o novo recurso deve ser dividido nos limites dos communespolígonos, e cada linha resultante deve ter alguns atributos do polígono do qual faz parte
  • atualização: o recurso existente deve continuar a existir, seu id não deve mudar, mas deve ter uma nova divisão geográfica nos limites dos polígonos e novos atributos, se necessário. Novos recursos devem ser inseridos na mesma tabela para cada novo polígono que a NEWgeometria intercepta.

Meu código abaixo atualiza o recurso existente de acordo, mas as linhas que deveriam ser inseridas na tabela não o são. Acho que tem algo a ver com a forma como lido com os retornos da função, mas tentei muitas combinações sem sucesso.

CREATE OR REPLACE FUNCTION split_segments_on_communes_and_update_collectivities_codes() 
RETURNS TRIGGER AS
$BODY$
DECLARE
rec RECORD;
BEGIN
    FOR rec IN
        SELECT c.insee_com AS insee_com,
               '0' || c.insee_dep AS insee_dep,
               c.insee_reg AS insee_reg,
               c.code_epci AS code_epci,
               (ST_Dump(ST_Intersection(NEW.geom, c.geom))).geom AS geom,
               (ST_Dump(ST_Intersection(NEW.geom, c.geom))).path[0] AS "path",
               NEW.id_statut,
               NEW.propriete,
               NEW.author,
               NEW.create_date
          FROM donnees_externes.commune c
          WHERE ST_Intersects(NEW.geom, c.geom)
    LOOP
        IF rec."path" = 1 OR rec."path" IS NULL THEN
            NEW.geom = rec.geom;
            NEW.code_reg = rec.insee_reg;
            NEW.code_dep = rec.insee_dep;
            NEW.code_epci = rec.code_epci;
            NEW.code_commune = rec.insee_com;
        ELSE
            INSERT INTO reference.segment (geom,
                                           id_statut,
                                           propriete,
                                           code_epci,
                                           code_dep,
                                           code_reg,
                                           author,
                                           create_date,
                                           code_commune)
            VALUES (
                   rec.geom,
                   NEW.id_statut,
                   NEW.propriete,
                   rec.code_epci,
                   rec.insee_dep,
                   rec.insee_reg,
                   NEW.author,
                   NEW.create_date,
                   rec.insee_com);
        END IF;
    END LOOP;
RETURN NEW;
END;
$BODY$
LANGUAGE plpgsql;


CREATE OR REPLACE TRIGGER split_segments_on_communes_and_update_collectivities_codes
 BEFORE INSERT
    OR UPDATE OF geom
    ON reference.segment
 FOR EACH ROW
 WHEN (pg_trigger_depth() < 1)
 EXECUTE PROCEDURE split_segments_on_communes_and_update_collectivities_codes();
postgresql
  • 1 1 respostas
  • 41 Views

1 respostas

  • Voted
  1. Best Answer
    idrizza
    2023-12-20T00:31:54+08:002023-12-20T00:31:54+08:00

    O problema estava na pathdefinição naquela linha: (ST_Dump(ST_Intersection(NEW.geom, c.geom))).path[0] AS "path".

    Aqui está a maneira correta de fazer isso:

    WITH parts AS (
                SELECT c.insee_com AS insee_com,
                       '0' || c.insee_dep AS insee_dep,
                       c.insee_reg AS insee_reg,
                       c.code_epci AS code_epci,
                       (ST_Dump(ST_Intersection(NEW.geom, c.geom))).geom AS geom,
                       NEW.id_statut,
                       NEW.propriete,
                       NEW.author,
                       NEW.create_date
                  FROM donnees_externes.commune c
                 WHERE ST_Intersects(NEW.geom, c.geom)
            )
            SELECT *,
                   ROW_NUMBER() OVER () AS "path"
              FROM parts
    

    E a função completa:

    CREATE OR REPLACE FUNCTION reference.split_segments_on_communes_and_update_collectivities_codes()
     RETURNS trigger
     LANGUAGE plpgsql
    AS $function$
    DECLARE
    rec RECORD;
    newgeom geometry;
    BEGIN
        FOR rec IN
            WITH parts AS (
                SELECT c.insee_com AS insee_com,
                       '0' || c.insee_dep,
                       c.insee_reg AS insee_reg,
                       c.code_epci AS code_epci,
                       (ST_Dump(ST_Intersection(NEW.geom, c.geom))).geom AS geom
                  FROM donnees_externes.commune c
                 WHERE ST_Intersects(NEW.geom, c.geom)
            )
            SELECT *,
                   ROW_NUMBER() OVER () AS "path"
              FROM parts
        LOOP 
            IF ST_GeometryType(rec.geom) != 'ST_Point' THEN
               IF rec."path" = 1 THEN
                  NEW.geom := rec.geom;
                  NEW.code_reg := rec.insee_reg;
                  NEW.code_dep := rec.insee_dep;
                  NEW.code_epci := rec.code_epci;
                  NEW.code_commune := rec.insee_com;
               ELSE
                  INSERT INTO reference.segment (geom,
                                                 id_statut,
                                                 propriete,
                                                 code_epci,
                                                 code_dep,
                                                 code_reg,
                                                 author,
                                                 create_date,
                                                 code_commune)
                  VALUES (
                         rec.geom,
                         NEW.id_statut,
                         NEW.propriete,
                         rec.code_epci,
                         rec.insee_dep,
                         rec.insee_reg,
                         NEW.author,
                         NEW.create_date,
                         rec.insee_com);
               END IF;
            END IF;
        END LOOP;
    RETURN NEW;
    END;
    $function$
    ;
    
    • 0

relate perguntas

  • Posso ativar o PITR depois que o banco de dados foi usado

  • Práticas recomendadas para executar a replicação atrasada do deslocamento de tempo

  • Os procedimentos armazenados impedem a injeção de SQL?

  • Sequências Biológicas do UniProt no PostgreSQL

  • Qual é a diferença entre a replicação do PostgreSQL 9.0 e o Slony-I?

Sidebar

Stats

  • Perguntas 205573
  • respostas 270741
  • best respostas 135370
  • utilizador 68524
  • Highest score
  • respostas
  • Marko Smith

    conectar ao servidor PostgreSQL: FATAL: nenhuma entrada pg_hba.conf para o host

    • 12 respostas
  • Marko Smith

    Como fazer a saída do sqlplus aparecer em uma linha?

    • 3 respostas
  • Marko Smith

    Selecione qual tem data máxima ou data mais recente

    • 3 respostas
  • Marko Smith

    Como faço para listar todos os esquemas no PostgreSQL?

    • 4 respostas
  • Marko Smith

    Listar todas as colunas de uma tabela especificada

    • 5 respostas
  • Marko Smith

    Como usar o sqlplus para se conectar a um banco de dados Oracle localizado em outro host sem modificar meu próprio tnsnames.ora

    • 4 respostas
  • Marko Smith

    Como você mysqldump tabela (s) específica (s)?

    • 4 respostas
  • Marko Smith

    Listar os privilégios do banco de dados usando o psql

    • 10 respostas
  • Marko Smith

    Como inserir valores em uma tabela de uma consulta de seleção no PostgreSQL?

    • 4 respostas
  • Marko Smith

    Como faço para listar todos os bancos de dados e tabelas usando o psql?

    • 7 respostas
  • Martin Hope
    Jin conectar ao servidor PostgreSQL: FATAL: nenhuma entrada pg_hba.conf para o host 2014-12-02 02:54:58 +0800 CST
  • Martin Hope
    Stéphane Como faço para listar todos os esquemas no PostgreSQL? 2013-04-16 11:19:16 +0800 CST
  • Martin Hope
    Mike Walsh Por que o log de transações continua crescendo ou fica sem espaço? 2012-12-05 18:11:22 +0800 CST
  • Martin Hope
    Stephane Rolland Listar todas as colunas de uma tabela especificada 2012-08-14 04:44:44 +0800 CST
  • Martin Hope
    haxney O MySQL pode realizar consultas razoavelmente em bilhões de linhas? 2012-07-03 11:36:13 +0800 CST
  • Martin Hope
    qazwsx Como posso monitorar o andamento de uma importação de um arquivo .sql grande? 2012-05-03 08:54:41 +0800 CST
  • Martin Hope
    markdorison Como você mysqldump tabela (s) específica (s)? 2011-12-17 12:39:37 +0800 CST
  • Martin Hope
    Jonas Como posso cronometrar consultas SQL usando psql? 2011-06-04 02:22:54 +0800 CST
  • Martin Hope
    Jonas Como inserir valores em uma tabela de uma consulta de seleção no PostgreSQL? 2011-05-28 00:33:05 +0800 CST
  • Martin Hope
    Jonas Como faço para listar todos os bancos de dados e tabelas usando o psql? 2011-02-18 00:45:49 +0800 CST

Hot tag

sql-server mysql postgresql sql-server-2014 sql-server-2016 oracle sql-server-2008 database-design query-performance sql-server-2017

Explore

  • Início
  • Perguntas
    • Recentes
    • Highest score
  • tag
  • help

Footer

AskOverflow.Dev

About Us

  • About Us
  • Contact Us

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve