Comment obtenir tous les accords lat / long qui définissent un fichier de formes

J'ai un Shapefile avec plusieurs formes à l'intérieur.

En utilisant MapWindow, j'ai pu déterminer que la forme dont j'ai besoin est celle avec shapeId 19.

Comme vous pouvez probablement le constater, j'ai très peu d'expérience avec le SIG en général, mais je pense que je suis venu au bon endroit pour demander de l'aide.

Ce dont j'ai besoin, c'est d'extraire toutes les coordonnées lat / long qui déterminent la forme.

Est-il possible de le faire en utilisant MapWindow ou dois-je utiliser un autre logiciel?

Si plus d'informations sont nécessaires pour vous aider, veuillez commenter et je mettrai à jour dès que possible.

Toute aide est appréciée car cela me rend fou!

QGIS peut vous aider. Vérifiez cette réponse (la partie WKT) à une question similaire: Conversion des polygones de carte fiscale de Shapefile en table de numéro de carte et de coordonnées de coin

• shp2text

• Geo :: ShapeLib

• ShpDump

• Plus...

Vous trouverez ci-dessous un moyen d'accéder aux coordonnées de latitude et de longitude du fichier de formes ESRI, entre autres bits d'informations tels que la référence spatiale, les attributs de champ, les valeurs de champ, etc., à l'aide de Python. Le code ci-dessous ne fonctionne que pour les polygones et les points (car je n'ai pas encore écrit le code pour les polylignes). J'ai essentiellement bricolé du code que j'ai trouvé éparpillé dans l'aide d'ArcGIS Desktop 9.3, ajouté quelques-unes des miennes et rassemblé dans une seule fonction. Il a été écrit avec ArcGIS 9.3. Vous devriez pouvoir passer un fichier de formes polygonal ou un fichier de formes ponctuelles et la logique le dirigera en conséquence.

 def get_shapefile( shape_file ):
    # Import native arcgisscripting module
    import arcgisscripting

    # Create the geoprocessor object
    gp = arcgisscripting.create(9.3)

    # Identify the geometry field
    desc = gp.Describe( shape_file )
    shapefieldname = desc.ShapeFieldName

    # Get shapefile Name
    print
    print 'Shapefile Name: ', desc.Name

    # Get the spatial reference
    spatial_reference = desc.SpatialReference.Name
    print 'Spatial Reference: ', spatial_reference
    print

    # Create search cursor
    rows = gp.SearchCursor( shape_file )
    row = rows.Next()

    # Enter while loop for each feature/row
    while row:

        # Create the geometry object
        feat = row.GetValue(shapefieldname)

        print '*' * 30
        print

        print 'Geometry related Information'
        print
        # Get Geometry Type
        geometry_Type = feat.Type
        print 'Geometry Type: ', geometry_Type

        # Get the area of the feature
        geometry_Area = feat.Area
        print 'geometry_Area; ', geometry_Area

        # Get the centroid for the feature
        geometry_Centroid = feat.Centroid
        print 'geometry_Centroid:', geometry_Centroid

        # Get the extent for the feature
        geometry_Extent = feat.Extent
        print 'geometry_extent: ', geometry_Extent

        print
        print 'Get Attribute Table Information'

        # Get all the fields for the feature class
        fields = desc.Fields

        total_number_of_fields = len( fields )
        print 'Total number of fields: ', total_number_of_fields
        print

        print 'List attribute table related information:'
        print

        field_num_cntr = 0

        # Loop through all the fields in the feature class
        for field in fields:

            print '*'*5, field_num_cntr, '*'*5
            print
            print 'field Type: ', field.Type
            print 'Scale: ', str(field.Scale)
            print 'Precision: ', str(field.Precision)
            print field.Name, '=> ', row.GetValue( field.Name )
            print

            field_num_cntr += 1


        if geometry_Type == 'polygon':

            # Variable to keep track of how many multipart polygons are in
            # featureclass
            partnum = 0 

            # Count the number of points in the current multipart feature
            partcount = feat.PartCount

            print
            print 'Number of polygons in feature class: ', partcount
            print

            # Enter while loop for each part in the feature (if a singlepart feature
            # this will occur only once)
            while partnum < partcount:

                # Print the part number
                print "Part ", str(partnum), "of", partcount, ":"
                print
                part = feat.GetPart(partnum)
                pnt = part.Next()

                pntcount = 0

                # Enter while loop for each vertex
                while pnt:

                    # Print x,y coordinates of current point
                    print 'X coord:', pnt.x, 'Y coord:', pnt.y
                    pnt = part.Next()
                    pntcount += 1

                    # If pnt is null, either the part is finished or there is an interior ring
                    if not pnt:
                        pnt = part.Next()
                        if pnt:
                            print "Interior Ring:"
                partnum += 1

                print
                print 'Number of coordinates in feature class: ', pntcount - 1
                print

        elif geometry_Type == 'point':

            feat = row.GetValue(shapefieldname)

            # Get coords
            pnt = feat.GetPart()

            # Print x,y coordinates of current point object
            print 'X coord:', pnt.x, 'Y coord:', pnt.y


        row = rows.Next()


 your_shapefile = 'Path\To\Your\Shapefile.shp'
 get_shapefile( your_shapefile )