Message posté par : T. Rossini
----------------------------------------
Un autre moyen sale :
-----------------
Code :
SELECT
('POINT' || replace(geoloc_etablissement, ',', ' '))::geometry('POINT')
FROM
sirene
;
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329433#p329433
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
----------------------------------------
-----------------
Nicolas Ribot écrit :
Bonjour,
Un moyen bien crad' (virer les parenthèses et fabriquer un tableau de valeurs en prenant ',' comme séparateur):
-----------------
Code :
with tmp as (
select string_to_array(replace(replace(geoloc_etablissement, '(', ''), ')', ''), ',') as coords
from sirene
) select coords[1] as latitude, coords[2] as longitude;
-----------------
Pour fabriquer un point, il suffira de faire: st_setSRID(st_makePoint(coords[2], coords[1]), 4326)
-----------------
Petite erreur: faut effectivement caster en float:
-----------------
Code :
st_setSRID(st_makePoint((coords[2])::float, (coords[1])::float), 4326)
-----------------
Nico
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329432#p329432
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
----------------------------------------
Avec une regexp:
-----------------
Code :
with tmp as (
select regexp_matches(geoloc_etablissement, '[-+]?[0-9]*\.?[0-9]+', 'g') as coords
from sirene
) select coords[1] as latitude, coords[2] as longitude;
-----------------
(je ne comprends rien à cette regexp, je l'ai trouvée sur le net :D )
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329431#p329431
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 : Djperou
----------------------------------------
Bonjour,
J'imagine que la colonne est en texte. Du coup tu peux utiliser split_part(<texte>,<délimiteur>,<colonne à renvoyer>) pour retrouver tes deux colonnes :
-----------------
Code :
SELECT split_part( "geolocalisation de établissement" , ',' , 1 ) as latitude,
split_part( "geolocalisation de établissement" , ',' , 2 ) as longitude
FROM sirene_v3
-----------------
Je ne sais pas comment tu veux utiliser cette info mais tu devras peut être utiliser la fonction CAST pour déclarer en real tes datas, ou encore utiliser ST_MakePoint() si tu veux les géométries PostGis de ces points.
J'espère que ça t'aidera
Bonne journée,
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329430#p329430
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,
Un moyen bien crad' (virer les parenthèses et fabriquer un tableau de valeur en prenant ',' comme séparateur):
-----------------
Code :
with tmp as (
select string_to_array(replace(replace(geoloc_etablissement, '(', ''), ')', ''), ',') as coords
from sirene
) select coords[1] as latitude, coords[2] as longitude;
-----------------
Pour fabriquer un point, il suffira de faire: st_setSRID(st_makePoint(coords[2], coords[1]), 4326)
Ca serait plus propre avec une expression régulière mais je n'y arrive pas :)
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329428#p329428
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 a tout le monde et une bonne santé a tous les membres.
Je cherche dans la base sirene V3 issu de open data soft a construire un point a partir du x et Y qui sont dans la même colonne "geolocalisation de établissement" (48.539364,2.653041). Je cherche a construire 2 champs distincts par un substring qui s'arreterait a la ',' car le nombre de caractères change a chaque fois...auriez vous une idée sous postgres.
Merci a vous
----------------------------------------
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=329421#p329421
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 : Mossy_D
----------------------------------------
Bonjour,
Merci pour votre participation.
Je reviens vers vous avec une bonne nouvelle puisqu'à nouveau j'ai rencontré ce problème d’auto-incrémentation partitionnée selon un certain critère... et cet fois je l'ai bel et bien résolu. La première solution était en fait la bonne. Toutefois, seule la requête SELECT fonctionne. Ensuite il est possible de générer une couche avec DB MANAGER. Si l'on fait UPDATE, puis SELECT, le code s'exécute mais le résultat est faussé. Là je ne saurais dire pourquoi... les "Windows fonctions" ne sont peut être pas conçues pour fonctionner avec UPDATE..??
Bref, voici le code donnant un résultat satisfaisant :
-----------------
Code :
SELECT *, ROW_NUMBER ()
OVER (PARTITION BY parcelle) AS no_bati
FROM bati
ORDER BY parcelle
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329262#p329262
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 : zakarya2626
----------------------------------------
Bonjour tout le monde ,
L’INSEE ne diffuse plus la données "natetab" concernant les personnes physiques . Les personnes physiques correspondent au code 1000 dans le champs 'categoriejuridiqueunitelegale' dans la nouvelle base de données SIRENE V3.
Ma question est la suivante : Y-a-t-il un moyen pour reconstituer le champs "NATETAB" et avoir en détail la nature des personnes physiques ( artisans , commerçant , profession libérale ..)?
Merci à vous
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329257#p329257
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 : Brindavoine
----------------------------------------
Aloha
A première vue votre count(*) porte sur les deux tables rpg et rivieres. Si vous faites un count(rpg.*) cela devrait marcher.
Vous pouvez aussi tester quelque chose comme
-----------------
Code :
select s.libelle_c, rpg.geom
from
(select libelle_c from rpg, riviers riv
where st_intersects(st_buffer(riv.geom,50), rpg.geom) and riv.layer = 'Ondaine_et_affluents') s inner join rpg on s.libelle_c=rpg.libelle_c
group by s.libelle_c, rpg.geom
order by s.libelle_c
-----------------
En passant par une sous requête ou mieux une CTE, vous devriez épurer les doublons. J'ai ajouté les geom dans les colonnes pour les récupérer.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329231#p329231
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,
Je souhaiterais déterminer le type des parcelles qui se font intersecter par des rivières. Avec la formule que j'applique, je me rend compte qu'il me duplique (ou triple) les géométrie. Lorsque j'exporte le résultat sur QGis, je suis obligé de "Supprimer les géométries dupliquées" pour avoir le bon résultat.
J'utilise cette formule :
select count(*)
from rpg, rivieres riv
where st_intersects(st_buffer(riv.geom,50), rpg.geom) and riv.layer = 'Ondaine_et_affluents'
Qu'ai-je oublié ?
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=329223#p329223
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
