增量维纳滤波法在波前探测解卷积中的应用

一个常规的自适应光学系统通常包含三个重要环节：波前探测、波前校正和波前重构。因此对系统的技术要求非常高,造成系统复杂,成本昂贵。基于哈特曼-夏克波前探测的图像解卷积处理就是“事后”处理的一种,它省去了波前校正环节,使常规的自适应光学系统得到简化,降低了系统成本。其基本原理为对瞬时波面进行短时间曝光探测,同时记录相应的短时间曝光图像,来进行解卷积处理。将增量维纳滤波法应用于基于哈特曼-夏克波前探测的解卷积中,并对室内模拟点源情况下的三组畸变光斑图像数据进行了解卷积恢复处理。结果表明,将增量维纳滤波法应用于基于波前探测的解卷积是完全可行的,在室内模拟点源情况下,恢复的图像可以达到衍射极限分辨力。与维纳滤波相比,它扩展了噪声抑制因子的选取范围,在噪声抑制因子选取不正确的情况下,仍能得到比维纳滤波更好的结果。     

Hybrid Deconvolution of Adaptive Optics Retinal Images from Wavefront Sensing

Adaptive optics can be used to compensate for the wave aberration of the human eyes to achieve high-resolution imaging in real time. However the correction ＆ partial due to the limitation of hardware. We propose a kind of hybrid image post-processing method, which uses the blind deconvolution combined with the residual data in wavefront sensor to restore the partially adaptive optics corrected retinal image. This method is applied in the image restoration of the vivid human retinal images. The results show that it is effective to improve the retinal image quality.     

基于目标和点扩展函数联合估计的点源目标图像近视解卷积

在传统的基于波前探测的解卷积方法中,由波前探测得到的点扩展函数被认为是精确的,并用维纳滤波进行复原,但是点扩展函数不可避免地存在误差,所以最终的复原目标图像质量不佳.为了解决该难题,提出了基于目标和点扩展函数联合估计的图像近视解卷积算法.它运用了点扩展函数和目标的先验信息,对点扩展函数和目标进行了规整和进一步约束,从而得到更优的恢复图像质量.对该方法的原理和实现过程进行了阐述,并将其运用于室内点源目标数据中.实验结果证明,与维纳滤波方法相比,该方法使图像恢复的效果得到明显改善.     

Estimation of the ocular point spread function by retina modeling

Retinal imaging often suffers from blurring aberrations. With knowledge of the blurring point spread function (PSF), better images can be reconstructed by deconvolution techniques. We demonstrate a method to enhance the contrast of retinal cells by estimating the ocular PSF. This is done by finding the cells' positions and their intensity distributions and using these as a model for the image. The feasibility of this method is demonstrated by Wiener deconvolution both for adaptively and nonadaptively corrected images.     

Adaptive optics retinal image registration from scale-invariant feature transform

With the use of adaptive optics(AO), high-resolution microscopic imaging of the living human retina in the single cell level has been achieved. In an AO retinal imaging system, with a small field size (about 1°, 300 μm) the motion of the eye severely affects the stabilization of the real-time video images results in significant distortions of the retina images. Scale-invariant feature transform (SIFT) algorithm is applied to automatically abstract corner points with subpixel resolution and match the points in two frames. With the matched corner points, we estimate and remove the motions of 20 frames of photoreceptor cells and capillary blood vessels, respectively. The maximal translational motion is about 30 and 44 pixels in the 20 frames whose size is 416 × 416 pixels. More general motions can be considered by the SIFT algorithm, but only simple translational motion can be considered by cross-correlation algorithm.     

Simulated annealing in ocular adaptive optics

We present what is to our knowledge a first hardware realization of a simulated annealing algorithm in an adaptive optics system designed to image the retina of the human eye. The algorithm is applied to the retinal image itself without the need for wavefront sensors in the system. We find that this optimization algorithm can be an alternative to the traditional Hartmann-Shack sensing. We also compare the simulated annealing algorithm to the stochastic parallel gradient descent algorithm.     

双树复数小波变换的视网膜图像半盲解卷积复原

针对导致自适应光学视网膜图像降质退化的原因,提出了一种结合双树复数小波变换(DTCWT)和图像半盲解卷积复原算法的方法。首先,对经过自适应光学实时校正技术得到的视网膜图像进行DT-CWT分解,得到低频和高频部分应的图像。将自适应光学成像系统中残余像差重建的光学传递函数作为图像复原模型的初始估计点扩散函数(PSF),并对低频部分图像进行条件约束的迭代半盲解卷积复原;对高频部分的图像进行去噪处理。最后,将处理后的高频和低频部分图像进行双树复数小波逆变换,获得复原图像。实验和结果表明:由该方法处理的视网膜细胞图像质量得到明显提高,图像客观质量评价参数相对于原始图像提高了5倍多;在视网膜细胞的空间频率范围内(70~90(°)~1),复原图像功率谱平均值提高了5倍左右,有助于对视网膜细胞的高分辨率观察。     

Phase diversity for calibrating noncommon path aberrations of adaptive optics system under nonideal measurement environment

When the defocus cannot be measured and the wavefront solution set is restricted by a multi-channel, some practical problems exist in the calibration of the noncommon path aberrations of the adaptive optics system. To solve these problems, an evaluation function of phase diversity algorithm is constructed in this paper. We use the method that the estimated aberration and the modulated deformable mirror iterate each other to make up the nonideal measurement environment. Then the ill-pose problem of the solution by phase diversity, produced as relaxing constraints of the diversity defocus on the wavefront solution set, is solved. We have adopted the proposed method to measure the noncommon path aberrations of the adaptive optics system on a 1.23 m telescope. Experimental results demonstrate that wavefront solution is more accurate and the whole imaging quality is improved effectively by using the deformable mirror to compensate the aberration measured.     

Roughness measurement of metals using a modified binary speckle image and adaptive optics

This paper proposes an integrated roughness measurement system that is based on adaptive optics (AO) and binary analysis of speckle pattern images. The aim of this study was to demonstrate the necessity for AO compensation in regions containing both heat and fluid flow turbulences. A speckle image was obtained by projecting a laser beam onto the specimen surface, and the laser pattern image reflected from the surface was binarized to experimentally correlate the intensity with the surface roughness. In the absence of the AO correction scheme, induced turbulences can severely increase the residual rms error from 0.14 to 1.4 μm. After a real-time closed-loop AO correction, we can reduce the wavefront root mean square (rms) error to 0.12 μm, which not only compensates for the aberration error from induced disturbances but also improves the overall performance of the optical system. In addition, an AO system having different gains was investigated, and a threshold gain value was found to be able to steadily compensate for the wavefront errors in less than 2 s. Measurement results of five steel samples having roughness ranging from 0.2 to 3.125 μm (0.3λ and 5λ, where λ is the diode laser wavelength) demonstrate an excellent correlation between the intensity distribution of binary images and average roughness with a correlation coefficient of 0.9982. Furthermore, the proposed AO-assisted system is in good agreement with the stylus method and less than 9.73% error values can be consistently obtained.     

Deconvolution of vibroacoustic images using a simulation model based on a three dimensional point spread function

Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalcifications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7–0.8 mm. Current VA systems based on confocal or linear array transducers generate C-scan images at the beam focal plane. Images on the axial plane are also possible, however the system resolution along depth worsens when compared to the lateral one. Typical axial resolution is about 1.0 cm. Furthermore, the elevation resolution of linear array systems is larger than that in lateral direction. This asymmetry degrades C-scan images obtained using linear arrays. The purpose of this article is to study VA image restoration based on a 3D point spread function (PSF) using classical deconvolution algorithms: Wiener, constrained least-squares (CLSs), and geometric mean filters. To assess the filters’ performance on the restored images, we use an image quality index that accounts for correlation loss, luminance and contrast distortion. Results for simulated VA images show that the quality index achieved with the Wiener filter is 0.9 (when the index is 1.0 this indicates perfect restoration). This filter yielded the best result in comparison with the other ones. Moreover, the deconvolution algorithms were applied to an experimental VA image of a phantom composed of three stretched 0.5 mm wires. Experiments were performed using transducer driven at two frequencies, 3075 kHz and 3125 kHz, which resulted in the difference-frequency of 50 kHz. Restorations with the theoretical line spread function (LSF) did not recover sufficient information to identify the wires in the images. However, using an estimated LSF the obtained results displayed enough information to spot the wires in the images. It is demonstrated that the phase of the theoretical and the experimental PSFs are dissimilar. This fact prevents VA image restoration with the current theoretical PSF. This study is a preliminary step towards understanding the restoration of VA images through the application of deconvolution filters.     

波前解卷积方法中的高频噪声抑制

 介绍了一种波前解卷积中噪声抑制规整化的新方法，并将此方法应用于室内模拟点源实验中。该方法通过在图像复原算法中增加针对图像高频部分的限制条件来抑制高频噪声，以达到对图像复原问题病态特性的规整化。实验结果表明：该规整化方法可以有效地抑制解卷积过程中高频噪声的影响，恢复出达到理论衍射极限分辨率的图像。对于噪声水平较高的降质图像，通过这种解卷积方法可以有效地提高信噪比。同维纳逆滤波方法相比，该方法可以在有效抑制导致病态的高频噪声的基础上充分保持图像的低频；与基于贝叶斯估计的近视解卷积算法相比，该方法不需要知道噪声水平或噪声类型等先验知识，只是从噪声本质出发，通过抑制降质图像高频部分，有效地解决了病态特性问题。     

A study on a fast measuring technique of wavefront using a Shack–Hartmann sensor

The measuring resolution and speed for wavefront are important to improve the performance of an adaptive optics system. In this paper, we propose a fast measuring technique with high resolution in the Shack–Hartmann wavefront sensor. The proposed measuring technique of wavefront combines the conventional center of mass algorithm with an estimated weighting factor. The estimated weighting factor is a real value anticipating the real center of mass in a wavefront spot image. This estimated weighting factor can be calculated at the initialization step of an adaptive optics system. We designed a robust adaptive optics system with this proposed measuring algorithm. The measuring accuracy and speed are investigated by using the proposed algorithm compared with the conventional center of mass algorithm in experiments.     

基于位相变更的非相干数字全息自适应成像

菲涅耳非相干数字全息作为一种非扫描的三维成像技术具有其独特的优势,但其成像过程中会受到各种像差的影响,导致成像分辨率、再现像的质量降低.为了解决这一问题,可以结合适当的自适应光学技术对波前像差进行探测和校正.位相变更是一种基于两幅具有已知位相差的强度图像实现波前探测和像差校正的技术.本文发展了基于位相变更的非相干数字全息自适应成像技术,不需要引入引导星,利用全息记录过程中的两幅相移全息图,实现波前像差的探测.本文给出了所发展技术的数值仿真和实验结果,结合位相变更算法求解出系统像差的位相分布,将像差的共轭位相加载到光瞳面上,在全息图记录的同时校正像差,从而提高重建像的质量.     

Estimation of speed of sound in dual-layered media using medical ultrasound image deconvolution

The speed of sound in soft tissues is usually assumed to be 1540 m/s in medical pulse-echo ultrasound imaging systems. When the true speed is different, the mismatch can lead to distortions in the acquired images, and so reduce their clinical value. Previously we reported a new method of sound-speed estimation in the context of image deconvolution. Unlike most other sound-speed estimation methods, this enables the use of unmodified ultrasound machines and a normal scanning pattern. Our approach was validated for largely homogeneous media with single sound speeds. In this article, we demonstrate that sound speeds of dual-layered media can also be estimated through image deconvolution. An ultrasound simulator has been developed for layered media assuming that, for moderate speed differences, the reflection at the interface may be neglected. We have applied our dual-layer algorithm to simulations and in vitro phantoms. The speed of the top layer is estimated by our aforesaid method for homogeneous media. Then, when the layer boundary position is known, a series of deconvolutions are carried out with dual-layered point-spread functions having different lower-layer speeds. The best restoration is selected using a correlation metric. The error level (e.g., a mean error of −9 m/s with a standard deviation of 16 m/s) for in vitro phantoms is found to be not as good as that of our single-speed algorithm, but is comparable to other local speed estimation methods where the data acquisition may not be as simple as in our proposed method.     

开环双脉冲液晶自适应光学视网膜成像系统

为了获得高分辨率视网膜图像,利用液晶空间光调制器作为波前校正器建立了一套开环液晶自适应光学视网膜成像系统。与闭环模式相比,采用开环模式后,系统的能量利用率提高了1倍。系统采用双脉冲照明方式,以减少人眼曝光量,保护人眼安全。在照明光学系统中加入了大小视场切换装置使成像视场由之前的0.8°增至1.7°。同时优化了系统的时序控制流程,对人眼像差连续校正的同时快速调节成像相机的前后位置至最佳像面。对于开环模式对动态人眼像差的校正精度进行了测量,实验测得,经开环校正后,残差波面的均方根值约为0.09λ;相应的斯特雷尔(Strehl)比高于0.70,系统分辨率接近光学衍射极限的分辨率。对两名志愿者进行了实验,获得了清晰的眼底视网膜细胞图像。     

探测器噪声对波前编码成像系统图像复原过程影响的理论研究

波前编码光学系统其成像过程由光学成像和图像复原两部分组成,光学成像过程探测器上所成的编码图像含有噪声,该噪声在图像复原过程将被放大,为了分析探测器噪声对波前编码系统在图像复原过程中产生的影响,通过添加高斯随机噪声来模拟探测器噪声,使用直接逆滤波、维纳滤波和Lucy—Richardson滤波算法,对波前编码成像系统焦点位...     

A weighted least squares algorithm for time-of-flight depth image denoising

In this paper, we present a new denoising method for the depth image of a time-of-flight (ToF) camera, based on weighted least squares (WLS) framework. The common method for ToF depth image denoising is to use bilateral filter. However, the ability of bilateral filter in edge preservation would be reduced while we attempt to smooth out larger spatial scale noise. In order to avoid this problem and preserve the edge information as much as possible, we introduce a new way to construct edge-preserving ToF depth image denoising based on WLS. We are to our knowledge the first to present a WLS-based method for ToF depth image denoising. Experimental results demonstrate that compared with bilateral filter, our proposed algorithm not only achieves better performance in edge preservation, but also improves the PSNR values of the denoised images by 0.5–2.6 dB.     

人眼视网膜成像自适应光学系统设计

设计一套基于液晶空间光调制器的人眼视网膜成像自适应光学系统,以获得高分辨率视网膜图像,并且使该系统实现体积小,功耗低,成本低等优点.采用夏克-哈特曼探测器和基于硅基板上的液晶器件分别作为波前探测器和波前校正器.系统采用双对准光源以主观方式来使人眼对准,近红外光探测成像以减小对人眼的刺激.使人眼对有限距离对焦,以减小离焦对成像的影响,使该系统既可用于正常眼,又可用于近视眼.用ZEMAX软件对系统进行了模拟分析,认为该系统可获得高于3 μm的视网膜分辨率,该系统设计是合理可行的.     

基于哈特曼-夏克波前传感器的模糊图像复原方法

离焦模糊图像的清晰度较低,因此必须对其进行复原。传统方法通常采用圆盘或高斯函数来近似离焦造成的点扩散函数,复原效果不够理想。为此,提出利用哈特曼-夏克波前传感器探测离焦波前,根据所得波前计算光学系统的点扩散函数,并采用Richardson-Lucy算法对模糊图像进行复原。搭建了实验用的光学系统,采集了离焦模糊图像以及相应的波前信息,获得了清晰的复原图像,并利用客观图像评价方法对退化图像和复原图像进行了评价,同时与传统方法得到的复原图像进行了比较。实验结果表明,该方法能精确重建点扩散函数,有效提高图像的质量。     

A piecewise local regularized Richardson-Lucy algorithm for remote sensing image deconvolution

As it is not possible to obtain an accurate point spread function (PSF) in remote sensing imaging, classic deconvolution methods such as Wiener filtering often introduce strong noise and ringing artifacts, which contaminate the restored images. In this paper, we modify the standard Richarson-Lucy (RL) algorithm with a piecewise local regularization term and combine it with residual deconvolution method. Experimental results show that it is effective in suppressing negative effects, and images with rich details and sharp edges are obtained.