Effect of intensity quantization level in parallel phase-shifting digital holography

Parallel phase-shifting digital holography (PPSDH) enables the instantaneous recording of three-dimensional fields. The three-dimensional field can be reconstructed using a computer by numerical scalar wave propagation. In PPSDH, we record a space-division multiplexed hologram that includes the required phase retardation in the spatial distribution of the reference wave and then interpolate the data at blanked pixels for each phase retardation to obtain the complex amplitude distribution of an object wave. The recorded quality of the multiplexed hologram influences the reconstruction quality. In this study, we investigate the effect of the intensity quantization of a multiplexed hologram on the reconstruction quality. We compare the influence of intensity quantization in PPSDH with that in the conventional phase-shifting method. Random noise is also added to the multiplexed hologram. The required intensity quantization level is helpful for selecting a digital image sensor.     

Phase deviation analysis and phase retrieval for partial intensity saturation in phase-shifting projected fringe profilometry

Intensity saturation may take place as a sinusoidal fringe pattern is projected onto an object which has a relative high reflective index on some regions of the surface. If a phase-shifting method is used, the illuming light intensities at the same point of the object may have different values for each projected phase-shifting fringe pattern. Therefore, when the intensity at a point of the object exceeds the saturation level for a certain phase step of projected fringe pattern, other intensities obtained at the same point with different phase steps of projected fringe patterns do not exceed the saturation level. This kind of saturation is defined as partial intensity saturation. In the case of partial intensity saturation, a relative larger phase error is introduced when a conventional phase reconstruction algorithm is used in the phase-shifting projected fringe profilometry (PSPFP). To get rid of such disadvantage, an improved algorithm for phase reconstruction is proposed in this paper. By using the new algorithm, the effect of partial intensity saturation for the phase reconstruction can be decreased greatly and a good quality of reconstructed phase map can be obtained. Phase deviation of the reconstructed phase is also analyzed. Finally, an experimental result with PSPFP is presented to validate the feasibility of the proposed algorithm.     

Terahertz imaging with compressed sensing and phase retrieval

We describe a novel, high-speed pulsed terahertz (THz) Fourier imaging system based on compressed sensing (CS), a new signal processing theory, which allows image reconstruction with fewer samples than traditionally required. Using CS, we successfully reconstruct a 64 x 64 image of an object with pixel size 1.4 mm using a randomly chosen subset of the 4096 pixels, which defines the image in the Fourier plane, and observe improved reconstruction quality when we apply phase correction. For our chosen image, only about 12% of the pixels are required for reassembling the image. In combination with phase retrieval, our system has the capability to reconstruct images with only a small subset of Fourier amplitude measurements and thus has potential application in THz imaging with cw sources.     

A modified discrete algebraic reconstruction technique for multiple grey image reconstruction for limited angle range tomography

The `missing wedge', which is due to a restricted rotation range, is a major challenge for quantitative analysis of an object using tomography. With prior knowledge of the grey levels, the discrete algebraic reconstruction technique (DART) is able to reconstruct objects accurately with projections in a limited angle range. However, the quality of the reconstructions declines as the number of grey levels increases. In this paper, a modified DART (MDART) was proposed, in which each independent region of homogeneous material was chosen as a research object, instead of the grey values. The grey values of each discrete region were estimated according to the solution of the linear projection equations. The iterative process of boundary pixels updating and correcting the grey values of each region was executed alternately. Simulation experiments of binary phantoms as well as multiple grey phantoms show that MDART is capable of achieving high‐quality reconstructions with projections in a limited angle range. The interesting advancement of MDART is that neither prior knowledge of the grey values nor the number of grey levels is necessary.     

一种三维物体相息图的迭代计算方法

在传统迭代傅里叶变换算法的基础上,提出了一种计算三维物体相息图的新方法.基于层析法将三维物体的多个分层物面作为衍射再现图像,在一个输入面(相息图)和多个输出面(再现像)之间进行迭代.通过在傅里叶迭代运算中引入距离相位因子,表示物体不同物面的深度,体现了物体的三维特征.实验结果证明了本文算法良好的收敛特性和再现性能.最后,分析了物面数量和间距对全息再现质量的影响,利用液晶空间光调制器采用时分复用的方法还原了三维物体的多个物面.     

基于Poisson-Markov分布最大后验概率的多通道超分辨率盲复原算法

在Poisson和Markov随机场分布假设基础上,改进了存在噪声时多通道互限制原则,提出了一种基于最大后验概率判据的多通道图像盲复原算法,可在无需已知通道内离散点扩散函数性质、类型和具体分布等先验知识的情况下,在图像幅度限制约束空间和离散点扩散函数能量恒定约束空间内,通过对图像和离散点扩散函数的交替迭代,使两者收敛到全局最优解,最终复原出超分辨率图像.对计算机模拟和实际存在微平移、微离焦情况的图像处理表明：算法在不同信噪比水平和非命中PSF尺寸下具有良好的图像复原效果和稳定性. 关键词： 泊松-马尔科夫分布 最大后验概率 多通道 盲复原     

Adaptive window iteration algorithm for enhancing 3D shape recovery from image focus

Depth from focus(DFF) is a technique for estimating the depth and three-dimensional(3D) shape of an object from a multi-focus image sequence. At present, focus evaluation algorithms based on DFF technology will always cause inaccuracies in deep map recovery from image focus. There are two main reasons behind this issue. The first is that the window size of the focus evaluation operator has been fixed. Therefore, for some pixels, enough neighbor information cannot be covered in a fixed window and is easily disturbed by noise, which results in distortion of the model. For other pixels, the fixed window is too large, which increases the computational burden. The second is the level of difficulty to get the full focus pixels, even though the focus evaluation calculation in the actual calculation process has been completed. In order to overcome these problems, an adaptive window iteration algorithm is proposed to enhance image focus for accurate depth estimation. This algorithm will automatically adjust the window size based on gray differences in a window that aims to solve the fixed window problem.Besides that, it will also iterate evaluation values to enhance the focus evaluation of each pixel. Comparative analysis of the evaluation indicators and model quality has shown the effectiveness of the proposed adaptive window iteration algorithm.     

An adaptive squared-distance-weighted interpolation for volume reconstruction in 3D freehand ultrasound

Volume reconstruction is a key procedure in 3D ultrasound imaging. An algorithm named as squared-distance-weighted (SDW) interpolation has been earlier proposed to reduce the blurring effect in the 3D ultrasonic images caused by the conventional distance weighted (DW) interpolation. However, the SDW parameter alpha, which controls the weight distribution, is a constant assigned by operators so that the interpolation effect is invariant for both sharp edges and speckle noises. In this paper, we introduced a new adaptive algorithm based on SDW interpolation for volume reconstruction of 3D freehand ultrasound. In the algorithm, the local statistics of pixels surrounding each voxel grid were used to adaptively adjust the parameter alpha in SDW. The voxel grids with a higher ratio of local variance and mean in their neighbourhoods would have a smaller alpha to make the image details sharper, while the voxel grids locating in regions with a lower ratio of local variance and mean would have a larger alpha to smooth image content in homogeneous regions, where speckle noise is usually observed and damages the image quality. By comparing the simulation results using the SDW and new adaptive algorithm, it was demonstrated that this new algorithm worked well in both edge preservation and speckle reduction.     

代数重建技术在板壳结构断层重建中的应用

针对板壳结构用锥束射线倾斜扫描时数据缺失等问题,利用代数重建技术(ART)进行了断层重建。计算机模拟结果表明,由ART重建出的切片图像质量远好于先前用滤波反投影(FBP)重建出的切片图像质量。对影响ART算法的各种因素进行了调整,即选取适当的初始图像和收敛因子就可获得最佳的图像质量。     

X-ray interferometric Fourier holography

An X-ray interferometric Fourier holography was proposed and theoretically investigated. X-ray Interferometric Young fringes and the reconstruction of an object image were investigated by the Fourier transform method. It was shown that on the output surface of the analyzer crystal (the third plate of the interferometer) the interference pattern of two slits gives X-ray interferometric Young fringes. An expression for the period of X-ray interferometric Young fringes was obtained. The subsequent reconstruction of the slit image as an object is conducted by means of the Fourier transform method of intensity distribution on the hologram. Three methods for reconstruction of the function of complex transmission of the object are presented: an analytical one–the approximate method, the iteration method and the step-by-step approach. As examples a recording of X-ray interferometric Fourier hologram and the reconstruction of the function of complex amplitude transmission of a beryllium circular cylinder are given.     

基于极大似然模型和期望最大化算法的闪光图像重建

针对闪光照相图像低信噪比的特点,研究了带约束的基于极大似然模型和期望最大化算法（ML-EM）的闪光照相图像重建算法。该算法在ML-EM迭代重建算法的基础上,根据闪光照相图像低信噪比的特点及被重建客体具有分段平滑的先验信息,在迭代重建的过程中进行了相应的噪声抑制,抑制了噪声对重建结果的影响,同时很好地保持了客体的边界特征。数值模拟结果表明,基于噪声约束的ML-EM重建算法能取得较好的重建效果。     

余弦二值编码纯相位全息图的数字微镜器件显示

分析了采用错误减算法的迭代过程,在已知物波函数傅里叶谱的振幅和物波函数振幅的情况下恢复出纯相位的物波函数,最大限度地保留物波的振幅及相位信息.提出采用余弦二值编码生成二值全息图,即全息图的透射率函数取0或1.二值全息图通过数字微镜器件全息显示系统进行了重构显示,重构效果很好.理论分析了数字微镜器件的衍射效率,表明其最大衍射效率仅和微镜之间的间隔尺寸与微镜边长之比有关.余弦二值编码方法从理论上消除了零级衍射,可以制作像素较多的全息图.     

四光束剪切相干成像目标重构算法研究

针对传统剪切光束成像技术的准实时性问题,提出用口字形排布的四束光代替传统L形三束剪切光照射目标,研究了四光束剪切相干成像目标重构算法.只需单次测量就能同时重构出四幅目标图像,减少了用于降低散斑噪声、获取高质量图像所需的测量次数,同时大大减少了多组发射时的光束切换次数,提高了成像效率.在算法实现中,通过最小二乘法恢复出四组波前相位,利用散斑幅值的简单代数运算恢复波前幅值,从而重构出目标图像.仿真结果表明,与传统方法相比,在图像质量相同的前提下,本文方法所需的数据采集时间减少了至少1/2,不但提高了目标重构效率,还可为远程运动目标的成像识别提供更好的手段.     

两步相移数字全息物光重建误差分析与校正

系统研究了两步相移数字全息干涉术中相移误差引起的波前再现误差的计算和校正方法. 基于衍射物光相位分布的随机性和振幅相位的相互独立性原理,介绍了相移数字全息中物光波前再现误差的表达形式,推导出步长为π/2的两步算法中物光重建误差的表达式. 通过进一步分析这一重建误差的结构和特点,结合物光表达式,给出了自动校正相移误差引起的波前重建误差的校正方法. 该方法无需增加测量,在未知相移误差大小的情况下,只对标准两步相移算法恢复的物光复振幅进行处理就可以实现对物光振幅和相位的同时校正. 计算机模拟结果表明,校正后可将 关键词： 相移干涉术 数字全息 物光重建 误差校正     

A comparison of iterative Fourier transform algorithms for image quality estimation

A comparison was established between two iterative Fourier transform algorithms (IFTAs), such as the original Gerchberg–Saxton (GS) and the mixed-region amplitude freedom (MRAF) algorithms, for the hologram reconstruction of different target images through the full reference image quality estimation (IQE) and pixel homogeneity in the Fourier plane presented theoretically and experimentally. The comparison was applied depending upon both algorithms based on a computer-generated hologram (CGH) implemented utilizing a reflective phase-modulated liquid-crystal spatial light modulator (LC-SLM) to obtain the digital kinoform holograms of the desired intensity distributions. These digital holograms were applied to reconstruct the intensity patterns for 852 nm, which represents a laser beam source. The theoretical and experimental results of the reconstructed patterns obtained using the MRAF algorithm were found to be smoother and better than the patterns obtained using the GS algorithm. Unmodulated light beam (dc term) is removed from the reconstructed patterns attributed to digital kinoform holograms of MRAF algorithm as an alternative to the theoretical and experimental results without using any additional optic equipment at the light path. Moreover, this paper discussed the full reference objective quality estimations, such as mean square error (MSE), peak signal-to-noise ratio (PSNR), structural content (SC), normalized absolute error (NAE), normalized cross correlation (NK), and homogeneity of pixels, through the contrast (Cont) and inverse difference moment (IDM) for numerical and experimental results. According to the two desired intensity distributions processed theoretically and experimentally, the results of MRAF algorithm were found to be in the highly accurate recovered phase, the quality of image was enhanced, and the dc term was decreased. Image quality estimation of full reference objective relay on the feedback algorithms experimental attestation has not been implemented yet.     

Speckle suppression and contrast enhancement in reconstruction of freehand 3D ultrasound images using an adaptive distance-weighted method

In three-dimensional (3D) freehand ultrasound (US), reconstructing a set of B-scans into a regular voxel array is the key procedure for consequent visualization and analysis. This paper presents a new adaptive interpolation algorithm for computing the voxel array to suppress speckle noises and enhance contrast. The local statistics of homogeneous regions including mean and variance were measured and the ratio of variance to mean was used as homogeneity criteria. For the computation of each voxel, the interpolation method was adaptively determined with respect to its local statistics. If the neighbouring pixels of a voxel satisfied the homogeneity criterion, its value was computed with an arithmetic mean filter. Otherwise, the voxel was probably locating in an inhomogeneous region and an adaptive distance-weighted (ADW) interpolation method was employed to compute its value. A resolution phantom and a subject’s forearm were reconstructed using the proposed algorithm and two other well-known methods - conventional distance-weighted (DW) and voxel nearest neighbourhood (VNN) interpolations. The comparison results demonstrated that the adaptive interpolation algorithm was able to suppress speckles, preserve edges and enhance contrast effectively for the volume reconstruction.     

动态背景下基于光流场分析的运动目标检测算法

针对现有动态背景下运动目标检测算法的不足,提出一种基于光流场分析的运动目标检测算法.首先根据前背景在光流梯度幅值和光流矢量方向上的差异确定目标的大致边界,然后通过点在多边形内部原理获得边界内部的稀疏像素点,最后以超像素为节点,利用混合高斯模型拟合的表观信息和超像素的时空邻域关系构建马尔可夫随机场模型的能量函数,并通过使目标函数能量最小化得到最终的运动目标检测结果.该算法不需要任何先验假设,能够同时处理动态背景和静态背景两种情况.多组实验结果表明,本文算法在检测的准确性和处理速度上均优于现有算法.     

Sub-Rayleigh imaging via undersampling scanning based on sparsity constraints

We demonstrate that, by undersampling scanning object with a reconstruction algorithm related to compressed sensing, an image with the resolution exceeding the finest resolution defined by the numerical aperture of the system can be obtained. Experimental results show that the measurements needed to achieve sub-Rayleigh resolution enhancement can be less than 10% of the pixels of the object. This method offers a general approach applicable to point-by-point illumination super-resolution techniques.     

一种用于在线三维测量的五步非等步相移算法

提出了一种五步非等步相移算法,并成功应用于在线三维测量中。仅投影一固定正弦光栅条纹到物体上,通过物体在线运动产生等效相移,在一个条纹周期范围内任意采集五帧变形条纹图,通过像素匹配使各帧条纹图中的物点一一对应并计算出相应的等效相移量。采用所提出的五步非等步相移算法,即可重构物体。计算机模拟与实验验证了该方法的可行性和有效...     

层析成像图像重建算法综述

介绍了层析成像技术的图像重建算法,并从正向问题数学模型的简化和反向问题数学模型的映射结构的角度比较了各种算法的特点和优劣。研究表明:用本质是线性算法的各种变换方法重建图像存在严重失真,而卷积滤波的引入可以使变换方法的重建效果有所改善;基于导数搜索的迭代算法对初始值依赖性强、收敛速度慢并且容易陷入局部最优解;基于Fourier变换的方法具有本质的局限性;小波变换则可以同时刻画图像时域和频域的细节特征;有限元法通过重建对象像素的智能划分可以简化正问题的复杂性;而具有物理背景的蒙特卡罗法、模拟退火法、遗传算法、粒子滤波法及神经网络法更适合于复杂且非线性的图像重建;智能化、仿生化、并行化以及各种算法的融合是层析成像图像重建算法的发展趋势。