Le portail francophone de la géomatique

Bonjour à tous !

Je suis actuellement et ce depuis trois jours en train de m’arracher les cheveux sur un script python. Je serai le plus heureux des étudiants sigistes si quelqu'un parvenait à m'éclairer, et ça laisserai mes cheveux tranquille.

J'utilise ArcGIS 10.0 et j'ai écrit un script pour automatiser la création de polygones de voronoï. Ce traitement découpe autour de points appelées  « germes » une mosaïque de polygones ayant ses frontières à la médiane des distances qu’il a avec les autres points frontière (http://help.arcgis.com/fr/arcgisdesktop … 001m000000).
Je travaille sur la modélisation de landmark, je dois donc appliquer ce traitement à une couche de landmark officiel, type parmi d’autres. J'ai mis au point un petit algorithme permettant de vérifier dans les autres couches de landmarks s’il y en aurait d'autres qui pourraient être inséré dans le fichier des landmark officiels.

J’ai modélisé le problème et ai programmé différentes fonction, celle qui compte les landmarks me pose un problème. Cette fonction a pour but de compter le nombre de landmarks identiques aux germes des polygones de voronoi. Elle devra donc ajouter +1 au compteur du landmark_germe autant de fois que le même landmark est repéré dans le polygone. Si elle détecte de nouveaux landmarks, la fonction appellera une autre fonction pour entrer ce point  dans une table de tri.


Le code :

def nb_lm(lm_dataset):

#Nb d’utilisation des landmarks des tables LM sugg et LM user, sous contraintes spatiales

    rows_voro=arcpy.UpdateCursor(lm_voro)
    for row_voro in rows_voro:
        dat = str(row_voro.Input_FID)
        arcpy.SelectLayerByAttribute_management(lm_voro,"NEW_SELECTION",'"Input_FID" ='+dat)
        arcpy.SelectLayerByLocation_management(lm_dataset, "WITHIN", lm_voro, autour, "NEW_SELECTION")
        print "fait"
        arcpy.CopyFeatures_management(lm_dataset,"temp")
        if int(arcpy.GetCount_management("temp").getOutput(0)) == int(arcpy.GetCount_management(lm_dataset).getOutput(0)) : #si la selection est nulle, landmark pas a supprmier
            print("rien dans les parages pour"+row_voro.LMNOM)
            row_voro.suppr = "non"
            rows_voro.updateRow(row_voro)
        
        else:
            print "trouvé"

#suite du traitement

J’ai cette erreurs: Runtime error <class 'arcgisscripting.ExecuteError'>: ERROR 000358: Expression non valide Échec de l’exécution de (SelectLayerByAttribute)

SQL donc.. Pourtant quand je lance cette ligne toute seul dans l’interpréteur en mettant a dat une valeur, ça marche. Et quand je print dat avant ce traitement, c’est bien un chiffre de type str. Je comprends pas pourquoi ça marche pas dans le code. J’ai essayé de tourner le truc dans tous les sens, et suis sûr que c’est un truc tout con. Ou c'est ArcGIS. Ou mes compétences restreintes en programmation.


Quelqu’un voit?


Merci pour votre aide!