Message posté par : Christophe Vergon
----------------------------------------
suite
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322123#p322123
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Christophe Vergon
----------------------------------------
Bonjour,
@yartostout : c'est pas plus efficace, mais plus proche de la réalité à mon sens. En comparant le résultat à une ortho HR on constate cela.
@Sylvain : capture16 : méthode 1 (tampon sur bâti puis enveloppe puis tampon)
capture 17 : joindre bâti adjacent + buffer + enveloppe
capture 18 (message suivant) la méthode 1 avec bâti original.
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322122#p322122
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : yartostout
----------------------------------------
Merci pour le partage !
J'ai privilégié la méthode classique par tampons pour notre suivi du SCoT sur notre territoire mais il faudrait que je teste ta méthode qui semble être plus efficace d'après tes dirs !
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322120#p322120
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Sylvain M.
----------------------------------------
Merci Christophe pour ce partage !
J'ai hâte de tester cette méthode sur mon territoire, mais ce ne sera pas dans les semaines à venir étant occupé par d'autres tâches (pas la tache urbaine ;-) ).
D'ici là, si jamais tu as quelques captures d'écran des résultats pour illustrer tes requêtes, ce serait très intéressant.
A+
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322115#p322115
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Christophe Vergon
----------------------------------------
De rien Renaud ;)
J'ai constaté que l'on peut faire varier la géométrie de l'enveloppe en fonction de l'ordre dans lequel on exécute les dilatation.
Version 1 : buffer sur bâti 10 m puis enveloppe orientée puis jonction puis buffer
Version 2 : jonction bâti puis buffer 10 m puis enveloppe puis buffer
..
A adapter selon les besoins.
Pour répondre à Marc Leobet : La tâche urbaine ci-définie à pour but l'estimation des changements cadastraux depuis la rénovation d'une commune sur les section rénovées par voie de mise à jour. Et ce pour une volumétrie de futur CCTP. Comme quoi une méthode peut s'appliquer à plusieurs définitions.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322071#p322071
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Renaud Mouche (renaud(a)georezo.net)
----------------------------------------
Merci Christophe
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322066#p322066
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Christophe Vergon
----------------------------------------
Bonjour,
Une variante qui permet de compter le nombre de bâtiments (2 bâtiments adjacent comptés comme 1) dans une grappe.
1) Fusion des bâtiments adjacents
2) Buffer 10 + BBox orientée
3) Fusion des éléments précédents adjacents
4) buffer 10
5) Fusion des éléments adjacent du 4
Résultat : nombre d'éléments du 4, nombre d'éléments du 2
-----------------
Code :
WITH
p as (SELECT st_collectionextract(unnest(st_clusterintersecting(st_accum(the_geom))),3) as geom FROM cadastre.batiment ),
t as (SELECT row_number() OVER() as id, geom FROM p),
t1 as (SELECT id,((st_dump(geom)).geom) as geom FROM t),
t2 as (SELECT id,count(id),st_orientedenvelope(st_buffer(st_union(st_makevalid(geom)),10)) as geom FROM t1 GROUP BY id),
t3 as (SELECT st_collectionextract(unnest(st_clusterintersecting(st_accum(geom))),3) as geom FROM t2),
t4 as (SELECT row_number() over() as id, geom FROM t3),
t5 as (SELECT id,((st_dump(geom)).geom) as geom FROM t4),
t6 as (SELECT id,count(id) as nbbat,st_union(geom) as geom FROM t5 GROUP BY id),
t7 as (SELECT st_collectionextract(unnest(st_clusterintersecting(st_accum(st_buffer(geom,10)))),3) as geom FROM t6),
t8 as (SELECT row_number() over() as id2,geom FROM t7),
t9 as (SELECT id2,(st_dump(geom)).geom as geom FROM t8),
tt as (SELECT id2,t9.geom,nbbat FROM t6,t9 WHERE st_within(t6.geom,t9.geom)),
t10 as (SELECT id2,count(id2) as nbgrappe,sum(nbbat) as nbbat, st_union(geom) as geom FROM tt GROUP BY id2)
SELECT *
into buff_bati_10
FROM t10
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322050#p322050
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : T. Rossini
----------------------------------------
Peut être pourriez vous poster la requête que vous éxecuter ?
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=322034#p322034
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Christophe Vergon
----------------------------------------
Bonjour,
Je vais détailler sur ce fil, comment à partir du couple PG10 et Postgis 2.5 et de la couche bâtiment du cadastre disponible sur data.gouv.fr et intégrée avec notre outil (cf. forum données), créer une tâche urbaine un peu plus jolie et paramétrable que la méthode dite des tampons. Ça doit marcher aussi avec la couche bâtiment fournie par data.gouv.fr (non testé je vous laisse adapter).
Première requête :
-----------------
Code :
WITH p as (SELECT batiment.*,ST_OrientedEnvelope(st_buffer(the_geom,10)) as obox_geom FROM cadastre.batiment),
p1 as (SELECT st_buffer((st_dump(st_union(st_clusterintersecting(st_accum(obox_geom))))).geom,10) as geom FROM p ),
p2 as (SELECT (st_dump(st_union(st_clusterintersecting(st_accum(geom))))).geom as resgeom FROM p1 )
SELECT row_number() over() , resgeom
INTO testgeom3
FROM p2
-----------------
Détails :
-----------------
Code :
st_buffer(the_geom,10)
-----------------
Classique un buffer de 10 mètres autour de la géométrie du bâtiment.
-----------------
Code :
ST_OrientedEnvelope()
-----------------
Nouvelle fonction postgis 2.5 qui fournit la BBox orientée de la géométrie
-----------------
Code :
(st_dump(st_union(st_clusterintersecting(st_accum(obox_geom))))).geom
-----------------
Création de multipolygons constitués des BBox orientées des bâtiments qui se touchent (st_clusterintersecting()) puis fragmentation en polygons élémentaires (st_dump()).geom
Application d'un buffer de 10 mètres sur les polygones résultants
Répétition du process sur les géométries résultantes.
Exercice : Utilises la fonction recherche du forum et trouve comment compter le nombre de bâtiments originaux qui sont compris dans un polygones final, but du jeu éliminer les bâtis isolés.
Requête numéro 2 qui permet d'éliminer les trous d'une certaine surface.
-----------------
Code :
WITH p as (SELECT row_number as a,unnest((st_dumprings((st_dump (resgeom)).geom)).path) as b,((st_dumprings((st_dump (resgeom)).geom)).geom) as c
FROM testgeom3 ), -- On transforme les multipolygones à trous en n polygones sans trous si path = 0 exterior ring si >0 trous
p1 as (SELECT a,b,st_exteriorring(c) as e,st_area(c) as d FROM p order by a,b), -- on transforme les polygones ext en polyline
p2 as (SELECT a, st_accum(e) as sg FROM p1 WHERE b<>0 AND d >=1000 GROUP BY a), les trous de plus de 1000 m²
p3 as (SELECT a, e FROM p1 WHERE b=0) les enveloppes ext (shells)
SELECT p2.a,CASE WHEN not sg is null THEN st_makepolygon(e,sg) ELSE st_makepolygon(e) END
into testpoly1 FROM p3 LEFT JOIN p2 ON p2.a=p3.a ORDER BY a -- On reconstitue les polygones soit avec les trous qui restent soit sans.
-----------------
Voilà
A+
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=321968#p321968
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Ericge
----------------------------------------
Bonjour,
J'ai bien testé la 3ième version et bien spécifié l'event_id de f.
Si je décompose, la requête f me crée bien une table avec pour chaque event_id, une colonne bounds contenant entre crochet, les valeurs min et max de trac_id.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=321960#p321960
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
