Discussion :

[rabeh.ram]

j'ai obtenu un bon résultat lorsque j'exécute ce code dans le cas de coupe une ficher shapefile.shp de type LIGNE par rectangle (xmin, ymin, xmax, ymax) mais je n'ai trouvé aucun résultat dans le cas où le type de shapefile.shp est zone fermé POLYGONE
pour cela veuillez m'aide pour rectifier ce code ou un autre code pour résoudre le problème

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
42
import shapefile
# Create a reader instance for our US Roads shapefile
 
# Create a writer instance copying the reader's shapefile type
save = "clip"
r  = shapefile.Reader("Myfile.shp")
w = shapefile.Writer(save, r.shapeType)
 
# Copy the database fields to the writer
w.fields = list(r.fields)
 
# Our selection box that contains Puerto Rico
#personal definition of coordinates (top left and bottom right coords)
 
xmin = -67.5
xmax = -65.0
ymin = 17.8
ymax = 18.6
 
# Iterate through the shapes and attributes at the same time
for road in r.iterShapeRecords():
    # Shape geometry
    geom = road.shape
    # Database attributes
    rec = road.record
    # Get the bounding box of the shape (a single road)
    sxmin, symin, sxmax, symax = geom.bbox
    # Compare it to our Puerto Rico bounding box.
    # go to the next road as soon as a coordinate is outside the box
    if sxmin < xmin:
        continue
    elif sxmax > xmax:
        continue
    elif symin < ymin:
        continue
    elif symax > ymax:
        continue
    # Road is inside our selection box.
    # Add it to the new shapefile
 
    w.shape(geom)
    w.record(*rec)

[plxpy]

Bonjour

nb: je ne reviens pas sur le détail du code et la façon d'organiser les choses...

Rien n'empêche ton code de fonctionner correctement sur un shapefile de type polygone. Si, à l'exécution, le shapefile créé (clip.shp, .shx et .dbf) est vide, c'est qu'aucun polygone en entrée n'est à l'intérieur de la bounding box que tu as définie.

Sur les copies d'écran, les polygones en entrée sont en violet. Une fois ton code exécuté, le polygone à l'intérieur de la bounding box (adapté ici à mes objets) est bien conservé dans le shapefile résultat (en rouge)

[rabeh.ram]

Bonjour

Merci plxpy de m'aider

Rien n'empêche ton code de fonctionner correctement sur un shapefile de type polygone. Si, à l'exécution, le shapefile créé (clip.shp, .shx et .dbf) est vide, c'est qu'aucun polygone en entrée n'est à l'intérieur de la bounding box que tu as définie.

j'ai testé mon code plusieurs fois sur des zones dans lequel le polygone coupeur est bien chevauché avec le shapefile qui être coupé mais le resultat elle est toujours un fichier .shp vide

-- vous pouvez partager avec moi les fichiers .shp qui vous avez utilisé dans le test de code ?
-- quel est le format des coordonnées a utilisé dans votre fichiers .shp (UTM / Latitude-Longitude) ?

[plxpy]

j'ai testé mon code plusieurs fois sur des zones dans lequel le polygone coupeur est bien chevauché avec le shapefile qui être coupé mais le resultat elle est toujours un fichier .shp vide

Ton code ne retient que les objets dont le rectangle englobant (bounding box) est complètement à l'intérieur de ta zone rectangulaire (autre bounding box) définie par (xmin, ymin, xmax, ymax) aux lignes 15 à 18.

Si ce n'est pas ce que tu cherches à faire ... Exécuter ton script sur tout un tas de shapefiles en entrée n'y changera rien : il continuera à faire la même chose.

Pour des sélections d'objets un peu plus pointues, je te conseille d'utiliser des packages comme Shapely.