ArcEngine对面进行抽稀(减少节点数)

public class clsSimplyfyPolygon
{

    //对指定图元,按指定buffer距离做抽稀
    public IPolygon doPolygonSimpled(IGeometry pGeometry, double dblDistance, int maxPointCount)
    {
        IPolygon pPolygon;
       
        ISegmentCollection pSegCol;
        IGeometryCollection pGeoCol;

        IPointCollection pntDoneCol = new MultipointClass();


        //  pGeometry = polygonFeature.Shape;
        pPolygon = (IPolygon)pGeometry;
        pSegCol = (ISegmentCollection)pPolygon;
        pGeoCol = (IGeometryCollection)pPolygon;

        for (int i = 0; i < pGeoCol.GeometryCount; i++)
        {
            IGeometry subGeometry;
            IPointCollection subPointCol;

            subGeometry = pGeoCol.get_Geometry(i);
            subPointCol = (IPointCollection)subGeometry;
            pntDoneCol.AddPointCollection(doSimpledPointCollection(ref subPointCol, dblDistance, maxPointCount));

        }
        return pPolygon;
    }
    //循环抽稀,ipointCollection
    public IPointCollection doSimpledPointCollection(ref IPointCollection pntCol, double dblDistance, int maxCount)
    {
        if (pntCol.PointCount > maxCount)
        {
            IPoint pntStart;
            IPoint pntMid;
            IPoint pntEnd;

            ArrayList arryDelIndexs = new ArrayList();
            arryDelIndexs.Clear();

            int modby3 = pntCol.PointCount % 3;
            int maxPointCount = pntCol.PointCount - modby3;

            for (int i = 0; i < maxPointCount; i = i + 3)
            {
                if ((i + 2) < pntCol.PointCount)
                {
                    pntStart = pntCol.get_Point(i);
                    pntMid = pntCol.get_Point(i + 1);
                    pntEnd = pntCol.get_Point(i + 2);

                    if (getMidDistanceFromThreePoints(pntStart, pntMid, pntEnd) < dblDistance)
                    {
                        arryDelIndexs.Add(i + 1);
                        //pntCol.RemovePoints(i + 1, 1);
                        //doSimpledPointCollection(ref pntCol, dblDistance, maxCount);
                    }
                }

            }

            if (arryDelIndexs.Count != 0)
            {
                for (int j = 0; j < arryDelIndexs.Count; j++)
                {

                    pntCol.RemovePoints(Convert.ToInt32(arryDelIndexs[j]) - j, 1);
                }
                doSimpledPointCollection(ref pntCol, dblDistance, maxCount);
            }

            if (pntCol.PointCount > maxCount)
            {
                doSimpledPointCollection(ref pntCol, dblDistance + dblDistance, maxCount);
            }
        }
        return pntCol;
    }

    //计算三点中,中间点到另外两点的直线距离
    public double getMidDistanceFromThreePoints(IPoint pntStart, IPoint pntMid, IPoint pntEnd)
    {
        double A = pntStart.Y - pntEnd.Y;
        double B = pntEnd.X - pntStart.X;
        double C = pntStart.X * pntEnd.Y - pntEnd.X * pntStart.Y;
        double distance = Math.Abs(A * pntMid.X + B * pntMid.Y + C) / Math.Sqrt(A * A + B * B);
        return distance;


    }

}

}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值