PreprintUnder double-blind reviewOAIZIB-CM · n=102

MBEDNet

A Mamba–CNN network with dual edge-guided decoders for segmentation of tibiofemoral bone and cartilage.

Anonymized authors · Anonymized affiliations
17.4 M parameters · 93.11% mean Dice · 10-page submission

Sagittal knee MRI slice.The same slice with MBEDNet's five-class segmentation painted over it.
Sagittal · 0.365 × 0.365 × 0.7 mm
Scroll
01 — Input

One sagittal slice

Volumetric knee MRI at 0.365 × 0.365 × 0.7 mm. The articular cartilage running between femur and tibia is two to three voxels thick and barely separable from the bone beneath it.

02 — Segmentation

Five classes, one pass

MBEDNet labels femoral and tibial bone plus femoral, medial and lateral tibial cartilage — a semantic decoder and a boundary decoder running together, so the bone–cartilage interface is predicted rather than inferred from region overlap.

03 — Reconstruction

Straight into surgical meshes

The predicted mask converts to triangular meshes. Eight femoral landmarks localised on those meshes land within 0.964 mm of the ground truth — the geometry a patient-specific implant is planned against.

Raw sagittal knee MRI slice.MBEDNet's predicted five-class mask over the same slice.
Reconstructed 3D mesh of the distal femur with cartilage.Reconstructed 3D mesh of the proximal tibia with cartilage.
Mean Dice · held-out test set
0.00%
Highest of the nine compared architectures, four of which are 1.3–5.5× larger.
Parameters
0.0M
23% below the nearest-scoring model.
HD95
0.00mm
Full model, all five classes.
Registration error
0.000mm
Tracked 3D-printed phantom.
Femoral boneFemoral cartilageTibial boneMedial tibial cartilageLateral tibial cartilage
00
Abstract

Preoperative 3D segmentation of the tibiofemoral joint guides image-based Total Knee Arthroplasty (TKA) and supports the design of patient-specific implants. Automated 3D segmentation in TKA still remains challenging due to thin, low contrast cartilage layers between the tibiofemoral bone structures.

We propose Mamba Edge-Net (MBEDNet), which integrates the local feature extraction of convolution with long dependency modelling of Mamba to effectively and efficiently segment tibiofemoral anatomies, with sharp contrast between bone and cartilage. MBEDNet is lightweight while achieving the highest average Dice similarity coefficient among compared state-of-the-art architectures, several of which require 1.3 to 5.5 times more parameters. The network is also optimized to handle severe osteophytes of Kellgren–Lawrence grade IV.

Evaluated on the publicly available OAIZIB-CM knee MRI dataset, MBEDNet reaches Dice scores of 98.46% and 98.64% on the femur and tibia, and 89.03% and 86.29% on femoral and tibial cartilage.

01
Background

Accurate TKA relies on precise segmentation of the femur, tibia and the articular cartilage around them — and manually tracing those joint lines across slices is laborious. Automating it runs into a split requirement: the cartilage–bone interface needs voxel-level acuity, while femur–tibia alignment needs the whole volume in view. No single backbone family gives both.

Where each family breaks down, as characterised in the paper
Backbone familyReceptive fieldComplexityWeak point
Convolutional U-NetsLocalLinearCannot model femur–tibia spatial relations
Vision TransformersGlobalO(N²) in tokensProhibitive on volumes; weak fine detail
Swin / windowed attentionWindowedReducedBuys efficiency by giving up long-range context
State space models (Mamba)GlobalLinearState averaging smooths out voxel-level structure

MBEDNet’s answer is not to pick one. Each encoder stage runs a convolutional branch and a state-space branch in parallel and fuses them, and a second decoder is dedicated entirely to boundaries — so the interface evidence the state-space path smooths away is reconstructed by a stream built to find it.

02 — Architecture

Two branches in, two decoders out

MBEDNet architecture: a Mamba-CNN encoder with parallel ResBlock3D+MPSK and TSMambaBlock branches fused by CBFFM, a two-block TSMamba bottleneck, and a dual decoder in which a PEE3D edge stream is injected into an ECA3D-gated mask stream at every level.
Whole networkEncoder, bottleneck, dual decoderA 7³ stride-2 stem, three parallel-branch encoder stages, a 256-channel state-space bottleneck, and two level-aligned decoders emitting a five-class mask and an edge map.
Left — encoderConvolution and state space, side by sideEach stage runs an MPSK branch and a tri-oriented TSMamba branch on the same input, then fuses them through a channel-wise bilateral gate. The fused tensor is both the skip and the next stage’s input.
Centre — bottleneckWhere femur meets tibiaTwo stacked TSMambaBlocks at 256 channels and 1/16 resolution model relations a local receptive field cannot reach — relative femur–tibia alignment above all.
Right — dual decoderEdges injected into the maskA PEE edge stream takes the raw, un-gated skip, and its boundary evidence is injected into the ECA-gated mask stream at every level — one direction only, edge to mask.
02
Components

Stem7³ · stride 2

A 7×7×7 convolution lifts the volume to 32 channels and halves every spatial dimension, giving (32, D/2, H/2, W/2) before the first encoder stage.

TSMambaglobal branch

The volume is refined three times, once along each axis. Each pass keeps selective state deciding what to retain or discard; the three scans merge through a 1×1×1 convolution and instance norm.

FTOM = IN( Conv1×1×1[ SD(X) SH(X) SW(X) ] )

Three separate passes keep the cost linear at O(D + H + W) each — against O((D × H × W)²) for full 3D self-attention.

MPSKlocal branch

Parallel 3D convolutions at 1³, 3³ and 5³ feed a learned per-channel selective-kernel mechanism, so the network uses small kernels for thin cartilage texture and large ones for bone. Pyramid pooling at 2×2, 4×4 and 8×8 adds multi-scale context.

FMPSK = Fin + FSK + Fpool
CBFFMfusion

The branches concatenate, pass through a 1×1×1 convolution, GELU and instance norm to give Z, then gate by a squeeze-and-excitation weight. The result serves both as the encoder skip and as the next stage’s input, after a stride-2 downsample doubling the channels (32 → 64 → 128).

g = σ( W2 GELU( W1 GAP(Z) ) )     Ffuse = g Z
Bottleneck256 ch · 1/16

Two stacked TSMambaBlocks model relative femur–tibia alignment. Fine boundary detail is reintroduced downstream by the skips and the edge decoder.

Mask streamECA-gated

At each of three upscaling stages, a stride-2 transposed convolution concatenates with an ECA-recalibrated skip and is refined by two 3×3×3 conv → GELU → IN blocks.

M(ℓ) = Conv3[ ConvT( M(ℓ−1) ) (ℓ) ]
Edge streamPEE

The boundary stream takes the raw, un-gated skip — ECA’s recalibration would attenuate exactly the high-frequency content the extractor needs. A Progressive Edge Extractor subtracts multi-scale low-pass responses (3D average pooling, k ∈ {3, 5}) from a squeezed feature to expose scale-specific unsharp detail.

E(ℓ) = PEE[ ConvT( E(ℓ−1) ) S(ℓ) ]
Cross-task fusionedge → mask

Boundary evidence is injected into the semantic stream in one direction only, at every decoder level, and the fused representation is what propagates onward. That matters most for cartilage spanning two to three voxels, where region overlap alone is ambiguous.

(ℓ) = GELU( Conv1[ M(ℓ) E(ℓ) ] )     DS = Conv1( M̃(ℓ) )
03
Loss & setup

L = αLDice + (1 − α)LCE + λELEdge + Ldeep
Dice, class-weighted
Weights fixed per class to penalise errors on the smaller cartilage classes more heavily.
Distance-weighted CE
Weights computed per sample from a Euclidean Distance Transform, exponentially upweighting voxels directly at the bone–cartilage interface. A schedule parameter α trades it off against Dice over training.
Edge loss λₑ = 0.1
Supervises the boundary decoder against ground-truth edges generated internally by morphological max pooling and erosion.
Deep supervision
Applied to intermediate mask-decoder outputs with decaying weights, preventing gradient dilution across the upscaling stages.

Training setup

Dataset
OAIZIB-CM — 507 knee MRIs, expert-annotated for five classes, each graded on the Kellgren–Lawrence scale
Voxel spacing
0.365 × 0.365 × 0.7 mm, native anisotropic
Training patches
64 × 192 × 192 sub-volumes, cropped on all three axes and biased toward bone and cartilage so thin classes are adequately represented
Sampling
KL-grade-aware weighted random sampler oversamples severe OA, guaranteeing exposure to osteophytic morphology
Inference
Full volume reconstructed from overlapping patch predictions by Gaussian-weighted sliding-window aggregation
Schedule
200 epochs, five-fold cross-validation
Optimizer
AdamW, lr 1×10−4 with cosine annealing to 1×10−6
Batch
2 per step, gradient accumulation over 4 steps → effective batch 8 (memory-bound)
Hardware
Single NVIDIA A100 80 GB
04
Results

MBEDNet matches the larger MSPF-VM-Unet with 23% fewer parameters, and its cartilage-focused training improves cartilage segmentation by roughly 1–2% over SwinUNETR. Boundary-aware Mamba encoders capture thin-structure geometry more effectively than pure transformer or heavier Mamba–CNN architectures.

Parameter count against mean Dice

MBEDNet in blue, baselines in grey. Up and to the left is better.

90919293025 M50 M75 M100 MPARAMETERSMEAN DSC (%)SwinUNETR62.2 M · 92.29MtRA-Unet96.5 M · 93.05MSPF-VM-Unet22.7 M · 93.05MBEDNet17.4 M · 93.11
Table 1 · Segmentation performance across the test set (DSC, %)
MethodParamsFBTBFCTCAvg
U-Net31.0 M96.4396.0584.7383.2190.11
V-Net71.0 M97.0496.5287.7584.0991.35
nnU-Net30.0 M97.8797.1788.5185.4992.26
I-Mask-RCNN44.0 M97.6097.5383.3280.6789.78
SwinUNETR62.2 M98.4698.4787.7784.4792.29
MtRA-Unet96.5 M98.5398.4489.1986.0293.05
VM-UNet27.0 M97.3596.3185.3383.1790.54
MSPF-VM-Unet22.7 M98.6198.3789.5785.6393.05
MBEDNet (ours)17.4 M98.4698.6489.0386.2993.11

Read this carefully. Baseline results are 5-fold cross-validation means reproduced from the MSPF-VM-Unet paper; MBEDNet and SwinUNETR were each trained on OAIZIB-CM under the same protocol with separate independent splits, and MBEDNet is evaluated on our held-out test set (n=102). No significance testing is reported for this table — without per-case predictions under a matched protocol, the paired, fold-dependence-aware inference a valid comparison requires is not possible. FB femoral bone · TB tibial bone · FC femoral cartilage · TC tibial cartilage, the mean of medial and lateral.

Ablation

Removing the Mamba branch costs the most, confirming that global context is essential for joint segmentation; removing the edge decoder degrades HD95 specifically, reflecting its role in preserving thin cartilage boundaries.

Dice lost when a component is removed

Change in five-class mean DSC against the full model (92.00%). Percentage points.

0.00.51.01.52.0DSC LOST vs FULL MODEL (POINTS)w/o Mamba branch2.20w/o Edge decoder1.30w/o Mamba bottleneck1.00w/o Edge loss0.40
Table 2 · Ablation across the five classes
ConfigurationDSC ↑ (%)HD95 ↓ (mm)ASSD ↓ (mm)
w/o Mamba branch89.801.440.400
w/o Edge decoder90.701.290.360
w/o Bottleneck (Mamba → conv)91.001.140.330
w/o Edge loss91.601.030.310
Full model92.000.950.289

HD95 — 95th percentile Hausdorff distance. ASSD — average symmetric surface distance. Table 2 averages over all five annotated classes, so its DSC is not directly comparable with the four-structure average in Table 1. Model differences were assessed on per-case Dice with a two-sided Wilcoxon signed-rank test paired across the 102 held-out cases; per-structure p-values were Holm–Bonferroni corrected.

Across disease severity

Scores hold up across KL grades, with a modest decline at grade IV where cartilage degradation is most pronounced and small spurious voxels occasionally appear near severely eroded interfaces.

Mean DSC by KL grade

Higher is better. Sample sizes below each grade.

9293940n=201n=122n=213n=274n=1592.18KELLGREN–LAWRENCE GRADE

HD95 by KL grade

Lower is better. Millimetres.

0.600.801.000n=201n=122n=213n=274n=151.07KELLGREN–LAWRENCE GRADE
Table 3 · Performance across KL grades (DSC, %)
KL gradenFBTBFCTCAvg DSCHD95ASSD
02098.7698.8288.6285.9993.050.730.149
11298.7498.7189.7488.0693.810.700.141
22198.5998.6288.4886.9293.150.790.151
32798.5898.6688.8884.8492.740.820.167
41598.6498.3988.7782.9392.181.070.238
All graded9598.6698.6488.9085.7592.990.950.289

Of the 102 test cases, 7 carry no KL grade and are excluded here. Grade 1 sits highest but rests on only 12 cases — the grade III and IV rows, at 27 and 15 cases, carry the weight of the robustness claim.

05
Clinical

Predicted masks were assessed surgically — surface geometry, landmark consistency, registration and tracking — on KL grade IV femoral segmentation, the hardest case the dataset offers. Predicted and ground-truth masks were converted to triangular meshes, and an in-house annotation module localised eight clinically relevant femoral landmarks across five grade IV test cases.

Mean landmark error
0.000mm
Euclidean, predicted vs ground-truth mesh (SD 0.528 mm), across eight femoral landmarks on five KL IV cases.
Target registration error
0.000mm
IR stereo camera and blunt probe, DRB space.
Tool tip to bone surface
0.000mm
Well inside clinical safety margins.
Per-structure STL distance
0.34 / 0.45mm
Femur / tibia, from the tracked phantom.
Fig. 2B–C · Phantom trackingIR stereo camera · blunt probe
Left: a 3D-printed femur phantom on a dynamic reference base touched with a tracked blunt probe, beside the registered STL showing a tool-tip distance of 0.34 millimetres. Right: the same setup for the tibia, at 0.45 millimetres.
Fig. 2B–C. A randomised test subject was 3D-printed as a physical phantom and evaluated under real-time optical tracking. Predicted landmarks were registered to manually annotated landmarks in DRB space; the segmented STL was then transformed by that registration and the tracked tool-tip distance to the bone surface measured.
06
Limits

No significance test vs baselines
Table 1 baselines are reproduced from prior work rather than recomputed on this test set, so a valid paired comparison is not available.
Grade IV is still the weak point
Tibial cartilage falls to 82.93% and HD95 rises to 1.07 mm at KL IV; small spurious voxels occasionally appear near severely eroded interfaces.
Single dataset, single modality
All results are on OAIZIB-CM knee MRI. Adapting the architecture to CT is named as future work, alongside further noise reduction and cartilage-specific improvements.

The conclusion the authors do draw: MBEDNet matches, and in places outperforms, existing state-of-the-art architectures at a lower computational cost, and stays diagnostically reliable in advanced disease — giving surgeons accurate preoperative masks regardless of severity.

07
Cite
@inproceedings{mbednet,
  title     = {MBEDNet: A Mamba-CNN Network with Dual Edge-Guided
               Decoders for Segmentation of TibioFemoral Bone
               and Cartilage},
  author    = {Anonymized Authors},
  note      = {Under double-blind review},
  url       = {https://anonymous.4open.science/r/MBEDNet-D6EC}
}

Author and venue fields are withheld while the submission is under double-blind review.

References

  1. Atabansi, C.C., et al.: DCM-Net: dual-encoder CNN-Mamba network with cross-branch fusion for robust medical image segmentation. BMC Medical Imaging 25(1), 395 (2025)
  2. Chen, J., et al.: TransUNet: Transformers make strong encoders for medical image segmentation. arXiv:2102.04306 (2021)
  3. Daydar, A., Pramanick, A., Sur, A., Kanagaraj, S.: Segmentation of tibiofemoral joint tissues from knee MRI using MtRA-Unet. arXiv:2401.12932 (2024)
  4. Du Toit, C., et al.: Automatic femoral articular cartilage segmentation using deep learning in 3D ultrasound images of the knee. Osteoarthritis and Cartilage Open 4(3), 100290 (2022)
  5. Felfeliyan, B., et al.: Improved-Mask R-CNN: towards an accurate generic MSK MRI instance segmentation platform. Comput. Med. Imaging Graph. 97, 102056 (2022)
  6. Hatamizadeh, A., et al.: Swin UNETR: Swin transformers for semantic segmentation of brain tumors in MRI images. arXiv:2201.01266 (2022)
  7. Humayun, A., Rehman, M., Liu, B.: A method framework of semi-automatic knee bone segmentation and reconstruction from CT images. Quant. Imaging Med. Surg. 14(10), 7151–7175 (2024)
  8. Isensee, F., et al.: nnU-Net: a self-configuring method for deep learning-based biomedical image segmentation. Nature Methods 18(2), 203–211 (2021)
  9. Kakavand, R., et al.: Integration of Swin UNETR and statistical shape modeling for semi-automated segmentation of the knee. Scientific Reports 14, 2748 (2024)
  10. Kakavand, R., et al.: Swin UNETR segmentation with automated geometry filtering for biomechanical modeling of knee joint cartilage. Ann. Biomed. Eng. 53(4), 908–922 (2025)
  11. Kompella, G., et al.: Segmentation of femoral cartilage from knee ultrasound images using Mask R-CNN. In: EMBC, pp. 966–969 (2019)
  12. Lee, C.Y., Xie, S., Gallagher, P., Zhang, Z., Tu, Z.: Deeply-supervised nets. In: AISTATS (2015)
  13. Liu, Z., et al.: Swin Transformer: hierarchical vision transformer using shifted windows. In: ICCV, pp. 10012–10022 (2021)
  14. Marsilio, L., et al.: CEL-UNet: a novel CNN architecture for 3D segmentation of knee bones affected by severe osteoarthritis. In: EMBC (2022)
  15. Milletari, F., Navab, N., Ahmadi, S.A.: V-Net: fully convolutional neural networks for volumetric medical image segmentation. In: 3DV, pp. 565–571 (2016)
  16. Murugesan, B., et al.: Recon-GLGAN: a global-local context based GAN for MRI reconstruction. In: MLMIR @ MICCAI, LNCS 11905, pp. 3–15 (2019)
  17. Polat, G., et al.: Class distance weighted cross-entropy loss for ulcerative colitis severity estimation. In: MIUA, pp. 157–171 (2022)
  18. Ronneberger, O., Fischer, P., Brox, T.: U-Net: convolutional networks for biomedical image segmentation. In: MICCAI, pp. 234–241 (2015)
  19. Ruan, J., Li, J., Xiang, S.: VM-UNet: vision Mamba UNet for medical image segmentation. arXiv:2402.02491 (2024)
  20. Sudre, C.H., et al.: Generalised Dice overlap as a deep learning loss function for highly unbalanced segmentations. In: DLMIA, LNCS 10553, pp. 240–248 (2017)
  21. Varacallo, M.A., et al.: Total knee arthroplasty techniques. In: StatPearls (2024)
  22. Wang, Q., et al.: ECA-Net: efficient channel attention for deep convolutional neural networks. In: CVPR, pp. 11534–11542 (2020)
  23. Wang, X., et al.: VM-UNet enhanced with multi-scale pyramid feature extraction for segmentation of tibiofemoral joint tissues from knee MRI. PLoS One 20(8), e0330740 (2025)
  24. Wang, Z., Zheng, J., Zhang, Y., et al.: Mamba-UNet: UNet-like pure visual Mamba for medical image segmentation. arXiv:2402.05079 (2024)
  25. Xing, Z., Ye, T., Yang, Y., et al.: SegMamba: long-range sequential modeling Mamba for 3D medical image segmentation. In: MICCAI, pp. 578–588 (2024)
  26. Yang, M., et al.: Automated knee cartilage segmentation for heterogeneous clinical MRI using GANs with transfer learning. Quant. Imaging Med. Surg. 12(5), 2620–2633 (2022)
  27. Yao, Y.: OAIZIB-CM: dataset from the CartiMorph project. Zenodo, v1.0 (2025)
  28. Zhou, Z., et al.: UNet++: a nested U-Net architecture for medical image segmentation. In: DLMIA, pp. 3–11 (2018)