使用Scipy对伽马光谱进行高斯拟合时拟合质量评估结果异常的问题咨询
Scipy对伽马光谱进行高斯拟合时拟合质量评估结果异常的问题咨询
我有伽马光谱数据,正在用Python结合Scipy对谱峰做高斯拟合,拟合过程本身没问题,但尝试用数值评估拟合质量(为了后续自动化处理)时,得到的结果却非常反常。
我使用的Scipy命令如下:
response = scipy.optimize.curve_fit(gauss, e, c, param0, full_output=True, absolute_sigma=True, method=FitMethod, bounds=bounds)
我通过下面的代码计算拟合质量:
fitQuality = np.linalg.cond(response[1]) # PCOV是response[1]
这里的拟合质量值范围是从0(拟合极佳)到无穷大(拟合极差)。举个例子,我对CsI探测器中的K-40伽马谱线做了拟合(绿色阴影部分),从数据的波动情况来看,我对这个拟合结果很满意,Scipy给出的拟合质量值是18,这也符合预期。
接下来是对高分辨率锗探测器中的Na-22谱线的拟合(红色阴影部分),这个拟合结果我同样非常满意,看起来几乎完美,但Scipy给出的拟合质量值却高达56,982,136,这意味着拟合质量极差。
这完全说不通啊,明明拟合效果近乎完美!
是不是我用来计算拟合质量的公式不合适?需要调整吗?还是我从根本上就用错了评估指标?请大家帮忙解答一下。
备注:内容来源于stack exchange,提问作者ullix
相关产品推荐
相关产品推荐

