Giacomo Caviglia

The linear sampling method and energy conservation

In this paper we explain the linear sampling method and its performances under various scattering conditions by means of an analysis of the far-field equation based on the principle of energy conservation. Specifically, we consider the conservation of energy along the flow strips of the Poynting vector associated with the scattered field whose far-field pattern is one of the two terms in the far-field equation. The behavior of these flow lines is numerically investigated and theoretically described. Appropriate assumptions on the flow lines, based on the numerical results, allow characterizing a set of approximate solutions of the far-field equation which can be used to visualize the boundary of the scatterer in the framework of the linear sampling method. In particular, under the same assumptions, we can show that Tikhonov regularized solutions belong to this set of approximate solutions for appropriate choices of the regularization parameter.

An optimisation approach to multiprobe cryosurgery planning

In cryosurgery operations, tumoural cells are killed by means of a freezing procedure realised with the insertion of cryoprobes in the diseased tissue. Cryosurgery planning aims at establishing the best values for operation parameters like number and position of the probes or temperature and duration of the freezing process. Here, we present an application of ant colony optimisation (ACO) to cryosurgery planning, whereby the ACO cost function is computed by numerically solving several direct Stefan problems in biological tissues. The method is validated in the case of a 2D phantom of a prostate cross section.

A novel description of FDG excretion in the renal system: application to metformin-treated models

This paper introduces a novel compartmental model describing the excretion of 18F-fluoro-deoxyglucose (FDG) in the renal system and a numerical method based on the maximum likelihood for its reduction. This approach accounts for variations in FDG concentration due to water re-absorption in renal tubules and the increase of the bladders volume during the FDG excretion process. From the computational viewpoint, the reconstruction of the tracer kinetic parameters is obtained by solving the maximum likelihood problem iteratively, using a non-stationary, steepest descent approach that explicitly accounts for the Poisson nature of nuclear medicine data. The reliability of the method is validated against two sets of synthetic data realized according to realistic conditions. Finally we applied this model to describe FDG excretion in the case of animal models treated with metformin. In particular we show that our approach allows the quantitative estimation of the reduction of FDG de-phosphorylation induced by metformin.

Estimate of FDG Excretion by means of Compartmental Analysis and Ant Colony Optimization of Nuclear Medicine Data

[18F]fluoro-2-deoxy-D-glucose (FDG) is one of the most utilized tracers for positron emission tomography (PET) applications in oncology. FDG-PET relies on higher glycolytic activity in tumors compared to normal structures as the basis of image contrast. As a glucose analog, FDG is transported into malignant cells which typically exhibit an increased radioactivity. However, different from glucose, FDG is not reabsorbed by the renal system and is excreted to the bladder. The present paper describes a novel computational method for the quantitative assessment of this excretion process. The method is based on a compartmental analysis of FDG-PET data in which the excretion process is explicitly accounted for by the bladder compartment and on the application of an ant colony optimization (ACO) algorithm for the determination of the tracer coefficients describing the FDG transport effectiveness. The validation of this approach is performed by means of both synthetic data and real measurements acquired by a PET device for small animals (micro-PET). Possible oncological applications of the results are discussed in the final section.

A new compartmental method for the analysis of liver FDG kinetics in small animal models

BackgroundCompartmental analysis is a standard method to quantify metabolic processes using fluorodeoxyglucose-positron emission tomography (FDG-PET). For liver studies, this analysis is complex due to the hepatocyte capability to dephosphorylate and release glucose and FDG into the blood. Moreover, a tracer is supplied to the liver by both the hepatic artery and the portal vein, which is not visible in PET images. This study developed an innovative computational approach accounting for the reversible nature of FDG in the liver and directly computing the portal vein tracer concentration by means of gut radioactivity measurements.MethodsTwenty-one mice were subdivided into three groups: the control group ‘CTR’ (n = 7) received no treatment, the short-term starvation group ‘STS’ (n = 7) was submitted to food deprivation with free access to water within 48 h before imaging, and the metformin group ‘MTF’ (n = 7) was treated with metformin (750 mg/Kg per day) for 1 month. All mice underwent a dynamic micro-PET study for 50 min after an 18F-FDG injection. The compartmental analysis considered two FDG pools (phosphorylated and free) in both the gut and liver. A tracer was carried into the liver by the hepatic artery and the portal vein, and tracer delivery from the gut was considered as the sole input for portal vein tracer concentration. Accordingly, both the liver and gut were characterized by two compartments and two exchange coefficients. Each one of the two two-compartment models was mathematically described by a system of differential equations, and data optimization was performed by applying a Newton algorithm to the inverse problems associated to these differential systems.ResultsAll rate constants were stable in each group. The tracer coefficient from the free to the metabolized compartment in the liver was increased by STS, while it was unaltered by MTF. By contrast, the tracer coefficient from the metabolized to the free compartment was reduced by MTF and increased by STS.ConclusionsData demonstrated that our method was able to analyze FDG kinetics under pharmacological or pathophysiological stimulation, quantifying the fraction of the tracer trapped in the liver or dephosphorylated and released into the bloodstream.

A physiology-based parametric imaging method for FDG–PET data

Parametric imaging is a compartmental approach that processes nuclear imaging data to estimate the spatial distribution of the kinetic parameters governing tracer flow. The present paper proposes a novel and efficient computational method for parametric imaging which is potentially applicable to several compartmental models of diverse complexity and which is effective in the determination of the parametric maps of all kinetic coefficients. We consider applications to [{18}F]-fluorodeoxyglucose Positron Emission Tomography (FDG-PET) data and analyze the two-compartment catenary model describing the standard FDG metabolization by an homogeneous tissue and the three-compartment non-catenary model representing the renal physiology. We show uniqueness theorems for both models. The proposed imaging method starts from the reconstructed FDG-PET images of tracer concentration and preliminarily applies image processing algorithms for noise reduction and image segmentation. The optimization procedure solves pixelwise the non-linear inverse problem of determining the kinetic parameters from dynamic concentration data through a regularized Gauss-Newton iterative algorithm. The reliability of the method is validated against synthetic data, for the two-compartment system, and experimental real data of murine models, for the renal three-compartment system.

The linear sampling method in a lossy background: an energy perspective

In this article we propose a physical approach to the linear sampling method (LSM) for a possibly inhomogeneous and lossy background, whereby the modified far-field equation at the basis of the method can be regarded as a constraint on the power fluxes carried by the Poynting vector associated with the scattered field. Under appropriate assumptions on the flow lines of this Poynting vector, the general theorem inspiring the LSM (and concerning the existence of approximate solutions to the modified far-field equation) can be reformulated in a different way, which is more appropriate to explain the performance of the method.

The Role of Point Sources and Their Power Fluxes in the Linear Sampling Method

In this paper we investigate the linear sampling method by focusing on energy conservation inside a lossless background, in the case of a three-dimensional, impenetrable, and acoustic scattering set-up. We analyze, from both a numerical and a theoretical viewpoint, how average power fluxes are carried throughout the host medium by the flow tubes of radiating fields. As a result, the far-field equation, which is the core of the linear sampling method, can be regarded as a physical constraint linking the power flux of the scattered wave with that of the field radiated by a point source. Then, we show that this constraint, together with appropriate assumptions on the flow tubes of the scattered field, gives rise to a physical framework whereby some theoretical flaws of the linear sampling method can be overcome.

Parametric Imaging of FDG-PET Data Using Physiology and Iterative Regularization: Application to the Hepatic and Renal Systems

The present paper proposes a novel computational method for parametric imaging of nuclear medicine data. The mathematical procedure is general enough to work for compartmental models of diverse complexity and is effective in the determination of the parametric maps of all kinetic parameters governing tracer flow. We consider applications to [18F]-fluorodeoxyglucose Positron Emission Tomography (FDG-PET) data and analyze the two-compartment catenary model describing the standard FDG metabolization by an homogeneous tissue, e.g. the liver, and the three-compartment non-catenary model representing the renal physiology. The proposed imaging method starts from the reconstructed FDG-PET images of tracer concentration and preliminarily applies image processing algorithms for noise reduction and image segmentation processes for selecting the region enclosing the organ of physiologic interest. The optimization scheme solves pixelwise the non-linear inverse problem of determining the kinetic parameters from dynamic concentration data through a Gauss-Newton iterative algorithm with a penalty term accounting for the ill-posedness of the problem. We tested our imaging approach on FDG-PET data of murine models obtained by means of a dedicated microPET system, and we analyzed different PET slices containing axial sections of the liver and axial sections of the kidneys. The reconstructed parametric images proved to be reliable and qualitatively effective in the description of the local FDG metabolism with respect to the different physiologies.

Energy streamlines for qualitative inverse scattering in fluids and solids

In time-harmonic acoustic fields, energy streamlines are defined as the integral curves of the power-flux density vector, averaged over a period. They provide a tool to visualize the details of propagation of energy. After reviewing the role of energy streamlines in the linear sampling method for acoustic inverse scattering, this work formulates a physical interpretation of the same qualitative method in the case of an isotropic homogeneous solid matrix. Specifically, it is shown that the linear sampling method results from conservation of energy along streamline tubes of energy flow associated with elastic waves.