การกำหนดการฉายภาพสำหรับหลายรูปร่างใน ArcMap?

ฉันมีไฟล์รูปร่างมากกว่า 100 ไฟล์ที่ไม่มีไฟล์. prj และเมื่อฉันนำมาไว้ใน ArcMap 10 พวกเขาแสดงระบบพิกัดว่าไม่ทราบ ฉันรู้ว่าระบบพิกัดไฟล์รูปร่างทั้งหมดคือ GCS WGS 1984 ฉันยังรู้ว่าฉันสามารถใช้เครื่องมือ Define Projection GP เพื่อกำหนดระบบพิกัดให้กับแต่ละไฟล์ แต่จะใช้เวลาตลอดไป

ฉันหวังว่าจะมีเครื่องมือ GP ในการกำหนดสิ่งเหล่านี้ แต่ฉันไม่เห็น ต่อไปฉันคิดว่าบางทีฉันสามารถใช้ python เพื่อทำสิ่งนี้ดังนั้นฉันจึงดูในเมนูช่วยเหลือและพบสคริปต์ แต่มันทำให้ฉันมีข้อผิดพลาด

นี่คือรหัสหลามที่ฉันลองใช้ (นี่เป็นไฟล์ shp ไฟล์เดียวดังนั้นฉันยังคงมีความเจ็บปวดในการพิมพ์ชื่อสำหรับแต่ละไฟล์:

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)

ฉันคิดว่าพวกคุณคิดมากเรื่องนี้ ...

1. คลิกขวาที่เครื่องมือ "กำหนดการฉาย" ในกล่องเครื่องมือ
2. เลือก "แบทช์"
3. ลากและวางเลเยอร์ของคุณลงในคอลัมน์ "ชุดข้อมูลที่ป้อน"
4. คลิกขวาในช่อง "Coordinate System" ช่องแรกเพื่อกรอกข้อมูลการฉายที่ถูกต้อง
5. จากนั้นคลิกขวาที่ภาพฉายที่คุณเพิ่งเลือกและเลือก "เติม" ซึ่งจะกรอกข้อมูลประมาณการที่เหลือทั้งหมดให้คุณ
6. กด "ตกลง" และคุณทำเสร็จแล้ว

หากอยู่ในไดเรกทอรีเดียวกันสิ่งนี้จะใช้งานได้ (เพียงแทนที่เส้นทางของคุณอาจเพิ่มในการจัดการข้อยกเว้นบางอย่าง):

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))

เพียงทำสำเนาของ. prj และเปลี่ยนชื่อ ดังนั้นเช่นหากคุณมี 3 รูปร่างไฟล์:

• one.shp,
• two.shp,
• three.shp

กำหนดการฉายภาพสำหรับ one.shp และคุณจะมี one.prj ในไดเรกทอรี คัดลอก one.prj ไปยังไดเรกทอรีของ two.shp และเปลี่ยนชื่อเป็น two.prj ทำซ้ำสำหรับรูปร่างทุกไฟล์ . prj เป็นเพียงไฟล์ข้อความ ตราบใดที่ยังมี. prj อยู่ในไดเรกทอรีเดียวกันกับ. shp และด้วยชื่อเดียวกันซอฟต์แวร์จะทำการเลือก อัตโนมัติด้วยเครื่องมืออะไรก็ตามที่คุณคุ้นเคยสำหรับการคัดลอกและเปลี่ยนชื่อไฟล์

แก้ไข:

อาจเป็นวิธีที่ง่ายที่สุด นำพวกเขาทั้งหมดเข้ามาและกำหนด crs ที่ถูกต้องให้กับเอกสาร พวกเขาจะ "ปราศจากโครงการ" แต่นั่งอยู่ในที่ที่ถูกต้อง จากนั้นเพียงส่งออก (หลายรายการ) ไปยังตำแหน่งใหม่

ฉันเห็นเครื่องมือ GP สำหรับการฉายภาพเป็นกลุ่ม

มันอาจไม่ถูกเปิดเผยในอาร์คี ??

ฉันคิดว่าฉันจะสร้าง pgdb, (ไฟล์ส่วนบุคคลหรือแม้กระทั่ง sde) แล้วสร้าง fds (ชุดข้อมูลคุณลักษณะ)

กำหนด crs ที่เหมาะสมให้กับ fds นั้น

นำเข้าไฟล์รูปร่าง (หลายรายการ) [fc ทั้งหมดใน fds สืบทอด crs ของ fds]

จากนั้นคุณสามารถส่งออกไปยังไฟล์รูปร่าง (หลายรายการ)

นี่จะเป็นการแทนการฉายแบบแบทช์และฉันแน่ใจว่าสามารถเขียนสคริปต์ได้

สคริปต์ shell และ ogr2ogr ง่ายแค่ไหน? ดูhttp://www.gdal.org/ogr2ogr.htmlสำหรับไวยากรณ์ที่ถูกต้อง

นี่คือสิ่งที่ฉันใช้ ... มันจะกำหนดเฉพาะการฉายภาพสำหรับไฟล์แรสเตอร์ที่ไม่มีการฉายภาพ หวังว่ามันจะช่วย นอกจากนี้ยังสร้างรายการไฟล์ที่ไม่มีการฉายภาพเพื่อจุดประสงค์ในการประกันคุณภาพ

สำหรับเวกเตอร์ของคุณเพียงแค่ต้องการ mod เล็กน้อย - 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)