As the title says I've several feature classes. I want to check whether they have same spatial reference using arcpy
asked Jun 4, 2012 at 7:09
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...
answered Jun 4, 2012 at 7:19
-
I'm just a programmer working on a GIS software. Do
spatialRef.Name
always identifies the projection? Is checkingspatialRef.Name
equals is enoguh?user– user2012年06月04日 07:22:42 +00:00Commented Jun 4, 2012 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
urcm– urcm2012年06月04日 08:59:32 +00:00Commented Jun 4, 2012 at 8:59
lang-py