Nested Radially Monotone Polar Occupancy Estimation: Clinically-Grounded Optic Disc and Cup Segmentation for Glaucoma Screening
Organizations: Graduate School of Advanced Integrated Studies in Human Survivability (GSAIS), Kyoto University, Kyoto 606-8306, Japan
Abstract
Valid segmentation of the optic disc (OD) and optic cup (OC) from fundus photographs is essential for glaucoma screening. Unfortunately, existing deep learning methods do not guarantee clinical validness including star-convexity and nested structure of OD and OC, resulting corruption in diagnostic metric, especially under cross-dataset domain shift. To adress this issue, this paper proposed NPS-Net (Nested Polar Shape Network), the first framework that formulates the OD/OC segmentation as nested radially monotone polar occupancy estimation.This output representation can guarantee the aforementioned clinical validness and achieve high accuracy. Evaluated across seven public datasets, NPS-Net shows strong zero-shot generalization. On RIM-ONE, it maintains 100% anatomical validity and improves Cup Dice by 12.8% absolute over the best baseline, reducing vCDR MAE by over 56%. On PAPILA, it achieves Disc Dice of 0.9438 and Disc HD95 of 2.78 px, an 83% reduction over the best competing method.
Figures & tables
| Stage | Channels | Norm | Convolution | Params |
|---|---|---|---|---|
| Enc 1 | GN (8) | CircPad | 31.04 M | |
| Enc 2 | GN (8) | CircPad | ||
| Enc 3 | BN | CircPad | ||
| Enc 4 | BN | CircPad | ||
| Bottleneck | BN | CircPad | ||
| Shape Prior | — | BN1d | Circ1D | 0.24 M |
| Loss | Supervises | Epoch | |
|---|---|---|---|
| Dice + BCE | Cartesian masks | 1.0 | 0 |
| Dice + BCE | Polar masks | 0.7 | 0 |
| SmoothL1 | Rim profile ( 16 ) | 0.5 | 0 |
| Cross-entropy | Shape distributions | 0.3 | 20 |
| SmoothL1 | Shape radii | 0.5 | 20 |
| SmoothL1 | Contour smoothness | 0.05 | 20 |
| Dataset | Ethnicity | Hardware | Country | Images |
| Internal Train / Validation / Test (70/15/15%) | ||||
| RIGA–BinRushed | Arab | Topcon TRC-NW7SF | Saudi Arabia | 233 |
| RIGA–Messidor | Caucasian | Topcon TRC-NW6 | France | 463 |
| ORIGA | Malay | Canon CR-DGi | Singapore | 654 |
| NETRA [ 36 ] | Nepalese | Canon | Nepal | 480 |
| Cross-Domain Benchmark (DRISHTI-GS + REFUGE Combined) | ||||
| Dice | HD95 | vCDR | Rim | Rim | |||
|---|---|---|---|---|---|---|---|
| Cup | Disc | Cup | Disc | ||||
| UNet | .8929 | .9643 | 4.65 | 0.90 | .0480 | .0258 | .6309 |
| Attn-UNet | .8946 | .9649 | 4.19 | 0.71 | .0481 | .0254 | .6450 |
| ResUNet | .8941 | .9650 | 4.70 | 0.86 | .0520 | .0253 | .6408 |
| Polar UNet | .8830 | .9719 | 5.08 | 0.76 | .0527 | .0268 | .6520 |
| TransUNet | .8852 | .9710 | 4.69 | 0.28 | .0534 | .0255 | .6194 |
| Dice | HD95 | vCDR | Rim | Rim | |||
|---|---|---|---|---|---|---|---|
| Cup | Disc | Cup | Disc | ||||
| UNet | .5162 | .5733 | 19.01 | 33.00 | .1781 | .1046 | .0828 |
| Attn-UNet | .5035 | .5399 | 20.37 | 31.26 | .1913 | .1114 | .0380 |
| ResUNet | .4625 | .5230 | 21.33 | 38.88 | .2171 | .1114 | .1338 |
| Polar UNet | .6094 | .7285 | 19.30 | 23.50 | .1546 | .0783 | .1999 |
| TransUNet | .6789 | .8539 | 18.53 | 13.51 | .1189 | .0646 | .2076 |
| Cup Dice | Disc Dice | vCDR MAE | ||||
|---|---|---|---|---|---|---|
| UNet | +.29 | ∗∗∗ | +.34 | ∗∗∗ | .11 | ∗∗ |
| Attn-UNet | +.30 | ∗∗∗ | +.37 | ∗∗∗ | .12 | ∗∗∗ |
| ResUNet | +.34 | ∗∗∗ | +.39 | ∗∗∗ | .15 | ∗∗∗ |
| Polar UNet | +.20 | ∗∗∗ | +.19 | ∗∗∗ | .09 | ∗∗∗ |
| TransUNet | +.13 | ∗∗ | +.06 | ∗∗∗ | .05 | ∗∗∗ |
| Dice | HD95 | vCDR | Rim | Rim | |||
|---|---|---|---|---|---|---|---|
| Cup | Disc | Cup | Disc | ||||
| UNet | .3786 | .4199 | 28.84 | 55.86 | .2404 | .2322 | .043 |
| Attn-UNet | .4023 | .4232 | 27.01 | 54.88 | .2694 | .2349 | .074 |
| ResUNet | .3222 | .3640 | 29.40 | 62.79 | .2689 | .2347 | .030 |
| Polar UNet | .4937 | .7316 | 33.69 | 24.93 | .2924 | .1847 | .085 |
| TransUNet | .6153 | .8185 | 26.71 | 18.15 | .2260 | .1441 | .307 |
| Dice | HD95 | vCDR | Rim | Rim | |||
|---|---|---|---|---|---|---|---|
| Cup | Disc | Cup | Disc | ||||
| UNet | .8603 | .9065 | 7.65 | 6.30 | .0858 | .0642 | .6130 |
| Attn-UNet | .8512 | .8997 | 8.21 | 6.78 | .0744 | .0579 | .6108 |
| ResUNet | .8497 | .9132 | 8.73 | 5.08 | .0692 | .0535 | .6372 |
| Polar UNet | .8166 | .9349 | 11.66 | 2.69 | .0859 | .0625 | .6595 |
| TransUNet | .8565 | .8536 | 8.14 | 11.82 | .1158 | .0780 | .5717 |
| Model | Cup | Disc | HD95 | vCDR | Rim | Rim |
|---|---|---|---|---|---|---|
| M1: Polar UNet | .609 | .729 | 19.30 | .155 | .078 | .200 |
| M2: +Monotone | .760 | .929 | 13.64 | .084 | .045 | .425 |
| M3: +Nesting † | .781 | .935 | 11.48 | .079 | .040 | .431 |
| M4: +Shape prior | .783 | .930 | 11.29 | .071 | .037 | .519 |
| M5: +Polar-TTA | .806 | .914 | 10.25 | .068 | .035 | .555 |
| † Violation rate = 0% from M3 onward ( everywhere). | ||||||