Set Vertical Datum using Arcpy

2618
1
01-21-2013 11:37 AM
VikramS
Occasional Contributor
Hello ,

Is it possible to set Z coordinate System using python ? I am creating the Spatial Reference object like below

sp = arcpy.CreateSpatialReference_management("GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],VERTCS['NAVD_1988',VDATUM['North_American_Vertical_Datum_1988'],PARAMETER['Vertical_Shift',0.0],PARAMETER['Direction',1.0],UNIT['Meter',1.0]];-400 -400 1000000000;-100000 10000;-100000 10000;8.98315284119521E-09;0.001;0.001;IsHighPrecision","#","#","#","#","#","0")

As you can see I have set the vertical datum . I have useed this spatial Reference object for creating the feature class . I have also enabled M and Z Values for the feature class ..

After creating the feature class , When I use the describe object on the feature class like below

test = arcpy.Describe("test123")
test.SpatialReference.exporttostring()

u"GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]];-400 -400 1000000000;-100000 10000;-100000 10000;8.98315284119521E-09;0.001;0.001;IsHighPrecision"


As you can see it is not returning the Z Cooridnate System .. Has the Z Coorindate System been set to feature class or not ? If yes , How to confirm it

Can anyone please help me in clearing this doubt

Thanks
Tags (2)
0 Kudos
1 Reply
ChrisFox3
Occasional Contributor III
I tested on my end and it is applied, it doesn't seem to be exposed through the describe object. However if you open the feature class properties you will see a Z Coordinate System tab.