ArcMap에서 여러 모양 파일에 대한 투영을 정의 하시겠습니까?

.prj 파일이없는 100 개 이상의 모양 파일이 있으므로 ArcMap 10으로 가져올 때 좌표계를 알 수없는 것으로 표시합니다. 모든 형상 파일 좌표계가 GCS WGS 1984라는 것을 알고 있습니다. 또한 투영 GP 정의 도구를 사용하여 좌표계를 각 파일에 개별적으로 할당 할 수 있지만 시간이 오래 걸릴 것입니다.

배치를 정의하는 GP 도구가 있기를 바랐지만 보이지 않습니다. 다음으로 파이썬을 사용 하여이 작업을 수행 할 수 있다고 생각했기 때문에 도움말 메뉴를보고 스크립트를 찾았지만 오류가 발생했습니다.

여기에 내가 시도한 파이썬 코드가 있습니다 (이것은 단일 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. 첫 번째 '좌표계'상자를 마우스 오른쪽 버튼으로 클릭하여 올바른 투영을 채 웁니다.
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 개의 shapefile이있는 경우 :

• one.shp,
• two.shp,
• 삼 .hp.

one.shp의 투영을 정의하면 디렉토리에 one.prj가 있습니다. one.prj를 two.shp의 디렉토리에 복사하고 이름을 two.prj로 바꾸고 모든 shapefile에 대해 반복하십시오. .prj는 텍스트 파일 일뿐입니다. .shp와 동일한 디렉토리에 동일한 이름을 가진 .prj가 있으면 소프트웨어가이를 선택합니다. 파일 복사 및 이름 바꾸기에 익숙한 도구로 자동화하십시오.

편집하다:

아마도 가장 쉬운 방법 일 것입니다. 그것들을 모두 가져 와서 올바른 crs를 문서에 할당하십시오. 그들은 "투영"될 것이지만 올바른 장소에 앉아있을 것입니다. 그런 다음 새 위치로 여러 항목을 내보내십시오.

배치 프로젝션을위한 GP 도구가 표시됩니다.

아크 피에 노출되지 않을 수 있습니다. ??

pgdb (File, personal 또는 sde)를 만든 다음 fds (feature data set)를 만들 것이라고 생각합니다.

해당 crd를 해당 fd에 지정하십시오.

모양 파일을 가져옵니다 (여러 개). [fds의 모든 fc는 fds의 crs를 상속합니다].

그런 다음 모양 파일 (여러)로 내보낼 수 있습니다.

이것은 배치 프로젝션 대신에 스크립트로 작성 될 수 있습니다.

간단한 쉘 스크립트와 ogr2ogr은 어떻습니까? 정확한 구문 은 http://www.gdal.org/ogr2ogr.html 을 참조 하십시오 .

여기에 내가 사용하는 것이 있습니다 ... 투영이없는 래스터 파일의 투영 만 정의합니다. 도움이 되길 바랍니다. 또한 품질 보증 목적으로 프로젝션없이 파일 목록을 작성합니다.

벡터의 경우 약간의 모드가 필요합니다-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)