Message posté par : preliator
----------------------------------------
Bonjour,
Disposant d'une couche multilinestring représentant les routes d'un département, je souhaiterais couper cette couche en lignes de 500m.
Voici un aperçu de ma couche de routes : https://zupimages.net/viewer.php?id=20/16/eh9e.png
Après quelques recherche sur internet, je suis parti sur cette démarche :
- Regrouper ma couche de route en un seul bloc (j'ai trouvé cela cohérent car la couche de base est constitués de plusieurs lignes de tailles variables)
- Transformer ce bloc en linestring pour la requête suivante
- Créer des points tous les 500m le long de ma linestring avec st_lineinterpolatepoints (pas possible de créer les points sur ma couche regroupée, car cette fonction n'accepte que du linestring)
- Couper la couche regroupée sur cette couche de points.
-----------------
Code :
-- REGROUPER MA COUCHE DE ROUTE EN UN SEUL BLOC MULTILINESTRING
drop table if exists ligne_regroup;
create table ligne_regroup as
select St_linemerge(ST_Union(geom)) as geom
from public.route_echantillon;
-- TRANSFORMER LE BLOC EN LINESTRING
drop table if exists ligne_regroup_linestring;
CREATE TABLE ligne_regroup_linestring AS
SELECT (ST_Dump(ligne_regroup.geom)).geom::Geometry(LineString,2154) AS geom
FROM ligne_regroup;
-- CREER UN IDENTIFIANT UNIQUE SUR LA LINESTRING
alter table ligne_regroup_linestring add column id serial;
-- CREER LES POINT TOUS LES 500 METRES
drop table if exists point_500;
create table point_500 as
with tmp as(
select *
from ligne_regroup_linestring
where st_length(ligne_regroup_linestring.geom) > 500)
SELECT id, st_lineinterpolatepoints((tmp.geom),500/st_length(tmp.geom)) as geom
from tmp;
-- ON COUPE LA ROUTE REGROUPEE SUR LES POINTS
drop table if exists route_cut_500_metres;
create table route_cut_500_metres as
select ST_SNAP((ST_DUMP(st_difference(ligne_regroup.geom,point))).geom,all_point,0.1) as geom
from
(select ST_Multi(ST_Union(st_expand(point_500.geom, 0.05))) as point from point_500 ) as t1,
ligne_regroup,(select ST_MULTI(ST_COLLECT(point_500.geom))as all_point from point_500) as t2;
-----------------
Malheureusement, j'obtiens des résultats aberrants. Les lignes semblent bien se découper, mais il manque de nombreuses lignes.
Comme on peut voir sur l'exemple suivant : https://zupimages.net/viewer.php?id=20/16/oaem.png
En vert : ma couche de route regroupée (ligne_regroup)
En rouge : le résultat de "route_cut_500_metres"
Les points : Résultat de la requête qui me crée des points tous les 500 mètres.
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332169#p332169
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 : preliator
----------------------------------------
Bonjour à tous,
J'aimerais attribuer une valeur à un champ de ma couche parcelle (nommée "parcelles") si son centroide se situe dans un autre polygone représentant des espaces protégés (nommée "protege"). J'ai un problème avec la syntaxe à avoir dans ce cas.
Voici mon code :
-----------------
Code :
update parcelles p set note_esp_protege =
case
when st_within(st_centroid(p.geom), protege.geom) then 0
else -1
end
from protege
-----------------
La requête tourne en boucle depuis plusieurs minutes, il y a donc très probablement une erreur.
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329980#p329980
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 : Nicolas Ribot
----------------------------------------
Vous pouvez aussi tenter "à l'ancienne", en faisant l'union de la couche et le dump du résultat (regardez dans la doc quelle fonction st_union ou st_memUnion est la plus adaptée).
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332011#p332011
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 : preliator
----------------------------------------
Bonjour à tous,
Je travaille actuellement sur une étude liée à la position des bâtiments dans un département. Pour avoir une base de données de bâtiments la plus complète possible, j'utilise deux couches de deux bdd différente : Les bâtiments d'OSM et de datagouv. Effectivement, certains bâtis d'OSM (environ 10%) ne sont pas présents dans les bâtis de Datagouv. Mais les 90% sont des doublons des bâtis de datagouv.
Je me retrouve dans une situation aussi jubilante qu'embêtante, car les bâtis sont scannés ... en décalé !
Voici un exemple : https://zupimages.net/viewer.php?id=20/18/71wt.png
En violet, nous avons les polygones de datagouv. En gris, ceux d'OSM.
Entre deux polygones qui représentent une même maison, la superficie varie entre 1 et 2 m² mais la forme est visiblement très homogène.
Je cherche donc à supprimer les bâtis d'OSM qui sont des doublons de ceux de datagouv. Mais comment reconnaître ces doublons ? Certains d'entre eux sont simplement superposés (donc un st_intersects suffira pour les selectionner et supprimer), mais dans cette exemple et dans d'autre, les bâtis sont très espacés (plusieurs dizaines de mètres).
Existe t-il d'autres manière de trouver ces doublons ? Je pensais à des formules qui se base sur la ressemblance de la forme de polygones.
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332557#p332557
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 : MAHE_Alain (amahe2(a)externe.generali.fr)
----------------------------------------
Bonjour,
J'ai installe et configure pgpool-II, mais j'ai l'erreur suivante:
[postgres@mfle1356p log]$ psql -h vippgpoolpoc -p 9999 -U postgres
psql: ERROR: SCRAM authentication failed
DETAIL: unable to decrypt password from pool_passwd
HINT: verify the valid pool_key exists
[postgres@mfle1356p pgpool-II]$ cat pool_hba.conf
# "local" is for Unix domain socket connections only
local all all trust
# IPv4 local connections:
host all all 127.0.0.1/32 trust
host all all ::1/128 trust
host all all 10.66.13.0/0 scram-sha-256
[postgres@mfle1356p pgpool-II]$ cat pool_passwd
pgpool:AESUlhhzCC3fyJ6JPRfMQd4bg==
postgres:AESmg984I6Thqc/kTSF83nUDg==
pgpool.conf
# - Authentication -
enable_pool_hba = on
pool_passwd = 'pool_passwd'
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=331882#p331882
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 : sporito
----------------------------------------
Bonjour,
mes données sont stockées sont Postgis et je les importe sous QGIS, je souhaiterais pouvoir paramétrer un champ sous forme d'une liste déroulante.
Sous Postgis, il existe le type énumération mais lorsque je l'utilise pour un champ, une fois fois ma couche importée sous QGIS, il n'est pas conçu comme une liste déroulante.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332395#p332395
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 : Theos2000
----------------------------------------
Bonsoir, j'aurais souhaité savoir comment cacher les schéma system sous pgadmin 4 par requete ou clic bouton ? Je sais que ces schémas sont important mais ne m'en servant pas j'aurais souhaité les masquer ? Est ce qu'on peut le faire directement via l'interface graphique car je n'ai pas trouvé de doc clair sur le sujet.
Merci
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332313#p332313
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 : Theos2000
----------------------------------------
Bonjour, j''ai une colonne que j'ai converti en tableau, mais suite a une restauration de ma base ,postgres le lit comme du texte
ex :{010010000} => Text .
Or je cherche une manipulation pour dans un update le transformer en tableau => {010010000} text []
Si je refais un array dessus ca ne va pas puisque je me retrouve avec ce genre de résultat {"{920640103}"}
Est ce que quelqu'un aurai une solution ou tombé sur un cas similaire ?
Merci
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332178#p332178
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 : Ronan
----------------------------------------
Bonjour,
je suis en train de travailler à l'implémentation d'une BDD postgis gérant des nuages de points. Certains sont de type raster, d'autres vecteur.
Les nuages de points raster peuvent faire de 1 à 1.5go (il y a des milliers de nuages de points de ce type).
Je ne sais pas ce qui serait le plus pertinent :
- Stocker uniquement l'URL du fichier dans la base PostGIS et stocker les fichiers dans une arborescence windows?
- Intégrer ces données en direct dans une table de type sonde_raster?
- Autres?
Jusqu'à maintenant, j'ai eu toujours à géré les données de type raster dans des fichiers organisés dans des répertoires, mais ici vu la volumétrie, cela me questionne, mais je n'ai pas d'expérience sur le fait de stocker en direct les raster dans postgis.
Merci d'avance
Ronan
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332286#p332286
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 : Nicolas Ribot
----------------------------------------
bonjour ,
Oui tout a fait:
Un trigger se fait en 2 étapes:
1) la creation d'une fonction de trigger
2) la creation du trigger en associant cette fonction à une table et des opérations sur la table (delete, insert, update)
La même fonction pourra etre utilisée plusieurs fois, mais il faudra créer autant de trigger sur les tables que nécessaire (create trigger...)
Le nom des colonnes peut etre différent: il suffit que la fonction de trigger le gère et sache quoi faire avec quelle colonne
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332246#p332246
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