Suite

Lire la carte vectorielle GRASS en python


Je sais que des questions similaires ont été posées plusieurs fois, mais aucune des solutions ne semble adaptée à mon problème ou ne fonctionne pas…

J'ai calculé les lignes de flux à l'aide de l'outil r.flow de GRASS GIS, puis je les exporte à l'aide de v.out.ogr qui me donne les fichiers flowlines.dbf, flowlines.prj, flowlines.shp et flowlines.shx. L'image montre la figure en GRASS de mes flowlines sur un DEM.

Maintenant, je dois lire ce fichier de formes en python et reproduire cette même image en python. Pour lire le fichier, j'ai téléchargé le module shapefile, qui semble fonctionner correctement, mais après cela, je suis perdu. Un ami m'a dit d'utiliser le module fiona mais il ne veut pas l'installer sur ma machine ubuntu, je ne sais pas quel est le problème, donc une solution sans fiona serait géniale…

Quelqu'un a-t-il de bons conseils ?


1) vous n'avez pas besoin d'exporter un fichier de formes. Avec le module Python GDAL (osgeo) vous pouvez lire directement la couche depuis la hiérarchie des dossiers GRASS (regardez la terminologie GRASS 6)

depuis osgeo import ogr # ouvrez le calque Grass shape1 ds = ogr.Open('/Users/grassdata/geol/MNT/vector/shape1/head') layer = ds.GetLayer(0)…

2) Si vous voulez lire le shapefile, vous pouvez utiliser différents modules -> regardez Est-il possible de regarder le contenu de Shapefile en utilisant Python sans licence ArcMap ?


Voir la vidéo: GRASS GIS: HOW TO CREATE A LOCATION (Octobre 2021).