Le portail francophone de la géomatique

Bonsoir FloSnow,
Merci de votre réponse.
Entre temps j'ai trouvé ma réponse qui se résume a ceci :

# Après avoir effectué la polygonisation via

gdal_polygonize.py input.tiff output.shp -b 1 -f "ESRI Shapefile" tileName DN

# Je répare la polygonisation via "Réparer les géométries" sous QGIS

# Je sélectionne la couche qui m'intéresse (0 = terre) pour le futur Export kml via :

ogr2ogr output.shp input.shp -where "DN=0"

# Je converti les polygones issus de la commande ci dessus en ligne via QGIS

# J'exporte en KML via :

ogr2ogr -f KML output.kml input.shp

Avec gdal/ogr je peux effectuer la quasi totalité des conversions en partant d'un png non géoréférencé pour en obtenir un KML des terres (512 tuiles de 32 x 32 tuiles de 256px)
Seuls les points "Réparer les géométries" et "Polygones en Lignes" sont effectués avec QGIS ce qui me permet sous Mac de faire un script Bash qui m'effectue une majorité des commandes, pour la partie QGIS j'utilise le traitement par lot mais j'ai toujours des soucis sur le nom des fichiers que je dois retravailler après.

Si des solutions existent pour ces 2 points, je suis preneur (la première je suis presque sur que ce n'est pas possible, la seconde en revanche doit être possible mais pas trouvé) et pour la question, non je le fais pas en Python parce que je ne connais pas du tout ce langage (à part des choses hyper basiques)

Sur 512 tuiles en KML, 4 sont en "vrac" je ne sais pas trop pourquoi sachant que les shp sont ok, j'ai une hypothèse sur le nombre de lignes trop nombreuses ...
Pour ça je vais "découper" mes tuiles "en vrac" de 32 x 32 en 4 x (16 x 16)

Dernière modification par GeGaX (Mon 13 January 2020 21:09)

Bonjour GeGaX.

Intéressant pour la limitation lors de la création des KML, comme quoi le volume des données compte en effet beaucoup. 


- Pour l'étape "Réparer les polygones" :

Après quelques recherches, il ne semble pas y avoir encore de moyen d'utiliser cet outil de traitement directement sur console. J'ai trouvé une discussion pouvant apporter une solution dans l'avenir avec les nouvelles versions de QGIS mais sans ogr/gdal : https://github.com/qgis/QGIS-Enhancemen … issues/140

Au cas où, il est toujours possible de partir d'un script .py (script réalisé à partir du topic : https://gis.stackexchange.com/questions … processing) :

import sys

from qgis.core import QgsApplication
from qgis.analysis import QgsNativeAlgorithms

# Chemin des modules python
sys.path.append(r"C:\Program Files\QGIS 3.10\apps\qgis\python\plugins")

# Initialisation de QGIS, sans l'interface
qgs = QgsApplication([], False)
qgs.initQgis()

# Import et initialisation des outils de traitement
import processing
from processing.core.Processing import Processing
Processing.initialize()

# Cette ligne permet de prendre en compte les outils "natifs" de QGIS, dont fait partie "Réparer les géométries"
QgsApplication.processingRegistry().addProvider(QgsNativeAlgorithms())

# Lancement du traitement "Réparer les géométries"
processing.run("native:fixgeometries", {'INPUT':'input.shp','OUTPUT':'output.shp'})

Puis de le faire tourner en insérant son chemin dans la console (en partant du OSGeo4W Shell) :

python-qgis script.py

- Pour l'étape "Polygones vers Lignes" :

Il est possible d'extraire uniquement les contours des géométries, et donc de créer un shape de lignes, en utilisant une requête SQL via ogr2ogr :

ogr2ogr output.shp input.shp -dialect sqlite -f "ESRI Shapefile" -sql "select ST_ExteriorRing(geometry) as geometry, * from input"

"-dialect sqlite" semble être obligatoire pour permettre à la requête de fonctionner.
Il risque cependant d'y avoir des erreurs de géométries s'il y a des entités multi-parties avec des trous.


Cordialement.