求一个matlab算光斑直径的算法

六月等风起

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

光杆司令

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

六月等风起

光杆司令 发表于 2017-10-31 13:10
将图像转换为灰度图，并沿直径方向取灰度轮廓，应该可以得到一个准高斯曲线，按如下高斯公式拟合，可以得到 ...

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

光杆司令

六月等风起 发表于 2017-10-31 15:41
不需要降噪、拟合什么的吗

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

六月等风起

光杆司令 发表于 2017-10-31 18:51
一般情况下，不需要降噪，因为拟合反映大的趋势；
下面贴一份高斯拟合的python代码，txt文件记录了沿直 ...

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