Suite

De featureclass (ou shp) à WKT e GeoJSON


J'utilise ArcGIS 10.1 sans extension Data Interoperability et je souhaite convertir certaines classes d'entités (points, lignes et polygones) contenues dans un fichier .mdb. J'ai lu quelques articles mais je n'arrive pas à convertir mon fichier.

Voici le code que j'ai (qui ne fonctionne pas) :

import arcpy cur=arcpy.da.SearchCursor("CESTINI",["[email protected]"]) pour la ligne dans cur : x,y = row[0] print ("{0},{1}".format(x ,y))

Quelqu'un peut m'aider?

J'ai vu ce lien :

Convertir la géométrie en WKT à l'aide d'ArcPy ?

http://resources.arcgis.com/en/help/main/10.1/index.html#//002z0000001t000000


J'essaye avec ça mais le résultat est :

import arcpy infc = arcpy.GetParameterAsText(0) pour la ligne dans arcpy.da.SearchCursor(infc, ["[email protected]", "[email protected]"]): print("feature {0}:".format(row[0 ])) partnum = 0 pour la partie dans la ligne[1] : print("Partie {0} :".format(partnum)) pour pnt en partie : si pnt : print("{0}, {1}".format (pnt.x, pnt.y)) else print("Anneau intérieur :") partnum+=1

Erreur d'analyse SyntaxError : syntaxe invalide (ligne 9)

J'ai essayé avec ce que tu m'as suggéré mais le résultat est :

importer arcpy

shp = r"C:Lavoroprova_arcpypunt.shp" avec arcpy.da.SearchCursor(shp,['[email protected]']) comme cur : pour la ligne dans cur : imprimer la ligne[0]

Erreur d'exécution Traceback (appel le plus récent en dernier) : Fichier "", ligne 2, dans RuntimeError : jeton de géométrie inconnue '[email protected]' Que recommandez-vous ?

Quelqu'un peut m'expliquer la vraie différence entre JSON et GeoJSON ?

J'ai pu utiliser ce script pour enregistrer dans une géométrie .txt mais mon problème est que je n'enregistre que la géométrie en 2D et je n'exporte pas les attributs liés à la géométrie.

importer arcpy

shp = r"C:Lavoroprova_arcpylin.shp" f= open ("C:/Lavoro/prova_arcpy/conv.wkt","w") avec arcpy.da.SearchCursor(shp,['[email protected] ']) as cur : pour la ligne dans cur : f.write("wkt_geom "+row[0]) f.close()

Comment puis-je faire?


Vous définissez votre curseur de manière incorrecte. Vous devez utiliser le mot-clé with pour travailler avec arcpy.da (module d'accès aux données) curseurs. (Notez qu'il y a un ancien curseur, que vous ne devriez pas utiliser, et un nouveau curseur, que vous devriez utiliser)

Voici comment mettre en place un da curseur et imprimez le WKT de chaque fonctionnalité :

import arcpy shp = r"D:PythonArcPyTestingairports_eur.shp" avec arcpy.da.SearchCursor(shp,['[email protected]']) as cur: for row in cur: print row[0]

J'utilise un fichier de formes de points, le résultat ressemble donc à ceci :

POINT (4.4338443496287647 51.949130189938153) POINT (8.7858617703132005 53.052287104156015) POINT (15.067460500705343 37.470072800341015) POINT (8.7129595350243676 45.627440514038092) POINT (11.07741797973909556 49.494505217034451)


Voir la vidéo: Creating Feature Classes in ArcGIS Pro (Octobre 2021).