基于压缩感知的差分关联成像方案研究

关联成像可提供一种运用常规手段难以获得清晰图像的方法, 能够解决一些常规成像技术不易解决的问题, 是近些年来量子光学领域的前沿和热点之一.本文提出一种基于压缩感知的差分关联成像方案(简称, 差分压缩关联成像方案), 将高斯分布的热光源强度分布作为压缩感知的测量矩阵, 差分物体信息作为被成像物体信息, 利用差分关联成像方案的高成像信噪比和压缩感知技术的低采样次数, 通过正交匹配追踪算法, 高质量地恢复出物体信息. 并以二灰度“双缝”、“NUPT”, 多灰度Lena图和Boats图为例, 数值仿真差分压缩关联成像过程; 同时将本方案350次测量的结果与差分关联成像方案30000次测量的结果进行对比, 研究结果表明针对不同的被成像物体(二灰度“双缝”、“NUPT”, 以及多灰度Lena图和Boats图), 10次成像的均方误差平均值分别降低了97.7%, 93.9%, 92.5%和71.4%; 与压缩鬼成像方案相比, 同样测量次数条件下均方误差值对于二灰度双缝和多灰度Lena图、Boats图等目标物 体分别有50.4%, 72.9%和66.8%的降低. 差分压缩关联成像方案极大地提高了成像信噪比, 降低了成像时间. 关键词： 关联成像 差分 压缩感知 均方误差     

合成孔径压缩感知超声成像中的高效能稀疏字典设计

针对目前常见的稀疏字典缺乏针对性,在合成孔径医学超声成像中的应用效果不佳,难以在低压缩率下保证重构图像质量的问题,本文设计了一种高效能的稀疏字典。根据超声回波信号是由发射脉冲信号经过不同延时衰减后叠加的特点,利用发射脉冲作为基函数构造稀疏字典,回波信号在该稀疏字典确定的变换域中具备很好的稀疏性,理论上能使其稀疏表示系数的稀疏度等于超声阵元接收到的反射回波数。通过FieldⅡ对简单点目标和复杂目标的仿真结果表明:在相同的重构算法和压缩率下该稀疏字典重构的平均绝对误差明显小于常见的稀疏字典,其值仅为DWT的几分之一,DFT和DCT的几十分之一,能让回波信号以更低的压缩率实现相同的恢复效果。本文最后使用体模的实际采集数据对算法的实际效果进行检测,实验结果也与仿真结果基本一致。基于该稀疏字典的压缩感知算法可以进一步减少合成孔径成像所需存储的数据量、降低系统的复杂度。     

压缩感知理论在光学成像中的应用

压缩感知以信号的稀疏性或可压缩性为条件,以远低于耐奎斯特采样频率对信号数据进行采样和编码。简要概括了压缩感知的基本理论,它采用非自适应线性投影来保持信号的原始结构, 能通过数值最优化问题精确或高概率地重构原始信号。详细介绍了其在光学成像系统中的应用,主要包括单像素相机、超薄成像、编码孔径成像、多路技术智能成像、多光谱成像和CMOS成像等成像系统。最后对该理论的应用前景进行了阐述。     

基于压缩感知的菲涅尔孔径编码无透镜成像（特邀）

为了降低菲涅尔孔径编码成像的传输带宽,改善小尺寸像感器拼接成像质量,提出菲涅尔孔径编码成像的压缩感知重建算法。基于压缩感知理论框架对部分采样的编码图像重建可行性进行分析,指出压缩感知重建误差随波带片常数缩小而降低。根据重建图像中孪生像和原始像在梯度域稀疏性的差异,引入全变差正则化实现抑制孪生像的效果。建立压缩感知理论框架下图像重建目标函数,利用交替方向乘子法进行求解。结合编码图像的能量分布和采样模式的可实现性,测试了矩形采样和辐射线采样两种模式,对不同采样数据量下重建图像质量进行分析。结果表明,辐射线采样模式相比矩形采样模式具有更高的图像采样效率,且仅通过7.3%的实验测量数据就可以获得质量良好的图像,为多块小尺寸像感器拼接成像提供了理论基础,有利于拓展编码掩膜无透镜成像的应用领域。     

一种数字微镜阵列分区控制和超分辨重建的压缩感知成像法

提出一种压缩感知成像框架结构.该结构采样端用新建的采样矩阵实现数字微镜阵列分区控制,可增强信息获取的准确性,测量得到与新数字微镜阵列对应的压缩采样值;重构端由采样值优化重构出低分辨率图像后,根据分区控制过程建立压缩感知理论框架下的超分辨重建模型,利用梯度稀疏约束优化算法进行求解,恢复出原高分辨率图像.实验结果表明:数字微镜阵列分区控制与超分辨重建相结合的方法可以明显降低压缩感知成像系统的计算量,缩短成像时间,并且具有较高的图像重构质量.     

最小均方误差解调算法在彩色超声血流成像中的研究

研究了最小均方误差正交解调算法在超声血流成像中的应用。首先由设计的期望信号与接收信号在最小均方误差原则下得到其迭代解调形式,获得回波信号同相分量和正交分量,然后由得到的正交解调信号通过自相关的方法对血流流速进行估计。解调算法的仿真结果表明,对高斯噪声信噪比为0.5—10dB的正弦波调制信号,解调输出平均信噪比与Hilbert变换法和I/Q解调法相比分别提高了15dB和4dB;血流成像的仿真结果表明,在流速估计性能相当的情况下,解调的乘法运算量仅分别为上述对比方法的18%和9%。因此在超声血流成像中应用最小均方误差正交解调算法,对于提高估计性能和降低运算量都有一定意义。     

基于空间广义相干系数加权的医学合成孔径超声成像

为解决广义相干系数用于合成孔径成像中所存在的运算量大,图像对比度改善有限等问题,提出空间广义相干系数加权成像方法。该算法根据单个孔径成像结果之间的相干性来计算相干系数,通过加权空间合成进行成像。采用对FieldⅡ仿真点目标和吸声斑目标的数据进行成像表明,算法不仅使运算量减少N(N为阵元数)倍,而且相对于传统的广义相干系数算法,算法对散射点成像可提高信噪比7 dB,对于吸声斑成像可提高对比度3.2 dB。可见算法可以有效地提高成像速度和成像质量。     

The application of compressed sensing in synthetic transmit aperture medical ultrasound imaging

In synthetic transmit aperture medical ultrasound imaging field,a compressed sensing ultrasound imaging method based on the sparsity in frequency domain is presented in order to reduce huge amount of data and large numbers of receiving channels.First,the sparsity in frequency domain is verified.Then the echo signal is compressively sampled in time-spatial domain based on compressed sensing and the echo signal is reconstructed by solving an optimization problem.Finally the image is made by using the synthetic transmit aperture approach.The experiments based on point target and fetus target are used to verify the proposed method.The MSE,resolution and image quality of reconstructed image and those of original image are compared and analyzed.The results show that only 30%amount of data and 50%of receiving channels were used to implement ultrasound imaging without reducing the quality of image in experiment.The amount of data and the complexity of system are reduced greatly by the proposed method based on compressed sensing.     

Modified synthetic transmit aperture algorithm for ultrasound imaging

The modified synthetic transmit aperture (STA) algorithm is described. The primary goal of this work was to assess the possibility to improve the image quality achievable using synthetic aperture (SA) approach and to evaluate the performance and the clinical applicability of the modified algorithm using phantoms. The modified algorithm is based on the coherent summation of back-scattered RF echo signals with weights calculated for each point in the image and for all possible combinations of the transmit-receive pairs. The weights are calculated using the angular directivity functions of the transmit-receive elements, which are approximated by a far-field radiation pattern of a narrow strip transducer element vibrating with uniform pressure amplitude over its width. In this way, the algorithm takes into account the finite aperture of each individual element in the imaging transducer array. The performance of the approach developed was tested using FIELD II simulated synthetic aperture data of the point reflectors, which allowed the visualization (penetration) depth and lateral resolution to be estimated. Also, both simulated and measured data of cyst phantom were used for qualitative assessment of the imaging contrast improvement. The experimental data were obtained using 128 elements, 4 MHz, linear transducer array of the Ultrasonix research platform. The comparison of the results obtained using the modified and conventional (unweighted) STA algorithms revealed that the modified STA exhibited an increase in the penetration depth accompanied by a minor, yet discernible upon the closer examination, degradation in lateral resolution, mainly in the proximity of the transducer aperture. Overall, however, a considerable (12 dB) improvement in the image quality, particularly in the immediate vicinity of the transducer’s surface was demonstrated. The modified STA method holds promise to be of clinical importance, especially in the applications where the quality of the “near-field” image, that is the image in the immediate vicinity of the scanhead is of critical importance such as for instance in skin- and breast-examinations.     

Optimization of sparse synthetic transmit aperture imaging with coded excitation and frequency division

An effective aperture approach is used for optimization of a sparse synthetic transmit aperture (STA) imaging system with coded excitation and frequency division. A new two-stage algorithm is proposed for optimization of both the positions of the transmit elements and the weights of the receive elements. In order to increase the signal-to-noise ratio in a synthetic aperture system, temporal encoding of the excitation signals is employed. When comparing the excitation by linear frequency modulation (LFM) signals and phase shift key modulation (PSKM) signals, the analysis shows that chirps are better for excitation, since at the output of a compression filter the sidelobes generated are much smaller than those produced by the binary PSKM signals. Here, an implementation of a fast STA imaging is studied by spatial encoding with frequency division of the LFM signals. The proposed system employs a 64-element array with only four active elements used during transmit. The two-dimensional point spread function (PSF) produced by such a sparse STA system is compared to the PSF produced by an equivalent phased array system, using the Field II simulation program. The analysis demonstrates the superiority of the new sparse STA imaging system while using coded excitation and frequency division. Compared to a conventional phased array imaging system, this system acquires images of equivalent quality 60 times faster, when the transmit elements are fired in pairs consecutively and the power level used during transmit is very low. The fastest acquisition time is achieved when all transmit elements are fired simultaneously, which improves detectability, but at the cost of a slight degradation of the axial resolution. In real-time implementation, however, it must be borne in mind that the frame rate of a STA imaging system depends not only on the acquisition time of the data but also on the processing time needed for image reconstruction. Comparing to phased array imaging, a significant increase in the frame rate of a STA imaging system is possible if and only if an equivalent time efficient algorithm is used for image reconstruction.     

基于压缩感知理论的逆合成孔径雷达抗欺骗干扰方法

针对逆合成孔径雷达（ISAR）成像对抗中的基于数字图像合成（DIS）技术欺骗干扰方法,在分析现有基于频率调制率微调的抗DIS欺骗干扰方法工程实现方面缺陷的基础上,通过将基于压缩感知理论的随机抽取采样和稀疏重构方法融入到信号接收和处理部分中,提出了基于压缩感知理论的改进抗DIS欺骗干扰方法,进一步提高基于频率调制率微调抗DIS欺骗干扰方法的工程实用性。仿真实验结果显示了基于压缩感知理论的改进抗DIS欺骗干扰方法的有效性。     

Compressed sensing in synthetic aperture photoacoustic tomography based on a linear-array ultrasound transducer

Photoacoustic tomography(PAT) has the unique capability of visualizing optical absorption inside several centimeters-deep biological tissue with a high spatial resolution. However, single linear-array transducer-based PAT suffers from the limited-view challenge, and thus the synthetic aperture configuration is designed that still requires multichannel data acquisition hardware. Herein, a feasible synthetic aperture PAT based on compressed sensing reconstruction is proposed. Both the simulation and experimental results tested the theoretical model and validated that this approach can improve the image resolution and address the limited-view problem while preserving the target information with a fewer number of measurements.     

压缩感知光谱成像技术的编码孔径与探测器匹配优化

编码孔径光谱成像仪在实际应用中存在着编码模板与探测器分辨率不匹配从而降低系统分辨率的问题。针对该问题进行了两种情况分析,并通过数学理论建模给出了相应的解决方案。对于编码模板分辨率高于探测器分辨率这一情况,提出引入邻域嵌入超分辨技术的方法,实现了基于压缩感知的超分辨光谱成像。对于编码模板分辨率低于探测器分辨率这一情况,提出区块阈值划分的编码孔径,将编码微元按照区块阈值重新划分并进行灰度分级,从而实现低分辨率编码模板的高分辨率编码孔径。利用梯度投影稀疏重构(GPSR)算法进行数据立方体重建,实验结果表明:运用基于超分辨理论的编码孔径快照光谱成像系统所测得的光谱图像更精准,内容更丰富;采用基于区块阈值划分的编码孔径的编码孔径快照光谱成像系统具有更高的空间分辨率和光谱分辨率。结果证实优化后的编码孔径快照光谱成像系统,其分辨率和成像质量大幅度提升,并实现了高分辨率元件的100%利用。     

Synthetic aperture ultrasound imaging

The paper describes the use of synthetic aperture (SA) imaging in medical ultrasound. SA imaging is a radical break with today's commercial systems, where the image is acquired sequentially one image line at a time. This puts a strict limit on the frame rate and the possibility of acquiring a sufficient amount of data for high precision flow estimation. These constrictions can be lifted by employing SA imaging. Here data is acquired simultaneously from all directions over a number of emissions, and the full image can be reconstructed from this data. The paper demonstrates the many benefits of SA imaging. Due to the complete data set, it is possible to have both dynamic transmit and receive focusing to improve contrast and resolution. It is also possible to improve penetration depth by employing codes during ultrasound transmission. Data sets for vector flow imaging can be acquired using short imaging sequences, whereby both the correct velocity magnitude and angle can be estimated. A number of examples of both phantom and in vivo SA images will be presented measured by the experimental ultrasound scanner RASMUS to demonstrate the many benefits of SA imaging.     

Three-dimensional synthetic aperture integral imaging

We propose synthetic aperture integral imaging, in which an effectively enlarged aperture (or field of view) is obtained by movement of small integral imaging system. This system substantially increases the field of view and the viewing resolution. The feasibility of our approach is experimentally demonstrated. To the best of our knowledge, this is the first time the synthetic aperture technique has been applied to three-dimensional integral imaging.