Random Wavelet nnU-Net: Prostate Lesion Masks Across bpMRI Sites

Fig. 1 of arXiv:2611.08535 by Heintz and colleagues, Weill Cornell Medicine and Cornell Tech: T2, ADC and DWI of a PI-CAI prostate MRI case as acquired and after the random sym4 wavelet filter used to train the open nnU-Net lesion model.

If you have ever tried to train a prostate lesion model on more than one public dataset, you know the problem. Every cohort was scanned differently and labeled by different rules, and the lesions themselves are small and low contrast on biparametric MRI (bpMRI). In arXiv:2611.08535, Thibault Heintz and colleagues at Weill Cornell Medicine and Cornell Tech pool five of those cohorts into one nnU-Net, add a deliberately destructive augmentation based on random wavelet filters, and release both the trainer code and the trained checkpoint, so you can download the model and run it on your own studies.

What the model does

You give it three co-registered sequences of a prostate bpMRI study: T2-weighted, the ADC map and the high b-value DWI, with ADC and DWI resampled onto the T2 grid. It returns a 3D mask with a lesion label. During training it also learns a prostate gland label, a trick called BackSplit (co-developed by one of the authors) that gives the network extra signal about the organ around a small lesion. At inference only the lesion label is used.

For a viewer or a reporting tool, the mask stands in for the contour a radiologist draws by hand today before an MRI-ultrasound fusion biopsy or a structured report. The paper cites only moderate agreement between readers for those contours, so a model-drafted first contour has a clear place in that workflow.

Five cohorts, five labeling rules

The training pool has 1,579 lesion-containing cases, which the authors describe as the largest collection used for this task. It combines a derived PROMIS set, PI-CAI (with the extra annotations by Pooch et al.), Prostate158, the UCLA Prostate-MRI-US-Biopsy collection (PMUB) and an internal Weill Cornell cohort. The cohorts differ in three ways that matter for training:

  • Scanners. PROMIS is entirely 1.5T, PI-CAI mixes 1.5T and 3T, the rest are 3T. High b-values range from 1000 to 1400 s/mm².
  • What counts as a lesion. Prostate158 annotates PI-RADS 4 and above, the internal cohort generally PI-RADS 3 and above, PMUB only the biopsy targets, and PI-CAI lesions of ISUP grade 2 or higher.
  • Lesion size. Median lesion volume runs from 0.182 mL in the internal data to 1.121 mL in PROMIS, mostly a consequence of those labeling rules.

Destructive augmentation with random wavelets

With probability 0.5 per training sample, on top of the standard nnU-Net augmentations, each sequence is replaced by the response of a randomly built wavelet filter. The filter decomposes the volume at four scales, mixes the detail bands of the three axes along a random direction, sums the scales with random weights and rescales the result to the original intensity range. A new filter is drawn for each sequence.

The result looks nothing like an MRI. The featured image of this post shows Fig. 1 of the paper: the top row is a PI-CAI case as acquired, the bottom row is the same case after the sym4 filter. Textures are scrambled and edges are displaced. The authors argue that this helps because the labels are fuzzy and the cohorts differ, so the network gains little from fitting exact image-to-mask alignment. They compare four wavelet families (Haar, db2, db4, sym4) against random convolutions and nnU-Net’s stronger DA5 preset, and sym4 comes out best on average Dice. The gain is largest on PROMIS, the only cohort scanned entirely at 1.5T and a minority of the training data.

The second change is the loss. Plain Dice plus cross-entropy is computed over all lesion voxels at once, so one large lesion can dominate a case with several small ones. They add a connected-component loss (CC-Dice-CE) that splits the volume into Voronoi cells around each reference lesion and scores each lesion in its own cell, and they show it works together with BackSplit.

How well it works

On the official Prostate158 test split, against Reader 1 labels, the final model reaches a Dice of 0.614, compared with the 0.450 baseline reported by the dataset authors. Across all five test sets, the full recipe raises mean Dice from 0.391 for a plain nnU-Net ResEnc-M to 0.422.

Both the wavelet augmentation and the CC loss raise lesion-level recall and lower precision. The final model’s lesion F1 across all test sets is 0.571, below the 0.596 of the BackSplit-only model, though still above plain nnU-Net. The authors say they prefer recall because a missed lesion costs more than a false positive. In a viewer, that means more candidate blobs for the reader to dismiss.

Where it can fail

  • No negative cases. Cases without an annotated lesion were excluded from training and testing, so there is no evidence on how often the model marks something in a lesion-free gland. Run it on your own negatives before you show its output to anyone.
  • Lenient detection rule. A predicted lesion counts as found if it overlaps a reference lesion with Dice above 0.1. The detection numbers are therefore generous about boundary quality.
  • Small public test set. The Prostate158 test split behind the headline number has 19 cases.
  • Single fold, final checkpoint. Every configuration was trained once on one fold, so small differences between augmentations may sit inside run-to-run noise.
  • Input preparation is on you. The model expects ADC and DWI already resampled to the T2 grid and stored as NIfTI in nnU-Net’s channel order. When you start from DICOM series, that conversion and resampling step is where most integration bugs will hide.
  • Not tested out of domain. The authors list a true out-of-domain evaluation as future work.

Code, weights and license

The trainers live on GitHub at tppm/prostate_lesion_segmentation under the MIT license. They plug into nnU-Net v2 as extra trainer classes driven by small YAML configs, and the repo includes the evaluation script and the train, validation and test split lists for the public cohorts. The final checkpoint (ResEnc-M, 3d_fullres, fold 0, about 400 MB as safetensors) is on Hugging Face at heintz/prostate_lesion_seg, also marked MIT. A few rough edges as of 9 November 2026: the README still shows a placeholder for the Hugging Face repo name, the trainer imports a bicc_loss module that is not in the repo, and the README says the models were trained on an nnU-Net fork without linking the commit. Expect some plumbing before the first prediction.

The MIT tag covers the code and the weights as published. The training data come from public cohorts with their own terms plus an internal cohort that is not shared, so check those terms before building a commercial product on the checkpoint. The model is research software and is not cleared as a medical device.

Paper: arXiv:2611.08535 · Code: GitHub (MIT) · Weights: Hugging Face (MIT)

We build custom medical imaging platforms — advanced DICOM viewers, AI segmentation, and the clinical systems around them.

Get in Touch

Copyright © 2026 PYCAD. All Rights Reserved.