Discussion :

[wissssam]

J'ai un shapefile linéaire de formes qui contient un seul élément comme vous le voyez sur l'image et je veux le diviser en plusieurs tronçon selon l'azimut et la direction de la ligne.

j'ai ce code :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
import os
import fiona
from shapely.geometry import mapping
from shapely.geometry import Point,LineString
import math
 
os.chdir('SHP')
 
def azimut(point1, point2):
       #Retourne l'azimuth de la ligne entre 2 points shapely
       angle = math.atan2(point2.x - point1.x, point2.y - point1.y)
       return math.degrees(angle)
       print angle
       if angle &>= 0: else math.degrees(angle) + 180
 
 
def paires(liste):
     '''parcourt une liste par paires'''
     for i in range(1, len(liste)):
         yield liste[i-1], liste[i]
 
with fiona.collection('BOM1.shp', 'r') as entree:
    # copie du schema de la couche et création d'un nouveau champ 'azimut'
    schema = entree.schema.copy()
    schema['properties']['azimut'] = 'int'
    # création d' une nouvelle couche avec le schéma résultant
    with fiona.collection('testligne_azim.shp', 'w', 'ESRI Shapefile', schema) as sortie:
        for ligne in entree:
        # utilisation de la fonction paire() pour extraire les segments de lignes
            for seg_start, seg_end in paires(ligne['geometry']['coordinates']):
                # création d'une ligne en fonction des points des segments
                line_start =Point(seg_start)
                line_end = Point(seg_end)
                segment = LineString([line_start.coords[0],line_end.coords[0]])
                # copie des attributs d'entrée et ajout de la valeur résultante
                # de la fonction azimut()
                elem = {}                  
                elem['properties'] = ligne['properties'] 
                elem['properties']['azimut'] = azimut(line_start, line_end)
                elem['geometry'] = mapping(segment)
                sortie.write(elem)

mais le problem du code est :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

if angle &>= 0: else math.degrees(angle) + 180

j'ai un syntax error angle &>= 0

et aussi :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
def paires(liste):
         '''parcourt une liste par paires'''
         for i in range(1, len(liste)):
             yield liste[i-1], liste[i]

j'ai une erreur de position 29: ordinal not in range(128)

[VinsS]

Salut,

L'erreur de syntaxe est évidente, non ?
&> ce n'est pas du Python mais semble tout droit sortit d'un HTML

De toutes façons ces lignes de code ne devraient pas s'exécuter puisqu'il y a un return avant

Pour la deuxième erreur, ça ne semble pas correspondre au code que tu montres, alors copie-colle plutôt le traceback complet de l'erreur.

[genemartin]

C'est apparemment ce qui arrive lorsqu'on copie du code sans comprendre. Le script provient de Python: traitement des couches vectorielles dans une perspective géologique, lecture et enregistrement des couches sous forme de dictionnaires avec le module Fiona

1. si vous connaissez Python, vous savez que >0 ne signifie rien et est une erreur (c'est aussi une erreur de transformation HTML sur le Portail qui varie suivant les navigateurs)
2. en HTML "&>" = >
3. si vous connaissez en plus les règle de géométrie angulaire (de 0 à 360°), il est facile de corriger le tir

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
def azimut(point1, point2):
       '''Retourne l'azimuth de la ligne entre 2 points shapely'''
       angle = math.atan2(point2.x - point1.x, point2.y - point1.y)
       return math.degrees(angle)if angle > 0 else math.degrees(angle) + 180

Vous pourriez tout aussi bien utiliser

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
 return math.degrees(angle) % 360

Le résultat de la deuxième partie du script est un simple générateur Python présenté sur ce Forum dans parcours d'une liste par paires, quelques solutions . Donc il ne peut pas être la cause de votre erreur.

Vous avez posé la même question sur Gis Stackexchange: Cut a linear feature according to the directions or azimuths/q/180098/2581

[wissssam]

Salut,

L'erreur de syntaxe est évidente, non ?
&> ce n'est pas du Python mais semble tout droit sortit d'un HTML

De toutes façons ces lignes de code ne devraient pas s'exécuter puisqu'il y a un return avant

Pour la deuxième erreur, ça ne semble pas correspondre au code que tu montres, alors copie-colle plutôt le traceback complet de l'erreur.

merci pour votre réponse oui c'est bien un code que j'ai trouvé , et qui n'est pas claire jarrive pas de comprendre la partie de la fonction paires
l'erreur compléte est :

Traceback (most recent call last):

File "C:\Anaconda\lib\site-packages\IPython\core\interactiveshell.py", line 3052, in run_code
self.showtraceback()

File "C:\Anaconda\lib\site-packages\IPython\core\interactiveshell.py", line 1851, in showtraceback
value, tb, tb_offset=tb_offset)

File "C:\Anaconda\lib\site-packages\IPython\core\ultratb.py", line 1240, in structured_traceback
self, etype, value, tb, tb_offset, number_of_lines_of_context)

File "C:\Anaconda\lib\site-packages\IPython\core\ultratb.py", line 1157, in structured_traceback
self, etype, value, elist, tb_offset, number_of_lines_of_context

File "C:\Anaconda\lib\site-packages\IPython\core\ultratb.py", line 511, in structured_traceback
lines = ''.join(self._format_exception_only(etype, value))

File "C:\Anaconda\lib\site-packages\IPython\core\ultratb.py", line 623, in _format_exception_only
Colors.Normal, s))

UnicodeDecodeError: 'ascii' codec can't decode byte 0xe9 in position 29: ordinal not in range(128)


In [14]:

j'ai remarqué que y'a quelque chose qui ne vas pas sur import fiona c'est a cause de sa j'ai cette erreur quand je supprime import fiona je n'est plus cette erreur