磁共振成像谱仪B0信号的高精度发生与测试

设计了一种基于现场可编程门阵列（FPGA）的磁共振成像（MRI）谱仪B₀信号的高精度发生方法,并对产生的B₀信号经高速采集卡采集之后进行测试和验证.FPGA从外部读取波形数据和参数,分别存储在双端口随机存取存储器（RAM）和参数寄存器中,根据预补偿算法实现B₀信号的发生,并通过对时间参数和幅度参数的控制,产生不同的B₀信号,时间分辨率为1 μs.对谱仪的B₀输出进行采集,再进行最小二乘拟合,以验证B₀信号发生的准确性.经实验证实,该设计可以产生正确、可控的高精度B₀信号.     

基于基追踪去噪的水声正交频分复用稀疏信道估计

针对传统的l₂-范数信道估计精度低的问题, 提出了一种基于基追踪去噪(BPDN)的水声正交频分复用稀疏信道估计方法, 该方法针对水声信道的稀疏特性, 利用少量的观测值即可以很高的精度估计出信道冲激响应. 与贪婪追踪类算法相比, 基于BPDN算法的稀疏信号估计具有全局最优解, 采用l₂-l₁范数准则估计信号, 同时考虑了观测值含噪情况, 通过调整正则化参数控制估计信号稀疏度和残余误差之间的平衡. 仿真分析了导频分布、正则化参数等对BPDN 算法的影响以及BPDN算法与最小平方(LS)、正交匹配追踪(OMP)信道估计算法的性能. 湖试结果表明, 在稀疏信道下, 基于BPDN的信道估计方法明显优于LS和OMP信道估计方法.     

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

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

块稀疏水声信道的改进压缩感知估计

压缩感知信道估计可利用信道稀疏特性提高估计性能,但对于具有典型块稀疏分布的水声信道,经典的l₀或l₁范数无法很好地描述块稀疏特性。利用水声信道块稀疏分布规律特性提出一种能够识别块稀疏结构的块稀疏似零范数,并在稀疏恢复信道估计算法中引入块稀疏似零范数约束项,进一步推导了复数域块稀疏似零范数恢复迭代算法,该算法通过对块稀疏似零范数进行梯度下降迭代并将梯度解投影至解空间来获得水声信道的块稀疏似零范数估计。数值仿真和海上水声通信实验结果表明该算法相对经典的稀疏信道估计算法有较明显的性能改善。通过算法推导、仿真和实验可获取结论:利用水声信道的块稀疏特性进行压缩感知重构可有效提高信道估计性能。     

基于多弛豫信号补偿的快速磁共振T1ρ散布成像

定量磁共振成像（MRI）可量化组织特性,是科学研究和临床研究的重要工具.旋转坐标系下的自旋-晶格弛豫时间（T_1ρ）能反映水与大分子之间的低频交互作用,在3 T及以上的高场环境下,T_1ρ受水和不稳定质子之间化学交换的影响较大,通过测量弛豫率随自旋锁定场强度的变化而得到其分布情况（T_1ρ散布）,可用于分析和量化质子的交换过程,因此T_1ρ散布是一种重要的定量MRI技术.然而,获得不同自旋锁定场强下T_1ρ加权图像的时间过长,限制了其应用范围.针对这一问题,本研究提出一种基于多弛豫信号补偿策略的快速T_1ρ散布成像方法.该方法将不同锁定频率下的T_1ρ加权图像补偿到同一信号强度水平,并结合低秩与稀疏建立重建模型.实验结果表明,该方法在加速倍数高达7倍时仍获得了较好的重建结果.     

双扩展水声信道的时延多普勒域稀疏估计

利用双扩展水声信道在时延-多普勒域存在的稀疏结构,将信道估计转化为压缩感知框架下的稀疏恢复问题可改善估计性能。但是,稀疏恢复经典方法如l_1范数、近似l₀范数无法适应水声信道时延-多普勒域稀疏度的动态变化,而匹配追踪(Matching pursuit,MP)、正交匹配追踪(Orthogonal Matching Pursuit,OMP)等贪婪类算法则存在着易进入局部最优解、二维搜索导致运算复杂度高等问题。提出在时延-多普勒域稀疏恢复的目标函数中引入非均匀范数约束(Non-uniform Norm Constraint,NNC),即在时延-多普勒域信道响应中根据每个时延-多普勒域位置的幅值分别分配为l₀或l₁范数约束,因而可通过不同范数约束组合的方式适应不同的时延-多普勒域稀疏度;同时,通过对非均匀范数代价函数进行梯度下降迭代求解并将梯度解投影至解空间推导了非均匀范数稀疏恢复的迭代求解方法,从而实现双扩展水声信道时延-多普勒估计。数值仿真和实验数据处理表明该算法相对经典方法有较明显的性能改善。通过仿真、海上水声通信实验结果可获取结论,利用时延-多普勒域稀疏特性的信道估计方法结合均衡器可有效提高双扩展信道条件下的水声通信性能。     

基于近似l0范数最小化的NMR波谱稀疏重建算法

在核磁共振（NMR）波谱中,过长的数据采集时间会使很多化学以及分子生物学领域的高分辨率多维谱应用难以实现. 传统的解决办法是使用随机非均匀采样代替奈奎斯特采样,但这样会使谱图质量受损. 压缩传感的出现为此提供了更好的解决办法,合适的压缩传感重建算法可以通过很少的随机非均匀采样将谱图高质量的重建出来. 该文先介绍了一种可用于谱图重建的压缩传感重建算法,名为“平滑l₀范数最小化法”,然后针对该算法对采样噪声鲁棒性较差的缺点进行了改进. 通过将改进后的算法与原算法在一维实数域信号以及NMR波谱信号重建实验中进行对比后表明,改进后的算法对噪声的鲁棒性明显提高,并能获得更好的重建性能.     

基于蒙特卡罗算法的低场核磁共振测量效率提高方法

核磁共振（NMR）的纵向弛豫时间（T₁）、横向弛豫时间（T₂）、自扩散系数（D₀）,以及T₂-T₁、T₂-D₀测量目前广泛应用于石油测井行业.在测量D₀的SGSE序列中,通过逐渐增大90°和180°脉冲之间的时间间隔（T_d）,可以对液体扩散行为产生的影响进行调节.然而T_d的"起点"、"步进数"和"终点"等参数必须设置得当才能准确测量T₁和D₀.目前参数的设置依赖多次的人工调整和测量人员的经验,耗时且使用门槛较高.本文用蒙特卡罗方法进行大量随机模拟,根据前面若干点的测量结果筛选出满足要求的随机值,预测下一个测量点的位置.该算法可以实时更新参数设置,实现自动化测量,达到降低测量门槛、缩短测量时间的目的.经验证,该算法可以适用于T₁、D₀的测量.     

基于非均匀螺旋线数据和布雷格曼迭代的快速磁共振成像方法

数据采集时间长是制约磁共振成像技术发展的重要瓶颈.为了解决这一问题, 本文基于压缩感知成像理论, 提出了一种结合非均匀螺旋线磁共振数据采集序列和布雷格曼迭代重建的快速磁共振成像方法, 通过欠采样缩短数据采集时间.欠采样引起混迭伪影则通过非均匀螺旋线欠采样特性和布雷格曼迭代重建去除.水模磁共振成像实验和在体磁共振成像实验结果表明: 欠采样情况下, 所提出的方法能有效去除欠采样导致的混迭伪影, 获得的图像结构信息完整的成像结果, 在缩短采样时间的同时, 具有较高的准确度. 关键词： 磁共振成像 非均匀螺旋线 全变分 布雷格曼迭代     

三维磁共振电阻抗成像的改进分层灵敏度算法

针对人体组织电导率的三维成像问题,提出一种改进的分层灵敏度磁共振电阻抗重建算法.利用单方向磁感应强度信息,对三维电导率图像实行分层重建,每层重建仅利用该层磁通密度分量测量数据,然后对单层重建结果进行修正以获得三维电导率重建图像.三介质长方体模型上的仿真实验证明,改进的分层重建算法改善了层间串扰现象,可以获得比一般分层算法甚至整体算法更高的图像分辨率,而且重建时间较整体算法显著减少;基于人体腿模型的仿真实验表明该算法对复杂模型三维重构的可行性;最后通过仿体实验验证算法的重建效果.改进的分层灵敏度重建算法降低了灵敏度矩阵法的计算机硬件需求,减少了重建时间,对MREIT的三维重建具有较高的成像精度和求解效率.     

Real-time dynamic vocal tract imaging using an accelerated spiral GRE sequence and low rank plus sparse reconstruction

PurposeTo develop a real-time dynamic vocal tract imaging method using an accelerated spiral GRE sequence and low rank plus sparse reconstruction.MethodsSpiral k-space sampling has high data acquisition efficiency and thus is suited for real-time dynamic imaging; further acceleration can be achieved by undersampling k-space and using a model-based reconstruction. Low rank plus sparse reconstruction is a promising method with fast computation and increased robustness to global signal changes and bulk motion, as the images are decomposed into low rank and sparse terms representing different dynamic components. However, the combination with spiral scanning has not been well studied. In this study an accelerated spiral GRE sequence was developed with an optimized low rank plus sparse reconstruction and compared with L1-SPIRiT and XD-GRASP methods. The off-resonance was also corrected using a Chebyshev approximation method to reduce blurring on a frame-by-frame basis.ResultsThe low rank plus sparse reconstruction method is sensitive to the weights of the low rank and sparse terms. The optimized reconstruction showed advantages over other methods with reduced aliasing and improved SNR. With the proposed method, spatial resolution of 1.3*1.3 mm² with 150 mm field-of-view (FOV) and temporal resolution of 30 frames-per-second (fps) was achieved with good image quality. Blurring was reduced using the Chebyshev approximation method.ConclusionThis work studies low rank plus sparse reconstruction using the spiral trajectory and demonstrates a new method for dynamic vocal tract imaging which can benefit studies of speech disorders.     

Temporal/spatial resolution improvement of in vivo DCE-MRI with compressed sensing-optimized FLASH

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides critical information regarding tumor perfusion and permeability by injecting a T(1) contrast agent, such as Gd-DTPA, and making a time-resolved measurement of signal increase. Both temporal and spatial resolutions are required to be high to achieve an accurate and reproducible estimation of tumor perfusion. However, the dynamic nature of the DCE experiment limits simultaneous improvement of temporal and spatial resolution by conventional methods. Compressed sensing (CS) has become an important tool for the acceleration of imaging times in MRI, which is achieved by enabling the reconstruction of subsampled data. Similarly, CS algorithms can be utilized to improve the temporal/spatial resolution of DCE-MRI, and several works describing retrospective simulations have demonstrated the feasibility of such improvements. In this study, the fast low angle shot sequence was modified to implement a Cartesian, CS-optimized, sub-Nyquist phase encoding acquisition/reconstruction with multiple two-dimensional slice selections and was tested on water phantoms and animal tumor models. The mean voxel-level concordance correlation coefficient for Ak(ep) values obtained from ×4 and ×8 accelerated and the fully sampled data was 0.87±0.11 and 0.83±0.11, respectively (n=6), with optimized CS parameters. In this case, the reduction of phase encoding steps made possible by CS reconstruction improved effectively the temporal/spatial resolution of DCE-MRI data using an in vivo animal tumor model (n=6) and may be useful for the investigation of accelerated acquisitions in preclinical and clinical DCE-MRI trials.     

Spatial and temporal resolution effects on dynamic contrast-enhanced magnetic resonance mammography

We tested the hypothesis that partial volume effects due to poor in-plane resolution and/or low temporal resolution used in clinical dynamic contrast-enhanced magnetic resonance imaging results in erroneous diagnostic information based on inaccurate estimates of tumor contrast agent extravasation and tested whether reduced encoding techniques can correct for dynamic data volume averaging. Image spatial resolution was reduced from 469 x 469 microm2 to those reported below by selecting a subset of k-space data. We then compared the top five K(trans)/V(T) "hot spots" obtained from the original data set, 469 x 469-microm in-plane spatial resolution and an 18-s temporal resolution processed by fast Fourier transform (FFT), with values obtained from data sets having in-plane spatial resolutions of 938 x 938, 1875 x 1875 and 2500 x 2500 microm2 and a temporal resolution of 18 s, or data sets with temporal resolutions of 36, 54 and 72 and a spatial resolution of 469 x 469 microm2, and found them to statistically differ from the parent data sets. We then tested four different post processing methods for improving the spatial resolution without sacrificing temporal resolution: zero-filled FFT, keyhole, reduced-encoding imaging by generalized-series reconstruction (RIGR) and two-reference RIGR (TRIGR). The top five values of K(trans)/V(T) obtained from data sets, the in-plane spatial resolutions of which were improved to 469 x 469 microm2 by zero-filling FFT, Keyhole and RIGR, statistically differed from those obtained from the original 469 x 469 microm2 FFT parent image data set. Only the 938 x 938 and 1875 x 1875 microm2 data sets reconstructed to 469 x 469 microm2 with TRIGR reconstruction method yielded values of the top five K(trans)/V(T) hot spots statistically the same as the original parent data set, 469 x 469 microm2 in-plane spatial and 18-s temporal-resolution FFT. That is, partial volume effects from data sets of different in-plane spatial resolution resulted in statistically different values of the top five K(trans)/V(T) hot spots relative to a high spatial and temporal resolution data set, and TRIGR reconstruction of these low resolution data sets to high resolution images provided statistically similar values with a savings in temporal resolution of 2 to 4 times.     

Interactive reduced FOV imaging for projection reconstruction and spiral acquisition.

MR fluoroscopy is likely to gain increasing importance for the visualization of dynamic processes such as cardiac function and for the guidance of interventional procedures. In many applications the dynamic processes are restricted to a part of the object under study making reduced field of view (rFOV) imaging desirable. The restriction to a smaller FOV can either be used to increase the spatial or the temporal resolution. In projection reconstruction (PR) and spiral imaging severe backfolding artifacts occur if a rFOV is used. In this paper efficient suppression schemes are proposed for PR- and spiral imaging to avoid backfolding artifacts. Evaluation of the proposed schemes was done on an interactive real-time MR-scanner. Cardiac function studies clearly showed the potential of this technique for PR- and spiral imaging.     

Retrospective cardiac gating: a review of technical aspects and future directions

The advent of short TR gradient-echo imaging has made it possible to acquire cine images of the heart with conventional whole body MRI scanners. In this paper, technical details of the data collection and image reconstruction process for cine MRI using retrospective cardiac gating are presented. Specifically, the following issues are discussed: data sorting and interpolation; time resolution; motion compensation and phase information; the type of steady state sequence including optimal flip angle; respiratory motion and correction; and the potential of 3D imaging.     

胶束型磁共振成像分子探针的设计与应用

飞速发展的分子影像学在肿瘤的早期诊断及检测中发挥着越来越重要的作用.磁共振成像(MRI)是分子影像学的重要分支，具有其他成像技术不可比拟的优越性和广阔的发展前景.它不需要放射性示踪剂，没有电离辐射，具有高的空间、时间分辨率和组织对比度.近年来，新型磁共振分子探针及成像序列取得了一系列进展，包括环境响应型分子探针、¹⁹F成像、¹²⁹Xe超极化成像以及化学交换饱和转移成像等，进一步拓展了MRI的应用范围.研究和开发靶向性好、弛豫效率高且安全性好的新型多模态MRI造影剂，进一步提高灵敏度是MRI领域的一项重要课题，例如将胶束的特性与一些MRI新方法结合，寻找合适的胶束体系，以提高MRI分子探针的灵敏度；或者引入多模态分子探针，弥补磁共振方法的不足.本文综述了胶束型MRI分子探针核心技术的研究进展与应用，并指出分子影像技术在生物医学工程研究和临床诊断中的重要性.     

高分辨率及高灵敏度小动物PET研究进展

作为一种高灵敏度且具有定量测量能力的功能分子影像技术,小动物PET越来越广泛地用于各种生物医学研究,例如疾病动物模型研究、新药物研发和新治疗方法评估等。首先回顾小动物PET成像系统的发展历史、效率和空间分辨率等性能的改进和产业化;其次,讨论了影响PET空间分辨率和效率的各种因素,包括晶体大小、探测器几何、正电子射程、光子非共线效应、图像重建算法和阻碍PET系统同时达到高空间分辨率和高效率的相互作用深度不确定效应;最后,介绍了小动物PET成像系统在以下几个方面的取得的最新进展:（1）高密度、小的光衰减常数和高光产额的闪烁晶体;（2）体积小、增益大、时间性能好、工作电压低和磁兼容的新型硅光电倍增管光探测器;（3）各种深度测量PET探测器,详细介绍了一个可达到分辨0.43 mm×0.43 mm×20 mm晶格和达到2.4 mm深度分辨率的双端读出探测器;（4）使用深度测量探测器的小动物PET成像系统,详细介绍了一个使用高分辨率双端读出探测器,全视野达到0.55 mm平均位置分辨率的小动物PET原型系统;（5）磁兼容插件式小动物PET成像系统和PET/MRI同时成像的优点;（6）小动物PET图像重建和数据校正的特点、传统的滤波反投影算法和新的迭代算法的优缺点和PET图像重建算法未来的发展方向。As the most sensitive and quantitative molecular imaging technique,small animal positron emission tomography (PET) has become a widely used tool in biomedical research such as in animal model of human disease,development of new drugs and the evaluation of new therapeutics.In this paper,first the history,the efforts to improve the spatial resolution and sensitivity as well as the commercialization process of small animal PET scanner are reviewed.Then the factors that affect the spatial resolution and sensitivity of PET scanner such as crystal size,detector geometry,positron range,photon noncollinearity and imaging reconstruction are discussed in detail.The depth of interaction effect which hinders the simultaneous achievement of PET spatial resolution and sensitivity are also discussed.Finally the recent progress made in the following areas of small animal PET instrumentation are introduced:(1) high density,short light decay constant and bright scintillator,(2) compact,high gain,good timing resolution,low bias voltage and MRI compatible silicon photomultiplier,(3) depth encoding detectors by using different methods,a detector using dual-ended readout,identifying 0.43 mm×0.43 mm×20 mm crystals and achieving a 2.4 mm depth of interaction resolution was introduced in detail,(4) small animal PET scanners using depth encoding detectors,a prototype scanner using high resolution dualended readout detectors and achieving an average of 0.55 mm spatial resolution in the whole field of view was introduced in detail.(5) MRI compatible small animal PET inserts and the advantage of simultaneous PET/MRI imaging,(6) image reconstruction and data correction of small animal PET,the filter back projection and iterative reconstruction algorithms are compared and a few key directions of PET image reconstruction will be presented.     

Referenceless PRFS MR Thermometry Using Partial Separability Model

Many areas of magnetic resonance (MR)-guided thermal therapy research would benefit from temperature maps with high spatial and temporal resolution. Conventional thermometry relies on the subtraction of baseline images, which makes it sensitive to tissue motion and frequency drift during the course of treatment. For another case is the limit of magnetic resonance imaging sampling speed, it is hard to accurately achieve MR thermometry with high spatiotemporal resolution especially for dynamic organs. To address these issues, a novel method for MR thermometry is presented by exploiting the data redundancy based on partial separability (PS) model and the referenceless thermometry. The PS model highly sparse sample two datasets in the (k–t) space for image reconstruction, which respectively determine the spatial and temporal resolutions. After the phase information is extracted from the images reconstructed by the PS model, the background phase outside the heated region from each acquired phase image through a polynomial fitting is estimated. Extrapolation of the polynomial to the heated region serves as the background phase estimate, which is then subtracted from the actual phase. The thermometry results showed that this method could accurately capture the dynamic change of MR thermometric images with 1.5 mm × 1.5 mm spatial resolution and 250 ms temporal resolution, respectively. The in vivo experiment of MR-guided high intensity focused ultrasound research and the cardiac dynamic MR thermometry are shown to demonstrate the benefits of the proposed method in high spatiotemporal resolution MR thermometry.     

时间分辨电子显微镜的数值模拟

对时间分辨电子显微镜进行了数值模拟。通过求解从样品透射出来的电子在静态磁场和动态电场的混合场中的运动,评价时间分辨电子显微镜的动态时空特性。根据该数值模拟,时间分辨电子显微镜能够在荧光屏上获得样品在不同时刻的6幅显微分幅图像。     

基于Hartley变换的磁共振成象

许多磁共振成象的应用场合需要利用正交双通道来采集一个具有时间和空间分辨力的对象序列。传统的基于Ｆｏｕｒｉｅｒ变换的成象方法，一方面，图象序列的重建时各帧图象是独立地进行重建的，因而图象序列的时间分辨力受到编码的限制；另一方面，来自两个通道之间的Ｆａｒｔｌｅｙ变换的磁共振成象技术。