Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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

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

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.