public static bool IsCircle(List<Point> points, out Point center, out double radius)
{
 int len = points.Count;
 center = new Point();
 radius = 0;

 // 判斷是否為封閉圖形
 if (!IsClosedFigure(points))
 return false;

 int judgePointNum = len * 50 / 100;
 if (len < judgePointNum)
 return false;

 // 取鏈表上三個(gè)點(diǎn)作為判斷圓的根據(jù)
 Point p1 = points[len / 6];
 Point p2 = points[len / 2];
 Point p3 = points[len * 5 / 6];
 if ((Math.Abs(p1.X - p2.X) < 100 && Math.Abs(p1.Y - p2.Y) < 100)
 || (Math.Abs(p1.X - p3.X) < 100 && Math.Abs(p1.Y - p3.Y) < 100)
 || (Math.Abs(p2.X - p3.X) < 100 && Math.Abs(p2.Y - p3.Y) < 100))
 return false;

 // 三個(gè)點(diǎn)確定圓的方程，獲取圓心坐標(biāo)及半徑
 GetCircle(p1, p2, p3, out center, out radius);

 // 獲取圓上平均分部的多個(gè)點(diǎn)，判斷其到圓心的距離與半徑之差是否在精度內(nèi)
 for (int i = 0; i < judgePointNum; ++i)
 {
 // 獲取圓上點(diǎn)
 Point p = points[len * i / judgePointNum];
 double deviation = Math.Abs(GetDistance(center, p) - radius);

 // 點(diǎn)在圓上的偏移量與半徑的比值若大于固定值，則不為圓
 if (deviation/radius > MaxRatio)
  return false;
 }

 return true;
}

/// <summary>
/// 最小二乘法求回歸直線(xiàn)方程
/// </summary>
/// <param name="points">輸入數(shù)據(jù)</param>
/// <param name="k">直線(xiàn)斜率</param>
/// <param name="b">直線(xiàn)截距</param>
/// <param name="type">直線(xiàn)類(lèi)型 1：水平線(xiàn) 2：垂直線(xiàn) 3：一般直線(xiàn)</param>
/// <returns></returns>
public static bool IsLine(List<Point> points, out double k, out double b, out int type)
{
 k = 0;
 b = 0;
 type = 0;

 if (points.Count < 2) return false;

 double averageX = 0, averageY = 0, n = 0;
 n = points.Count;
 foreach (Point p in points)
 {
 averageX += p.X;
 averageY += p.Y;
 }
 averageX /= n;
 averageY /= n;

 double numerator = 0, denominator = 0;
 foreach (Point p in points)
 {
 numerator += (p.X - averageX) * (p.Y - averageY);
 denominator += (p.X - averageX) * (p.X - averageX);
 }

 if (numerator == 0) //平行于X軸為水平線(xiàn)，返回縱坐標(biāo)平均值
 {
 b = averageY;
 type = 1;
 }
 else if (denominator == 0)//平行于Y軸為垂直線(xiàn)，返回橫坐標(biāo)平均值
 {
 b = averageX;
 type = 2;
 }
 else
 {
 type = 3;
 }

 k = numerator / denominator;
 b = averageY - k * averageX;

 foreach (Point p in points)
 {
 dis = GetPoint2LineDistance(p, k, b, type);
 if (dis > MAX_POINT_LINE_DIS) return false; //點(diǎn)到擬合直線(xiàn)距離過(guò)大
 }

 return true;
}

/// <summary>
/// 計(jì)算點(diǎn)到直線(xiàn)的距離
/// </summary>
/// <param name="p">待計(jì)算點(diǎn)</param>
/// <param name="k">直線(xiàn)斜率</param>
/// <param name="b">直線(xiàn)截距</param>
/// <param name="type">直線(xiàn)類(lèi)型 1：水平線(xiàn) 2：垂直線(xiàn) 3：一般直線(xiàn)</param>
/// <returns>距離</returns>
private static double GetPoint2LineDistance(Point p, double k, double b, int type)
{
 if (type == 1)
 {
 return Math.Abs(p.Y - b);
 }
 else if (type == 2)
 {
 return Math.Abs(p.X - b);
 }
 else
 {
 double numerator = 0, denominator = 0;
 numerator = Math.Abs(k * p.X - p.Y + b);
 denominator = Math.Sqrt(k * k + 1);
 return numerator / denominator;
 }
}