点光源哈特曼最优阈值估计方法研究

针对夏克-哈特曼波前传感器探测系统中噪声随时间及空间变化频率较快的特点,为了准确估计系统的最优阈值,根据高斯光斑与噪声的分布特性,提出一种以滑动窗口内像素均值及图像信号的局部梯度作为参数,构造关于噪声权重函数的方法,由此获得子孔径阈值的最优估计值,并详细分析了算法的基本原理和实现过程.以典型处理方法获取的阈值与理论最优阈值的误差作为评价标准,仿真和实验结果表明本文提出的阈值估计方法在不同信噪比、不同光斑大小的条件下,均能取得优于典型阈值处理方法获得的结果,且与理论最优阈值的误差小于10%.

Shack-Hartmann optimum threshold estimation for the point source

The Shack-Hartmann wavefront sensor (SHWFS) is an optical detection device based on the measurements of wavefront slopes. It is widely used in an adaptive optics system due to its simple structure and strong environment adaptability. The measuring accuracy of the SHWFS depends mainly on the accuracy of the spot image centroid in each sub-aperture. There are many centroid algorithms including the center of gravity algorithm, Gauss fitting algorithm, and correlation algorithm. As to the simplicity, robustness, high accuracy and stability, the center of gravity algorithm is more widely used. However, the accuracy of gravity algorithm is sensitive to the noise including discretization, aliasing, photon noise, readout noise, stray light, and direct current bias. To improve the accuracy of centroid, the output signals of SHWFS must be pre-processed to suppress the noise effect by using the method of thresholding in general. Many threshold methods have been presented to reduce the error of centroid and there theoretically exists an optimum threshold which causes the minimum error of centroid based on the characteristics of SHWFS and noise. However, it is difficult to separate the signals from the noises, and the optimum threshold cannot be estimated accurately in real time in the SHWFS systems. In this paper aiming at noises in SHWFS, which vary with time and space rapidly, a method based on the noise weighted function of the mean value of pixels and the local gradient direction of image signals in the moving windows is presented according to the characteristics of the Gaussian spot and noise distributions. Moreover, the theory and parameters determination of the method are analyzed. The method utilizes the probability that the pixels in the moving windows belong to the noise, and the probability is inversely proportional to the mean value of pixels and the local gradient direction of image signals, and so the monotonically reducing probability function of pixels is constructed. Finally, the standard deviation and mean value of noise can be obtained, and the estimation value of optimum threshold is equal to the mean value of noise plus three times the standard deviation of noise. To investigate the effects of the optimum threshold estimation with the different spot sizes, spot strengths and noise levels, the proposed algorithm is compared with traditional methods. The simulation and experimental results show that the proposed method could achieve higher accuracy, and the error between the threshold obtained by the method presented in this paper and theoretical optimum threshold is less than 10%, which is less than those from the traditional methods.