Recherche minière: D'un fichier CSV vers polygones (permis miniers) par PyQGis
par
, 07/03/2021 à 10h34 (1158 Affichages)
Bonjour,
Dans la gestion minière, il arrive qu'on cherche à traiter un lot de permis miniers (carrés et/ou polygones concaves) fourni sous forme de tableau et comportant des paramètres géométriques, géographiques et administratifs. Dès lors, pour traiter le tableau il faut avoir à sa disposition un Système d'Information Géographique (SIG) et un programme conséquent pour se faciliter la tâche. Le lot ou le tableau peut contenir des dizaines voir des centaines de permis.
Le programme ou script (en fichier lié) construit sous PyQGis pour QGis 3. et postérieur va aider toute personne intéressée.
Procédure à suivre
Les données tout d'abord doivent être dans un fichier genre tableur (Excel de Microsoft ou Calc du Libre Office).
- Filtrez le fichier dans le classeur, et créez des feuilles par type de permis (permis de recherche 'PR', licence d'exploitation 'LE' et autorisation d'exploration 'AE') et par zone de projection.
Le 'PR' ici est un permis carré (4000x4000m) à positionner à l'aide d'un point pivot fourni par l'administration de tutelle (cas dans la gestion minière au Maroc). Les 'PE' et 'AE' se traitent de la même manière (à mettre ensemble dans leur projection respective).
Les tables de chaque feuille doivent respecter un ordre précis.
pour les PR:
- Colonne 1: Numéro_permis
- Colonne 2: Xpivot
- Colonne 3: Ypivot
- Colonne 4: px (est ou ouest par rapport au point pivot)
- Colonne 5: dx (distance x, centre du carré/point pivot)
- Colonne 6: py (nord ou sud par rapport au point pivot)
- Colonne 7: dy (distance y, centre du carré/point pivot)
- les colonnes suivantes sont au libre-choix.
Pour les LE et AE:
- Colonne 1: Numéro_permis
- Colonne 2: xangle (coordonnée x de chaque angle du polygone)
- Colonne 3: yangle (Coordonnée y de chaque angle du polygone)
- les colonnes suivantes sont au libre choix.
Les angles de chaque polygone doivent être ordonnés dans le sens horaire ou anti-horaire.- A la fin de chaque polygone ajoutez une ligne remplie de zéro (2ème et 3ème colonne) .
Le script en y arrivant va sauter cette ligne pour passer au traitement du polygone suivant.- Une fois le classeur reconstitué et les feuilles renommées, sauvegardez chaque feuille sous le format csv (séparateur ";" dans Excel).
- Téléchargez le script (fichier lié '*.py') qui traitera les fichiers csv.
- Dans QGis 3. ou postérieur, ouvrez la console python et chargez le fichier *.py.
Lancez le script et suivez les instructions des boîtes de dialogue.
Des exemples de fichiers csv sont à disposition (fichiers liés) pour voir comment ils sont construits. Les têtes des colonnes sont les mêmes qui se présenteront par défaut dans les dialogues du script. Il est possible de les reprendre à son compte pour faciliter l'exécution.
Fichiers liés:
- Script python (*.py)
- Fichiers csv compressés