Definieren von Projektionen für mehrere Shapefiles in ArcMap?

Ich habe über 100 Shape-Dateien, die keine .prj-Datei haben. Wenn ich sie in ArcMap 10 bringe, wird das Koordinatensystem daher als unbekannt angezeigt. Ich weiß, dass das Koordinatensystem aller Formdateien GCS WGS 1984 ist. Ich weiß auch, dass ich mit dem GP-Tool "Projektion definieren" das Koordinatensystem jeder Datei einzeln zuweisen kann, aber das wird ewig dauern.

Ich hatte gehofft, dass es ein GP-Tool gibt, mit dem sich diese stapelweise definieren lassen, aber ich sehe keines. Als nächstes dachte ich, ich könnte vielleicht Python verwenden, also habe ich im Hilfemenü ein Skript gefunden, aber es gibt mir einen Fehler.

Hier ist der Python-Code, den ich ausprobiert habe (dies ist für eine einzelne shp-Datei, so dass ich immer noch die Mühe habe, den Namen für jede Datei einzugeben:

import arcpy
infc = r"C:\Documents and Settings\User\My Documents\ArcGIS\shpfiles\Site_2.shp"
prjfile = r"<install directory>\Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj"
arcpy.DefineProjection_management(infc, prjfile)

Ich denke ihr überdenkt diesen einen ...

1. Klicken Sie mit der rechten Maustaste auf das Tool "Define Projection" in der Toolbox.
2. Wählen Sie "Charge",
3. Ziehen Sie Ihre Ebenen per Drag & Drop in die Spalte "Input Dataset".
4. Klicken Sie mit der rechten Maustaste in das erste Feld "Koordinatensystem", um die richtige Projektion auszufüllen.
5. Klicken Sie dann mit der rechten Maustaste auf die gerade ausgewählte Projektion und wählen Sie "Ausfüllen", um alle übrigen Projektionen für Sie auszufüllen.
6. Klicken Sie auf "OK" und Sie sind fertig.

Wenn sie sich im selben Verzeichnis befinden, würde so etwas funktionieren (ersetzen Sie einfach Ihre Pfade, fügen Sie möglicherweise eine Ausnahmebehandlung hinzu):

import os, shutil

wgs84prjpath='c:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Geographic Coordinate Systems/World/WGS 1984.prj'

def definewgs84(shpfilepath):
    tgtpath = os.path.splitext(shpfilepath)[0] + '.prj'
    shutil.copyfile(wgs84prjpath,tgtpath)


yourshapefiledirectory='C:/temp/'
extension='shp'

shpfileslist = [file for file in os.listdir(yourshapefiledirectory) if file.lower().endswith(extension)]

for filename in shpfileslist:
    definewgs84(os.path.join(yourshapefiledirectory,filename))

Einfach Kopien der .prj erstellen und umbenennen. Also zB wenn Sie 3 Shapefiles haben:

• one.shp,
• two.shp,
• drei.shp.

Definieren Sie die Projektion für one.shp und Sie haben one.prj im Verzeichnis. Kopieren Sie one.prj in das Verzeichnis two.shp und benennen Sie es in two.prj um. Wiederholen Sie diesen Vorgang für alle Shapefiles. Die .prj ist nur eine Textdatei. Solange sich eine .prj-Datei im selben Verzeichnis wie eine .shp-Datei mit demselben Namen befindet, wird sie von der Software abgerufen. Automatisieren Sie mit allen Tools, die Sie zum Kopieren und Umbenennen von Dateien kennen.

Bearbeiten:

Die wahrscheinlich einfachste Methode wäre. Bringen Sie sie alle ein und weisen Sie dem Dokument die richtigen crs zu. Sie würden "unprojektiert" sein, aber am richtigen Ort sitzen. Dann exportieren Sie einfach (mehrere) an einen neuen Speicherort.

Ich sehe ein GP-Tool für die Stapelprojektion.

Es darf nicht bogenförmig belichtet werden. ??

Ich denke, ich würde eine pgdb erstellen (Datei, persönlich oder sogar sde) und dann einen fds (Feature-Datensatz) erstellen.

Weisen Sie diesem fds das entsprechende crs zu.

Importieren Sie die Formdateien (mehrere). [alle fc in einem fds erben die crs des fds].

Dann können Sie in eine Formdatei exportieren (mehrere).

Dies würde die Stapelprojektion ersetzen, und ich bin sicher, dass Skripte erstellt werden könnten.

Wie wäre es mit einem einfachen Shell-Skript und ogr2ogr? Die genaue Syntax finden Sie unter http://www.gdal.org/ogr2ogr.html .

Hier ist, was ich benutze ... es wird nur die Projektion für Rasterdateien definieren, die keine Projektion haben. Ich hoffe es hilft. Es wird auch eine Liste der Dateien ohne Projektion für Qualitätssicherungszwecke erstellt.

Für deine Vektoren wird nur ein kleiner Mod benötigt - FileList = arcpy.ListFeatureClasses ()

# Defines projection for all rasters with undefined projection
# CAUTION - make sure you know that the projection you are defining is the correct one
# for all files in the current directory.

#Licence: Creative Commons
#Created by: George Corea; georgec@atgis.com.au, coreagc@gmail.com

import arcpy, glob, os, sys, arcgisscripting
from arcpy import env, mapping
path = os.getcwd()
env.workspace = path
env.overwriteOutput = True

print 'Reading files from ' + path
os.chdir(path)
prjFile=r'C:\\Python26\\GDA_1994_MGA_Zone_55.prj'

x=0
z=x+1
NoProjCount=0
FileList= arcpy.ListRasters()
for File in FileList:
    desc = arcpy.Describe(File)
    SR = desc.spatialReference
    if SR.name == "Unknown":
        print "Projection of " + str(File) + " is " + SR.name + " so defining projection."
        f = open('NoProjection.txt', 'a')
        f.write(str(File)+"\n")
        f.close()
        arcpy.DefineProjection_management(File, prjFile) 
        NoProjCount=NoProjCount+1
    else:    
        print File + " is projected " + str(SR.name)