一个完备的要素裁剪函数
要素裁剪的问题，即如何裁剪要素的SHAPE，还有许多网友在询问，其实这并没有什么特别之处。除了AO自带的Clip函数以外，我一般都是使用ITopologicalOperator接口提供的拓扑方法来自己搞定，原因很简单，每一条要素的裁剪情况都能自己掌握。
要素裁剪无非两种情况，一是不处理裁剪面边缘上的要素，统统删除，这个没有什么问题，空间关系选择Intersect即可；另一个是裁剪边缘上的要素，如Polygon如果与裁剪面相交，则保留不想交的部分，如果Polyline与裁剪面相交，也是保留面外面的一部分线段。
这几个拓扑方法在很多时候会爆出异常，其实解决的方法在网上早就有了，关键是两个Geoemtry的SpatialReference需要一致才行。下面给出一个完备的函数：
public static IGeometry getUpdateOuterGeometry2(IGeometry srcGeo, IGeometry refPolygon)
{
if (srcGeo.SpatialReference != null && refPolygon.SpatialReference == null)
refPolygon.Project(srcGeo.SpatialReference);
if (srcGeo.SpatialReference == null && refPolygon.SpatialReference != null)
srcGeo.Project(refPolygon.SpatialReference);
IGeometry pGeo = null;
ITopologicalOperator2 pTopo;
ITopologicalOperator2 pTopo2;
switch (srcGeo.GeometryType)
{
case esriGeometryType.esriGeometryPolyline:
pTopo = (ITopologicalOperator2)refPolygon ;
pTopo.IsKnownSimple_2 = false;
pTopo.Simplify();
pTopo2 = (ITopologicalOperator2)srcGeo;
pTopo2.IsKnownSimple_2 = false;
pTopo2.Simplify();
pGeo = pTopo2.Difference(pTopo.Intersect(srcGeo, esriGeometryDimension.esriGeometry1Dimension));
break;
case esriGeometryType.esriGeometryPolygon:
pTopo = (ITopologicalOperator2)srcGeo;
pTopo.IsKnownSimple_2 = false;
pTopo.Simplify();
pGeo = pTopo.Difference(refPolygon);
break;
default:
pGeo = srcGeo;
break;
}
return pGeo;
}
Related Items
Or, take a look at Archives and Categories
面裁减的算法，我们写过，文章在我的博客里
http://hi.baidu.com/geochenyj/blog/item/97d0d7fbcb1eff2a4f4aea59.html