数据适合多元高斯曲线

Question

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

光束点

放大表面情节

我想用高斯模型对z数据进行建模，然后找到＆＃34;半高全宽＃34;的曲线。如何将多元高斯拟合到我的数据中？是否有这样的Matlab / Python包或任何其他方式，代码等。？

Answer 1

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

如果您熟悉ROOT，则只需要创建TH2直方图并填充您的点，其中8位彩色图像的强度是直方图计数。然后将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