Le portail francophone de la géomatique

Bonjour !
Merci pour cette réponse motivante mais voyez vous une solution pour générer les couches Shp pour chaque commune (au nombre de 285) et sans faire 285 manipulations ?????
Merci

Dernière modification par remidici (Mon 18 December 2017 16:16)

Bonjour !
Oui je suis très attaché à ma 2CV; c'est un sujet de stage pour moi ! ;-)

Evidement ce n'est pas mon but, je désire avoir des couches avec les données attributaires le plus propre possible; je devrais trouver la solution dans une requete SQL pour bien choisir mes champs mais en attendant je n'ai rien comme idée concernant mon but de créer des couches avec le moins de manipulations possibles.
Je reste dans l'attente........et je continue à chercher ;-)

Une idée?

Bonjour,

Je félicite pas le prof ou le maître de stage qui vous demande de faire un truc pareil.

Car vous ne voyez pas la différence qui existe entre un vulgaire shp et une base de donnée.

SELECT * INTO public.shp1_com034 FROM cadastre.commune WHERE insee='034'

Puis on affiche le résultat sous Qgis et "enregistrer sous". Ou plus mieux utiliser OGR.

Bonjour et merci !
J'ai bien récupéré le code, je cherche maintenant à le réécrire afin de le faire fonctionner sur ma machine.
Pour l'instant le script m'a bien créé des dossiers et les shape correspondants.

Je poste mon code pour savoir ce que vous en pensez:
@echo on
setlocal enableextensions
set WORKSPACE=U:\Cadastre\com-054
cd /d %WORKSPACE%
for /f "tokens=*" %%A in ('dir /b /ad "%WORKSPACE%"') do (
    for /f "tokens=*" %%G in ('dir "%%~dpnA" /b /a-d') do (
        MKDIR %WORKSPACE%\%%A\TAR
        call "C:\Program Files\7-Zip\7z.exe" e %%~dpnA\%%G -o%WORKSPACE%\%%A\TAR\ -r
        MKDIR %WORKSPACE%\%%A\FEUILLES
        MKDIR %WORKSPACE%\%%A\FEUILLES\%%~nG
        call "C:\Program Files\7-Zip\7z.exe" e %WORKSPACE%\%%A\TAR\%%~nG -o%WORKSPACE%\%%A\FEUILLES\%%~nG\ -r
        MKDIR %WORKSPACE%\%%A\FEUILLES_SHP
        MKDIR %WORKSPACE%\%%A\FEUILLES_SHP\%%~nG
        call "C:\Program Files\PostgreSQL\9.5\bin\ogr2ogr" -f "ESRI Shapefile" %WORKSPACE%\%%A\FEUILLES_SHP\%%~nG %WORKSPACE%\%%A\FEUILLES\%%~nG\*.THF
        MKDIR %WORKSPACE%\%%A\COMMUNE
        call "C:\Program Files\PostgreSQL\9.5\bin\ogr2ogr" -f "ESRI shapefile" -append %WORKSPACE%\%%A\COMMUNE %WORKSPACE%\%%A\FEUILLES_SHP\%%~nG       
      )
    (
      MKDIR %WORKSPACE%\FUSION
    call "C:\Program Files\PostgreSQL\9.5\bin\ogr2ogr" -f "ESRI shapefile" -append %WORKSPACE%\FUSION %WORKSPACE%\%%A\COMMUNE
    )
)
pause
Endlocal
exit

Il me faut maintenant trouver comment faire pour sélectionner seulement les couches qui m'interessent. (Commune, parcelles, batis, sections, subdivisions fiscales, tronçons fluviaux et route)
Il me faut aussi trouver comment renommer les couches en ajoutant à la fin des noms le numéro insee de la commune ainsi que "_2017"
Et enfin, je cherche à étendre ce travail sur plusieurs dossiers en un seul script: où dois-je mettre un "/r" ?

ET problème suplémentaire, ce script ne génère pas de prj !!!

Quelqu'un a une/des idée(s)?

Dernière modification par remidici (Wed 20 December 2017 10:51)

Pour information le boulot est déjà fait en partie par Etalab, qui publie justement une version GeoJSON (facile ensuite à convertir en SHP).
C'est ici : https://cadastre.data.gouv.fr/datasets/cadastre-etalab

Par ailleurs les scripts sont ouverts : https://github.com/etalab/cadastre
Pour lire les fichiers EDIGÉO nous utilisons GDAL + un développement spécifique.

@ChristopheV vous allez très (trop) loin dans vos suppositions

Si ces données ont été produites ce n'est ni pour remplacer les données de référence, ni pour consacrer la BD Topo.
Nous avons fourni un premier travail pour faciliter la réutilisation, et vous conviendrez que le caractère juridiquement opposable ou la position des cimetières ne rentre pas dans notre définition de l'usage général.
Néanmoins ce travail n'est qu'une première étape, on poursuivra tout à long de l'année.

Oui nous utilisons la partie "OGR" de GDAL, mais pas "ogr2ogr".
De toute façon cela n'a pas d'importance puisque la prochaine version en préparation n'utilise ni l'un ni l'autre.

Enfin si les subdivisions de section sont des multi polygones, c'est parce que la documentation du standard le prévoit (page 72 => Géométrie => multi-faces).

À ta disposition pour poursuivre l'échange !

Dernière modification par jdesboeufs (Tue 02 January 2018 12:52)

Le PCI Vecteur est la donnée de référence.
Convertir toutes les couches du PCI dans des formats usuels, ça reste le PCI.

L'autre produit qu'on diffuse est en construction. Il ne peut que s'améliorer, mais pour ça on a besoin de retours et de temps.
On aurait pu attendre le produit parfait et le sortir en 2022. Mais ce n'est pas notre philosophie

Bonjour et encore merci pour vos réponses,
comme je l'ai dit plus haut, je suis un débutant, j'ai donc bien du mal à mettre en application ce que vous m'avez donné mais j'y travaille.

Ma difficulté actuelle est de convertir le fichier Json des parcelles sur un département entier (+ de 800Mo) en un fichier shape car Qgis ne veut pas le charger. J'ai aussi essayé de le faire transformer par le site "mygeodata-converter" mais il n'en finissait jamais....alors je m'inspire de ça:
https://www.geomatick.com/2017/03/02/po … a-la-volee
J'ai besoin de temps pour acquérir toutes ces compétences qui ne sont pas d'un niveau basique pour moi.

Dernière modification par remidici (Sun 07 January 2018 18:27)

Bonjour,

Le format Json possède les défauts de ses qualités. C'est un fichier texte facilement lisible par une machine ou par un humain.
Donc une géométrie si on choisi de la représenter dans un tableau d'octets pour être standard avec le Well Known Text et le Shape :
En JSON chaque valeur d'octet sera une chaîne de 1 à 3 caractères dans le fichier texte, alors un polygone avec 10 sommets cela fait tout de suite beaucoup de caractères, à lire, à convertir ... En Shape la même valeur est numérique et on code 4 valeurs pour un entier (et maintenant 8), en revanche difficile pour un humain de lire le Shape dans le texte.

En résumé la place prise et le temps de traitement d'une information (géométrique) Json sont beaucoup plus important que ceux pris pour une donnée shp.

Sans vouloir être négatif un simple constat :

Option 1 :

J'ai un format texte complexe et riche difficile à traiter : Edigéo
Je le traduit et simplifie (j'enlève de l'information) dans un format texte plus simple : Json
Je le traduit en Shp
Option : j'intègre le shp dans un bd postgis

Option 2 : je prend mes fichiers edigéo sur data.gouv.fr, avec ogr (ou le plugin cadastre hein !) j'intègre dans postgis.

Question aux chefs d'entreprises c'est la quelle qu'est économiquement la bonne ? (et ne raisonnons pas uniquement en terme de productivité,mais aussi en terme de qualité ).

Réponse cela dépend du territoire d'action de l'entreprise. ie si je suis un BE qui m'occupe d'une commune et que j'ai besoin d'un fond parcellaire uniquement pour le "graphisme" l'option 1 est la bonne.

si je veux faire de l'analyse de données à l'échelle d'un département ou d'une communauté d'agglo (volume) c'est la 2 qui est la bonne.

Bonjour,

Je crois vous avoir répondu :
Edigeo sur data.gouv
Plugin cadastre
Postgis vs Qgis

La limite étant une base de données par département.

Bonjour,

Houlala !!!!!! C'est à signaler à la communauté plugin-cadastre QGis. De façon à renouveler les tests et les comparatifs entre les diverses solutions et LA REFERENCE : cadastre.gouv.fr

Si c'est avéré ça craint ! (et ça conforte mes propos ici ou là)

Avec un peu de code en python ou autre, avec une lecture séquentielle du fichier texte et conversion unitaire dans la BD.

Bonjour,

Ok pour l'intention de réduire les manipulations lors de la création des données mais avec autant de couches créées, leur utilisation, et donc la manipulation au quotidien sera-t-elle plus aisée? Pourquoi autant saucissonner les données?
J'avoue que je lis le(s) discussion(s) sur ce sujet depuis le début et peine toujours à cerner la finalité (je dois manquer des infos de contexte peut-être)... mais bon, ce n'est peut-être que moi.

Bizarre quand même que QGIS galère autant à ouvrir et afficher un GeoJSON...

De mon côté j'arrive afficher les 1 261 344 parcelles du département 12 en une minute sur mon navigateur (Chrome 63).
https://cadastre.data.gouv.fr/data/etal … es.json.gz

À votre décharge j'ai une machine assez récente : Core i5 2,7 GHz et suffisamment de RAM.