Targeted-ROI imaging in electron paramagnetic resonance imaging

Electron paramagnetic resonance imaging (EPRI) is a technique that has been used for in vivo oxygen imaging of small animals. In continuous wave (CW) EPRI, the measurement can be interpreted as a sampled 4D Radon transform of the image function. The conventional filtered-backprojection (FBP) algorithm has been used widely for reconstructing images from full knowledge of the Radon transform acquired in CW EPRI. In practical applications of CW EPRI, one often is interested in information only in a region of interest (ROI) within the imaged subject. It is desirable to accurately reconstruct an ROI image only from partial knowledge of the Radon transform because acquisition of the partial data set can lead to considerable reduction of imaging time. The conventional FBP algorithm cannot, however, reconstruct accurate ROI images from partial knowledge of the Radon transform of even dimension. In this work, we describe two new algorithms, which are referred to as the backprojection filtration (BPF) and minimum-data filtered-backprojection (MDFBP) algorithms, for accurate ROI-image reconstruction from a partial Radon transform (or, truncated Radon transform) in CW EPRI. We have also performed numerical studies in the context of ROI-image reconstruction of a synthetic 2D image with density similar to that found in a small animal EPRI. This demonstrates both the inadequacy of the conventional FBP algorithm and the success of BPF and MDFBP algorithms in ROI reconstruction. The proposed ROI imaging approach promises a means to substantially reduce image acquisition time in CW EPRI.     

Exploiting rank deficiency and transform domain sparsity for MR image reconstruction

The reconstruction of magnetic resonance (MR) images from the partial samples of their k-space data using compressed sensing (CS)-based methods has generated a lot of interest in recent years. To reconstruct the MR images, these techniques exploit the sparsity of the image in a transform domain (wavelets, total variation, etc.). In a recent work, it has been shown that it is also possible to reconstruct MR images by exploiting their rank deficiency. In this work, it will be shown that, instead of exploiting the sparsity of the image or rank deficiency alone, better reconstruction results can be achieved by combining transform domain sparsity with rank deficiency.To reconstruct an MR image using its transform domain sparsity and its rank deficiency, this work proposes a combined l₁-norm (of the transform coefficients) and nuclear norm (of the MR image matrix) minimization problem. Since such an optimization problem has not been encountered before, this work proposes and derives a first-order algorithm to solve it.The reconstruction results show that the proposed approach yields significant improvements, in terms of both visual quality as well as the signal to noise ratio, over previous works that reconstruct MR images either by exploiting rank deficiency or by the standard CS-based technique popularly known as the ‘Sparse MRI.’     

双重分数傅里叶变换计算全息

提出了双重分数傅里叶变换计算全息，在这种方法中，将两个图像的信息分别经不同阶的分数傅里叶变换后，记录在同一张分数傅里叶变换计算全息图上，它需要两个特定的分数傅里叶变换系统才能再现出所记录的图像信息，利用其再现方式的特殊性，可制成一种新的安全认证系统。     

基于空域非均匀傅里叶变换的傅里叶望远镜

为了在稀疏发射阵列下清晰重构目标图像,提出了一种基于空域非均匀傅里叶变换(NDFT)的傅里叶望远镜信号处理方法。依据傅里叶望远镜的发射器位置与抽取的目标空间频率关系,结合MATLAB程序特点,完成了空域非均匀傅里叶逆变换,重构了目标图像。稀疏发射阵列配置方式为:T型阵列单臂放置11个发射望远镜,连续抽取目标的8个低频信息,再抽取3个高频分量。选择不同形状和灰度分布的4个卫星作为成像目标。与补零均匀快速傅里叶变换(FFT)方法重构的图像对比发现:信噪比为100 dB时,相比补零均匀FFT方法, NDFT方法重构图像的Strehl比都有所提升,最高提升了0.159 8。     

Synthetic Aperture Radar image nonlinear enhancement algorithm based on NSCT transform

Due to the imaging mechanism, Synthetic Aperture Radar (SAR) images are susceptible to speckle noise, which affects radar image interpretation. So image denoising and enhancement are important topics of improving SAR image performance. A nonlinear image enhancement algorithm based on nonsubsampled contourlet transform (NSCT) is proposed in this paper. The image is decomposed into coefficients of different scales and directions through nonsubsampled contourlet transform. It is denoised by the threshold method of the multi-scale product of NSCT coefficients. Then thresholds of the nonlinear enhancement function are determined according to the coefficients of each scale. The two parameters of the function, among which one is used to control the range of enhancement and the other can determine the strength of enhancement, are obtained by solving nonlinear equations. The coefficients processed by the enhancement function are used to reconstruct the image. The simulation results on the Matlab platform show that the algorithm has a good effect of enhancing details of images and suppressing noise signals meanwhile.     

螺旋MRI的网格化数据重建算法比较

螺旋MRI的原始数据是在不均匀的 k- 空间螺旋轨迹上采样得到的,需要通过网格化算法等手段将数据变成等间距的网格数据后,才能采用FFT进行重建,最后得到供临床使用的图像. 本文对Jackson网格化算法和Claudia大矩阵算法的重建速度和图像结果进行了比较, 并得出以下结论：1）在获得相近图像质量的情况下,Claudia大矩阵重采样算法比Jackson算法要快且更方便在仪器上实现. 2）在Jackson双倍细网格算法的实现方式中,数据驱动插值比网格驱动插值更有效率. 3）在Claudia大矩阵重采样算法中,对冗余比大于310∶1的数据进行图像重建的时候,网格点的幅值不平均化比平均化后的效果还要好. 这几个结论都将有利于MRI图像重建技术的进一步提高.      

Optical image compression based on adaptive directional prediction discrete wavelet transform

The traditional lifting wavelet transform cannot effectively reconstruct the nonhorizontal and nonvertical high-frequency information of an image. In this paper, we present a new image compression method based on adaptive directional prediction discrete wavelet transform (ADP-DWT). We first design a directional prediction model to obtain the optimal transform direction of the lifting wavelet. Then, we execute the directional lifting transform along the optimal transform direction. The edge and texture energy can be reduced in the nonhorizontal and nonvertical directions of the high-frequency sub-bands. Finally, the wavelet coefficients are coded with the set partitioning in hierarchical trees (SPIHT) algorithm. The new method holds the advantages of both adaptive directional lifting (ADL) and direction-adaptive discrete wavelet transform (DA-DWT), and the computational complexity is far lower than that in these methods. For the images containing regular and fine textures or edges, the coding preformance of ADP-DWT is better than that of ADL and DA-DWT.     

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

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

总变差约束的数据分离最小图像重建模型及其Chambolle-Pock求解算法

基于优化的迭代法,可以结合压缩感知和低秩矩阵等稀疏优化技术高精度地重建图像.其中,总变差最小（total variation minimization,TV）模型是一种简单有效的优化模型.传统的约束TV模型,使用数据保真项为约束项,TV正则项为目标函数.本文研究TV约束的、数据分离最小（TV constrained,data divergence minimization,TVcDM）新型TV模型及其求解算法.详细推导了TVcDM模型的Chambolle-Pock（CP）算法,验证了模型及算法的正确性;分析了算法的收敛行为;评估了模型的稀疏重建能力;分析了模型参数的选择对重建的影响及算法参数对收敛速率的影响.研究表明,TVcDM模型有高精度稀疏重建能力;TVcDM-CP算法确保收敛,但迭代过程中有振荡现象;TV限对重建有重要影响,参数值过大会引入噪声而过小会模糊图像细节;算法参数的不同选取会导致不同的收敛速率.     

复杂刚性运动物体的高精度三维测量算法

基于相移法的三维形貌重建精度高,对环境噪声和阴影等不敏感,但由于多幅条纹解相位,难以应用于动态物体的三维测量中,为此,提出了一种新的算法.基于Harris算法提取刚性运动棋盘格的角点,确定相邻两帧采集条纹图像之间的像素偏差并校正采集条纹图像;根据投影仪和摄像机的标定参数建立投影图像和采集图像之间的空间变换矩阵,并根据变...     

Time-domain surface profile imaging via a hyperspectral Fourier transform spectrometer

We describe a method for time-domain surface profile measurements via white-light reflection spectroscopy using a hyperspectral Fourier transform spectrometer (HS-FTS). This technique measures the frequency of the spectral modulation of reflected light from a multilayer optical surface and reports the spatially resolved optical thickness. Owing to the Fourier relationship, the Fourier transform spectrometer manifests this spectral modulation as temporal satellites in interferogram space. We show that measurement of the positions of these satellites can be used to reconstruct the optical thickness profile over a surface using the HS-FTS.     

被动式太赫兹图像目标检测研究

研究提出一种被动式太赫兹扫描成像的目标检测方法。用0.22 THz被动式扫描成像系统采集原始图像,在滤波去噪后用针孔像分析法获得太赫兹扫描成像系统的点扩展函数,进而用Lucy-Richardson算法重构图像,通过灰度变换和边缘检测增强图像对比度和目标分辨力。实验结果表明,算法能有效改善被动式太赫兹图像质量,提高成像系统探测隐藏可疑物的能力。     

An exact, linear solution to the problem of imaging through turbulence

We show how, in principle, to solve the ‘blind deconvolution' problem. This is in the context of the problem of imaging through atmospheric turbulence. The approach is digital but not iterative, and requires as input data but two short-exposure intensity images, without the need for reference point sources. By taking the Fourier transform of each image and dividing, a set of linear equations is generated whose unknowns are sampled values of the two random point spread functions that degraded the images. An oversampling by 50% in Fourier space equalizes the number of unknowns and independent equations. With some prior knowledge of spread function support, and in the absence of added noise of image detection, the inverted equations give exact solutions. The two observed images are then inverse filtered to reconstruct the object.     

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

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

Multiple linear feature detection through top-hat transform by using multi linear structuring elements

A multiple linear feature detection algorithm through top-hat transform using the constructed multiple linear structuring elements is proposed in this paper. The desired linear features are treated as a set. And, the set is divided into different subsets. Multi linear structuring elements corresponding to different subsets are constructed. After that, top-hat transform is performed by using the constructed linear structuring elements, and the results are combined to reconstruct the desired linear features. Then, the extracted linear features are binarized and processed to form the final detected binary linear features. Because of the effective performance of the top-hat transform using the constructed multi linear structuring elements, the linear features of different images from different applications could be well detected. The analysis and experimental results show that, the proposed algorithm could be well used for multiple linear feature detection in different applications.     

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

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

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.     

基于峰值信噪比和小波方向特性的图像奇异值去噪技术

提出一种利用小波变换子图像不同的方向特性和峰值信噪比进行奇异值分解的图像去噪算法。由于图像经过小波变换后,低频子图像集中了原图像的大部分能量噪声,故仅作简单维纳滤波;而噪声则主要集中在小波域中的三个不同方向的高频子图中,且系数较小,因此可以利用奇异值分解进行去噪处理,即用较大的奇异值和对应的特征向量重构出去噪图像,然而由于奇异值分解固有的行列方向性,对于高频对角线子图重构出的图像去噪效果不理想,故采取旋转至行列方向后再进行常用的奇异值滤波;最后将去噪后的低频和高频子图进行小波反变换重构出最终的去噪图像,其中重构所需的奇异值个数由图像的峰值信噪比确定。 实验结果表明,该方法在有效去噪的同时较好的保留了原有的高频细节信息。     

工业CT图像轮廓提取与优化系统设计

为了有效地提高面向逆向工程的工业CT图像轮廓处理的精准性和计算效率,开发了工业CT图像轮廓提取与优化系统。该系统主要包括图像处理模块、轮廓识别与提取模块、轮廓拟合与优化模块、图像处理经典算法库和数据接口等功能模块,实现了从工业CT图像到轮廓数据的三维点云与轮廓线的标准格式输出,并与逆向工程软件对接。同时,在轮廓拟合与优化模块中,提出了基于全局约束模型的轮廓拟合策略,并采用有向图分解和数值代数相结合的计算方法,提高了轮廓拟合的计算精度和效率。最后通过实例验证了系统的有效性和正确性。     

Improvement of reconstruction in 3D computer-generated holograms using 1D Fourier transform operations

Improvement of images reconstructed from 3D computer-generated holograms using 1D Fourier transform operations is presented. Required computing time for making 3D holograms is decreased when 1D Fourier transform method is used instead of 2D Fourier transform method, but reconstructed images from the holograms created by the 1D Fourier transform method are somewhat inferior in quality to the images by 2D Fourier transform method. In this paper, we propose a method that the disadvantageous properties of the holograms made by 1D Fourier transform are improved. We have performed two step improvements in this experiment. From the numerical experiments, the images that are reconstructed from the holograms created from 1D Fourier transform method show marked improvement.