在ArcEngine中创建带高程Z值的点和线图层

管线和其附属物的坐标数据都是带有Z值的

而且有些情况下，一个管段的两个端点的x,y值一模一样(垂直的管段)

这样的线，在直接生成shape图层的时候，就会产生问题，特别是

使用ArcSDE的C API直接创建到表中的时候你会发现，这样的数据是生成不了的

解决的方法就是为图层添加Z值，一个图层是否带有高程值，可以在ArcMap中通过查看

图层的属性表得知，带有Z值的图元，在Shape字段中，其类型后面会有个ZM字样

比如point ZM，Polyline ZM

创建带高程Z值的图层时，只设置Point的Z属性是没有用的，默认情况下

ArcEngine会忽略Z值，

两步走，解决问题：

第一步：创建Shape图层，使用IFeatureWorkspace.CreateFeatureClass

此函数的第二个参数是IFields，字段定义，在字段定义中对Shape字段的类型

使用IGeometryDefEdit.HasZ_2 = true指定其包含Z值，

IGeometryDef pGeometryDef = new GeometryDef(); // 为esriFieldTypeGeometry类型的字段创建几何定义，包括类型和空间参照 
                    IGeometryDefEdit pGeometryDefEdit = pGeometryDef as IGeometryDefEdit;
                    pGeometryDefEdit.GeometryType_2 = geometryType;
                    pGeometryDefEdit.HasZ_2 = true;//图层是有高程值的
                    pGeometryDefEdit.SpatialReference_2 = axmapcontrol.SpatialReference;                    
                    pFieldEdit.GeometryDef_2 = pGeometryDef;

第二步：添加图元，需要使用IZAware接口指定其ZAware为true

IPoint pFromPoint = new PointClass();
pFromPoint.PutCoords(fromX, fromY);                                                           
                                pFromPoint.Z = formZ;
                                IZAware fromZAware = pFromPoint as IZAware;
                                fromZAware.ZAware = true;
                                //IZ iFromZ = (IZ)pFromPoint;
                                IPoint pToPoint = new PointClass();
                                pToPoint.PutCoords(toX, toY);                                
                                pToPoint.Z = toZ;
                                IZAware toZAware = pToPoint as IZAware;
                                toZAware.ZAware = true;
                                //IZ iToZ = (IZ)pToPoint;
                                
                                
                                IPolyline pPolyline = new PolylineClass();
                                IZAware iPolylineAware = (IZAware)pPolyline;
                                iPolylineAware.ZAware = true;
                                pPolyline.FromPoint = pFromPoint;
                                pPolyline.ToPoint = pToPoint;

顺便提一下ArcSDE C API中的此种问题，

一开始我以为比较简单，使用SE_layerinfo_set_3D就能解决问题，

但老有几条记录导入不了，使用field calculator查看一下图元的Z值才发现，

高程都变成整数了，我可以确定传入的数据是double型的，但是调用

SE_shape_generate_point和SE_shape_generate_line之后，生成的数据就是高程为整数

查了查，网上有代码说使用SE_coordref_set_precision设置坐标为高精度的，

我在创建layerinfo和插入shape的时候都进行了设置，却没有任何效果，只好暂时作罢..

20100723更新：

上面高程精度丢失问题解决了，

本来以为很简单，直接生成shape再用ArcCatalog导入到ArcSDE中，

用的时候发现，有几个图层死活导入不了，

提示错误是ORA-01426:numericoverflow,

于是还是得写代码直接导入，今天折腾了将近一天，想法是用C++直接调用SDE的C API来做，试试看值传入的有没有问题，

搞到下午，折腾的差不多了的时候，突然发现这个函数SE_coordref_set_z_by_range

使用高程的最大值最小值加0.001后传入这个函数，效果是立竿见影啊，立即搞定了，

囧的是，为了这个问题，午饭都没吃