Source author record

Fatma Terzioglu

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

ResearcherUnclaimed source record

math.NA math.FA Numerical Analysis physics.data-an physics.med-ph

Catalog footprint

What is connected

5works

5topics

3close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

An inversion algorithm for P-functions with applications to Multi-energy CT

Multi-energy computed tomography (ME-CT) is an x-ray transmission imaging technique that uses the energy dependence of x-ray photon attenuation to determine the elemental composition of an object of interest. Mathematically, forward ME-CT measurements are modeled by a nonlinear integral transform. In this paper, local conditions for global invertibility of the ME-CT transform are studied, and explicit stability estimates quantifying the error propagation from measurements to reconstructions are provided. Motivated from the inverse problem of image reconstruction in ME-CT, an iterative inversion algorithm for the so-called P-functions is proposed. Numerical simulations for ME-CT, in two and three materials settings with an equal number of energy measurements, confirm the theoretical predictions.

preprint2020arXiv

Uniqueness criteria in multi-energy CT

Multi-Energy Computed Tomography (ME-CT) is a medical imaging modality aiming to reconstruct the spatial density of materials from the attenuation properties of probing x-rays. For each line in two- or three-dimensional space, ME-CT measurements may be written as a nonlinear mapping from the integrals of the unknown densities of a finite number of materials along said line to an equal or larger number of energy-weighted integrals corresponding to different x-ray source energy spectra. ME-CT reconstructions may thus be decomposed as a two-step process: (i) reconstruct line integrals of the material densities from the available energy measurements; and (ii) reconstruct densities from their line integrals. Step (ii) is the standard linear x-ray CT problem whose invertibility is well-known, so this paper focuses on step (i). We show that ME-CT admits stable, global inversion provided that (a well-chosen linear transform of) the differential of the transform in step (i) satisfies appropriate orientation constraints that makes it a P-matrix. We introduce a notion of quantitative P-function that allows us to derive global stability results for ME-CT in the determined as well as over-determined (with more source energy spectra than the number of materials) cases. Numerical simulations based on standard material properties in imaging applications (of bone, water, contrast agents) and well accepted models of source energy spectra show that ME-CT is often (always in our simulations) either (i) non-globally injective because it is non-injective locally (differential not of full rank), or (ii) globally injective as soon as it is locally injective (differentials satisfy our proposed constraints).

preprint2016arXiv

3D Image Reconstruction from Compton camera data

In this paper, we address analytically and numerically the inversion of the integral transform (\emph{cone} or \emph{Compton} transform) that maps a function on $\mathbb{R}^3$ to its integrals over conical surfaces. It arises in a variety of imaging techniques, e.g. in astronomy, optical imaging, and homeland security imaging, especially when the so called Compton cameras are involved. Several inversion formulas are developed and implemented numerically in $3D$ (the much simpler $2D$ case was considered in a previous publication). An admissibility condition on detectors geometry is formulated, under which all these inversion techniques will work.

preprint2016arXiv

Inversion of Weighted Divergent Beam and Cone Transforms

In this paper, we investigate the relations between the Radon and weighted divergent beam and cone transforms. Novel inversion formulas are derived for the latter two. The weighted cone transform arises, for instance, in image reconstruction from the data obtained by Compton cameras, which have promising applications in various fields, including biomedical and homeland security imaging and gamma ray astronomy. The inversion formulas are applicable for a wide variety of detector geometries in any dimension. The results of numerical implementation of some of the formulas in dimensions two and three are also provided.

preprint2015arXiv

Some Inversion Formulas for the Cone Transform

Several novel imaging applications have lead recently to a variety of Radon type transforms, where integration is done over a family of conical surfaces. We call them \emph{cone transforms} (in 2D they are also called \emph{V-line} or \emph{broken ray} transforms). Most prominently, they are present in the so called Compton camera imaging that arises in medical diagnostics, astronomy, and lately in homeland security applications. Several specific incarnations of the cone transform have been considered separately. In this paper, we address the most general (and overdetermined) cone transform, obtain integral relations between cone and Radon transforms in $\mathbb{R}^n$, and a variety of inversion formulas. In many applications (e.g., in homeland security), the signal to noise ratio is very low. So, if overdetermined data is collected (as in the case of Compton imaging), attempts to reduce the dimensionality might lead to essential elimination of the signal. Thus, our main concentration is on obtaining formulas involving overdetermined data.