Diviser un élément linéaire selon les directions ou azimuts de la ligne

Version imprimable

1 pièce(s) jointe(s)

Citation:

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.

Pièce jointe 200405

Citation:

j'ai ce code :

Code:

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 &amp;&gt;= 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)

Citation:

mais le problem du code est :

Code:

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

Citation:

j'ai un syntax error angle &>= 0

Citation:

et aussi :

Code:

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]

Citation:

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

10/02/2016, 08h17
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.

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

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

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:

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

10/02/2016, 20h50
wissssam

Citation:

Envoyé par 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.

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 :

Citation:

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]:

Citation:

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