Discussion :

[preliator]

Bonjour à tous,

Sur Python, je voudrais créer un nouveau Shape à partir d'un autre, en appliquant un certains nombre de filtres. Pour cela, j'utilise la library pyshp.

Voici mon code de départ :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
import shapefile
 
 
sf = shapefile.Reader("Parcelles/parcelle.shp")
 
w = shapefile.Writer("result/test.shp", shapeType = 5, encoding="utf8")
w.fields = sf.fields[1:]
 
for shaperec in sf.iterShapeRecords():
	w.record(*shaperec.record)
	w.shape(shaperec.shape)
w.close()

De cette manière, je crée simplement une copie du shape d'origine. Après avoir testé plusieurs combinaison, je ne sais pas où placer mes conditions de filtre (par exemple, champ1 = valeur ET champ2 = valeur, etc).

Merci à vous.

[plxpy]

Bonjour

si tu veux "intervenir" sur la valeur des attributs des objets (géographiques, pas de POO ici) shapefile, c'est entre la lecture d'un objet (dans sf) et la réécriture de ses attributs(dans w) ie entre la ligne 10 et la ligne 11 que tu dois faire quelque chose.

exemple sur un shapefile de communes où je mutiplie par deux le code INSEE de la commune (aucun intérêt, juste pour l'exemple !), entier codé sous forme de chaine, premier attribut (indice 0) de chaque objet :

Code :

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

1
2
3
4
5
6
7
8
9
>>> sfin = Reader('communes.shp')
>>> sfout = Writer('new.shp', shapeType=5, encoding='utf8')        
>>> sfout.fields = sfin.fields[1:]
>>> for feature in sfin.iterShapeRecords(): 
...     feature.record[0] = str(int(feature.record[0])*2) 
...     sfout.record(*feature.record) 
...     sfout.shape(feature.shape)
...
>>> sfout.close()

affichage des deux shapefiles et comparaison des attributs d'un même objet

[preliator]

Un grand merci à toi pour ces explications