AreEngine 求最小面积的外接矩形,非IEnvelope,表达不清楚了
来自:https://www.cnblogs.com/clgis/p/9831915.html
1,总是会得到一些奇奇怪怪的要求,求一个面对象的外接最小面积的矩形,和ArcToolBox中的Mininum Bounding Geometry功能下的RECTANGLE_BY_AREA想似。具体看下图:
区别如上图所示:IEnvelope 得到的是下图所示,需要的是第一种
(只是记录一下,可以解决问题,但不是最优方法,代码冗余量大,待解决)20200312已更新
基本思路(旧):获取所有边,以其中的某一条边入手,遍历所有点到这条边的距离,取最大距离的点(pointA),以此点为准画一条平行于第一条线的线,作为第二条边,再以此点(pointA)画一条垂直于第一条边的线作为辅助线,再次遍历所有点到此辅助线的距离,得到最大距离的点(pointB),以pointB画一条垂直第一条边的线作为第三条边,再次遍历所有点到第三条边的距离,取得最大距离的点pointC,以pointC画一条垂直第一条边的线,作为第四条边。然后相邻两条线求交点,得到四个点,用四个点构造矩形。计算面积,以此,将所有的矩形都得到然后,比较面积求最小的就好了,好鸡儿啰嗦我。
更改:旧方法遍历第一次边的时候,后续的操作即为计算该polygon以 各个边为基准的Envelope,后续操作可谓傻瓜至极。新操作看代码不多BB
private IGeometry MinArea(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<IGeometry> pGeos = new List<IGeometry>(); for (int i = 0; i < pSegments.SegmentCount; i++) { ISegment pSegment = pSegments.get_Segment(i); ILine pLine = new LineClass() { FromPoint = pSegment.FromPoint, ToPoint = pSegment.ToPoint }; axMapControl1.FlashShape(pLine.FromPoint); pTran.Rotate(pPointCent, pLine.Angle * -1); pIndex.Add((pGeometry.Envelope as IArea).Area); pGeos.Add(pGeometry.Envelope); pTran.Rotate(pPointCent, pLine.Angle); } IGeometry pGeoRe = pGeos[pIndex.IndexOf(pIndex.Min())]; return pGeoRe; }
可复制粘贴直接使用。