最小二乗法で円フィッティングを考える2

こちらの続きです。
最小二乗法で円フィッティングを考える1 - C#ATIA
・・・結局、探しても見つからない為、numpy使う事にしました。

こちらを参考にさせてもらいました。
円の最小二乗近似 - morikomorou’s blog

    # param: list(Point3D)
    def circleFitting(points: list):

        x = np.array([p.x for p in points]) 
        y = np.array([p.y for p in points]) 

        A = np.vstack((x, y, np.ones((len(x))))).T
        v = -(x ** 2 + y ** 2)
        u, residuals, rank, s = np.linalg.lstsq(A, v, rcond=None)

        cx_pred = u[0] / (-2)
        cy_pred = u[1] / (-2)
        r_pred = np.sqrt(cx_pred ** 2 + cy_pred ** 2 - u[2])

        return cx_pred, cy_pred, r_pred

    # param: list(Point3D)
    def getFittingCircle(points: list) -> adsk.core.Circle3D:

        x, y, r = circleFitting(points)
        return adsk.core.Circle3D.createByCenter(
            adsk.core.Point3D.create(x, y, 0),
            adsk.core.Vector3D.create(0, 0, 1),
            r
        )