Explainable anatomical shape analysis through deep hierarchical generative models.

Quantification of anatomical shape changes currently relies on scalar global indexes which are largely insensitive to regional or asymmetric modifications. Accurate assessment of pathology-driven anatomical remodeling is a crucial step for the diagnosis and treatment of many conditions. Deep learning approaches have recently achieved wide success in the analysis of medical images, but they lack interpretability in the feature extraction and decision processes. In this work, we propose a new interpretable deep learning model for shape analysis. In particular, we exploit deep generative networks to model a population of anatomical segmentations through a hierarchy of conditional latent variables. At the highest level of this hierarchy, a two-dimensional latent space is simultaneously optimised to discriminate distinct clinical conditions, enabling the direct visualisation of the classification space. Moreover, the anatomical variability encoded by this discriminative latent space can be visualised in the segmentation space thanks to the generative properties of the model, making the classification task transparent. This approach yielded high accuracy in the categorisation of healthy and remodelled left ventricles when tested on unseen segmentations from our own multi-centre dataset as well as in an external validation set, and on hippocampi from healthy controls and patients with Alzheimer's disease when tested on ADNI data. More importantly, it enabled the visualisation in three-dimensions of both global and regional anatomical features which better discriminate between the conditions under exam. The proposed approach scales effectively to large populations, facilitating highthroughput analysis of normal anatomy and pathology in largescale studies of volumetric imaging

Chromatic Illumination Discrimination Ability Reveals that Human Colour Constancy Is Optimised for Blue Daylight Illuminations

The phenomenon of colour constancy in human visual perception keeps surface colours constant, despite changes in their reflected light due to changing illumination. Although colour constancy has evolved under a constrained subset of illuminations, it is unknown whether its underlying mechanisms, thought to involve multiple components from retina to cortex, are optimised for particular environmental variations. Here we demonstrate a new method for investigating colour constancy using illumination matching in real scenes which, unlike previous methods using surface matching and simulated scenes, allows testing of multiple, real illuminations. We use real scenes consisting of solid familiar or unfamiliar objects against uniform or variegated backgrounds and compare discrimination performance for typical illuminations from the daylight chromaticity locus (approximately blue-yellow) and atypical spectra from an orthogonal locus (approximately red-green, at correlated colour temperature 6700 K), all produced in real time by a 10-channel LED illuminator. We find that discrimination of illumination changes is poorer along the daylight locus than the atypical locus, and is poorest particularly for bluer illumination changes, demonstrating conversely that surface colour constancy is best for blue daylight illuminations. Illumination discrimination is also enhanced, and therefore colour constancy diminished, for uniform backgrounds, irrespective of the object type. These results are not explained by statistical properties of the scene signal changes at the retinal level. We conclude that high-level mechanisms of colour constancy are biased for the blue daylight illuminations and variegated backgrounds to which the human visual system has typically been exposed

Natural images from the birthplace of the human eye

Here we introduce a database of calibrated natural images publicly available through an easy-to-use web interface. Using a Nikon D70 digital SLR camera, we acquired about 5000 six-megapixel images of Okavango Delta of Botswana, a tropical savanna habitat similar to where the human eye is thought to have evolved. Some sequences of images were captured unsystematically while following a baboon troop, while others were designed to vary a single parameter such as aperture, object distance, time of day or position on the horizon. Images are available in the raw RGB format and in grayscale. Images are also available in units relevant to the physiology of human cone photoreceptors, where pixel values represent the expected number of photoisomerizations per second for cones sensitive to long (L), medium (M) and short (S) wavelengths. This database is distributed under a Creative Commons Attribution-Noncommercial Unported license to facilitate research in computer vision, psychophysics of perception, and visual neuroscience.Comment: Submitted to PLoS ON

Charged-Lepton Flavour Physics

This writeup of a talk at the 2011 Lepton-Photon symposium in Mumbai, India, summarises recent results in the charged-lepton flavour sector. I review searches for charged-lepton flavour violation, lepton electric dipole moments and flavour-conserving CP violation. I also discuss recent progress in tau-lepton physics and in the Standard Model prediction of the muon anomalous magnetic moment.Comment: Presented at Lepton-Photon 2011, Mumbai, India; 23 pages, 14 figure

Turbulence and galactic structure

Interstellar turbulence is driven over a wide range of scales by processes including spiral arm instabilities and supernovae, and it affects the rate and morphology of star formation, energy dissipation, and angular momentum transfer in galaxy disks. Star formation is initiated on large scales by gravitational instabilities which control the overall rate through the long dynamical time corresponding to the average ISM density. Stars form at much higher densities than average, however, and at much faster rates locally, so the slow average rate arises because the fraction of the gas mass that forms stars at any one time is low, ~10^{-4}. This low fraction is determined by turbulence compression, and is apparently independent of specific cloud formation processes which all operate at lower densities. Turbulence compression also accounts for the formation of most stars in clusters, along with the cluster mass spectrum, and it gives a hierarchical distribution to the positions of these clusters and to star-forming regions in general. Turbulent motions appear to be very fast in irregular galaxies at high redshift, possibly having speeds equal to several tenths of the rotation speed in view of the morphology of chain galaxies and their face-on counterparts. The origin of this turbulence is not evident, but some of it could come from accretion onto the disk. Such high turbulence could help drive an early epoch of gas inflow through viscous torques in galaxies where spiral arms and bars are weak. Such evolution may lead to bulge or bar formation, or to bar re-formation if a previous bar dissolved. We show evidence that the bar fraction is about constant with redshift out to z~1, and model the formation and destruction rates of bars required to achieve this constancy.Comment: in: Penetrating Bars through Masks of Cosmic Dust: The Hubble Tuning Fork strikes a New Note, Eds., K. Freeman, D. Block, I. Puerari, R. Groess, Dordrecht: Kluwer, in press (presented at a conference in South Africa, June 7-12, 2004). 19 pgs, 5 figure

Single-Scale Natural SUSY

We consider the prospects for natural SUSY models consistent with current data. Recent constraints make the standard paradigm unnatural so we consider what could be a minimal extension consistent with what we now know. The most promising such scenarios extend the MSSM with new tree-level Higgs interactions that can lift its mass to at least 125 GeV and also allow for flavor-dependent soft terms so that the third generation squarks are lighter than current bounds on the first and second generation squarks. We argue that a common feature of almost all such models is the need for a new scale near 10 TeV, such as a scale of Higgsing or confinement of a new gauge group. We consider the question whether such a model can naturally derive from a single mass scale associated with supersymmetry breaking. Most such models simply postulate new scales, leaving their proximity to the scale of MSSM soft terms a mystery. This coincidence problem may be thought of as a mild tuning, analogous to the usual mu problem. We find that a single mass scale origin is challenging, but suggest that a more natural origin for such a new dynamical scale is the gravitino mass, m_{3/2}, in theories where the MSSM soft terms are a loop factor below m_{3/2}. As an example, we build a variant of the NMSSM where the singlet S is composite, and the strong dynamics leading to compositeness is triggered by masses of order m_{3/2} for some fields. Our focus is the Higgs sector, but our model is compatible with a light stop (with the other generation squarks heavy, or with R-parity violation or another mechanism to hide them from current searches). All the interesting low-energy mass scales, including linear terms for S playing a key role in EWSB, arise dynamically from the single scale m_{3/2}. However, numerical coefficients from RG effects and wavefunction factors in an extra dimension complicate the otherwise simple story.Comment: 32 pages, 3 figures; version accepted by JHE

Isomeric and β-decay spectroscopy of 173,174Ho

β-decay spectroscopy of 173,174Ho (Z =67, N =106,107) was conducted at Radioactive Isotope Beam Factory at RIKEN by using in-flight fission of a 345-MeV/u 238U primary beam. A previously unreported isomeric state at 405 keV with half-life of 3.7(12) μs and a spin and parity of (3/2+) is identified in 173 Ho. Moreover, a new state with a spin and parity of 9− was discovered in 174 Er. The experimental log ft values of 5.84(20) and 5.25(18) suggest an allowed-hindered β decay from the ground state of 174Ho to the Kπ=8−isomeric state in 174Er. Configuration-constrained potential energy surface (PES) calculations were performed and the predictions are in reasonable agreement with the experimental results.Part of the WAS3ABi was supported by the Rare Isotopes Science Project which is funded by MSIP and NRF of Korea. This work was supported by JSPS KAKENHI Grants No. 24740188, No. 25247045, and No. 25287065, STFC (UK authors), the UK National Measurement Office (P.H.R.), the U.S. Department of Energy, Office of Science, Office of Nuclear Physics under Contract No. DE-AC02-06CH11357 (F.G.K.), NRF Korea Grants No. 2016R1D1A1A09917463, No. 2017M2A2A6A02071071 (C.S.L.), No. 2019R1F1A1058370, and No. 2016R1A5A1013277 (K.Y.C.), and Science Foundation Ireland under Grant No. 12/IP/1288 (O.J.R.)

β-γ and isomeric decay spectroscopy of 168Dy

This contribution will report on the experimental work on the level structure of 168Dy. The experimental data have been taken as part of the EURICA decay spectroscopy campaign at RIBF, RIKEN in November 2014. In the experiment, a 238U primary beam is accelerated up to 345 MeV/u with an average intensity of 12 pnA. The nuclei of interest are produced by in-flight fission of 238U impinging on Be target with a thickness of 5 mm. The excited states of 168Dy have been populated through the decay from a newly identified isomeric state and via the β decay from 168Tb. In this contribution, scientific motivations, experimental procedure and some preliminary results for this study are presented