PhD thesis · Université Claude Bernard Lyon 1 · Defended September 17, 2024

Automatic Detection of Pathologies in Conventional Scanner Data

Détection automatique d’anomalies sur examens tomodensitométriques standards

Aissam Djahnine

Philips Research France · CREATIS (INSA Lyon, Université Claude Bernard Lyon 1) · Hospices Civils de Lyon

3
levels of supervision, from segmentation to none
1st
place, JFR 2022 data challenge (PE)
0.952
AUROC, unsupervised 3D flow
88.8%
recall of weakly supervised localization
CT-3DFlow pipeline: a 3D normalizing flow trained on normal patches scores overlapping patches of a CT scan into a log-likelihood map that is aggregated and thresholded into a normal or abnormal prediction
The final contribution of the thesis: a 3D normalizing flow trained only on normal lung patches flags any deviation from normal anatomy, without annotation.

Abstract

Computed tomography (CT) is a common medical imaging examination used to diagnose or monitor many conditions, particularly in the thoracic region. These CT scans produce large amounts of high-quality volumetric data. However, the increasing number of scans and images produced by acquisitions result in a significant workload for radiologists. Radiologists primarily focus on suspected anomalies based on the initial clinical indication, which may lead to underdiagnosis of other atypical or unrelated anomalies. This underscores the importance of developing automatic alert systems based on AI to assist in identifying and locating incidental pathologies. Our thesis work focused on developing this type of system, applied to the detection and characterization of pulmonary pathologies in CT scans, using machine learning methods that require minimal annotations and are independent of clinical indications.

In this thesis, we explore different supervision methods to gradually reduce the dependency on the number of annotated data needed to train anomaly detection algorithms in CT scans. Through supervised learning, we proposed a method for detecting pulmonary embolisms in scans based on the segmentation of blood clots responsible for the embolism. We compared this segmentation-based supervised approach to a classification approach supervised only by the indication of the presence of a pulmonary embolism, illustrating the impact of the type of supervision on performance. We also proposed methods for automatically estimating clinical biomarkers such as the Qanadli score and the VD/VG ratio.

Although this supervised approach works very well and is often preferred for developing clinical decision support systems specific to a pathology, an anomaly detection system must support all pathologies and anatomical conditions encountered, and we do not currently have a representative annotated database. Therefore, we opted for a semiological approach that identifies specific patterns and morphologies present in most pulmonary pathologies (emphysema, cysts, nodules, consolidation, and interstitial patterns) and proposed a weakly supervised approach based on contrastive self-supervised pre-training. This approach allows detecting and locating a wide range of anomalies without requiring exhaustive and detailed data annotation.

Finally, we proposed an unsupervised approach using a generative method. Generative methods allow learning and representing complex distributions ; we learned the distribution of "normal" patients to then identify anomalies by deviation from this distribution. Training this model only requires a set of examinations from patients without anomalies, thus overcoming the difficulties in creating a representative database of all observed pathologies.

Keywords: Anomaly detection · Computed Tomography · Supervised learning · Weakly supervised learning · Unsupervised learning

Objectives

The thesis focuses on chest CT: a widely used modality with many examinations and highly reproducible volumetric data. The goal is a reading assistant that detects and localizes thoracic anomalies on any 3D CT scan with little annotation, by progressively reducing the level of supervision, from supervised to unsupervised methods.

Table 1: Three levels of supervision explored in the thesis, from most to least annotation.

PartSupervisionTarget anomaliesAnnotation
ISupervisedPulmonary embolismPixel level (lesion segmentation)
IISelf- and weakly supervised5 common semiological signs: emphysema and cysts, nodules, consolidation, interstitial pattern, normalPatch level (image regions)
IIIUnsupervisedAny deviation from normal lung anatomyNone (normal scans only)

I. Supervised: pulmonary embolism

Pulmonary embolism (PE) is a frequent and serious condition in which a clot blocks the pulmonary arteries. Diagnosis relies on CT angiography, the RV/LV diameter ratio and the Qanadli obstruction score. We compared two levels of supervision: a 3D segmentation network trained on fine clot segmentations from 54 PE-positive patients, and a 2D CNN + RNN classifier trained on image-level labels from more than 7,000 patients.

  • Segmentation beats labels: learning from a few precise clot segmentations outperformed classification trained on thousands of image-level labels.
  • Severity, not just presence: the same pipeline estimates the Qanadli score and the RV/LV ratio automatically.
  • First place at the 2022 JFR data challenge (Journées Francophones de Radiologie) organized by the French Society of Radiology.

Full write-up: Detection and severity quantification of pulmonary embolism.

Supervised learning works very well for a specific pathology, but an anomaly detection system must cover every pathology and anatomical condition, and no representative annotated database exists. The next two parts reduce the need for annotation.

II. Self- and weakly supervised: semiological signs

Instead of one detector per disease, a semiological approach targets the patterns shared by most lung diseases: emphysema and cysts, nodules, consolidation and interstitial patterns. A 3D patch classifier is pre-trained with SimCLR contrastive learning on unlabeled patches from 1,337 patients, then fine-tuned on 4,000 visually selected patches from 1,330 patients. At inference, overlapping patches produce per-sign probability maps that are aggregated into a patient-level prediction.

Two-stage method: SimCLR pre-training of a patch classifier on unlabeled patches, supervised fine-tuning on labeled patches, then per-label probability maps aggregated into a patient prediction
Figure 1: Detection method. Stage 1: self-supervised SimCLR pre-training, then supervised fine-tuning of the patch classifier. Stage 2: per-label probability maps on overlapping 50³ patches, aggregated at patient level.
Latent space projections colored by sign (normal, emphysema, consolidation, nodules, interstitial) for random initialization, ImageNet initialization and SimCLR pre-training with fine-tuning
Figure 2: Structure of the latent space for random initialization (left), ImageNet weights (right) and SimCLR pre-training with fine-tuning (bottom). Pre-training separates the five classes much more clearly.
Axial CT slices with probability maps for emphysema, consolidation, nodules and interstitial pattern, and a 3D view, for four patients
Figure 3: Localization. CT slices with the prediction maps of the patch classifier for each sign (EMPH, COND, NOD, INTER) and a 3D view; warmer colors mean higher probability.
  • Self-supervision helps: SimCLR pre-training gave a clear boost in patient-level classification over random and ImageNet initialization, on a 251-patient test set.
  • Simultaneous localization: the maps localize several abnormal regions at once, consistent with the ground-truth annotations.
  • 88.8% recall and 88.3% precision of localization in a qualitative analysis by an expert radiologist.

Full write-ups: Weakly-supervised pathology detection and localization and Tailored 3D CT contrastive pretraining.

III. Unsupervised: normalizing flows

Normalizing flows are generative networks that learn the exact probability distribution of the data by transforming it into a simpler one. Trained on masked lung images from normal patients only, the model assigns low likelihood to anything that departs from normal anatomy. Thresholds and filters are selected on a validation set, and slice- or patch-level scores are aggregated into a patient-level decision.

Three variants: a 2D flow on axial slices, a 2.5D flow with adjacent slices, and a multi-plane flow combining axial, coronal and sagittal flows into a 3D binary map
Figure 4: 2D, 2.5D and multi-plane methods. The 2.5D variant adds spatial context from adjacent slices; the multi-plane variant runs one flow per plane and aggregates them at patient level.
Grid of real and generated lung CT images in axial, coronal, sagittal and 2.5D views
Figure 5: Real (left three columns) and generated (right three columns) CT images in (a) axial, (b) coronal and (c) sagittal planes, and (d) the 2.5D axial representation. The flow preserves vessels, bronchi and lung contours.
Images generated by the flow in axial, coronal and sagittal planes, conditioned on slice position from 0.09 to 0.99
Figure 6: Conditioned generation. The flow is conditioned on the slice position in the volume (0 at the top of the patient, 1 at the bottom), or on the X, Y, Z coordinates of a 3D patch.

The 3D variant, CT-3DFlow, works on 3D patches to catch pathologies that are hard to see in 2D or 2.5D and to localize them in 3D.

Table 2: Patient-level performance on a test set of 251 patients with various combinations of anomalies (defense slides). Best in bold; best unsupervised F1 and accuracy underlined.

MethodFamilyAUROCF1Accuracy
AEReconstruction0.6080.5380.556
VAEReconstruction0.6470.6070.598
One-Class SVMOne-class0.8720.8530.867
DRAEMReconstruction0.6510.6210.628
f-AnoGANGAN0.8210.7830.765
GANomalyGAN0.7720.740.694
FastFlowNormalizing flow0.9420.9120.903
CFLOW-ADNormalizing flow0.7320.7390.706
DifferNetNormalizing flow0.9010.8610.843
AnoDDPMDiffusion0.9240.8910.873
SimCLR + fine-tuning (Part II)Self-supervised0.9630.9740.952
NF 2DUnsupervised (ours)0.9330.9090.891
NF 2D, conditionedUnsupervised (ours)0.9710.9310.909
NF 2.5DUnsupervised (ours)0.9560.9030.916
NF multi-planeUnsupervised (ours)0.9180.90.895
NF 3D (CT-3DFlow)Unsupervised (ours)0.9520.9390.924
NF 3D, conditionedUnsupervised (ours)0.9270.9020.9
Axial CT slices next to anomaly probability maps for interstitial pattern, emphysema with mass, mass, and consolidation
Figure 7: Unsupervised localization: (a) interstitial pattern, (b) emphysema and mass, (c) mass, (d) consolidation. Warmer colors mean higher anomaly probability.
  • Beats other unsupervised families: our flows outperform one-class SVM, autoencoders, GANs, diffusion models and other normalizing flows.
  • 3D over 2D: 3D models perform better than the 2D versions.
  • Close to self-supervision: without any annotation, the flows are competitive with the SimCLR model of Part II.
  • Localization without supervision: the likelihood maps highlight the abnormal regions.

Full write-up: CT-3DFlow.

Conclusion & perspectives

  • Less annotation: exploring several levels of supervision reduced the dependence on annotations while keeping good performance.
  • Self-supervision and semiology: a promising route to generic reading-assistance systems.
  • Unsupervised methods: an interesting alternative for anomaly detection, outperforming existing state-of-the-art methods.
  • Current limits: performance does not yet allow routine clinical use, but opens promising perspectives.

Next steps include multimodal methods that combine imaging modalities (MRI, ultrasound, PET), text reports and biomarkers, and spectral CT, which adds information on tissue composition by separating materials by their absorption at different energies.

Jury

  • Pr. Isabelle Bloch, Sorbonne Université · Reviewer
  • Pr. Alexis Jacquier, Aix-Marseille Université · Reviewer
  • Pr. Adrian Basarab, Université Claude Bernard Lyon 1 · Examiner
  • Dr. Constance de Margerie-Mellon, Université Paris Cité · Examiner
  • Dr. Olivier Nempont, Philips Research France · Co-supervisor
  • Dr. Alexandre Popoff, Philips Research France · Co-supervisor
  • Pr. Loïc Boussel, UCBL, Hospices Civils de Lyon · Thesis director

Publications from the thesis

  1. Detection and severity quantification of pulmonary embolism with 3D CT data using an automated deep learning-based artificial solution, Diagnostic and Interventional Imaging, 2024.
  2. Weakly-supervised learning-based pathology detection and localization in 3D chest CT scans, Medical Physics, 2024.
  3. BERT-based natural language processing analysis of French CT reports: application to the measurement of the positivity rate for pulmonary embolism, Research in Diagnostic and Interventional Imaging, 2023.
  4. eCAPRI: automated TAVI risk assessment with imaging biomarkers, European Radiology, 2025.
  5. CT-3DFlow: Leveraging 3D Normalizing Flows for Unsupervised Detection of Pathological Pulmonary CT scans, arXiv, 2024.
  6. Tailored 3D CT contrastive pretraining to improve pulmonary pathology classification, IEEE ICSP, 2022.

BibTeX

@phdthesis{djahnine2024thesis,
  title  = {D{\'e}tection automatique d'anomalies sur examens
            tomodensitom{\'e}triques standards},
  author = {Djahnine, Aissam},
  school = {Universit{\'e} Claude Bernard Lyon 1},
  year   = {2024},
  number = {2024LYO10167},
  url    = {https://theses.hal.science/tel-05263492}
}

More research