-
Instance-level quantitative saliency in multiple sclerosis lesion segmentation
Authors:
Federico Spagnolo,
Nataliia Molchanova,
Roger Schaer,
Meritxell Bach Cuadra,
Mario Ocampo Pineda,
Lester Melie-Garcia,
Cristina Granziera,
Vincent Andrearczyk,
Adrien Depeursinge
Abstract:
In recent years, explainable methods for artificial intelligence (XAI) have tried to reveal and describe models' decision mechanisms in the case of classification tasks. However, XAI for semantic segmentation and in particular for single instances has been little studied to date. Understanding the process underlying automatic segmentation of single instances is crucial to reveal what information w…
▽ More
In recent years, explainable methods for artificial intelligence (XAI) have tried to reveal and describe models' decision mechanisms in the case of classification tasks. However, XAI for semantic segmentation and in particular for single instances has been little studied to date. Understanding the process underlying automatic segmentation of single instances is crucial to reveal what information was used to detect and segment a given object of interest. In this study, we proposed two instance-level explanation maps for semantic segmentation based on SmoothGrad and Grad-CAM++ methods. Then, we investigated their relevance for the detection and segmentation of white matter lesions (WML), a magnetic resonance imaging (MRI) biomarker in multiple sclerosis (MS). 687 patients diagnosed with MS for a total of 4043 FLAIR and MPRAGE MRI scans were collected at the University Hospital of Basel, Switzerland. Data were randomly split into training, validation and test sets to train a 3D U-Net for MS lesion segmentation. We observed 3050 true positive (TP), 1818 false positive (FP), and 789 false negative (FN) cases. We generated instance-level explanation maps for semantic segmentation, by developing two XAI methods based on SmoothGrad and Grad-CAM++. We investigated: 1) the distribution of gradients in saliency maps with respect to both input MRI sequences; 2) the model's response in the case of synthetic lesions; 3) the amount of perilesional tissue needed by the model to segment a lesion. Saliency maps (based on SmoothGrad) in FLAIR showed positive values inside a lesion and negative in its neighborhood. Peak values of saliency maps generated for these four groups of volumes presented distributions that differ significantly from one another, suggesting a quantitative nature of the proposed saliency. Contextual information of 7mm around the lesion border was required for their segmentation.
△ Less
Submitted 25 June, 2024; v1 submitted 13 June, 2024;
originally announced June 2024.
-
Assessing radiomics feature stability with simulated CT acquisitions
Authors:
Kyriakos Flouris,
Oscar Jimenez-del-Toro,
Christoph Aberle,
Michael Bach,
Roger Schaer,
Markus Obmann,
Bram Stieltjes,
Henning Mueller,
Adrien Depeursinge,
Ender Konukoglu
Abstract:
Medical imaging quantitative features had once disputable usefulness in clinical studies. Nowadays, advancements in analysis techniques, for instance through machine learning, have enabled quantitative features to be progressively useful in diagnosis and research. Tissue characterisation is improved via the 'radiomics' features, whose extraction can be automated. Despite the advances, stability of…
▽ More
Medical imaging quantitative features had once disputable usefulness in clinical studies. Nowadays, advancements in analysis techniques, for instance through machine learning, have enabled quantitative features to be progressively useful in diagnosis and research. Tissue characterisation is improved via the 'radiomics' features, whose extraction can be automated. Despite the advances, stability of quantitative features remains an important open problem. As features can be highly sensitive to variations of acquisition details, it is not trivial to quantify stability and efficiently select stable features. In this work, we develop and validate a Computed Tomography (CT) simulator environment based on the publicly available ASTRA toolbox (www.astra-toolbox.com). We show that the variability, stability and discriminative power of the radiomics features extracted from the virtual phantom images generated by the simulator are similar to those observed in a tandem phantom study. Additionally, we show that the variability is matched between a multi-center phantom study and simulated results. Consequently, we demonstrate that the simulator can be utilised to assess radiomics features' stability and discriminative power.
△ Less
Submitted 6 October, 2022;
originally announced October 2022.
-
The Sizes and Albedos of Centaurs 2014 YY $_{49}$ and 2013 NL $_{24}$ from Stellar Occultation Measurements by RECON
Authors:
Ryder H. Strauss,
Rodrigo Leiva,
John M. Keller,
Elizabeth Wilde,
Marc W. Buie,
Robert J. Weryk,
JJ Kavelaars,
Terry Bridges,
Lawrence H. Wasserman,
David E. Trilling,
Deanna Ainsworth,
Seth Anthony,
Robert Baker,
Jerry Bardecker,
James K Bean Jr.,
Stephen Bock,
Stefani Chase,
Bryan Dean,
Chessa Frei,
Tony George,
Harnoorat Gill,
H. Wm. Gimple,
Rima Givot,
Samuel E. Hopfe,
Juan M. Cota Jr.
, et al. (24 additional authors not shown)
Abstract:
In 2019, the Research and Education Collaborative Occultation Network (RECON) obtained multiple-chord occultation measurements of two centaur objects: 2014 YY$_{49}$ on 2019 January 28 and 2013 NL$_{24}$ on 2019 September 4. RECON is a citizen-science telescope network designed to observe high-uncertainty occultations by outer solar system objects. Adopting circular models for the object profiles,…
▽ More
In 2019, the Research and Education Collaborative Occultation Network (RECON) obtained multiple-chord occultation measurements of two centaur objects: 2014 YY$_{49}$ on 2019 January 28 and 2013 NL$_{24}$ on 2019 September 4. RECON is a citizen-science telescope network designed to observe high-uncertainty occultations by outer solar system objects. Adopting circular models for the object profiles, we derive a radius $r=16^{+2}_{-1}$km and a geometric albedo $p_V=0.13^{+0.015}_{-0.024}$ for 2014 YY$_{49}$, and a radius $r=66 ^{+5}_{-5}$km and geometric albedo $p_V = 0.045^{+0.006}_{-0.008}$ for 2013 NL$_{24}$. To the precision of these measurements, no atmosphere or rings are detected for either object. The two objects measured here are among the smallest distant objects measured with the stellar occultation technique. In addition to these geometric constraints, the occultation measurements provide astrometric constraints for these two centaurs at a higher precision than has been feasible by direct imaging. To supplement the occultation results, we also present an analysis of color photometry from the Pan-STARRS surveys to constrain the rotational light curve amplitudes and spectral colors of these two centaurs. We recommend that future work focus on photometry to more deliberately constrain the objects' colors and light curve amplitudes, and on follow-on occultation efforts informed by this astrometry.
△ Less
Submitted 5 February, 2021;
originally announced February 2021.
-
A Single-chord Stellar Occultation by the Extreme Trans-Neptunian Object (541132) Leleākūhonua
Authors:
Marc W. Buie,
Rodrigo Leiva,
John M. Keller,
Josselin Desmars,
Bruno Sicardy,
JJ Kavelaars,
Terry Bridges,
Robert Weryk,
Dave Herald,
Sean L. Haley,
Ryder Strauss,
Elizabeth Wilde,
Robert Baker,
Ken Conway,
Bryan Dean,
Mackenzie Dunham,
James J Estes,
Naemi Fiechter,
Rima Givot,
Cameron Glibbery,
Bruce Gowe,
Jennifer N. Hayman,
Olivia L Ireland,
Matthew Kehrli,
Erik M. Moore
, et al. (10 additional authors not shown)
Abstract:
A stellar occultation by the extreme large-perihelion trans-Neptunian object (541132) Leleākūhonua (also known by the provisional designation of 2015 TG387) was predicted by the Lucky Star project and observed with the Research and Education Collaborative Occultation Network on 2018 October 20 UT. A single detection and a nearby nondetection provide constraints for the size and albedo. When a circ…
▽ More
A stellar occultation by the extreme large-perihelion trans-Neptunian object (541132) Leleākūhonua (also known by the provisional designation of 2015 TG387) was predicted by the Lucky Star project and observed with the Research and Education Collaborative Occultation Network on 2018 October 20 UT. A single detection and a nearby nondetection provide constraints for the size and albedo. When a circular profile is assumed, the radius is $r={110}_{-10}^{+14}$ km, corresponding to a geometric albedo ${p}_{V}={0.21}_{-0.05}^{+0.03}$, for an adopted absolute magnitude of H V = 5.6, typical of other objects in dynamically similar orbits. The occultation also provides a high-precision astrometric constraint.
△ Less
Submitted 7 November, 2020;
originally announced November 2020.
-
Standardised convolutional filtering for radiomics
Authors:
Adrien Depeursinge,
Vincent Andrearczyk,
Philip Whybra,
Joost van Griethuysen,
Henning Müller,
Roger Schaer,
Martin Vallières,
Alex Zwanenburg
Abstract:
The Image Biomarker Standardisation Initiative (IBSI) aims to improve reproducibility of radiomics studies by standardising the computational process of extracting image biomarkers (features) from images. We have previously established reference values for 169 commonly used features, created a standard radiomics image processing scheme, and developed reporting guidelines for radiomic studies. Howe…
▽ More
The Image Biomarker Standardisation Initiative (IBSI) aims to improve reproducibility of radiomics studies by standardising the computational process of extracting image biomarkers (features) from images. We have previously established reference values for 169 commonly used features, created a standard radiomics image processing scheme, and developed reporting guidelines for radiomic studies. However, several aspects are not standardised. Here we present a complete version of a reference manual on the use of convolutional filters in radiomics and quantitative image analysis. Filters, such as wavelets or Laplacian of Gaussian filters, play an important part in emphasising specific image characteristics such as edges and blobs. Features derived from filter response maps were found to be poorly reproducible. This reference manual provides definitions for convolutional filters, parameters that should be reported, reference feature values, and tests to verify software compliance with the reference standard.
△ Less
Submitted 7 February, 2024; v1 submitted 9 June, 2020;
originally announced June 2020.
-
The Parallel Distributed Image Search Engine (ParaDISE)
Authors:
Dimitrios Markonis,
Roger Schaer,
Alba García Seco de Herrera,
Henning Müller
Abstract:
Image retrieval is a complex task that differs according to the context and the user requirements in any specific field, for example in a medical environment. Search by text is often not possible or optimal and retrieval by the visual content does not always succeed in modelling high-level concepts that a user is looking for. Modern image retrieval techniques consist of multiple steps and aim to r…
▽ More
Image retrieval is a complex task that differs according to the context and the user requirements in any specific field, for example in a medical environment. Search by text is often not possible or optimal and retrieval by the visual content does not always succeed in modelling high-level concepts that a user is looking for. Modern image retrieval techniques consist of multiple steps and aim to retrieve information from large--scale datasets and not only based on global image appearance but local features and if possible in a connection between visual features and text or semantics. This paper presents the Parallel Distributed Image Search Engine (ParaDISE), an image retrieval system that combines visual search with text--based retrieval and that is available as open source and free of charge. The main design concepts of ParaDISE are flexibility, expandability, scalability and interoperability. These concepts constitute the system, able to be used both in real-world applications and as an image retrieval research platform. Apart from the architecture and the implementation of the system, two use cases are described, an application of ParaDISE in retrieval of images from the medical literature and a visual feature evaluation for medical image retrieval. Future steps include the creation of an open source community that will contribute and expand this platform based on the existing parts.
△ Less
Submitted 19 January, 2017;
originally announced January 2017.
-
Using MapReduce for Large-scale Medical Image Analysis
Authors:
Dimitrios Markonis,
Roger Schaer,
Ivan Eggel,
Henning Müller,
Adrien Depeursinge
Abstract:
The growth of the amount of medical image data produced on a daily basis in modern hospitals forces the adaptation of traditional medical image analysis and indexing approaches towards scalable solutions. The number of images and their dimensionality increased dramatically during the past 20 years. We propose solutions for large-scale medical image analysis based on parallel computing and algorith…
▽ More
The growth of the amount of medical image data produced on a daily basis in modern hospitals forces the adaptation of traditional medical image analysis and indexing approaches towards scalable solutions. The number of images and their dimensionality increased dramatically during the past 20 years. We propose solutions for large-scale medical image analysis based on parallel computing and algorithm optimization. The MapReduce framework is used to speed up and make possible three large-scale medical image processing use-cases: (i) parameter optimization for lung texture segmentation using support vector machines, (ii) content-based medical image indexing, and (iii) three-dimensional directional wavelet analysis for solid texture classification. A cluster of heterogeneous computing nodes was set up in our institution using Hadoop allowing for a maximum of 42 concurrent map tasks. The majority of the machines used are desktop computers that are also used for regular office work. The cluster showed to be minimally invasive and stable. The runtimes of each of the three use-case have been significantly reduced when compared to a sequential execution. Hadoop provides an easy-to-employ framework for data analysis tasks that scales well for many tasks but requires optimization for specific tasks.
△ Less
Submitted 23 October, 2015;
originally announced October 2015.
-
An Efficient and Optimal Filter for Identifying Point Sources in Millimeter/Sub-Millimeter Wavelength Sky Maps
Authors:
T. A. Perera,
G. W. Wilson,
K. S. Scott,
J. E. Austermann,
J. R. Schaar,
A. Mancera
Abstract:
A new technique for reliably identifying point sources in millimeter/sub-millimeter wavelength maps is presented. This method accounts for the frequency dependence of noise in the Fourier domain as well as non-uniformities in the coverage of a field. This optimal filter is an improvement over commonly-used matched filters that ignore coverage gradients. Treating noise variations in the Fourier dom…
▽ More
A new technique for reliably identifying point sources in millimeter/sub-millimeter wavelength maps is presented. This method accounts for the frequency dependence of noise in the Fourier domain as well as non-uniformities in the coverage of a field. This optimal filter is an improvement over commonly-used matched filters that ignore coverage gradients. Treating noise variations in the Fourier domain as well as map space is traditionally viewed as a computationally intensive problem. We show that the penalty incurred in terms of computing time is quite small due to casting many of the calculations in terms of FFTs and exploiting the absence of sharp features in the noise spectra of observations. Practical aspects of implementing the optimal filter are presented in the context of data from the AzTEC bolometer camera. The advantages of using the new filter over the standard matched filter are also addressed in terms of a typical AzTEC map.
△ Less
Submitted 2 April, 2013; v1 submitted 1 April, 2013;
originally announced April 2013.