Source author record

Ben Cox

Ben Cox appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Phase-insensitive versus phase-sensitive ultrasound absorption tomography in the frequency domain

The sensitivity of phase-sensitive detectors, such as piezoelectric detectors, becomes increasingly directional as the detector element size increases. In contrast, pyroelectric sensors, which are phase-insensitive, retain their omni-directionality even for large element sizes, although they have significantly poorer temporal resolution. This study uses numerical models to examine whether phase-insensitive detectors can be used advantageously in ultrasound tomography, specifically absorption tomography, when the number of detectors is sparse. We present measurement models for phase-sensitive and phase-insensitive sensors and compare the quality of the absorption reconstructions between these sensor types based on image contrast metrics. We perform the inversion using synthetic data with a Jacobian-based linearised matrix inversion approach.

preprint2020arXiv

Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast

Ultrasound Tomography has seen a revival of interest in the past decade, especially for breast imaging, due to improvements in both ultrasound and computing hardware. In particular, three-dimensional ultrasound tomography, a fully tomographic method in which the medium to be imaged is surrounded by ultrasound transducers, has become feasible. In this paper, a comprehensive derivation and study of a robust framework for large-scale bent-ray ultrasound tomography in 3D for a hemispherical detector array is presented. Two ray-tracing approaches are derived and compared. More significantly, the problem of linking the rays between emitters and receivers, which is challenging in 3D due to the high number of degrees of freedom for the trajectory of rays, is analysed both as a minimisation and as a root-finding problem. The ray-linking problem is parameterised for a convex detection surface and three robust, accurate, and efficient ray-linking algorithms are formulated and demonstrated. To stabilise these methods, novel adaptive-smoothing approaches are proposed that control the conditioning of the update matrices to ensure accurate linking. The nonlinear UST problem of estimating the sound speed was recast as a series of linearised subproblems, each solved using the above algorithms and within a steepest descent scheme. The whole imaging algorithm was demonstrated to be robust and accurate on realistic data simulated using a full-wave acoustic model and an anatomical breast phantom, and incorporating the errors due to time-of-flight picking that would be present with measured data. This method can used to provide a low-artefact, quantitatively accurate, 3D sound speed maps. In addition to being useful in their own right, such 3D sound speed maps can be used to initialise full-wave inversion methods, or as an input to photoacoustic tomography reconstructions.

preprint2020arXiv

Towards accurate quantitative photoacoustic imaging: learning vascular blood oxygen saturation in 3D

Significance: 2D fully convolutional neural networks have been shown capable of producing maps of sO$_2$ from 2D simulated images of simple tissue models. However, their potential to produce accurate estimates in vivo is uncertain as they are limited by the 2D nature of the training data when the problem is inherently 3D, and they have not been tested with realistic images. Aim: To demonstrate the capability of deep neural networks to process whole 3D images and output 3D maps of vascular sO$_2$ from realistic tissue models/images. Approach: Two separate fully convolutional neural networks were trained to produce 3D maps of vascular blood oxygen saturation and vessel positions from multiwavelength simulated images of tissue models. Results: The mean of the absolute difference between the true mean vessel sO$_2$ and the network output for 40 examples was 4.4% and the standard deviation was 4.5%. Conclusions: 3D fully convolutional networks were shown capable of producing accurate sO$_2$ maps using the full extent of spatial information contained within 3D images generated under conditions mimicking real imaging scenarios. This work demonstrates that networks can cope with some of the confounding effects present in real images such as limited-view artefacts, and have the potential to produce accurate estimates in vivo.

preprint2018arXiv

Approximate k-space models and Deep Learning for fast photoacoustic reconstruction

We present a framework for accelerated iterative reconstructions using a fast and approximate forward model that is based on k-space methods for photoacoustic tomography. The approximate model introduces aliasing artefacts in the gradient information for the iterative reconstruction, but these artefacts are highly structured and we can train a CNN that can use the approximate information to perform an iterative reconstruction. We show feasibility of the method for human in-vivo measurements in a limited-view geometry. The proposed method is able to produce superior results to total variation reconstructions with a speed-up of 32 times.

preprint2018arXiv

Enhancing Compressed Sensing 4D Photoacoustic Tomography by Simultaneous Motion Estimation

A crucial limitation of current high-resolution 3D photoacoustic tomography (PAT) devices that employ sequential scanning is their long acquisition time. In previous work, we demonstrated how to use compressed sensing techniques to improve upon this: images with good spatial resolution and contrast can be obtained from suitably sub-sampled PAT data acquired by novel acoustic scanning systems if sparsity-constrained image reconstruction techniques such as total variation regularization are used. Now, we show how a further increase of image quality can be achieved for imaging dynamic processes in living tissue (4D PAT). The key idea is to exploit the additional temporal redundancy of the data by coupling the previously used spatial image reconstruction models with sparsity-constrained motion estimation models. While simulated data from a two-dimensional numerical phantom will be used to illustrate the main properties of this recently developed joint-image-reconstruction-and-motion-estimation framework, measured data from a dynamic experimental phantom will also be used to demonstrate their potential for challenging, large-scale, real-world, three-dimensional scenarios. The latter only becomes feasible if a carefully designed combination of tailored optimization schemes is employed, which we describe and examine in more detail.

preprint2018arXiv

Model based learning for accelerated, limited-view 3D photoacoustic tomography

Recent advances in deep learning for tomographic reconstructions have shown great potential to create accurate and high quality images with a considerable speed-up. In this work we present a deep neural network that is specifically designed to provide high resolution 3D images from restricted photoacoustic measurements. The network is designed to represent an iterative scheme and incorporates gradient information of the data fit to compensate for limited view artefacts. Due to the high complexity of the photoacoustic forward operator, we separate training and computation of the gradient information. A suitable prior for the desired image structures is learned as part of the training. The resulting network is trained and tested on a set of segmented vessels from lung CT scans and then applied to in-vivo photoacoustic measurement data.

preprint2016arXiv

Accelerated High-Resolution Photoacoustic Tomography via Compressed Sensing

Current 3D photoacoustic tomography (PAT) systems offer either high image quality or high frame rates but are not able to deliver high spatial and temporal resolution simultaneously, which limits their ability to image dynamic processes in living tissue. A particular example is the planar Fabry-Perot (FP) scanner, which yields high-resolution images but takes several minutes to sequentially map the photoacoustic field on the sensor plane, point-by-point. However, as the spatio-temporal complexity of many absorbing tissue structures is rather low, the data recorded in such a conventional, regularly sampled fashion is often highly redundant. We demonstrate that combining variational image reconstruction methods using spatial sparsity constraints with the development of novel PAT acquisition systems capable of sub-sampling the acoustic wave field can dramatically increase the acquisition speed while maintaining a good spatial resolution: First, we describe and model two general spatial sub-sampling schemes. Then, we discuss how to implement them using the FP scanner and demonstrate the potential of these novel compressed sensing PAT devices through simulated data from a realistic numerical phantom and through measured data from a dynamic experimental phantom as well as from in-vivo experiments. Our results show that images with good spatial resolution and contrast can be obtained from highly sub-sampled PAT data if variational image reconstruction methods that describe the tissues structures with suitable sparsity-constraints are used. In particular, we examine the use of total variation regularization enhanced by Bregman iterations. These novel reconstruction strategies offer new opportunities to dramatically increase the acquisition speed of PAT scanners that employ point-by-point sequential scanning as well as reducing the channel count of parallelized schemes that use detector arrays.

preprint2016arXiv

Certain families of Polynomials arising in the study of hyperelliptic Lie algebras

The associative ring $R(P(t))=\mathbb C[t^{\pm1},u \,|\, u^2=P(t)]$, where $P(t)=\sum_{i=0}^na_it^i=\prod_{k=1}^n(t-α_i)$ with $α_i\in\mathbb C$ pairwise distinct, is the coordinate ring of a hyperelliptic curve. The Lie algebra $\mathcal{R}(P(t))=\text{Der}(R(P(t)))$ of derivations is called the hyperelliptic Lie algebra associated to $P(t)$. In this paper we describe the universal central extension of $\text{Der}(R(P(t)))$ in terms of certain families of polynomials which in a particular case are associated Legendre polynomials. Moreover we describe certain families of polynomials that arise in the study of the group of units for the ring $R(P(t))$ where $P(t)=t^4-2bt^2+1$. In this study pairs of Chebychev polynomials $(U_n,T_n)$ arise as particular cases of a pairs $(r_n,s_n)$ with $r_n+s_n\sqrt{P(t)}$ a unit in $R(P(t))$. We explicitly describe these polynomial pairs as coefficients of certain generating functions and show certain of these polynomials satisfy particular second order linear differential equations.

preprint2016arXiv

Quantitative photoacoustic tomography using forward and adjoint Monte Carlo models of radiance

Forward and adjoint Monte Carlo (MC) models of radiance are proposed for use in model-based quantitative photoacoustic tomography. A 2D radiance MC model using a harmonic angular basis is introduced and validated against analytic solutions for the radiance in heterogeneous media. A gradient-based optimisation scheme is then used to recover 2D absorption and scattering coefficients distributions from simulated photoacoustic measurements. It is shown that the functional gradients, which are a challenge to compute efficiently using MC models, can be calculated directly from the coefficients of the harmonic angular basis used in the forward and adjoint models. This work establishes a framework for transport-based quantitative photoacoustic tomography that can fully exploit emerging highly parallel computing architectures.

preprint2015arXiv

Imaginary Verma Modules for $U_q(\widehat{\mathfrak{sl}(2)})$ and Crystal-like bases

We consider imaginary Verma modules for quantum affine algebraU_q(\widehat{\mathfrak{sl}(2)}) and define a crystal-like base which we call an imaginary crystal basis using the Kashiwara algebra K_q constructed in earlier work of the authors. In particular, we prove the existence of imaginary like bases for a suitable category of reduced imaginary Verma modules for U_q(\widehat{\mathfrak{sl}(2)}).

preprint2013arXiv

$N$-point locality for vertex operators: normal ordered products, operator product expansions, twisted vertex algebras

In this paper we study fields satisfying $N$-point locality and their properties. We obtain residue formulae for $N$-point local fields in terms of derivatives of delta functions and Bell polynomials. We introduce the notion of the space of descendants of $N$-point local fields which includes normal ordered products and coefficients of operator product expansions. We show that examples of $N$-point local fields include the vertex operators generating the boson-fermion correspondences of type B, C and D-A. We apply the normal ordered products of these vertex operators to the setting of the representation theory of the double-infinite rank Lie algebras $b_{\infty}, c_{\infty}, d_{\infty}$. Finally, we show that the field theory generated by $N$-point local fields and their descendants has a structure of a twisted vertex algebra.

preprint2013arXiv

$N$-point Virasoro Algebras and Their Modules of Densities

In this paper we introduce and study $n$-point Virasoro algebras, $\tilde{\W_a}$, which are natural generalizations of the classical Virasoro algebra and have as quotients multipoint genus zero Krichever-Novikov type algebras. We determine necessary and sufficient conditions for the latter two such Lie algebras to be isomorphic. Moreover we determine their automorphisms, their derivation algebras, their universal central extensions, and some other properties. The list of automorphism groups that occur is $C_n$, $D_n$, $A_4$, $S_4$ and $A_5$. We also construct a large class of modules which we call modules of densities, and determine necessary and sufficient conditions for them to be irreducible.

preprint2013arXiv

Representations of $a_{\infty}$ and $d_{\infty}$ with central charge 1 on the single neutral fermion Fock space $\mathit{F^{\otimes \frac{1}{2}}}$

We construct a new representation of the infinite rank Lie algebra $a_{\infty}$ with central charge $c=1$ on the Fock space $\mathit{F^{\otimes \frac{1}{2}}}$ of a single neutral fermion. We show that $\mathit{F^{\otimes \frac{1}{2}}}$ is a direct sum of irreducible integrable highest weight modules for $a_{\infty}$ with central charge $c=1$. We prove that as $a_{\infty}$ modules $\mathit{F^{\otimes \frac{1}{2}}}$ is isomorphic to the Fock space $\mathit{F^{\otimes 1}}$ of the charged free fermions. As a corollary we obtain the decompositions of certain irreducible highest weight modules for $d_{\infty}$ with central charge $c=\frac{1}{2}$ into irreducible highest weight modules for $d_{\infty}$ with central charge $c=1$.