Set Vertical Datum using Arcpy

VikramS — Mon, 21 Jan 2013 19:37:51 GMT

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

Re: Set Vertical Datum using Arcpy

ChrisFox3 — Thu, 31 Jan 2013 22:36:08 GMT

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.

topic Set Vertical Datum using Arcpy in Python Questions

Set Vertical Datum using Arcpy

Re: Set Vertical Datum using Arcpy