Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

As the title says I've several feature classes. I want to check whether they have same spatial reference using arcpy

share|improve this question
add comment

1 Answer 1

you can do it with this way:

import arcpy

env.workspace = "D:/test/data.gdb"
fcs = arcpy.ListFeatureClasses()

ary = []

for ft in fcs:    
    desc = arcpy.Describe(ft)
    spatialRef = desc.SpatialReference

    ary.append([spatialRef.Name, ft])
    print spatialRef.Name

# Finding the same Spatial Reference Object
for a in range(len(ary)):
    if ary[a][0] == 1:
       print ary[a]

# Finding different form other Spatial Reference Object
for a in range(len(ary)):
    if ary[a][0] > 1:
       print ary[a]

i hope it helps you...

share|improve this answer
    
I'm just a programmer working on a GIS software. Do spatialRef.Name always identifies the projection? Is checking spatialRef.Name equals is enoguh? –  user Jun 4 '12 at 7:22
    
spatialRef.Name will give you Map Projections as Mercator, Gauss-Kruger or Lambert Conformal Conic. if you want to learn spatial reference code you can use desc.SpatialReference.factoryCode –  Aragon Jun 4 '12 at 8:59
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.