calcul du plus grand nombre

Version imprimable

17/02/2011, 12h00
bernards111

calcul du plus grand nombre

Bonjour,

je débute un peu en python et je cherche un moyen de calculer un rang de valeur(1,2,3..) par rapport à des pourcentage dédié à un critère.
j'ai plusieurs pourcentage pour un critère et j'aimerais insérer dans un nouveau champ(mon calcul se base sur celui ci) des valeurs qui corresponde à

le % le plus élévé pour ce critère est rang1 , le second rang 2 et ainsi de suite.

Comment incrémenter cette valeur de rang pour cess différents critère?

merci de votre aide
17/02/2011, 12h13
plxpy

Je n'ai strictement rien compris à ce que tu veux faire !
17/02/2011, 12h20
bernards111

j'ai deux liste,
-une liste de métier (....)
-une liste de pourcentage
chaque métier correspond à plusieurs %
j'aimerais attribuer une valeur de rang (1,2,3..) en fonction de l'importance de ce %.
métier1 >>>(5%,10%)
10% correspond à rang : 1
5%corrspond à rang: 2
17/02/2011, 12h32
plxpy

toujours pas, même en y mettant de la bonne volonté !

C'est quoi ces "rangs" ? indice du pourcentage dans la liste des pourcentages ?
Où est l'info "tel métier est lié à tels pourcentages" ?
etc ...

les infos que tu donnes sont trop vagues et incomplètes
17/02/2011, 13h12
bernards111

rang(classement des valeurs de la plus grande à la plus petite ici ces valeurs sont les pourcentages)

ces infos sont rattachés dans une table attributaire
Une colonne métier et une colonne %
dans une 3ème colonne (rang) je cherche à insérer automatiquement ces valeurs de rang (1,2,3..)

metier % region rang
pêcheur 10 ..... 1
pêcheur 20 ...... 2
agriculteur 15 ....... 1
agriculteur 10 ..... 2
17/02/2011, 14h47
pfeuh

Salut,

Citation:

Envoyé par plxpy

Je n'ai strictement rien compris à ce que tu veux faire !

Moi non plus. C'est dommage, un problème bien énoncé est quasiment résolu.

A+

Pfeuh
17/02/2011, 14h55
valAa

edit : au vu du message #5 que je n'avais pas vu, c'est pas ce que je pensais. J'enlève mon message.
C'est vraiment pas clair.

Si tu veux classer des valeurs, regarde la fonction sorted() et ou la méthode de liste sort()
http://http://wiki.python.org/moin/HowTo/Sorting/
17/02/2011, 15h32
bernards111

1 pièce(s) jointe(s)

je devrais parler du logiciel que j'utilise avant tout: c'est arcgis
c'est un logiciel de cartographie où est inclus un module python
j'ai des tables attributaires reliés à une géométrie (cf pièces jointes)

j'ai regardé ton lien valAa, c'est pas tout à fait cela en effet
effort_maille>> c'est mes %
j'ai plusieurs % au sein d'un même code maille
j'aimerais classer ces pourcentages pour chaque code_maille par ordre de grandeur (de 1 à...)

merci
17/02/2011, 15h58
josmiley

metier % region rang
pêcheur 10 ..... 1
pêcheur 20 ...... 2
agriculteur 15 ....... 1
agriculteur 10 ..... 2

c'est quoi les critères ?
17/02/2011, 16h03
bernards111

le code_maille sur la pj et ce code correspond à une région
17/02/2011, 16h29
josmiley

posons la question autrement ...
si mériter vaut X et région vaut Y, combien vaut rang ?
17/02/2011, 17h29
bernards111

le rang se définit par rapport au pourcentage(effort_maille).
pour une région >>j'ai un métier, un pourcentage

j'ai région=[palangre,tremail,nasse]
"palangre" (un métier)= [10,5,2,7] (valeur de pourcentage)

le rang de "palangre"[10]=1
"palangre"[7]=2
....
ce n'est pas l'index que je veux obtenir, mais trier les valeurs de pourcentage et incrémenter dans le champ "rang" une valeur qui me permettrait de les classer du plus grand au plus petit(de 1 à 4 pour le palangre).

Je crois comprendre, mais pour une réponse plus complète, j'ai une petite question.

Pour une entrée (cf ta pièce jointe) qui a cette allure :

Code:

1
2
3
4
5
6
 
Métier    Effort maille    Maille    Rang
..........                 101
palangre   0.9             102        0
tremail    0.5             102        0
..........                 103

tu souhaites obtenir ça (les rangs sont calculés pour une maille) :

Code:

1
2
3
4
5
6
 
Métier    Effort maille    Maille    Rang
..........                 101
palangre   0.9             102        1
tremail    0.5             102        2
..........                 103

ou ça (les rangs sont calculés pour les couples (maille,métier)) :

Code:

1
2
3
4
5
6
 
Métier    Effort maille    Maille    Rang
..........                 101
palangre   0.9             102        1
tremail    0.5             102        1
..........                 103

17/02/2011, 20h08
bernards111

oui c'est la première proposition
"tu souhaites obtenir ça (les rangs sont calculés pour une maille) :"

suis désolé ,j'était tellement dans les tuto python toute la journée que j'ai le cerveau qu a cramé
en me relisant8O (aprés une bonne pause) je comprend a quelle point je n'était pas très clair:oops:

merci
18/02/2011, 06h40
plxpy

Il y a 2 aspects dans ce que tu veux faire.

Tu dois :

1 - causer avec ArcGIS : lire/récupérer tes objets (objets géographiques, pas objets au sens POO) et les "réécrire" une fois que tu auras déterminé ces fameux "rangs"

2 - mettre à jour les objets géographiques en calculant leur "rang"

Côté ArcGIS, tu dois te référer à la documentation fournie pour savoir comment et sous quelle forme tu récupères tes objets.

Je connais un peu ArcGIS, pas trop l'interface Python proposée, mais je ne pense pas prendre beaucoup de risque en imaginant qu'une fonction te permettra d'effectuer une requête (en précisant la classe d'objets voulue, telle ou telle condition sur tel ou tel attribut, une zone géographique, etc...) et de récupérer tes objets (une liste d'objets).

Quelle forme (Python) auront-ils ? des instances d'une classe "ArcGISFeature" ? des tuples ou des dictionnaires regroupant la géométrie et les attributs de chaque objet ? autre chose ? Là aussi, réfère-toi à la doc.

Côté "Python pur", à partir de ta liste (unique) d'objets géographiques, il faut que :

- tu construises autant de listes qu'il y a de mailles (tu classes les "métiers" dans chacune des mailles) avec les objets géographiques localisés dans une même maille
- tu appliques un tri basé la valeur de "effort maille" (tri inversé car tu veux en premier la plus forte valeur) voir cette discussion récente
- dans chaque liste, tu affectes le rang de chaque objet en fonction de sa position dans la liste triée (1 pour le premier, 2 pour le deuxième, etc...)

Après tu dois "faire prendre en compte" les changements effectués sur tes objets à ton "projet" ArcGIS (écriture ? sorte de "commit" ? voir la doc).

Sans connaître finement l'interface Python proposée par ArcGIS (les modules, les classes, les fonctions), je ne peux pas être plus précis.

ps : à lire ce billet concernant les SIG (dont ArcGIS) et la façon dont ils (mal)traitent Python
18/02/2011, 10h57
bernards111

merci de ta réponse!
Oui effectivement les différences sur python pour les sig est dommage car je ne sais pas quel tuto aborder.
Mais bon j'ai peut être trouvé deux fonctions qui peuvent répondre à mon problème

http://help.arcgis.com/fr/arcgisdesk...003m000000.htm
et
http://help.arcgis.com/fr/arcgisdesk...0000039000000/

Autant tirer parti du fait que l'on peut demander un curseur trié.
Quelque chose comme ceci conviendrait, il me semble:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 # ouvrir un UpdateCursor, reprenant les champs code_maille, effort_maille et rang, trié sur code_maille et sur effort_maille descendant rows = arcpy.UpdateCursor(...) maille_courante = None for row in rows: if maille_courante == row.code_maille: rank += 1 else: rank = 1 row.rang = rank rows.updateRow(row) del row del rows
Cela vaudrait sans doute la peine d'implémenter un context_manager pour rendre un peu plus propre la gestion du curseur...

21/02/2011, 09h05
bernards111

ok merci!

tu peux m'expliquer ce qu'est un context manager? stp

Attention il faut un peu s'accrocher ;)

Le problème est que le code ci-dessus n'est pas sûr. En me référant à la doc, je lis qu'un verrou exclusif est mis sur la table quand un UpdateCursor est ouvert. Si une exception se produit, par exemple lors de la mise à jour d'une ligne, le "del row; del rows" n'est pas exécuté et le verrou n'est pas libéré, ce qui rend la table inutilisable (jusqu'à la fin du processus).
La seule façon que je vois pour justifier ce que dit la doc est qu'une méthode __del__ est définie sur les classes Cursor et Row. En soi c'est déjà un choix discutable...

Un meilleur code est celui-ci:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 # ouvrir un UpdateCursor, reprenant les champs code_maille, effort_maille et rang, trié sur code_maille et sur effort_maille descendant rows = arcpy.UpdateCursor(...) try: row = None maille_courante = None for row in rows: if maille_courante == row.code_maille: rank += 1 else: rank = 1 row.rang = rank rows.updateRow(row) finally: if row: del row del rows
Le try/finally permet de s'assurer que le code dans la partie 'finally' est toujours exécuté, même si l'exécution est interrompue par une exception.

Mais écrire ce try/finally à chaque fois que l'on veut manipuler la DB est un peu barbant et facile à oublier. Un context manager permet de simplifier cela:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 from contextlib import contextmanager @contextmanager def update_cursor(*args): r = arcpy.UpdateCursor(*args) yield r del r # ouvrir un UpdateCursor, reprenant les champs code_maille, effort_maille et rang, trié sur code_maille et sur effort_maille descendant with update_cursor(...) as rows: row = None maille_courante = None try: for row in rows: if maille_courante == row.code_maille: rank += 1 else: rank = 1 row.rang = rank rows.updateRow(row) finally: if row: del row
On a remplacé l'appel direct à UpdateCursor par un appel à update_cursor. Dans update_cursor, la partie avant le yield est exécutée avant le corps du "with", le yield passe la main au corps du "with" et la partie après le yield correspond au "finally".

Mais on n'a que peu gagné car ici c'est relativement complexe; d'après la doc, il faut s'assurer de libérer à la fois le curseur (rows) et la ligne actuelle (row), et le context manager ne s'occupe que du curseur lui-même. Il ne peut en effet pas s'occuper de "row" vu que celle-ci n'est instanciée que plus tard et n'est pas encore connue. Il reste donc un try/finally pour s'occuper de "row".

Pour améliorer les choses, il faut donc se passer du décorateur (@contextmanager) et écrire un gestionnaire de contexte "à la main", qui sera aussi un proxy pour le curseur afin de protéger les "rows":
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 class SafeUpdateCursor(object): def __init__(self, *args): self.cur = arcpy.UpdateCursor(*args) self.row = None def __iter__(self): return self def __enter__(self): return self def __exit__(self, exc_type, exc_val, exc_tb): if self.row: del self.row del self.cur return False def next(self): self.row = self.cur.next() return self.row def __getattr__(self, name): # on délègue les autres méthodes au curseur return getattr(self.cur, name) # ouvrir un UpdateCursor, reprenant les champs code_maille, effort_maille et rang, trié sur code_maille et sur effort_maille descendant with SafeUpdateCursor(...) as rows: maille_courante = None for row in rows: if maille_courante == row.code_maille: rank += 1 else: rank = 1 row.rang = rank rows.updateRow(row)
Ce n'est pas simple je l'admets... La classe SafeUpdateCursor remplace la fonction UpdateCursor. Les méthodes __enter__ et __exit__ font de cette classe un gestionnaire de contexte qui peut être utilisé avec "with". Je ne tenterais pas d'expliquer plus en détail son fonctionnement, d'autant plus que je ne suis pas certain qu'il n'y ait pas d'erreur vu que je n'ai pas testé ce code (je n'ai pas arcgis). Une fois que tout est en place, tu peux remarquer que le code qui manipule le curseur est simplifié vu qu'il ne faut plus s'occuper explicitement ni des exceptions ni de libérer le curseur ou la ligne.

22/02/2011, 18h04
bernards111

Je te remercie pour le temps que tu m'a accordé, d'une simple lecture je ne comprend pas tout (je debute8O ) mais je vais essayer de le décortiquer et d'y voir plus clair demain.

merci:ccool:
28/02/2011, 09h42
bernards111

Re

j'ai essayé ton code et il m'indique une erreur sur la fonction (<type 'exceptions.RuntimeError'>: ERROR 999999: Erreur lors de l’exécution de la fonction.
Échec de l’exécution de (Script).)

dans ton titre tu mettais "trié sur effort maille et code maille mais je ne comprend pas la partie du code qui fait cela ? sinon je continue de parcourir la doc arcgis...

merci pour votre aide

j'avance petit a petit mais je n'obtiens pas ce que je veux

le tri des champs ce fait sur
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 import arcpy rows=arcpy.UpdateCursor"C:/stations/Nouveaudossier/Export_Output.shp","","","ID_POLY;RATIO;CALCUL","ID_POLY D;RATIO D") #D pour tri descendant >>> maille= None >>> for row in rows : ... while maille==row.ID_POLY: ... rank+=1 ... else: ... rank=1 ... row.CALCUL=rank ... rows.updateRow (row) ... del row ... del rows
avec ce code je n'obtiens que des 1 sur toute la colonne, je n'arrive pas à incrémenter la valeur de 1 tant que les codes maille sont identiques?

merci de votre aide

04/03/2011, 00h36
dividee

Citation:

Envoyé par bernards111

dans ton titre tu mettais "trié sur effort maille et code maille mais je ne comprend pas la partie du code qui fait cela ?

Ben j'avais écrit (...) à remplacer par les bons arguments... Comme je l'ai dit je n'ai pas arcgis et je n'ai aucune intention de l'acheter. La documentation seule m'indique que ce que j'ai mis en commentaire doit être possible en passant les bons arguments, mais est insuffisante pour déterminer la syntaxe exacte de ces arguments sans un peu d'expérimentation.

Le nom des champs que laisse entrevoir ton code ne correspond pas avec le nom des champs cités précédemment, difficile pour moi de faire le rapprochement.

Est-ce que ID_POLY correspond à code_maille ? Si tu dis que le résultat est toujours 1, cela indiquerait que les ID_POLY sont tous différents, ce qui ne correspond pas à code_maille dans tes posts précédents.

Ah et c'est un "if", pas un "while", sinon tu vas entrer dans une boucle infinie.

je reviens pour mettre le code qui marche

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
# tri  descendant code_maille a et effort_maille d
import arcpy

def rang ():

rows= arcpy.UpdateCursor("C:/benoit/Export_Output.shp","","","ID_MAILLE; Ratio ; RANG","ID_MAILLE D; RATIO D")
maille_courante=None
for row in rows:
    if maille_courante == row.ID_MAILLE:
        rank+=1
    else:
        rank=1
        maille_courante=row.ID_MAILLE
    row.rang = rank
    rows.updateRow(row)

__esri_field_calculator_splitter__
def rang ()

merci à tous pour votre aide et en particulier à dividee