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 : thomas.bagot (thomas.bagot(a)atmo-bfc.org)
----------------------------------------
Bonjour,
Je me pose une question dans la gestion des base de données. J'ai plusieurs base de données thématiques qui utilisent toutes des données géographiques parfois les mêmes ou parfois différentes données géographiques. Je me demande si je dois créer un schéma de données géographiques par base de données, ou alors une base de données géographiques pour toutes les bases de données ?
En vous remerciant d'avance, Thomas Bagot
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=333047#p333047
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 : Corossol Wali
----------------------------------------
Bonjour à tous,
J'ai commandé à une boîte privée de technicien et de géomaticien une "migration de ma base de donnée logée sous Access". La base de donnée a été créée à partir d'ACCESS. Pour des problms de bugs, je souhaite maintenant la migré vers une extension directe de bas de donnée de QGis afin de pouvoir effectuer les mise à jour, modifications...sous QGis et sa base de donnée directe.
C'est une base de donnée qui est alimentée en permanence avec les infrastructures publics de notre collectivité. Est-ce que la migration peut prendre du temps? parce-que j'ai fourni ma bdonnée depuis maintenant 6 ou 7 mois.
Est-ce un travail fastidieux?
Dans l'attente...merci de m'aiderrrrr!!!!
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=333190#p333190
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 : aguenec
----------------------------------------
Bonjour,
Pour pouvoir les extraire au format shape, il faut d'abord restaurer votre "backup" au sein d'une base PostGRES.
Le site https://datafoncier.cerema.fr/ressources/installer-une-base-donnees-fichier… est plutôt bien fait et décrit pas à pas la restauration.
Vous pourrez ensuite les exploiter ou les extraire au format souhaité depuis QGIS.
Cordialement
Antoine
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332933#p332933
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 : Elise Henrot
----------------------------------------
Bonjour,
Je reviens sur la question de François. J'ai également demandé les fichiers fonciers via le site "démarche simplifiée", et les données transmises sont au format backup (avec les DDT et DREAL on obtenait simplement les shp et tout allait bien). L'interlocuteur m'a précisé avoir "testé avec succès l’intégration des .backup dans une version PostgreSQL 9.5.10". Mais je suis géographe urbaniste, pas géomaticienne ou informaticienne. J'ai bien installé PostgreSQL 9.5.10, mais je n'y comprend rien. Y'aurait-il un tutoriel pas à pas pour pouvoir extraire les fichiers shp à partir de ces données "backup".
Bien cordialement.
Elise
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332925#p332925
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 dispose d'une table points "PointA" (500 000 points), et d'une seconde table de points "PointB" (73 000 points). Je voudrais faire la somme d'une colonne quantitative située dans "PointB" des points qui se situent à une distance fixe de chacun des points de "PointA". Pour cela, utiliser la clause St_dwithin me paraît être cohérent.
-----------------
Code :
select pta.id, sum(ptb.ind) somme_pop
from pointa pta, pointb ptb
where st_dwithin(pta.geom, ptb.geom, 10000)
group by pta.id
-----------------
Malheureusement, la requête dure plus d'une heure pour s'achever (malgré les index spatiaux).
Existe t-il une syntaxe qui me permettrait d'aller plus vite ?
Merci.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=333019#p333019
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 : uterpandragon (alain.terieur0806(a)gmail.com)
----------------------------------------
Bonjour,
Débutant sous PGSQL/Postgis, je cherche à sélectionner des communes traversées par des routes.
J'ai donc essayé ceci :
-----------------
Code :
SELECT st_intersection(communes.geom, routes.geom)
FROM public.communes, public.routes;
-----------------
Cela me renvoie une liste de géométrie, par ailleurs plus importante que je ne devrais avoir de communes sélectionnées.
Je souhaiterais avoir la liste des communes avec leur nom et pas une colonne avec la géométrie en hexadécimale apparemment.
Comment m'y prendre svp ?
Merci par avance, Cordialement.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332962#p332962
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 : white-shadow90
----------------------------------------
Bonjour à tous,
Je cherche à créer un polygone (sans trou et d'un seul tenant) qui regroupe X pourcents de points issus d'une même table.
Le polygone à créer doit regrouper les points les plus proches ENTRE EUX. En conséquence, je ne vois pas comment utiliser les fonctions dérivées de st_cluster telles que ST_ClusterDBSCAN ou ST_ClusterKMeans car ma contrainte n'est pas une distance en tant que telle ou un nombre de clusters.
Cela me fait penser à une classe issue d'une classification ascendante hiérarchique mais je sèche complètement pour la requête à effectuer.
Auriez-vous des pistes à me donner svp ? Je vous remercie par avance pour votre aide.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332832#p332832
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.
----------------------------------------
Re bonjour,
Suite à ce sujet (https://georezo.net/forum/viewtopic.php?pid=332590) qui ne concernait que la problématique des dates EXIF, je me permet de faire remonter mon appel à l'aide sur la communauté des utilisateurs de l'excellent programme ExifTool : https://exiftool.org
Cet outil permet de manipuler les métadonnées (EXIF, XMP, ...) des photos, pour les extraire, les mettre à jour, ...
Il m'a d'ors-et-déjà permis de générer une table de base de 120000 photos géolocalisées dans PostGre, et ceci en un temps tout à fait raisonnable (45 minutes de scan).
Cependant, je n'ai pour l'instant récupéré que quelques métadonnées de base (date, modèle d'appareil photo, chemin, coordonnées GPS, ...), via cette commande, qui pourra peut-être en intéresser certains :
-----------------
Code :
exiftool -T -r -Directory -filename -DateTimeOriginal -model -GPSLatitude# -GPSLongitude# -ext jpg -charset FileName=Latin "M:\Photos" > "M:\Photos\dbphotos.txt"
-----------------
Phil Harvey, le créateur d'Exiftool, m'a suggéré une méthode qui créé une base SQLite, mais ses scripts sont en Perl, et je n'aurai pas la compétence pour les adapter à mes besoins. Sa solution est en réponse à ma question sur ce fil (https://exiftool.org/forum/index.php?topic=11087.0) .
Peut-être qu'un autre géomaticien photographe plus compétent que moi serait en mesure de créer un script pour récupérer toutes les métadonnées d'une photothèque dans une base PostGreSQL ?
C'est en tout cas le but de l'ouverture de ce fil !
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=332591#p332591
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