为什么我的曲线拟合达到最大FEV=10000?

2022-07-06 12:58:42 发布

您现在位置:Python中文网/ 问答频道 /正文

我试图做适合,但我每个maxfev即使增加到10000!。我想我做错了什么

我将通过4个相同大小(30X240X240)的图像堆栈循环,并在4个不同的图像上获得相同体素的值。因此,作为参考,我的循环将使用n=30,d=240


v1, v2, v3, v4 = 0.000, 0.000, 0.000, 0.000

xdata = np.array([50,300,600,1000])

D_map = np.zeros((n,d,d))
K_map = np.zeros((n,d,d))

def kurtosis(x, D, K):
    return (-x * D) + (1/6 * (x **2)* (D **2) * K)

for i in range(n):
    for r in range(d):
        for c in range(d):
            v1 = b50_data[i,r,c]
            if v1 < 28:
                D_map[i,r,c] = 0
                K_map[i,r,c] = 0
            else:
                v2 = b300_data[i,r,c]
                v3 = b600_data[i,r,c]
                v4 = b1000_data[i,r,c]
                ydata = np.array([v1,v2, v3, v4])
                ydata = np.log(ydata/v1)

                popt, pcov = curve_fit(kurtosis, xdata, ydata, maxfev= 10000)
                D_map[i,r,c] = popt[0] * 1000
                K_map[i,r,c] = popt[0]


Tags: in图像mapfordatanprangev3v2v4v1xdataydatapoptmaxfev