C#求點集的最小包圍矩形
更新時間:2020年09月18日 11:37:03 作者:俠客行江湖
這篇文章主要為大家詳細介紹了C#求點集的最小包圍矩形,文中示例代碼介紹的非常詳細,具有一定的參考價值,感興趣的小伙伴們可以參考一下
C# 求點集的最小包圍矩形,供大家參考,具體內(nèi)容如下
思路:
1、求點集的中心點
2、將點集繞矩形進行一系列角度的旋轉,并求記錄旋轉點集的包圍矩形的面積和旋轉角度;
3、將面積最小的矩形繞點集中心點旋轉回去。
// 1.尋找多邊形的中心 public XYZ GetCenter(List<XYZ> pts) { double sumx = 0; double sumy = 0; foreach (var p in pts) { sumx = sumx + p.X; sumy = sumy + p.Y; } var pt = new XYZ(sumx/pts.Count(),sumy/pts.Count(),0); return pt; } // 2.旋轉多邊形,針對每個點實現(xiàn)繞中心點旋轉 public XYZ RotatePt(XYZ inpt ,XYZ centerPt ,double theta) { double ix = inpt.X; double iy = inpt.Y; double cx = centerPt.X; double cy = centerPt.Y; double Q = theta / 180 * 3.1415926; //角度 double ox, oy; ox = (ix - cx) * Math.Cos(Q) - (iy - cy) * Math.Sin(Q) + cx; //旋轉公式 oy = (ix - cx) * Math.Sin(Q) + (iy - cy) * Math.Cos(Q) + cy; var outpt = new XYZ(ox,oy,0); return outpt; } // 3.多邊形旋轉后求簡單外接矩形 public List<XYZ> GetRect(List<XYZ> inpts) { var outpts =new List<XYZ>(); int size = inpts.Count(); if (size == 0) return null; else { var tempx = new List<double>(); var tempy = new List<double>(); for (int i = 0; i < size; i++) { tempx.Add(inpts[i].X); tempy.Add(inpts[i].Y); } XYZ endpoint0 = new XYZ(tempx.Min(), tempy.Max(), 0); XYZ endpoint1 = new XYZ(tempx.Max(), tempy.Max(), 0); XYZ endpoint2 = new XYZ(tempx.Max(), tempy.Min(), 0); XYZ endpoint3 = new XYZ(tempx.Min(), tempy.Min(), 0); outpts.Add(endpoint0); outpts.Add(endpoint1); outpts.Add(endpoint2); outpts.Add(endpoint3); return outpts; } } // 4.存儲每個旋轉角度下多邊形的外接矩形,記錄外接矩形的頂點坐標、面積和此時多邊形的旋轉角度 public class RectData { public List<XYZ> boundary { get;set;} public XYZ center { get; set; } public double theta { get; set; } public double area { get; set; } } public RectData GetRotateRectDatas(List<XYZ> inpts, double theta) { XYZ center = GetCenter(inpts); var tempvertices = new List<XYZ>(); for (int i=0; i<inpts.Count();i++) { XYZ temp = RotatePt(inpts[i], center, theta); tempvertices.Add(temp); } List<XYZ> vertices = GetRect(tempvertices); double deltaX, deltaY; //求每個外接矩形的面積 deltaX = vertices[0].X - vertices[2].X; deltaY = vertices[0].Y - vertices[2].Y; var polygen = new RectData { area=Math.Abs(deltaY * deltaX), center= center, theta = theta, boundary= vertices }; return polygen; } //獲取所有新的矩形 public List<RectData> GetAllNewRectDatas(List<XYZ> inpts) { var polygens =new List<RectData>(); for (int theta = 0; theta <= 90;) { polygens.Add(GetRotateRectDatas(inpts, theta)); theta = theta + 5; } return polygens; } //獲取新的矩形 public RectData GetMinAreaRect(List<RectData> polygons) { double minarea = 100000000; int N =0; for ( int i=0; i< polygons.Count(); i++) { if (minarea > polygons[i].area) { minarea = polygons[i].area; N = i; } } var polygon = new RectData(); polygon = polygons[N]; //旋轉到最小面積的方向 XYZ centerPt = GetCenter(polygon.boundary); var boundary = new List<XYZ>(); foreach(var bound in polygon.boundary) { XYZ pt = RotatePt(bound, polygon.center, -polygon.theta); boundary.Add(pt); } var outpolygon = new RectData { center= polygon.center, area = polygon.area, theta = polygon.theta, boundary = boundary }; return outpolygon; }
以上就是本文的全部內(nèi)容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。
您可能感興趣的文章:
相關文章
淺拷貝和深拷貝深入理解(shallow copy VS deep copy)
淺拷貝和深拷貝深入理解(shallow copy VS deep copy) 本文重點討論引用類型變量的拷貝機制和實現(xiàn)2014-01-01C#中如何將MongoDB->RunCommand結果映射到業(yè)務類的方法總結
這篇文章主要給大家總結介紹了關于C#中如何將MongoDB->RunCommand結果映射到業(yè)務類的方法,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面來一起看看吧。2018-04-04