[求助] 求一个matlab算光斑直径的算法

如题，我现在拿ccd拍到了激光的光斑，求一个matlab算法可以算出1/e^-2的直径。

将图像转换为灰度图，并沿直径方向取灰度轮廓，应该可以得到一个准高斯曲线，按如下高斯公式拟合，可以得到半高全宽(FWHM)。1/e^2全宽的0.84大致等于FWHM。

不需要降噪、拟合什么的吗

本帖最后由 光杆司令 于 2017-10-31 20:20 编辑

一般情况下，不需要降噪，因为拟合反映大的趋势；
下面贴一份高斯拟合的python代码，txt文件记录了沿直径方向的灰度值，png是运行代码后生成的图像。

首先感谢您的帮助。
不过，这个python没用过，电脑里也没安装。
从算法上讲，您的这个只用了一个直径方向拟合而没有二维圆拟合，假设这个误差不太大，我在matlab里用同样的算法也可以了吧。
其次，假设由于透镜像差的原因导致高斯光束通过系统后CCD上不再是高斯分布，我可不可以找到灰度最大值后，圈出1/e^2点拟合出圆，求出半径。