Le portail francophone de la géomatique

Salut,

vous pouvez trouver la liste des tables de votre base dans la relation pg_catalog.pg_tables.

SELECT format($$
           SELECT ST_ConcaveHull( ST_Collect(geom),0.70,true ) 
           INTO %1$s.%2$s_concavhull 
           FROM %1$s.%2$s 
           $$, schemaname, tablename)
FROM pg_catalog.pg_tables
WHERE schemaname = 'mon_schema';

Dernière modification par tumasgiu (Thu 04 October 2018 18:20)

Bonjour,

Comme le dit le message, un SELECT en plpgsql doit etre stocké qq part (SELECT ... INTO), ou renvoyé par la fonction.
Pour exécuter des requetes SQL, utilisez:

EXECUTE format(....);

Nicolas

Merci pour votre aide.
J'ai légèrement modifié le code pour boucler sur le code de département de ma vue de départ.

CREATE OR REPLACE FUNCTION voronoi_pt_cuivre_dep() 
RETURNS VARCHAR AS $BODY$
DECLARE
        lv VARCHAR;
BEGIN
        FOR lv IN SELECT DISTINCT code_dep
               FROM cuivre.v_france_extract_dep
               LOOP
EXECUTE format($$
            -- Calcule l'enveloppe concave pour les points d'entrée
            -- les points sur l'enveloppe concave sont ceux qui ont des polygones non fermés dans le diagramme de Voronoï en sortie.
           DROP TABLE IF EXISTS voronoi.%2$s;
           CREATE TABLE voronoi.%2$s AS
        (SELECT a.id, a.geom  
         FROM
           (WITH conv_hull(geom) AS
            (SELECT ST_ConcaveHull( ST_Collect(geom),0.70,true ) 
           --INTO %1$s.%2$s
               FROM %1$s.%2$s meshpoints), 
           -- Trouve les points d'origine de l'enveloppe concave
           edge_points(points) AS 
               (SELECT mp.geom FROM %1$s.%2$s mp, conv_hull ch 
               WHERE ST_Touches(ch.geom, mp.geom)),
           -- Buffer autour de l'enveloppe (distance arbitraire)    
             buffered_points(geom) AS             
            (SELECT ST_Buffer(geom, 100) as geom FROM conv_hull),
           -- Trouve les points les plus proches du buffer de l'enveloppe la plus proche des points de l'étape 2                     
          closest_points(points) AS     
            (SELECT ST_Closestpoint(ST_Exteriorring(bp.geom), ep.points) as points,ep.points as epoints FROM buffered_points bp, edge_points ep),
          -- Ajouter ces points au jeu de données d'origine        
          combined_points(points) AS     
        (SELECT points FROM closest_points UNION SELECT geom FROM %1$s.%2$s meshpoints),
         -- Calcule les diagrammes de voronoi à partir de ce jeu de données combiné. Tous les points ont maintenant des polygones fermés                     
          voronoi (vor) AS 
       (SELECT ST_Dump(ST_VoronoiPolygons(ST_Collect(points))) as geom FROM combined_points)    
     SELECT 
     (vor).path[1] as id, 
     (vor).geom 
 --   INTO voronoi.%2$s                                  
     FROM voronoi)a)                              
    $$, schemaname, tablename)
FROM pg_catalog.pg_tables
WHERE schemaname = 'pt_cuivre'; 
END LOOP;                                
END
$BODY$ LANGUAGE plpgsql;

Là postgres me dit que la requête a renvoyé plus d'une ligne..

J'ai corrigé ce qui ma semblé être des erreurs,
et indenter correctement.

Notez que vous ne vous servez à aucun moment de la variable lv.

CREATE OR REPLACE PROCEDURE voronoi_pt_cuivre_dep()
RETURNS VARCHAR AS $BODY$
DECLARE
        lv VARCHAR;
        r RECORD;
BEGIN
        FOR lv IN SELECT DISTINCT code_dep
               FROM cuivre.v_france_extract_dep
        LOOP
            FOR r IN SELECT schemaname, tablename
                    FROM pg_catalog.pg_tables
                    WHERE schemaname = 'pt_cuivre'
            LOOP
                --TO DO utiliser lv
                EXECUTE $$
                            -- Calcule l'enveloppe concave pour les points d'entrée
                            -- les points sur l'enveloppe concave sont ceux qui ont des polygones non fermés dans le diagramme de Voronoï en sortie.
                            DROP TABLE IF EXISTS voronoi.$2;
                            CREATE TABLE voronoi.$2 AS
                            SELECT a.id, a.geom 
                            FROM
                               (WITH conv_hull(geom) AS
                                    (SELECT ST_ConcaveHull( ST_Collect(geom),0.70,true )
                                    --INTO $1.$2
                                    FROM $1.$2 meshpoints),
                               -- Trouve les points d'origine de l'enveloppe concave
                                edge_points(points) AS
                                    (SELECT mp.geom FROM $1.$2 mp, conv_hull ch
                                    WHERE ST_Touches(ch.geom, mp.geom)),
                               -- Buffer autour de l'enveloppe (distance arbitraire)   
                                buffered_points(geom) AS             
                                    (SELECT ST_Buffer(geom, 100) as geom FROM conv_hull),
                               -- Trouve les points les plus proches du buffer de l'enveloppe la plus proche des points de l'étape 2                     
                                closest_points(points) AS     
                                    (SELECT ST_Closestpoint(ST_Exteriorring(bp.geom), ep.points) as points,ep.points as epoints FROM buffered_points bp, edge_points ep),
                              -- Ajouter ces points au jeu de données d'origine       
                                combined_points(points) AS     
                                    (SELECT points FROM closest_points UNION SELECT geom FROM $1.$2 meshpoints),
                             -- Calcule les diagrammes de voronoi à partir de ce jeu de données combiné. Tous les points ont maintenant des polygones fermés                     
                                voronoi (vor) AS
                                    (SELECT ST_Dump(ST_VoronoiPolygons(ST_Collect(points))) as geom FROM combined_points)
                                 SELECT
                                 (vor).path[1] as id,
                                 (vor).geom
                                 FROM voronoi) a                 
                        $$ USING r.schemaname, r.tablename, lv;
            END LOOP;
        END LOOP;                               
END
$BODY$ LANGUAGE plpgsql;

Dernière modification par tumasgiu (Wed 24 October 2018 15:36)