Archived Forums

# gp_prj_assign.py
#
# Batch assignment of coordinate reference systems to feature classes
'''Batch assing coordinate systems'''

import os, sys, string, win32com.client

# start reporting
import time
me=os.path.splitext(os.path.basename(sys.argv[0]))[0]
print "\n",__doc__,"\n",me,"started at: " + str(time.asctime())
timestamp=time.clock()

# Create the Geoprocessor object
gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")

try:
    path=sys.argv[1]					   #1
except:
    path="G:\\beta_test\\Tests\\wf\\test_b_def\\"  #1

oplist=[]
gp.workspace = path
fclist = gp.ListFeatureClasses("*")                #2
fclist.reset()
fc=fclist.next()
while fc:
    desc=gp.describe(fc)
    SR = desc.SpatialReference
    if SR.Name == "Unknown":                       #3
        print fc,"has unknown CRS."
        oplist.append(fc)
    fc=fclist.next()
if len(oplist) == 0:
    print "All CRS are known."
    sys.exit(0)
    
crs="27700" # British National Grid

print "assigning",crs,"to",len(oplist),"feature classes..."
for fc in oplist:
    print fc
    try:
       gp.DefineProjection_management (fc, crs)    #4 
       print "CRS",crs,"assigend"
    except:
       print "could not assing CRS",crs,"to", fc
       gp.GetMessages()
       
delta=time.clock()-timestamp
timestamp=time.clock()
print me,"successfully finished:",str(time.asctime()),"\nTotal time in sec:",str(delta),"\n"    

# cpprj.py
# high speed copy projection files to new shapefiles
# because this not available in ArcGIS 9
# emulating AML conventions
# avoid Dispatch object
# but still add as a Toolbox tool
# 
# Kim Ollivier kimo@ollivier.co.nz
# 9 August 2004
"""
   Batch Projection Definition
   Usage cpprj <target_dir> <src_proj_file> 
"""

import os
import sys
import shutil
import glob

# optional, only for log messages
# from win32com.client import Dispatch
# gp = Dispatch("esriGeoprocessing.GpDispatch.1")
			
msgExit = "  Usage: " + os.path.basename(sys.argv[0]) + " <target_dir> {NZMG|srcprj_file} "

if len(sys.argv) == 1 :
	print msgExit
	sys.exit()

dirprj = sys.argv[1]


if not os.path.exists(dirprj) :
	print "  Dir path not found", dirprj
	print msgExit
	sys.exit()

# default projection file for me
default = os.path.normcase("c:\\arcgis\\nzmg.prj")

if len(sys.argv) == 2 :
	srcprj = os.path.normcase(default)
elif len(sys.argv) == 3 :
	srcprj = os.path.normcase(sys.argv[2])
else :
	print "  Too many parameters"
	print msgExit
	sys.exit()
	
if not os.path.isfile(srcprj) :
	print "  Projection file not found for: ", os.path.basename(srcprj), srcprj
	print msgExit
	print os.path.isfile(srcprj)
	sys.exit()

# change to directory 
os.chdir(dirprj)
# get list of shapefiles here
shps = glob.glob('*.shp')
# could use a walk of subdirectories

for shp in shps:
	parts = os.path.splitext(shp)
	newprj = parts[0] + '.prj'
	if not os.path.exists(newprj):
		shutil.copyfile(srcprj,newprj)
		print shp, " projection added"
		# gp.AddMessage(shp + " projection added")
	else:
		print shp, " already defined"
		# gp.AddMessage(shp + " projection already defined")

# command line example using a <code> not a <file.prj>
workspace d:/teararoa/motatapu/aerial
DefineProjection_management f40a.tif 43040

# dialog window:
Executing (DefineProjection_7): DefineProjection D:\teararoa\motatapu\aerial\f40a.tif GEOGCS['',DATUM['',SPHEROID['',1.0,0.0]],PRIMEM['',0.0],UNIT['',1.0]] D:\teararoa\motatapu\aerial\f40a.tif
Start Time: Tue Aug 10 09:22:18 2004
Executed (DefineProjection_7) successfully.
End Time: Tue Aug 10 09:24:02 2004 (Elapsed Time: 104.00 secs)

Public Sub MultiSpatRef()
  
  Dim pApp As esriCatalogUI.IGxApplication
  Dim pGxSelection As esriCatalog.IGxSelection
  Dim pEnumGX As esriCatalog.IEnumGxObject
  Dim pGxObject As esriCatalog.IGxObject
  Dim pName As esriSystem.IName
  Dim pDS As esriGeoDatabase.IDataset
  Dim pGDSE As esriGeoDatabase.IGeoDatasetSchemaEdit
  Dim pSpatRef As esriGeometry.ISpatialReference
  
  Set pApp = esriArcCatalog.Application
  Set pGxSelection = pApp.Selection
  Set pEnumGX = pGxSelection.SelectedObjects

  Set pGxObject = pEnumGX.Next
  
  'Exit if no dataset is selected
  If pGxObject Is Nothing Then
    MsgBox "Select a dataset", vbCritical, "No dataset selected"
    Exit Sub
  End If
    
  Set pSpatRef = SelectSpatRef
  
  'Exit out of sub if you "Cancelled" the Coordinate System dialog box
  If pSpatRef Is Nothing Then
    Exit Sub
  End If

  Do Until pGxObject Is Nothing
    If TypeOf pGxObject Is esriCatalog.IGxDataset Then
    
      Set pName = pGxObject.InternalObjectName
      Set pDS = pName.Open
      Set pGDSE = pDS
      
      With pGDSE
        If .CanAlterSpatialReference Then
          .AlterSpatialReference pSpatRef
        End If
      End With
      
    End If
    Set pGxObject = pEnumGX.Next
  Loop

End Sub


Private Function SelectSpatRef() As esriGeometry.ISpatialReference

  Dim pSpatRefDialog As esriCatalogUI.ISpatialReferenceDialog
  Set pSpatRefDialog = New esriCatalogUI.SpatialReferenceDialog
  
  Set SelectSpatRef = pSpatRefDialog.DoModalCreate(False, False, False, 0)

End Function

# ---------------------------------------------------------------------------
# Define projection for all feature classes in a main directory and subdirs
# Created: Dec 6 2004
# Author: Knut Jetlund (knut.jetlund@forsvarsbygg.no)
# ---------------------------------------------------------------------------

# Import system modules
import sys, string, os, win32com.client

# Create the Geoprocessor object
gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")

# Path to the main directory
ws = "D:\Test"
gp.workspace = ws
gp.AddMessage("Main directory: " + ws)

# Path and name for the projection file 
cs = "C:/Programfiler/ArcGIS/Coordinate Systems/Projected Coordinate Systems/Utm/Wgs 1984/WGS 1984 UTM Zone 32N.prj"

# Create a list of subdirectories (workspaces) under the main directory
workspaces = gp.listworkspaces()

# Loop through all workspaces, starting with the main directory
while gp.Exists(ws):
  gp.workspace = ws
  gp.AddMessage("Workspace: " + ws)

  #Find feature classes in the current workspace  
  datasets = gp.listfeatureclasses()
  dataset = datasets.Next()
  
  # Loop through all feature classes and set projection
  while dataset:
    gp.AddMessage("Defining projection for " + dataset)
    gp.DefineProjection_management(dataset, cs)
    dataset = datasets.Next()
    
  # Find next workspace
  ws = workspaces.Next()

# ---------------------------------------------------------------------------
# Define projection for all shapefiles in a directory
# Created on: to sep 16 2004 02:10:09 
# ---------------------------------------------------------------------------

# Import system modules
import sys, string, os, win32com.client

# Create the Geoprocessor object
gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")

# Load required toolboxes...
gp.AddToolbox("C:/Programfiler/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")

# Local variables...
# Path to the directory
gp.workspace = "D:/Kart/grunnkart/Vektor/K50lib/UTM_33_SHAPE_NY"
gp.toolbox = "management"

# extension
type = "shp"
# Path and name for the projection file 
cs = "C:/Programfiler/ArcGIS/Coordinate Systems/Projected Coordinate Systems/Utm/Wgs 1984/WGS 1984 UTM Zone 33N.prj"

# Find files in the directory
files = os.listdir(gp.workspace) 

# Loop through all files looking for shapefiles
for f in files :
  ext = string.lower(os.path.splitext(f)[1][1:])
  if ext == type :
    print "Definerer projeksjon for " + f + "... "
    # Process: Define Projection...
    gp.DefineProjection_management(f, cs)
    

Sub BatchDefineProject()

    Dim pGP As Object
    Set pGP = CreateObject("esriGeoprocessing.GpDispatch.1")
    
    'Set the workspace with the featureclasses
    pGP.Workspace = "C:\uc2004\shapefiles"
    
    'Set projection file
    prjfile = "C:\ArcGIS\Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj"
    'NOTE: This could be a path to an existing featureclass with the projecton I want.
    
    'Get list of shapefiles or featureclasses
    Dim fcs As Object
    Set fcs = pGP.ListFeatureClasses("*")
    
    'Get the first featureclass from list
    fc = fcs.Next()
    
    'Loop through and define a projection
    Do While Not fc = ""
       pGP.DefineProjection_management fc, prjfile
       fc = fcs.Next()
    Loop

End Sub