Source author record

Philipp Roser

Philipp Roser 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

16works
7topics
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

16 published item(s)

preprint2021arXiv

X-ray Scatter Estimation Using Deep Splines

Algorithmic X-ray scatter compensation is a desirable technique in flat-panel X-ray imaging and cone-beam computed tomography. State-of-the-art U-net based image translation approaches yielded promising results. As there are no physics constraints applied to the output of the U-Net, it cannot be ruled out that it yields spurious results. Unfortunately, those may be misleading in the context of medical imaging. To overcome this problem, we propose to embed B-splines as a known operator into neural networks. This inherently limits their predictions to well-behaved and smooth functions. In a study using synthetic head and thorax data as well as real thorax phantom data, we found that our approach performed on par with U-net when comparing both algorithms based on quantitative performance metrics. However, our approach not only reduces runtime and parameter complexity, but we also found it much more robust to unseen noise levels. While the U-net responded with visible artifacts, our approach preserved the X-ray signal's frequency characteristics.

preprint2020arXiv

Appearance Learning for Image-based Motion Estimation in Tomography

In tomographic imaging, anatomical structures are reconstructed by applying a pseudo-inverse forward model to acquired signals. Geometric information within this process is usually depending on the system setting only, i. e., the scanner position or readout direction. Patient motion therefore corrupts the geometry alignment in the reconstruction process resulting in motion artifacts. We propose an appearance learning approach recognizing the structures of rigid motion independently from the scanned object. To this end, we train a siamese triplet network to predict the reprojection error (RPE) for the complete acquisition as well as an approximate distribution of the RPE along the single views from the reconstructed volume in a multi-task learning approach. The RPE measures the motioninduced geometric deviations independent of the object based on virtual marker positions, which are available during training. We train our network using 27 patients and deploy a 21-4-2 split for training, validation and testing. In average, we achieve a residual mean RPE of 0.013mm with an inter-patient standard deviation of 0.022 mm. This is twice the accuracy compared to previously published results. In a motion estimation benchmark the proposed approach achieves superior results in comparison with two state-of-the-art measures in nine out of twelve experiments. The clinical applicability of the proposed method is demonstrated on a motion-affected clinical dataset.

preprint2020arXiv

Deep Scatter Splines: Learning-Based Medical X-ray Scatter Estimation Using B-splines

The idea of replacing hardware by software to compensate for scattered radiation in flat-panel X-ray imaging is well established in the literature. Recently, deep-learningbased image translation approaches, most notably the U-Net, have emerged for scatter estimation. These yield considerable improvements over model-based methods. Such networks, however, involve potential drawbacks that need to be considered. First, they are trained in a data-driven fashion without making use of prior knowledge and X-ray physics. Second, due to their high parameter complexity, the validity of deep neural networks is difficult to assess. To circumvent these issues, we introduce here a surrogate function to model X-ray scatter distributions that can be expressed by few parameters. We could show empirically that cubic B-splines are well-suited to model X-ray scatter in the diagnostic energy regime. Based on these findings, we propose a lean convolutional encoder architecture that extracts local scatter characteristics from X-ray projection images. These characteristics are embedded into a global context using a constrained weighting matrix yielding spline coefficients that model the scatter distribution. In a first simulation study with 17 thorax data sets, we could show that our method and the U-Net-based state of the art reach about the same accuracy. However, we could achieve these comparable outcomes with orders of magnitude fewer parameters while ensuring that not high-frequency information gets manipulated.

preprint2020arXiv

Effects of Tissue Material Properties on X-Ray Image, Scatter and Patient Dose Determined using Monte Carlo Simulations

With increasing patient and staff X-ray radiation awareness, many efforts have been made to develop accurate patient dose estimation methods. To date, Monte Carlo (MC) simulations are considered golden standard to simulate the interaction of X-ray radiation with matter. However, sensitivity of MC simulation results to variations in the experimental or clinical setup of image guided interventional procedures are only limited studied. In particular, the impact of patient material compositions is poorly investigated. This is mainly due to the fact, that these methods are commonly validated in phantom studies utilizing a single anthropomorphic phantom. In this study, we therefore investigate the impact of patient material parameters mapping on the outcome of MC X-ray dose simulations. A computation phantom geometry is constructed and three different commonly used material composition mappings are applied. We used the MC toolkit Geant4 to simulate X-ray radiation in an interventional setup and compared the differences in dose deposition, scatter distributions and resulting X-ray images. The evaluation shows a discrepancy between different material composition mapping up to 20 % concerning directly irradiated organs. These results highlight the need for standardization of material composition mapping for MC simulations in a clinical setup.

preprint2020arXiv

Fully-automatic CT data preparation for interventional X-ray skin dose simulation

Recently, deep learning (DL) found its way to interventional X-ray skin dose estimation. While its performance was found to be acceptable, even more accurate results could be achieved if more data sets were available for training. One possibility is to turn to computed tomography (CT) data sets. Typically, computed tomography (CT) scans can be mapped to tissue labels and mass densities to obtain training data. However, care has to be taken to make sure that the different clinical settings are properly accounted for. First, the interventional environment is characterized by wide variety of table setups that are significantly different from the typical patient tables used in conventional CT. This cannot be ignored, since tables play a crucial role in sound skin dose estimation in an interventional setup, e. g., when the X-ray source is directly underneath a patient (posterior-anterior view). Second, due to interpolation errors, most CT scans do not facilitate a clean segmentation of the skin border. As a solution to these problems, we applied connected component labeling (CCL) and Canny edge detection to (a) robustly separate the patient from the table and (b) to identify the outermost skin layer. Our results show that these extensions enable fully-automatic, generalized pre-processing of CT scans for further simulation of both skin dose and corresponding X-ray projections.

preprint2020arXiv

Pitfalls in interventional X-ray organ dose assessment -- combined experimental and computational phantom study: application to prostatic artery embolization

With X-ray radiation protection and dose management constantly gaining interest in interventional radiology, novel procedures often undergo prospective dose studies using anthropomorphic phantoms to determine expected dose values. Due to inherent uncertainties, these dose values might not allow for patient-specific risk assessment. Therefore, first the aim of this study is to quantify the influence of these parameters on local X-ray dose to evaluate their relevance in the assessment of patient-specific organ doses. Second, this knowledge further enables validating a simulation approach, which allows employing physiological material models and patient-specific geometries. Phantom dosimetry experiments using MOSFET dosimeters were conducted reproducing imaging scenarios in prostatic arterial embolization (PAE). Associated organ-equivalent dose of prostate, bladder, colon and skin was determined. Dose deviation induced by possible small displacements of the patient was reproduced by moving the X-ray source. Dose deviation induced by material differences was investigated by analyzing two different commonly used phantoms. We reconstructed the experiments using Monte Carlo (MC) simulations, a reference male geometry, and different material properties. Overall, MC simulated organ dose values are in accordance with the measured ones for the majority of cases. Marginal displacements of X-ray source relative to the phantoms lead to deviations of 6% to 135% in organ dose values, while skin dose remains relatively constant. Regarding the impact of phantom material composition, underestimation of internal organ dose values by 12% to 20% is prevalent in all simulated phantoms. Skin dose, however, can be estimated with low deviation of 1% to 8% at least for two materials. Prospective reference dose studies might not extend to precise patient-specific dose assessment.

preprint2020arXiv

Simultaneous Estimation of X-ray Back-Scatter and Forward-Scatter using Multi-Task Learning

Scattered radiation is a major concern impacting X-ray image-guided procedures in two ways. First, back-scatter significantly contributes to patient (skin) dose during complicated interventions. Second, forward-scattered radiation reduces contrast in projection images and introduces artifacts in 3-D reconstructions. While conventionally employed anti-scatter grids improve image quality by blocking X-rays, the additional attenuation due to the anti-scatter grid at the detector needs to be compensated for by a higher patient entrance dose. This also increases the room dose affecting the staff caring for the patient. For skin dose quantification, back-scatter is usually accounted for by applying pre-determined scalar back-scatter factors or linear point spread functions to a primary kerma forward projection onto a patient surface point. However, as patients come in different shapes, the generalization of conventional methods is limited. Here, we propose a novel approach combining conventional techniques with learning-based methods to simultaneously estimate the forward-scatter reaching the detector as well as the back-scatter affecting the patient skin dose. Knowing the forward-scatter, we can correct X-ray projections, while a good estimate of the back-scatter component facilitates an improved skin dose assessment. To simultaneously estimate forward-scatter as well as back-scatter, we propose a multi-task approach for joint back- and forward-scatter estimation by combining X-ray physics with neural networks. We show that, in theory, highly accurate scatter estimation in both cases is possible. In addition, we identify research directions for our multi-task framework and learning-based scatter estimation in general.

preprint2020arXiv

Tenfold your photons -- a physically-sound approach to filtering-based variance reduction of Monte-Carlo-simulated dose distributions

X-ray dose constantly gains interest in the interventional suite. With dose being generally difficult to monitor reliably, fast computational methods are desirable. A major drawback of the gold standard based on Monte Carlo (MC) methods is its computational complexity. Besides common variance reduction techniques, filter approaches are often applied to achieve conclusive results within a fraction of time. Inspired by these methods, we propose a novel approach. We down-sample the target volume based on the fraction of mass, simulate the imaging situation, and then revert the down-sampling. To this end, the dose is weighted by the mass energy absorption, up-sampled, and distributed using a guided filter. Eventually, the weighting is inverted resulting in accurate high resolution dose distributions. The approach has the potential to considerably speed-up MC simulations since less photons and boundary checks are necessary. First experiments substantiate these assumptions. We achieve a median accuracy of 96.7 % to 97.4 % of the dose estimation with the proposed method and a down-sampling factor of 8 and 4, respectively. While maintaining a high accuracy, the proposed method provides for a tenfold speed-up. The overall findings suggest the conclusion that the proposed method has the potential to allow for further efficiency.

preprint2016arXiv

Gravitation and cosmology with York time

[Shortened abstract:] In this thesis we investigate a solution to the `problem of time' in canonical quantum gravity by splitting spacetime into surfaces of constant mean curvature parameterised by York time. We argue that there are reasons to consider York time a viable candidate for a physically meaningful notion of time. We investigate a number York-time Hamiltonian-reduced cosmological models and explore some technical aspects, such as the non-canonical Poisson structure. We develop York-time Hamiltonian-reduced cosmological perturbation theory by solving the Hamiltonian constraint perturbatively around a homogeneous background for the physical (non-vanishing) Hamiltonian that is the momentum conjugate to the York time parameter. We proceed to canonically quantise the cosmological models and the perturbation theory and discuss a number of conceptual and technical points, such as volume eigenfunctions and the absence of a momentum representation due to the non-standard commutator structure. We propose an alternative, wavefunction-free method of quantisation based on an ensemble of trajectories and a dynamical configuration-space geometry and discuss its application to gravity. We conclude by placing the York-time theories explored in this thesis in the wider context of a search for a satisfactory theory of quantum gravity.

preprint2015arXiv

An extension of cosmological dynamics with York time

It has been suggested that the York parameter $T$ (effectively the scalar extrinsic curvature of a spatial hypersurface) may play the role of a fundamental time parameter. In a flat, forever expanding cosmology the York parameter remains always negative, taking values $T=-\infty$ at the big bang and approaching some finite non-positive value as $t\rightarrow\infty$, $t$ being the usual cosmological time coordinate. Based on previous results concerning a simple, spatially flat cosmological model with a scalar field, we provide a temporal extension of this model to include `times' $T>0$, an epoch not covered by the cosmological time coordinate $t$, and discuss the dynamics of this `other side' and its significance. We argue that the extension is necessary if a consistent quantisation scheme is to exist. Furthermore, we investigate which types of potentials lead to smooth transitions, paying particular attention to currently favoured inflaton potentials.

preprint2015arXiv

Cosmological perturbation theory with York time

One method to overcome the notorious problem of time in the quantisation of gravity is the identification of a physically preferred time parameter, a promising candidate being so-called `York time'. The dynamical equations for matter and spatial geometry in York time may be obtained via Hamiltonian reduction, that is, by solving the Hamiltonian constraint for the physical, non-vanishing Hamiltonian density identified as the variable conjugate to the chosen time parameter. Yet in general this equation cannot be solved algebraically. Here we show how in a cosmological scenario, where one may treat geometric and matter inhomogeneities as small perturbations, one is able to obtain the physical Hamiltonian density by solving the constraint equation perturbatively. By construction the Hamiltonian density is quadratic in the perturbation variables, which makes it easily quantisable, although subtleties arise due to the non-canonical form of the Poisson brackets and the time-dependent coefficients. The latter are determined by the evolution of the background variables.

preprint2015arXiv

Quantum mechanics as the dynamical geometry of trajectories

We illustrate how non-relativistic quantum mechanics may be recovered from a dynamical Weyl geometry on configuration space and an `ensemble' of trajectories (or `worlds'). The theory, which is free of a physical wavefunction, is presented starting from a classical `many-systems' action to which a curvature term is added. In this manner the equations of equilibrium de~Broglie-Bohm theory are recovered. However, naïvely the set of solution precludes solutions with non-zero angular momentum (a version of a problem raised by Wallstrom). This is remedied by a slight extension of the action, leaving the equations of motion unchanged.

preprint2015arXiv

Quantum mechanics on York slices

For some time the York time parameter has been identified as a candidate for a physically meaningful time in cosmology. An associated Hamiltonian may be found by solving the Hamiltonian constraint for the momentum conjugate to the York time variable, although an explicit solution can only be found in highly symmetric cases. The Poisson structure of the remaining variables is not canonical. Here we quantise this dynamics in an anisotropic minisuperspace model via a natural extension of canonical quantisation. The resulting quantum theory has no momentum representation. Instead the position basis takes a fundamental role. We illustrate how the quantum theory and the modified representation of its momentum operators lead to a consistent theory in the presence of the constraints that arose during the Hamiltonian reduction. We are able to solve for the eigenspectrum of the Hamiltonian. Finally we discuss how far the results of this model extend to the general non-homogeneous case, in particular perturbation theory with York time.

preprint2014arXiv

Classical and quantum cosmology with York time

We consider a solution to the problem of time in quantum gravity by deparameterisation of the ADM action in terms of York time, a parameter proportional to the extrinsic curvature of a spatial hypersurface. We study a minisuperspace model together with a homogeneous scalar field, for which we can solve the Hamiltonian constraint exactly and arrive at an explicit expression for the physical (non-vanishing) Hamiltonian. The scale factor and associated momentum cease to be dynamical variables, leaving the scalar field as the only physical degree of freedom. We investigate the resulting classical theory, showing how the dynamics of the scale factor can be recovered via an appropriate interpretation of the Hamiltonian as a volume. We then quantise the system in the Schrödinger picture. In the quantum theory we recover the dynamics of the scale factor by interpreting the spectrum and expectation value of the Hamiltonian as being associated with volume rather than energy. If trajectories in the sense of de~Broglie-Bohm are introduced in the quantum theory, these are found to match those of the classical theory. We suggest that these trajectories may provide the basis for a perturbation theory in which both background and perturbations are quantised.

preprint2012arXiv

Quantum computation from a de Broglie-Bohm perspective

We consider an example of a quantum algorithm from the point of view of the de Broglie-Bohm formulation of quantum mechanics. For concreteness we look at two particular implementations: one using spin-1/2 particles as described by a simple model due to Bell and the other as energy states in an infinite-well potential. We extend the analysis to encompass a complete set of quantum gates. We conclude by discussing the relevance of our investigations for the debate about the origin of the quantum-computational speed-up.

preprint2010arXiv

Quantum Computing in the de Broglie-Bohm Pilot-Wave Picture

Much attention has been drawn to quantum computing and the exponential speed-up in computation the technology would be able to provide. Various claims have been made about what aspect of quantum mechanics causes this speed-up. Formulations of quantum computing have traditionally been made in orthodox (Copenhagen) and sometimes many-worlds quantum mechanics. We will aim to understand quantum computing in terms of de Broglie-Bohm Pilot-Wave Theory by considering different simple systems that may function as a basic quantum computer. We will provide a careful discussion of Pilot-Wave Theory and evaluate criticisms of the theory. We will assess claims regarding what causes the exponential speed-up in the light of our analysis and the fact that Pilot-Wave Theory is perfectly able to account for the phenomena involved in quantum computing.