Sparsity-constrained SENSE reconstruction: An efficient implementation using a fast composite splitting algorithm

Parallel imaging and compressed sensing have been arguably the most successful and widely used techniques for fast magnetic resonance imaging (MRI). Recent studies have shown that the combination of these two techniques is useful for solving the inverse problem of recovering the image from highly under-sampled k-space data. In sparsity-enforced sensitivity encoding (SENSE) reconstruction, the optimization problem involves data fidelity (L2-norm) constraint and a number of L1-norm regularization terms (i.e. total variation or TV, and L1 norm). This makes the optimization problem difficult to solve due to the non-smooth nature of the regularization terms. In this paper, to effectively solve the sparsity-regularized SENSE reconstruction, we utilize a new optimization method, called fast composite splitting algorithm (FCSA), which was developed for compressed sensing MRI. By using a combination of variable splitting and operator splitting techniques, the FCSA algorithm decouples the large optimization problem into TV and L1 sub-problems, which are then, solved efficiently using existing fast methods. The operator splitting separates the smooth terms from the non-smooth terms, so that both terms are treated in an efficient manner. The final solution to the SENSE reconstruction is obtained by weighted solutions to the sub-problems through an iterative optimization procedure. The FCSA-based parallel MRI technique is tested on MR brain image reconstructions at various acceleration rates and with different sampling trajectories. The results indicate that, for sparsity-regularized SENSE reconstruction, the FCSA-based method is capable of achieving significant improvements in reconstruction accuracy when compared with the state-of-the-art reconstruction method.     

Temporal breakdown and Borel resummation in the complex Langevin method

We reexamine the Parisi–Klauder conjecture for complex e^iθ/2?⁴$e^{i\theta /2}{?}^4$ measures with a Wick rotation angle 0≤θ/2≤π/2$0\le \theta /2\le \pi /2$ interpolating between Euclidean signature and Lorentzian signature. Our main result is that the asymptotics for short stochastic times t$t$ encapsulates information also about the equilibrium aspects. The moments evaluated with the complex measure and with the real measure defined by the stochastic Langevin equation have the same t→0$t\to 0$ asymptotic expansion which is shown to be Borel summable. The Borel transform correctly reproduces the time dependent moments of the complex measure for all t$t$, including their t→∞$t\to \infty$ equilibrium values. On the other hand the results of a direct numerical simulation of the Langevin moments are found to disagree from the ‘correct’ result for t$t$ larger than a finite t_c$t_c$. The breakdown time t_c$t_c$ increases powerlike for decreasing strength of the noise’s imaginary part but cannot be excluded to be finite for purely real noise. To ascertain the discrepancy we also compute the real equilibrium distribution for complex noise explicitly and verify that its moments differ from those obtained with the complex measure.     

太赫兹电磁波大气吸收衰减逐线积分计算

大气对太赫兹辐射传输存在一定的非协作性,即存在吸收衰减。为了实现太赫兹辐射的有效应用必须细致地了解太赫兹辐射大气传输的窗口位置、宽度及大气透过率。选取了处理大气非均匀路径、吸收带重叠等大气辐射传输问题的最精确方法逐线积分法,发展计算程序,并基于HITRAN分子谱线数据,对水汽、氧气、臭氧、氮气、二氧化碳等单组分气体分子对太赫兹辐射传输的吸收衰减情况进行了计算与分析,并给出了在太赫兹电磁波大气传输衰减中占主要因素的水汽和氧气的衰减峰位置。     

Volterra算子及其伴随算子线性组合的数值域

经典Volterra算子$V$及其伴随算子$V^*$在复空间$L^2[0,1]$中起着关键作用. 关于$V$和$V^*$的线性组合的性质, 我们给出了确保$z_1V+z_2V^*(z_1,z_2\in\mathbb{C})$满足增生性质的等价条件. 本文描述了$(u+iv)I+mV+nV^*(u,v,m,n\in\mathbb{R},m+n\geq0)$数值域的精确表示.     

带非线性源项的Riesz回火分数阶扩散方程的预估校正方法

本文针对带非线性源项的Riesz回火分数阶扩散方程,利用预估校正方法离散时间偏导数,并用修正的二阶Lubich回火差分算子逼近Riesz空间回火的分数阶偏导数,构造出一类新的数值格式.给出了数值格式在一定条件下的稳定性与收敛性分析,且该格式的时间与空间收敛阶均为二阶.数值试验表明数值方法是有效的.     

$p$-双调和算子的第一特征值的等周上界和 Reilly 型不等式

在本文中,我们给出了嵌入到欧氏空间中的$n$维闭超曲面上$p$-双调和算子的第一特征值的一些等周上界.我们也给出了浸入到高维流形如欧氏空间,球面和射影空间中的闭子流形上$p$-双调和算子的第一特征值的一些Reilly-型不等式.     

一个具最佳常数的复合核积分不等式

引入Γ-函数等特殊函数,利用权函数方法和实分析技巧,建立一个具最佳常数的复合核积分不等式.考虑了它的等价式,证明了它们的常数因子是最佳的,并通过取特殊的参数值,得到一些有意义的结果.     

参数型Marcinkiewicz积分交换子在Hardy空间的有界性(英文)

研究了参数型Marcinkiewicz积分高阶交换子并证明了其从H1b(Rn)到L1(Rn)有界,其中b∈BMO.     

带有两个重夸克的重子谱计算中几个困难问题的研究

由两个重夸克和一个轻夸克组成的重子可以看作是一个两体系统.它的两个重夸克组成一个玻色型的双夸克团.利用B–S方程导出了它的轻夸克和重的双夸克之间的等效相互作用势.在利用这种势计算重子质量的过程中,发现有几个困难问题需要深入探讨.它们是:(1)算符排序,(2)由非相对论展开带来的误差,(3)自旋–自旋耦合,(4)在标量双夸克组成的重子态和矢量双夸克组成的重子态之间的混合.本文详细地讨论并适当地处理了这些问题.     

超快速COSY类2D NMR实验弛豫假峰的抑制

在脉冲梯度场基础上研究了进一步缩短COSY类2D NMR实验数据采集时间的新方法.利用主动屏蔽梯度场线圈的梯度场恢复时间仅为100μs数量级这一特点,提出一类基于脉冲梯度场消除弛豫假峰的脉冲序列.并以DQF-COSY脉冲序列为例加以验证.实验结果表明:采用文中的方法,取序列重复时间远小于T₁和T₂是切实可行的;梯度场技术能很好地抑制弛豫假蜂以达到2D NMK数据的快速采集.