求兩圓交點演算法代碼(golang實現)

這是一個建立於 的文章，其中的資訊可能已經有所發展或是發生改變。

數學求解原理請看:求兩圓相交的交點的方法

實現代碼：

package mainimport (    "math")//代表一個點，包含橫縱座標type Point struct {    X, Y float64}//代表一個圓，包含橫縱座標及半徑type Circle struct {    Point    R float64}//建立圓對象func NewCircle(x, y, r float64) *Circle {    return &Circle{Point{x, y}, r}}//求兩圓相交的交點，交點個數可能有0,1,2func Intersect(a *Circle, b *Circle) (p []Point) {    dx, dy := b.X - a.X, b.Y - a.Y    lr := a.R + b.R //半徑和    dr := math.Abs(a.R - b.R) //半徑差    ab := math.Sqrt(dx * dx + dy * dy) //圓心距    if ab <= lr && ab > dr {        theta1 := math.Atan(dy / dx)        ef := lr - ab        ao := a.R - ef / 2        theta2 := math.Acos(ao / a.R)        theta := theta1 + theta2        xc := a.X + a.R * math.Cos(theta)        yc := a.Y + a.R * math.Sin(theta)        p = append(p, Point{xc, yc})        if ab < lr { //兩個交點            theta3 := math.Acos(ao / a.R)            theta = theta3 - theta1            xd := a.X + a.R * math.Cos(theta)            yd := a.Y - a.R * math.Sin(theta)            p = append(p, Point{xd, yd})        }    }    return}