Ultra-limited-angle CT image reconstruction algorithm based on reweighting and edge-preserving

BACKGROUND: Ultra-limited-angle image reconstruction problem with a limited-angle scanning range less than or equal to π 2 is severely ill-posed. Due to the considerably large condition number of a linear system for image reconstruction, it is extremely challenging to generate a valid reconstructed image by traditional iterative reconstruction algorithms. OBJECTIVE: To develop and test a valid ultra-limited-angle CT image reconstruction algorithm. METHODS: We propose a new optimized reconstruction model and Reweighted Alternating Edge-preserving Diffusion and Smoothing algorithm in which a reweighted method of improving the condition number is incorporated into the idea of AEDS image reconstruction algorithm. The AEDS algorithm utilizes the property of image sparsity to improve partially the results. In experiments, the different algorithms (the Pre-Landweber, AEDS algorithms and our algorithm) are used to reconstruct the Shepp-Logan phantom from the simulated projection data with noises and the flat object with a large ratio between length and width from the real projection data. PSNR and SSIM are used as the quantitative indices to evaluate quality of reconstructed images. RESULTS: Experiment results showed that for simulated projection data, our algorithm improves PSNR and SSIM from 22.46db to 39.38db and from 0.71 to 0.96, respectively. For real projection data, our algorithm yields the highest PSNR and SSIM of 30.89db and 0.88, which obtains a valid reconstructed result. CONCLUSIONS: Our algorithm successfully combines the merits of several image processing and reconstruction algorithms. Thus, our new algorithm outperforms significantly other two algorithms and is valid for ultra-limited-angle CT image reconstruction.

Two new Hager-Zhang iterative schemes with improved parameter choices for monotone nonlinear systems and their applications in compressed sensing

Notwithstanding its efficiency and nice attributes, most research on the iterative scheme by Hager and Zhang [Pac. J. Optim. 2(1) (2006) 35-58] are focused on unconstrained minimization problems. Inspired by this and recent works by Waziri et al. [Appl. Math. Comput. 361(2019) 645-660], Sabi’u et al. [Appl. Numer. Math. 153(2020) 217-233], and Sabi’u et al. [Int. J. Comput. Meth, doi:10.1142/S0219876220500437], this paper extends the Hager-Zhang (HZ) approach to nonlinear monotone systems with convex constraint. Two new HZ-type iterative methods are developed by combining the prominent projection method by Solodov and Svaiter [Springer, pp 355-369, 1998] with HZ-type search directions, which are obtained by developing two new parameter choices for the Hager-Zhang scheme. The first choice, is obtained by minimizing the condition number of a modified HZ direction matrix, while the second choice is realized using singular value analysis and minimizing the spectral condition number of the nonsingular HZ search direction matrix. Interesting properties of the schemes include solving non-smooth functions and generating descent directions. Using standard assumptions, the methods’ global convergence are obtained and numerical experiments with recent methods in the literature, indicate that the methods proposed are promising. The schemes effectiveness are further demonstrated by their applications to sparse signal and image reconstruction problems, where they outperform some recent schemes in the literature.

Five-Year Prospective Roentgen Stereophotogrammetric and Clinical Outcomes of the BioPro MTP-1 Hemiarthroplasty

Background: Mixed results for functional outcomes and long-term fixation have been reported for first metatarsophalangeal arthroplasty. This prospective study was designed to evaluate the migration of the BioPro metatarsophalangeal-1 (MTP-1) joint hemiprosthesis with Roentgen stereophotogrammetric analysis (RSA). Migration patterns of the prosthesis, prosthesis-induced erosion of the metatarsal bone, and clinical outcomes were evaluated sequentially to 5 years postoperation (PO). Methods: Eleven female patients received the BioPro-1 hemiprosthesis. Prosthesis translation and metatarsal erosion were measured with RSA at immediately PO, 6 weeks, and 3, 6, 12, 36, and 60 months postoperatively. Clinical assessment was done by patient questionnaires. Results: RSA data of 9 patients were available for analysis. Median (range) number of markers used in RSA analysis, condition number, and mean error of markers around the prosthesis were 4 (3-7), 320 (208-862), and 0.13 (0.02-0.28), respectively. Progressive subsidence was seen up to 3 years PO (mean 2.1 mm, SE 0.32). Progressive metatarsal erosion was found from 1 year PO (mean 0.49 mm, SE 0.15). Pain, function, and quality scores improved after surgery and did not deteriorate at later follow-up moments. Conclusion: Model-based RSA of the BioPro-1 prosthesis shows nonstabilizing medial and distal translation and metatarsal erosion. Despite the measured migration and erosion, clinical outcomes improved and remained similar up to 5 years postoperation.

Large-deviation asymptotics of condition numbers of random matrices

Let $\mathbf{X}$ be a $p\times n$ random matrix whose entries are independent and identically distributed real random variables with zero mean and unit variance. We study the limiting behaviors of the 2-normal condition number k(p,n) of $\mathbf{X}$ in terms of large deviations for large n, with p being fixed or $p=p(n)\rightarrow\infty$ with $p(n)=o(n)$ . We propose two main ingredients: (i) to relate the large-deviation probabilities of k(p,n) to those involving n independent and identically distributed random variables, which enables us to consider a quite general distribution of the entries (namely the sub-Gaussian distribution), and (ii) to control, for standard normal entries, the upper tail of k(p,n) using the upper tails of ratios of two independent $\chi^2$ random variables, which enables us to establish an application in statistical inference.

On solving classes of positive-definite quantum linear systems with quadratically improved runtime in the condition number

Quantum algorithms for solving the Quantum Linear System (QLS) problem are among the most investigated quantum algorithms of recent times, with potential applications including the solution of computationally intractable differential equations and speed-ups in machine learning. A fundamental parameter governing the efficiency of QLS solvers is κ, the condition number of the coefficient matrix A, as it has been known since the inception of the QLS problem that for worst-case instances the runtime scales at least linearly in κ [Harrow, Hassidim and Lloyd, PRL 103, 150502 (2009)]. However, for the case of positive-definite matrices classical algorithms can solve linear systems with a runtime scaling as κ, a quadratic improvement compared to the the indefinite case. It is then natural to ask whether QLS solvers may hold an analogous improvement. In this work we answer the question in the negative, showing that solving a QLS entails a runtime linear in κ also when A is positive definite. We then identify broad classes of positive-definite QLS where this lower bound can be circumvented and present two new quantum algorithms featuring a quadratic speed-up in κ: the first is based on efficiently implementing a matrix-block-encoding of A−1, the second constructs a decomposition of the form A=LL† to precondition the system. These methods are widely applicable and both allow to efficiently solve BQP-complete problems.

Random Toeplitz matrices: The condition number under high stochastic dependence

In this paper, we study the condition number of a random Toeplitz matrix. As a Toeplitz matrix is a diagonal constant matrix, its rows or columns cannot be stochastically independent. This situation does not permit us to use the classic strategies to analyze its minimum singular value when all the entries of a random matrix are stochastically independent. Using a circulant embedding as a decoupling technique, we break the stochastic dependence of the structure of the Toeplitz matrix and reduce the problem to analyze the extreme singular values of a random circulant matrix. A circulant matrix is, in fact, a particular case of a Toeplitz matrix, but with a more specific structure, where it is possible to obtain explicit formulas for its eigenvalues and also for its singular values. Among our results, we show the condition number of a non-symmetric random circulant matrix [Formula: see text] of dimension [Formula: see text] under the existence of the moment generating function of the random entries is [Formula: see text] with probability [Formula: see text] for any [Formula: see text], [Formula: see text]. Moreover, if the random entries only have the second moment, the condition number satisfies [Formula: see text] with probability [Formula: see text]. Also, we analyze the condition number of a random symmetric circulant matrix [Formula: see text]. For the condition number of a random (non-symmetric or symmetric) Toeplitz matrix [Formula: see text] we establish [Formula: see text], where [Formula: see text] is the minimum singular value of the matrix [Formula: see text]. The matrix [Formula: see text] is a random circulant matrix and [Formula: see text], where [Formula: see text] are deterministic matrices, [Formula: see text] indicates the conjugate transpose of [Formula: see text] and [Formula: see text] are random diagonal matrices. From random experiments, we conjecture that [Formula: see text] is well-conditioned if the moment generating function of the random entries of [Formula: see text] exists.

cdev: a ground-truth based measure to evaluate RNA-seq normalization performance

Normalization of RNA-seq data has been an active area of research since the problem was first recognized a decade ago. Despite the active development of new normalizers, their performance measures have been given little attention. To evaluate normalizers, researchers have been relying on ad hoc measures, most of which are either qualitative, potentially biased, or easily confounded by parametric choices of downstream analysis. We propose a metric called condition-number based deviation, or cdev, to quantify normalization success. cdev measures how much an expression matrix differs from another. If a ground truth normalization is given, cdev can then be used to evaluate the performance of normalizers. To establish experimental ground truth, we compiled an extensive set of public RNA-seq assays with external spike-ins. This data collection, together with cdev, provides a valuable toolset for benchmarking new and existing normalization methods.