数据拟合多元高斯曲线

Question

我有一个质子束撞击闪烁片的图像。 该图像被过滤为灰度的8位bmp文件。 我可以使用python / matlab基本上获取大量数据点列表，其中每个点都包含(x,y)位置和(z)幅度。

束斑

放大表面图

我想用高斯建模我的z数据，然后找到曲线的“半峰全宽”。 如何使多元高斯拟合数据？ 是否有这样的Matlab / Python包或任何其他方式，代码等？

Answer 1

这似乎是一项非常适合ROOT的工作，并且ROOT提供python绑定。 https://root.cern.ch/pyroot http://www.rootpy.org/

如果您熟悉ROOT，则只需创建TH2直方图并填充您的点，其中8bit彩色图像的强度即为直方图计数。 然后，将TF2公式定义为二维高斯并将其与FIT方法（minuit算法）拟合。

一些链接： https : //root.cern.ch/root/html534/guides/users-guide/FittingHistograms.html https://root.cern.ch/root/html/tutorials/fit/fit2.C.html