Regularized reconstruction in quantitative SPECT using CT side information from hybrid imaging

A penalized-likelihood (PL) SPECT reconstruction method using a modified regularizer that accounts for anatomical boundary side information was implemented to achieve accurate estimates of both the total target activity and the activity distribution within targets. In both simulations and experimental I-131 phantom studies, reconstructions from (1) penalized likelihood employing CT-side information-based regularization (PL-CT), (2) penalized likelihood with edge preserving regularization (no CT) and (3) penalized likelihood with conventional spatially invariant quadratic regularization (no CT) were compared with (4) ordered subset expectation maximization (OSEM), which is the iterative algorithm conventionally used in clinics for quantitative SPECT. Evaluations included phantom studies with perfect and imperfect side information and studies with uniform and non-uniform activity distributions in the target. For targets with uniform activity, the PL-CT images and profiles were closest to the 'truth', avoided the edge offshoots evident with OSEM and minimized the blurring across boundaries evident with regularization without CT information. Apart from visual comparison, reconstruction accuracy was evaluated using the bias and standard deviation (STD) of the total target activity estimate and the root mean square error (RMSE) of the activity distribution within the target. PL-CT reconstruction reduced both bias and RMSE compared with regularization without side information. When compared with unregularized OSEM, PL-CT reduced RMSE and STD while bias was comparable. For targets with non-uniform activity, these improvements with PL-CT were observed only when the change in activity was matched by a change in the anatomical image and the corresponding inner boundary was also used to control the regularization. In summary, the present work demonstrates the potential of using CT side information to obtain improved estimates of the activity distribution in targets without sacrificing the accuracy of total target activity estimation. The method is best suited for data acquired on hybrid systems where SPECT-CT misregistration is minimized. To demonstrate clinical application, the PL reconstruction with CT-based regularization was applied to data from a patient who underwent SPECT/CT imaging for tumor dosimetry following I-131 radioimmunotherapy.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85409/1/pmb10_9_007.pd

Treatment of Advanced Emphysema with Emphysematous Lung Sealant (AeriSeal (R))

Background: This report summarizes initial tests of an emphysematous lung synthetic polymer sealant (ELS) designed to reduce lung volume in patients with advanced emphysema. Objectives: The primary study objective was to define a therapeutic strategy to optimize treatment safety and effectiveness. Methods: ELS therapy was administered bronchoscopically to 25 patients with heterogeneous emphysema in an open-label, noncontrolled study at 6 centers in Germany. Treatment was performed initially at 2-4 subsegments. After 12 weeks, patients were eligible for repeat therapy to a total of 6 sites. Safety and efficacy were assessed after 6 months. Responses were evaluated in terms of changes from baseline in lung physiology, functional capacity, and health-related quality of life. Follow-up is available for 21 of 25 patients. Results: Treatment was well tolerated. There were no treatment-related deaths (i.e. within 90 days of treatment), and an acceptable short-and long-term safety profile. Physiological and clinical benefits were observed at 24 weeks. Efficacy responses were better among Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage III patients {[}n = 14; change in residual volume/total lung capacity (Delta RV/TLC) = -7.4 +/- 10.3%; Delta forced expiratory volume in 1 s (Delta FEV(1)) = +15.9 +/- 22.6%; change in forced vital capacity (Delta FVC) = +24.1 +/- 22.7%; change in carbon monoxide lung diffusion capacity (Delta DLCO) = +19.3 +/- 34.8%; change in 6-min walk test (Delta 6MWD) = +28.7 +/- 59.6 m; change in Medical Research Council Dyspnea (Delta MRCD) score = -1.0 +/- 1.04 units; change in St. George's Respiratory Questionnaire (Delta SGRQ) score = -9.9 +/- 15.3 units] than for GOLD stage IV patients (n = 7; Delta RV/TLC = -0.5 +/- 6.4%; Delta FEV 1 = +2.3 +/- 12.3%; Delta FVC = +2.6 +/- 21.1%; Delta DLCO = -2.8 +/- 17.2%; Delta 6MWD = +28.3 +/- 58.4 m; Delta MRCD = 0.3 +/- 0.81 units; Delta SGRQ = -6.7 +/- 7.0 units). Conclusions: ELS therapy shows promise for treating patients with advanced heterogeneous emphysema. Additional studies to assess responses in a larger cohort with a longer follow-up are warranted. Copyright (C) 2011 S. Karger AG, Base

Acceleration of heavy and light particles in turbulence: comparison between experiments and direct numerical simulations

We compare experimental data and numerical simulations for the dynamics of inertial particles with finite density in turbulence. In the experiment, bubbles and solid particles are optically tracked in a turbulent flow of water using an Extended Laser Doppler Velocimetry technique. The probability density functions (PDF) of particle accelerations and their auto-correlation in time are computed. Numerical results are obtained from a direct numerical simulation in which a suspension of passive pointwise particles is tracked, with the same finite density and the same response time as in the experiment. We observe a good agreement for both the variance of acceleration and the autocorrelation timescale of the dynamics; small discrepancies on the shape of the acceleration PDF are observed. We discuss the effects induced by the finite size of the particles, not taken into account in the present numerical simulations.Comment: 7 pages, 4 figure

Large scale cosmic-ray anisotropy with KASCADE

The results of an analysis of the large scale anisotropy of cosmic rays in the PeV range are presented. The Rayleigh formalism is applied to the right ascension distribution of extensive air showers measured by the KASCADE experiment.The data set contains about 10^8 extensive air showers in the energy range from 0.7 to 6 PeV. No hints for anisotropy are visible in the right ascension distributions in this energy range. This accounts for all showers as well as for subsets containing showers induced by predominantly light respectively heavy primary particles. Upper flux limits for Rayleigh amplitudes are determined to be between 10^-3 at 0.7 PeV and 10^-2 at 6 PeV primary energy.Comment: accepted by The Astrophysical Journa

Mathematical Properties of a New Levin-Type Sequence Transformation Introduced by \v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la. I. Algebraic Theory

\v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la [J. Math. Phys. \textbf{44}, 962 - 968 (2003)] introduced in connection with the summation of the divergent perturbation expansion of the hydrogen atom in an external magnetic field a new sequence transformation which uses as input data not only the elements of a sequence $\{s_n \}_{n=0}^{\infty}$ of partial sums, but also explicit estimates $\{\omega_n \}_{n=0}^{\infty}$ for the truncation errors. The explicit incorporation of the information contained in the truncation error estimates makes this and related transformations potentially much more powerful than for instance Pad\'{e} approximants. Special cases of the new transformation are sequence transformations introduced by Levin [Int. J. Comput. Math. B \textbf{3}, 371 - 388 (1973)] and Weniger [Comput. Phys. Rep. \textbf{10}, 189 - 371 (1989), Sections 7 -9; Numer. Algor. \textbf{3}, 477 - 486 (1992)] and also a variant of Richardson extrapolation [Phil. Trans. Roy. Soc. London A \textbf{226}, 299 - 349 (1927)]. The algebraic theory of these transformations - explicit expressions, recurrence formulas, explicit expressions in the case of special remainder estimates, and asymptotic order estimates satisfied by rational approximants to power series - is formulated in terms of hitherto unknown mathematical properties of the new transformation introduced by \v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la. This leads to a considerable formal simplification and unification.Comment: 41 + ii pages, LaTeX2e, 0 figures. Submitted to Journal of Mathematical Physic

The ZEUS Forward Plug Calorimeter with Lead-Scintillator Plates and WLS Fiber Readout

A Forward Plug Calorimeter (FPC) for the ZEUS detector at HERA has been built as a shashlik lead-scintillator calorimeter with wave length shifter fiber readout. Before installation it was tested and calibrated using the X5 test beam facility of the SPS accelerator at CERN. Electron, muon and pion beams in the momentum range of 10 to 100 GeV/c were used. Results of these measurements are presented as well as a calibration monitoring system based on a $^{60}$Co source.Comment: 38 pages (Latex); 26 figures (ps

Decision-Directed Channel Estimation Implementation for Spectral Efficiency Improvement in Mobile MIMO-OFDM

Channel estimation algorithms and their implementations for mobile receivers are considered in this paper. The 3GPP long term evolution (LTE) based pilot structure is used as a benchmark in a multiple-input multiple-output (MIMO) orthogonal frequency division multiplexing (OFDM) receiver. The decision directed (DD) space alternating generalized expectation-maximization (SAGE) algorithm is used to improve the performance from that of the pilot symbol based least-squares (LS) channel estimator. The performance is improved with high user velocities, where the pilot symbol density is not sufficient. Minimum mean square error (MMSE) filtering is also used in estimating the channel in between pilot symbols. The pilot overhead can be reduced to a third of the LTE pilot overhead with DD channel estimation, obtaining a ten percent increase in data throughput. Complexity reduction and latency issues are considered in the architecture design. The pilot based LS, MMSE and the SAGE channel estimators are implemented with a high level synthesis tool, synthesized with the UMC 0.18 μm CMOS technology and the performance-complexity trade-offs are studied. The MMSE estimator improves the performance from the simple LS estimator with LTE pilot structure and has low power consumption. The SAGE estimator has high power consumption but can be used with reduced pilot density to increase the data rate.National Science FoundationTekesElektrobitRenesas Mobile EuropeAcademy of FinlandNokia Siemens NetworksXilin

Sparse Matrix-Based HPC Tomography

Tomographic imaging has benefited from advances in X-ray sources, detectors and optics to enable novel observations in science, engineering and medicine. These advances have come with a dramatic increase of input data in the form of faster frame rates, larger fields of view or higher resolution, so high performance solutions are currently widely used for analysis. Tomographic instruments can vary significantly from one to another, including the hardware employed for reconstruction: from single CPU workstations to large scale hybrid CPU/GPU supercomputers. Flexibility on the software interfaces and reconstruction engines are also highly valued to allow for easy development and prototyping. This paper presents a novel software framework for tomographic analysis that tackles all aforementioned requirements. The proposed solution capitalizes on the increased performance of sparse matrix-vector multiplication and exploits multi-CPU and GPU reconstruction over MPI. The solution is implemented in Python and relies on CuPy for fast GPU operators and CUDA kernel integration, and on SciPy for CPU sparse matrix computation. As opposed to previous tomography solutions that are tailor-made for specific use cases or hardware, the proposed software is designed to provide flexible, portable and high-performance operators that can be used for continuous integration at different production environments, but also for prototyping new experimental settings or for algorithmic development. The experimental results demonstrate how our implementation can even outperform state-of-the-art software packages used at advanced X-ray sources worldwide

Career changers and fast-track induction: teacher perspectives on their early professional development

In this article, the early professional development of mature, early career teachers who entered the profession via an employment-based route to teaching in England is presented and explored from the teachers’ own perspectives. From a larger sample in a longitudinal study, the development of four career changers is traced in detail, using a model of professional learning which highlights teaching and professional skills. There is some evidence that those who undertake the programme experience a smoother induction into teaching than those from more traditional routes, although the internship year itself is challenging and demands a high level of commitment and resilience on the part of the pre-service teachers. Contextual factors, such as school support and learning from experienced teachers, are vital in enabling early professional development and the absence of these can hinder learning and self-confidence. Findings are analysed with reference to the model, as well as research on teachers’ development, within a socio-cultural learning framework. The study makes a valuable contribution to knowledge about the professional development of mature entrants to teaching

Spatial extent and isolation of marine artificial structures mediate fish density

Installations of artificial structures in coastal oceans create de facto habitat for marine life. These structures encompass wide varieties of physical characteristics, reflecting their multiple, diverse purposes and creating a need to understand which characteristics maximize fish habitat. Here, we test how physical characteristics – horizontal area, vertical relief, and spatial isolation – relate to fish density from echosounder surveys over artificial structures like concrete pipes, train boxcars, and ships purposely sunk to function as reefs. Echosounder mapping of 31 artificial reef structures and associated fish across a 200 km linear length of the continental shelf of North Carolina, USA, revealed that structures with greater horizontal area and vertical relief host higher fish densities than smaller, shorter structures. Artificial structure spatial arrangement also relates to fish density, as isolated structures are generally associated with greater localized fish densities than structures closer to one another. Patterns in the relationships between fish density and reef characteristics differed for schooling fish, as there was some evidence that reefs of intermediate area exhibited higher schooling fish density. These results suggest that intentional design and spatial arrangement of marine built structures like artificial reefs relates to and can be deliberately incorporated into siting and deployment decisions to enhance their role as fish habitat