AreEngine 求最小面积的外接矩形，非IEnvelope，表达不清楚了

1，总是会得到一些奇奇怪怪的要求，求一个面对象的外接最小面积的矩形，和ArcToolBox中的Mininum Bounding Geometry功能下的RECTANGLE_BY_AREA想似。具体看下图：

 

区别如上图所示：IEnvelope 得到的是下图所示，需要的是第一种

（只是记录一下，可以解决问题，但不是最优方法，代码冗余量大，待解决）20200312已更新

基本思路（旧）：获取所有边，以其中的某一条边入手，遍历所有点到这条边的距离，取最大距离的点（pointA），以此点为准画一条平行于第一条线的线，作为第二条边，再以此点（pointA）画一条垂直于第一条边的线作为辅助线，再次遍历所有点到此辅助线的距离，得到最大距离的点（pointB），以pointB画一条垂直第一条边的线作为第三条边，再次遍历所有点到第三条边的距离，取得最大距离的点pointC,以pointC画一条垂直第一条边的线，作为第四条边。然后相邻两条线求交点，得到四个点，用四个点构造矩形。计算面积，以此，将所有的矩形都得到然后，比较面积求最小的就好了，好鸡儿啰嗦我。

更改：旧方法遍历第一次边的时候，后续的操作即为计算该polygon以 各个边为基准的Envelope，后续操作可谓傻瓜至极。新操作看代码不多BB

private IGeometry MinAreaGeometry(IPolygon pPolygon)
        {
            ITopologicalOperator pTopo = pPolygon as ITopologicalOperator;
            IGeometry pGeometry = pTopo.ConvexHull();
            IPoint pPointCent = (pGeometry as IArea).Centroid;

            ITransform2D pTran = pGeometry as ITransform2D;
            ISegmentCollection pSegments = pGeometry as ISegmentCollection;
            List<double> pIndex = new List<double>();
            List<IEnvelope> pGeos = new List<IEnvelope>();
            List<double> listangle = new List<double>();
            for (int i = 0; i < pSegments.SegmentCount; i++)
            {
                ISegment pSegment = pSegments.get_Segment(i);
                ILine pLine = new LineClass() { FromPoint = pSegment.FromPoint, ToPoint = pSegment.ToPoint };
                pTran.Rotate(pPointCent, pLine.Angle * -1);
                pIndex.Add((pGeometry.Envelope as IArea).Area);
                pGeos.Add(pGeometry.Envelope);
                listangle.Add(pLine.Angle);
                pTran.Rotate(pPointCent, pLine.Angle);
            }

            IEnvelope pEnvRe = pGeos[pIndex.IndexOf(pIndex.Min())];

            IPointCollection pPoints = new MultipointClass();
            pPoints.AddPoint(pEnvRe.UpperLeft);
            pPoints.AddPoint(pEnvRe.UpperRight);
            pPoints.AddPoint(pEnvRe.LowerRight);
            pPoints.AddPoint(pEnvRe.LowerLeft);

            IGeometry pGeoEnvMin = ConstructPolygon(pPoints);

            ITransform2D pTranEnv = pGeoEnvMin as ITransform2D;
            pTranEnv.Rotate(pPointCent, listangle[pIndex.IndexOf(pIndex.Min())]);

            return pGeoEnvMin;
        }

 

可复制粘贴直接使用。

 

分类: ArcEngine开发