Single image non-uniformity correction using compressive sensing

A non-uniformity correction (NUC) method for an infrared focal plane array imaging system was proposed. The algorithm, based on compressive sensing (CS) of single image, overcame the disadvantages of “ghost artifacts” and bulk calculating costs in traditional NUC algorithms. A point-sampling matrix was designed to validate the measurements of CS on the time domain. The measurements were corrected using the midway infrared equalization algorithm, and the missing pixels were solved with the regularized orthogonal matching pursuit algorithm. Experimental results showed that the proposed method can reconstruct the entire image with only 25% pixels. A small difference was found between the correction results using 100% pixels and the reconstruction results using 40% pixels. Evaluation of the proposed method on the basis of the root-mean-square error, peak signal-to-noise ratio, and roughness index (ρ) proved the method to be robust and highly applicable.     

基于Hadamard矩阵优化排序的快速单像素成像

为提升单像素成像速度,提出了基于Hadamard矩阵优化排序的压缩采样解决方案.利用数值仿真和室外实验对提出的5种排序方法进行了对比分析.研究结果表明:按Haar小波变换系数绝对值排序时单像素成像效果最优,排序对应到Walsh序后可利用快速变换重建图像,速度达300帧/秒@64×64像素;最优排序下,采样率25%仍可重建图像,采样速度可提升4倍.针对排序方法与成像信噪比关系,从关联成像角度给出了其物理解释:测量基矩阵元邻域数值相等的区域面积等效于光场二阶相干面积,当光场二阶相干面积随测量基由大到小排序时成像效果最优.本文研究成果可用于提升单像素成像速度,具有实用价值.     

Accelerated electron paramagnetic resonance imaging using partial Fourier compressed sensing reconstruction

PurposeElectron paramagnetic resonance (EPR) imaging has evolved as a promising tool to provide non-invasive assessment of tissue oxygenation levels. Due to the extremely short T₂ relaxation time of electrons, single point imaging (SPI) is used in EPRI, limiting achievable spatial and temporal resolution. This presents a problem when attempting to measure changes in hypoxic state. In order to capture oxygen variation in hypoxic tissues and localize cycling hypoxia regions, an accelerated EPRI imaging method with minimal loss of information is needed.MethodsWe present an image acceleration technique, partial Fourier compressed sensing (PFCS), that combines compressed sensing (CS) and partial Fourier reconstruction. PFCS augments the original CS equation using conjugate symmetry information for missing measurements. To further improve image quality in order to reconstruct low-resolution EPRI images, a projection onto convex sets (POCS)-based phase map and a spherical-sampling mask are used in the reconstruction process. The PFCS technique was used in phantoms and in vivo SCC7 tumor mice to evaluate image quality and accuracy in estimating O₂ concentration.ResultsIn both phantom and in vivo experiments, PFCS demonstrated the ability to reconstruct images more accurately with at least a 4-fold acceleration compared to traditional CS. Meanwhile, PFCS is able to better preserve the distinct spatial pattern in a phantom with a spatial resolution of 0.6 mm. On phantoms containing Oxo63 solution with different oxygen concentrations, PFCS reconstructed linewidth maps that were discriminative of different O₂ concentrations. Moreover, PFCS reconstruction of partially sampled data provided a better discrimination of hypoxic and oxygenated regions in a leg tumor compared to traditional CS reconstructed images.ConclusionsEPR images with an acceleration factor of four are feasible using PFCS with reasonable assessment of tissue oxygenation. The technique can greatly enhance EPR applications and improve our understanding cycling hypoxia. Moreover this technique can be easily extended to various MRI applications.     

三种降采样成像策略的性能优化以及与传统傅里叶望远镜的比较

为了获得成像质量较好且成像时间较少的新型傅里叶望远镜成像策略,本文比较了三种降采样成像策略(压缩感知方法 (CS)、低频全采样方法 (LF)和变密度随机采样方法 (VD))与传统傅里叶望远镜(FT)在图像质量和成像时间上的差异.分析方法如下:利用传统FT外场实验所获得的目标频谱数据作为基础,三种降采样方法 (LF, VD和CS)分别按照各自的采样模式和重构方法实现目标图像的重构;通过直观观察和Strehl比两种方法比较三种降采样方法与传统FT在图像质量上的差异;通过分析成像时间的组成要素,初步比较三种降采样方法与传统FT在成像时间上的差异.分析表明:1)压缩感知方法的图像质量优于其他两种降采样方法 (LF和VD),但略低于传统成像结果; 2)压缩感知方法在成像质量上略低于传统FT,但在成像时间上却明显小于传统FT; 3)分析中采用的外场数据均含噪声,这说明上述三种降采样重构过程对噪声有较好的鲁棒性.综合上述分析结果可以看出,基于压缩感知的傅里叶望远镜(CS-FT)是在实际含噪情况下可大幅减少成像时间的优良成像策略.     

Improved simultaneous multislice magnetic resonance imaging using total variation regularization

Controlled aliasing in parallel imaging results in higher acceleration(CAIPIRINHA) for simultaneous multislice imaging has been proposed recently, which combines multiband excitation and phase cycling techniques to reduce scan time and improve subsequent imaging reconstruction. In this work, the total variation(TV) regularization method is used to further improve CAIPIRINHA. The TV regularization uses an edge-preserving prior, which establishes a relationship between neighboring pixels for image reconstruction. It reduces artifacts and suppresses noise amplification simultaneously.The results are presented with a standard eight-channel head coil with an acceleration factor of 4, where the TV-regularized CAIPIRINHA generates an improved reconstruction as compared with a typical nonregularized CAIPIRINHA.     

一种结合并行成像和压缩感知的快速磁共振成像新方法

压缩感知（CS）技术和并行成像技术（主要是SENSE技术、GRAPPA技术等）都能通过减少k空间数据的采集量来加快磁共振成像速度,目前已有一些将两种方法相结合进一步加速磁共振成像速度的方法（例如CS-GRAPPA）.本文针对数据采集和重建这两方面对现有CS-GRAPPA方法进行了改进,采集方式上采用了局部等间隔采集模板以满足GRAPPA重建的要求,并对采集模板进行随机放置以满足CS重建的要求;数据重建时,根据自动校正数据估算GRAPPA算法中欠采行的重建误差,并利用误差的大小确定在CS算法中保真的程度.不同磁共振图像重建实验的结果表明：与现有方法相比,本文方法能够更好地保留原有图像细节并有效减少伪影.     

傅里叶望远镜重构图像虚像分析

研究了采用T型激光发射阵列情况下傅里叶望远镜重构图像的虚像问题, 基于系统成像基本原理, 明确了虚像的来源. 分析认为, 发射光束在x和y轴扫描时, 在俯仰角误差的作用下, 抽取目标的空间频率和设定值之间存在一定的偏差; 当进行轴向和象限相位闭合、计算目标的单一傅里叶分量时, 该偏差会对频谱的频移造成随机的相位影响, 反映在空域上使得重构图像存在虚影现象. 在不同扫描方式情况下, 采用下一步外场实验参数, 通过计算机模拟证实了上述分析的正确性.     

A simple application of compressed sensing to further accelerate partially parallel imaging

Compressed sensing (CS) and partially parallel imaging (PPI) enable fast magnetic resonance (MR) imaging by reducing the amount of k-space data required for reconstruction. Past attempts to combine these two have been limited by the incoherent sampling requirement of CS since PPI routines typically sample on a regular (coherent) grid. Here, we developed a new method, “CS+GRAPPA,” to overcome this limitation. We decomposed sets of equidistant samples into multiple random subsets. Then, we reconstructed each subset using CS and averaged the results to get a final CS k-space reconstruction. We used both a standard CS and an edge- and joint-sparsity-guided CS reconstruction. We tested these intermediate results on both synthetic and real MR phantom data and performed a human observer experiment to determine the effectiveness of decomposition and to optimize the number of subsets. We then used these CS reconstructions to calibrate the generalized autocalibrating partially parallel acquisitions (GRAPPA) complex coil weights. In vivo parallel MR brain and heart data sets were used. An objective image quality evaluation metric, Case-PDM, was used to quantify image quality. Coherent aliasing and noise artifacts were significantly reduced using two decompositions. More decompositions further reduced coherent aliasing and noise artifacts but introduced blurring. However, the blurring was effectively minimized using our new edge- and joint-sparsity-guided CS using two decompositions. Numerical results on parallel data demonstrated that the combined method greatly improved image quality as compared to standard GRAPPA, on average halving Case-PDM scores across a range of sampling rates. The proposed technique allowed the same Case-PDM scores as standard GRAPPA using about half the number of samples. We conclude that the new method augments GRAPPA by combining it with CS, allowing CS to work even when the k-space sampling pattern is equidistant.     

小波变换在太赫兹三维成像探测内部缺陷中的应用

采用Syn View Head 300对内部有胶和空气孔的样件进行了太赫兹二维扫描(xy轴方向),系统通过线性调频连续波技术得到样件内部的三维信息.检测薄层时,由于太赫兹源的波长在亚毫米量级,薄层的上下表面反射峰相距太近而难以辨别.为了提高太赫兹探测的纵向分辨率,采用小波变换对探测信号进行处理,对小波系数进行三维重构,获得的三维小波系数图像比原始三维探测信号更加精确.该方法有效提高了太赫兹成像的纵向检测精度,纵向分辨率可达1 mm.     

鬼成像中一些数学问题

鬼成像是一种与传统成像方式不同的通过光场涨落的高阶关联获得图像信息的新型成像方式。近年来,相比传统成像方式,鬼成像所拥有的一些优点如高灵敏度、超分辨能力、抗散射等,使其在遥感、多光谱成像、热X射线衍射成像等领域得到广泛研究。随着对鬼成像的广泛研究,数学理论和方法在其中发挥的作用愈显突出。例如,基于压缩感知理论,可以进行鬼成像系统采样方式优化、图像重构算法设计及图像重构质量分析等研究工作。本文旨在探索鬼成像中的一些有趣的数学问题,主要包括:系统预处理方法、光场优化及相位恢复问题。对这些问题的研究既可以丰富鬼成像理论,又能推动它在实际应用中的发展。     

Compressed sensing sparse reconstruction for coherent field imaging

Return signal processing and reconstruction plays a pivotal role in coherent field imaging, having a significant influence on the quality of the reconstructed image. To reduce the required samples and accelerate the sampling process,we propose a genuine sparse reconstruction scheme based on compressed sensing theory. By analyzing the sparsity of the received signal in the Fourier spectrum domain, we accomplish an effective random projection and then reconstruct the return signal from as little as 10% of traditional samples, finally acquiring the target image precisely. The results of the numerical simulations and practical experiments verify the correctness of the proposed method, providing an efficient processing approach for imaging fast-moving targets in the future.     

一种基于时空变换和压缩感知的磁共振螺旋采样的图像重建算法

螺旋采样磁共振快速成像在功能性成像、并行成像和动态成像等领域发挥着越来越重要的作用.螺旋采样图像重建的传统算法是用核函数将螺旋状分布的k空间数据插值到均匀网格中,再利用傅里叶变换和最小二乘法进行重建.但是基于网格化的算法对核函数过于依赖,在网格化过程中产生难以避免的误差.该文提出了基于时空变换和压缩感知的l1范数的最优化模型和重建算法.时空变换矩阵描述了空间上的磁共振图像与采集到的时域信号间的关系,使得算法直接使用采集到的数据作为保真约束项,避免了网格化过程产生的误差.此外,基于图像处理单元的并行计算被用来提高时空变换矩阵的运算速度,使得算法具有较强的应用价值.     

Nuclear norm-regularized SENSE reconstruction

SENSitivity Encoding (SENSE) is a mathematically optimal parallel magnetic resonance (MRI) imaging technique when the coil sensitivities are known. In recent times, compressed sensing (CS)-based techniques are incorporated within the SENSE reconstruction framework to recover the underlying MR image. CS-based techniques exploit the fact that the MR images are sparse in a transform domain (e.g., wavelets). Mathematically, this leads to an l(1)-norm-regularized SENSE reconstruction. In this work, we show that instead of reconstructing the image by exploiting its transform domain sparsity, we can exploit its rank deficiency to reconstruct it. This leads to a nuclear norm-regularized SENSE problem. The reconstruction accuracy from our proposed method is the same as the l(1)-norm-regularized SENSE, but the advantage of our method is that it is about an order of magnitude faster.     

相控阵超声缺陷检测数据压缩感知

为了解决超声相控阵信号采集、存储和传输中数据量大的问题,研究了压缩感知在相控阵无损检测信号和图像压缩重构中应用的可行性。首先使用5种贪婪算法对相控阵仿真信号进行压缩重构,根据百分比均方误差选取最优算法并考虑了噪声对精度的影响,结果表明压缩感知可以用低于奈奎斯特极限的测量点数准确重构原始图像;其次用人工缺陷回波信号进行实验验证,通过稀疏度计算选择适用相控阵信号的最优稀疏基,并通过5种传感矩阵的优化选择进一步提高了重构精度。实验结果尽管达不到仿真中的理想效果,但是能以少量测量值准确恢复图像,并能保证缺陷的识别,说明压缩感知算法可以有效提高相控阵缺陷检测效率。此外,在保持测量点数相同的情况下,仿真和实验都研究了不同采样率对重构精度的影响,当测量点数超过一定值时,证实了压缩感知实际与采样率无关。      

Computational ghost imaging with deep compressed sensing

Computational ghost imaging (CGI) provides an elegant framework for indirect imaging, but its application has been restricted by low imaging performance. Herein, we propose a novel approach that significantly improves the imaging performance of CGI. In this scheme, we optimize the conventional CGI data processing algorithm by using a novel compressed sensing (CS) algorithm based on a deep convolution generative adversarial network (DCGAN). CS is used to process the data output by a conventional CGI device. The processed data are trained by a DCGAN to reconstruct the image. Qualitative and quantitative results show that this method significantly improves the quality of reconstructed images by jointly training a generator and the optimization process for reconstruction via meta-learning. Moreover, the background noise can be eliminated well by this method.     

Radial stereo imaging system for three-dimensional reconstruction

In this paper, we propose a radial stereo imaging system for three-dimensional reconstruction of an object. This system uses a novel optical device which consists of a camera and two reflected mirrors. It can reconstruct three-dimensional model of an object by capturing one image. Depending on the reflection of the mirrors, each world point in the scene has a pair of points which sit in the inner and outer of the captured image, respectively. We firstly analyze the imaging constraint of the proposed system, and then derive the geometrical relationship between world point and image point pair to calculate the three-dimensional coordinate values of an arbitrary point on the object, which means that the three-dimensional model of the object can be recovered. Experiments are carried out to prove our design, and the results show that our radial stereo imaging system is convenient and effective for the three-dimensional reconstruction of an object. This system can be used in variety of applications, such as virtual reality and commodity show.     

Continuous-wave terahertz digital holography by use of a pyroelectric array camera

Terahertz (THz) digital holography is realized based on a 2.52?THz far-IR gas laser and a commercial 124 × 124 pyroelectric array camera. Off-axis THz holograms are obtained by recording interference patterns between light passing through the sample and the reference wave. A numerical reconstruction process is performed to obtain the field distribution at the object surface. Different targets were imaged to test the system's imaging capability. Compared with THz focal plane images, the image quality of the reconstructed images are improved a lot. The results show that the system's imaging resolution can reach at least 0.4 mm. The system also has the potential for real-time imaging application. This study confirms that digital holography is a promising technique for real-time, high-resolution THz imaging, which has extensive application prospects.     

二维透射式Terahertz波时域谱成像研究

通过THz(Terahertz)波的相位信息和振幅信息来解析出样品某些可用于成像的参量;并且对THz波相位信息成像和振幅信息成像进行了比较.结果显示：当样品对THz波的吸收比较小时,利用THz波的相位信息进行成像比通过振幅信息进行成像质量要好一些;时域谱峰值位置的漂移会对相位信息成像产生较大的影响.实验中发现输出光谱中出现了劈裂现象,即由输入的单峰变成了双峰.分析表明,这一现象是由样品内部厚度或折射率突变导致的,该现象可用于THz波轮廓成像研究.     

Calculating depth maps from digital holograms using stereo disparity

Depth extraction is an important aspect of three-dimensional (3D) image processing with digital holograms and an essential step in extended focus imaging and metrology. All available depth extraction techniques with macroscopic objects are based on variance; however, the effectiveness of this is object dependent. We propose to use disparity between corresponding points in intensity reconstructions to determine depth. Our method requires a single hologram of a scene, from which we reconstruct two different perspectives. In the reconstruction the phase information is not needed, which makes this method useful for in-line digital holography. To our knowledge disparity based 3D image processing has never been proposed before for digital holography.     

Visually weighted reconstruction of compressive sensing MRI

Compressive sensing (CS) enables the reconstruction of a magnetic resonance (MR) image from undersampled data in k-space with relatively low-quality distortion when compared to the original image. In addition, CS allows the scan time to be significantly reduced. Along with a reduction in the computational overhead, we investigate an effective way to improve visual quality through the use of a weighted optimization algorithm for reconstruction after variable density random undersampling in the phase encoding direction over k-space. In contrast to conventional magnetic resonance imaging (MRI) reconstruction methods, the visual weight, in particular, the region of interest (ROI), is investigated here for quality improvement. In addition, we employ a wavelet transform to analyze the reconstructed image in the space domain and fully utilize data sparsity over the spatial and frequency domains. The visual weight is constructed by reflecting the perceptual characteristics of the human visual system (HVS), and then applied to ?₁ norm minimization, which gives priority to each coefficient during the reconstruction process. Using objective quality assessment metrics, it was found that an image reconstructed using the visual weight has higher local and global quality than those processed by conventional methods.