Le portail francophone de la géomatique

jonathan · Tue 15 July 2008 17:05

Bonjour,

Je travail sur postgis/winXP.
Je souhaite créer une table basée sur deux autres (couches de points) en associant aux objets de la première une propriété de l'objet le plus proche dans la deuxième.

J'utilise une requête :

Code:

select a.*, b.alti from tmp_pts a,mnt b
where Distance(a.the_geom,b.the_geom)=(select min(foo.valeur) from (select d.alti as z,distance(d.the_geom,c.the_geom) 
as valeur from tmp_pts c, mnt d)foo);

basée sur les publications de Jean David Techer (merci à lui ).

http://www.davidgis.fr/documentation/wi … #id2508149

Cependant le traitement est très long (une couche >1.5millions de points/1>50000points), en attendant le résultat pouvais vous me faire vos remarques éventuelles sur la pertinence (ou la bonne adaptation) de cette requête.

Merci

(desole pour le titre...)

En fait il s'agit de faire l'élévation d'un couche par rapport à un MNT (c'est vrai que c'est plus clair comme ca).

Je n'ai pour l'instant pas eu de résultat probant...

Dernière modification par jonathan (Thu 17 July 2008 14:55)

jonathan · Thu 17 July 2008 16:29

J'ai créé des index gist et adapté ma requête. Cependant j'obtiens toujours une erreur.

Code:

create table result as select a.*, b."X", b."Y", b."Z" from tmp_pts_aep a,mnt_lg b
where Distance(a.the_geom,b.the_geom)=(select min(distance(d.the_geom,c.the_geom)) from tmp_pts_aep c, mnt_lg d
where d.the_geom && 'BOX3D(x(c.the_geom)-30 y(c.the_geom)-30, x(c.the_geom)+30 y(c.the_geom)+30)'::box3d);

erreur :

Code:

ERREUR:  BOX3D parser - couldnt parse.  It should look like: BOX3D(xmin ymin zmin,xmax ymax zmax) or BOX3D(xmin ymin,xmax ymax)

********** Erreur **********

ERREUR: BOX3D parser - couldnt parse.  It should look like: BOX3D(xmin ymin zmin,xmax ymax zmax) or BOX3D(xmin ymin,xmax ymax)
État SQL :XX000

Merci d'avance pour votre aide.

Dernière modification par jonathan (Thu 17 July 2008 16:33)

ranalisolofo.larrys · Fri 18 July 2008 08:33

Au lieu d'utiliser BOX3D utilise plutôt BOX2D parce que tu n'as que le x et y.

jonathan · Fri 18 July 2008 08:48

Bonjour,

Je me suis aussi fais cette réflexion mais il semble qu'on puisse utiliser la box3d de cette façon.

Sinon :

Code:

....'BOX2D(x(c.the_geom)-25 y(c.the_geom)-25, x(c.the_geom)+25 y(c.the_geom)+25)'::box2d);

erreur :

Code:

ERREUR:  BOX2DFLOAT4 parser - doesnt start with BOX(

et

Code:

....'BOX(x(c.the_geom)-25 y(c.the_geom)-25, x(c.the_geom)+25 y(c.the_geom)+25)'::box2d);

erreur :

Code:

ERREUR:  BOX2DFLOAT4 parser - couldnt parse.  It should look like: BOX(xmin ymin,xmax ymax)

Merci.

Dernière modification par jonathan (Fri 18 July 2008 08:49)

jonathan · Fri 18 July 2008 10:14

Le problème, c'est ce qu'il y a dans les box3D : x(a.the_geom)+30...

Peut on utiliser une fonction à cette endroit ? Si oui comment ?

Merci

Guillaume Sueur · Fri 18 July 2008 10:28

et pourquoi pas ST_expand(geometry, float) ?
cela permet de définir une BBOX à partir de l'extent d'une géométrie augmentée dans tous les sens de la valeur du deuxième paramètre.

Guillaume

jonathan · Fri 18 July 2008 11:03

Super !!!

La requête est partie. J'attends le résultat.

Merci.

jonathan · Mon 21 July 2008 08:50

La requête ne m'a retourné que deux lignes : celles pour lesquelles les points d'une couche avaient les coordonnées x et y égales (sur la partie entière) aux points de l'autre couche.
2 points sur environ 20000 !

jonathan · Wed 23 July 2008 08:44

Code:

SELECT a.*, b.gid, b."X",b."Y",b."Z"
FROM tmp_pts_aep a,mnt_lg b
WHERE b.gid=(SELECT  d.gid FROM tmp_pts_aep c, mnt_lg d
WHERE d.the_geom && buffer(c.the_geom,30) 
ORDER BY distance(a.the_geom,b.the_geom) LIMIT 1)
;

La requête tourne depuis plus de 20 heures sans résultat...

Au secours

Yves · Wed 23 July 2008 09:27

Bonjour,

Quels sont les index créés sur cette requête ?

Y.

jonathan · Wed 23 July 2008 10:58

Des index gist sont créés sur chacune des tables.

Yves · Wed 23 July 2008 11:37

Bonjour,

Je pense pas que la manière dont tu t'y prennes soit la meilleur. Les Buffer, ORDER BY, distance() sont des fonctions/commandes qui demandent pas mal de ressource !

Voici les améliorations que je te propose :
1/ créer un champ buffer dans ta table tmp_pts_aep
2/ Créer une table temporaire :

Code:

CREATE TABLE tmp AS SELECT  d.gid FROM tmp_pts_aep c, mnt_lg d
WHERE d.the_geom && c.buffer 
ORDER BY distance(a.the_geom,b.the_geom) LIMIT 1

3/ faire ta requête finale :

Code:

SELECT a.*, b.gid, b."X",b."Y",b."Z"
FROM tmp_pts_aep a,mnt_lg b
INNER JOIN tmp USING(gid);

Y.

jonathan · Wed 23 July 2008 13:52

Merci,

Le problème, c'est que cette requête retourne une table d'une seule ligne.

Code:

CREATE TABLE tmp AS SELECT b.gid FROM tmp_pts_aep a, mnt_lg b
WHERE b.the_geom && a.buffer 
ORDER BY distance(a.the_geom,b.the_geom) LIMIT 1

Yves · Wed 23 July 2008 15:13

Bien évidement, avec un LIMIT 1

La clause LIMIT est constituée de deux sous-clauses indépendantes :

LIMIT { nombre | ALL }
OFFSET début

nombre spécifie le nombre maximum de lignes à renvoyer alors que début spécifie le nombre de lignes à passer avant de commencer à renvoyer des lignes. Lorsque les deux clauses sont spécifiées, début lignes sont passées avant de commencer à compter les nombre lignes à renvoyer.

[source : http://docs.postgresqlfr.org/8.2/sql-se … #sql-limit]

Y.

Dernière modification par Yves (Wed 23 July 2008 15:14)

jonathan · Wed 23 July 2008 16:27

Je voyais plustot ca comme un seul point retourné par point de reference. (le point le plus pres de chaque point d'une autre couche (20000pts)

GRASS sait-il faire cette élévation.

Merci

jonathan · Mon 28 July 2008 08:42

Bonjour,

J'ai résolu mon problème sans fonctionnalité spatiale (en partie seulement).

-création d'une table (ID_DIST) avec les identifiants des points des deux tables et la distance qu'il y a entre chacun deux (dans la limite d'un buffer).
-auto jointure qui retourne l'occurrence ayant la distance minimale pour chaque identifiant d'une des couches.
-mise à jour d'une des tables en fonction de ID_DIST.

(L'élévation que je souhaite obtenir sera possible grâce à une macro fournie par le logiciel que j'utilise.)

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Tue 15 July 2008 17:05

mise à jour en fonction

Code:

#2 Thu 17 July 2008 16:29

Re: mise à jour en fonction

Code:

Code:

#3 Fri 18 July 2008 08:33

Re: mise à jour en fonction

#4 Fri 18 July 2008 08:48

Re: mise à jour en fonction

Code:

Code:

Code:

Code:

#5 Fri 18 July 2008 10:14

Re: mise à jour en fonction

#6 Fri 18 July 2008 10:28

Re: mise à jour en fonction

#7 Fri 18 July 2008 11:03

Re: mise à jour en fonction

#8 Mon 21 July 2008 08:50

Re: mise à jour en fonction

#9 Wed 23 July 2008 08:44

Re: mise à jour en fonction

Code:

#10 Wed 23 July 2008 09:27

Re: mise à jour en fonction

#11 Wed 23 July 2008 10:58

Re: mise à jour en fonction

#12 Wed 23 July 2008 11:37

Re: mise à jour en fonction

Code:

Code:

#13 Wed 23 July 2008 13:52

Re: mise à jour en fonction

Code:

#14 Wed 23 July 2008 15:13

Re: mise à jour en fonction

#15 Wed 23 July 2008 16:27

Re: mise à jour en fonction

#16 Mon 28 July 2008 08:42

Re: mise à jour en fonction

Pied de page des forums