python实现根据一个斜椭圆上的点计算出椭圆的长短轴的代码

以下是 Python 代码实现：

import numpy as np
from scipy.optimize import minimize

def fit_ellipse(x, y):
    # 定义误差函数
    def f_ellipse(c):
        a, b, cx, cy = c
        return ((x - cx) ** 2 / a ** 2 + (y - cy) ** 2 / b ** 2 - 1) ** 2
    
    # 初始化参数
    x_mean = np.mean(x)
    y_mean = np.mean(y)
    x_diff = x - x_mean
    y_diff = y - y_mean
    init = [1, 1, x_mean, y_mean]
    
    # 最小化误差函数
    res = minimize(f_ellipse, init, method='Powell')
    a, b, cx, cy = res.x
    return a, b

# 示例
x = np.array([0.0, 1.0, 0.0, -1.0])
y = np.array([1.0, 0.0, -1.0, 0.0])
a, b = fit_ellipse(x, y)
print("长轴：", max(a, b))
print("短轴：", min(a, b))

这里使用了 Scipy 库的 minimize 函数来最小化误差函数，求解椭圆的长短轴。需要传入椭圆上的点的坐标 x 和 y，返回长轴 a 和短轴 b。