---
title: Skeletal Muscle Area (SMA) in Medical Imaging
url: https://www.emergentmind.com/topics/skeletal-muscle-area-sma
type: topic
---

# Skeletal Muscle Area (SMA) in Medical Imaging

Searching arXiv for the provided SMA-related papers to ground the article in current literature.
Skeletal Muscle Area (SMA) is the cross-sectional area of muscle tissue on a single image plane. In radiologic practice, recent work defines SMA as the cross-sectional area of all muscle tissue, measured in cm², on a single axial CT slice; related MRI and ultrasound workflows compute the same quantity from segmented muscle pixels and physical pixel area. SMA is used in longitudinal monitoring of muscle atrophy and recovery in intensive care, in sarcopenia assessment, and in identification and tracking of cancer cachexia. A major standardization convention is the mid-third lumbar vertebral level (L3), where prior studies cited in recent work report a high correlation (\(r\approx0.95\)) between muscle area and whole-body muscle volume, whereas longitudinal ultrasound studies emphasize retrieval of the same anatomical view across visits rather than a vertebral landmark [2503.16556] [2306.04739].

## 1. Definition, anatomical scope, and derived indices

In the CT-based oncology literature, SMA is defined as the cross-sectional area of all muscle tissue on a single axial CT slice, measured in cm² [2503.16556]. In one CT sarcopenia workflow, manual annotation at L3 explicitly included the psoas, erector spinae, quadratus lumborum, and abdominal wall muscles, and the traced contour was rasterized into a binary mask whose foreground pixels corresponded to the expert’s SMA measurement [2508.17275]. MRI and ultrasound studies use the same geometric principle: a muscle mask is estimated on a selected slice or view, and area is computed from the count of muscle pixels and the physical pixel area [2409.06942] [2306.04739].

The L3 convention is central in abdominal CT and MRI body-composition analysis. Recent work cites Shen et al. (2004) and Mourtzakis et al. (2008) for the observation that muscle area at the L3 level correlates highly with whole-body muscle volume, with \(r\approx0.95\), and notes that L3 is readily identified on routine abdominal CT scans as a standardized, reproducible anatomical landmark [2503.16556]. By contrast, longitudinal ultrasound of the rectus femoris in ICU patients requires repeated acquisition of the same anatomical “slice” over time, because the principal source of error is not only boundary delineation but also inconsistency in view matching across days and operators [2306.04739].

A derived index used in oncology is the Skeletal Muscle Index (SMI), defined as
$$
\mathrm{SMI} \;=\; \frac{\mathrm{SMA}}{\bigl(\text{height in m}\bigr)^2}
\quad \bigl[\mathrm{cm}^2/\mathrm{m}^2\bigr].
$$
BMI is used alongside SMI to adjust for overall body habitus in downstream predictive models [2503.16556]. A common misconception is to treat SMA and SMI as interchangeable; the cited formulation distinguishes absolute cross-sectional muscle area from a height-normalized index.

## 2. Image acquisition and slice or view localization

Accurate SMA estimation begins with correct slice or view selection. In the SMAART-AI CT pipeline, all DICOM CT series are retrieved from PACS/VNA, axial series are identified via the `ImageOrientationPatient` header, and TotalSegmentator is applied to label 117 structures; slices with vertebra index \(=29\) are then selected as L3 [2503.16556]. This makes slice localization an explicit computational stage rather than an informal precondition.

MRI-based L3 workflows can still depend on expert localization. In DAFS Express, all In-Phase Dixon MRI scans are imported, metadata such as pixel spacing and slice thickness are read, and a trained rater visually identifies and extracts the true axial slice at the L3 vertebral level before automated segmentation proceeds [2409.06942]. The MRI workflow also includes resampling to a consistent in-plane resolution and definition of a loose abdominal ROI to accelerate segmentation [2409.06942].

A separate CT sarcopenia workflow used manual L3 identification by a trained radiologist, who selected the single transverse slice passing through the inferior endplate of the third lumbar vertebra and exported it as a 2-D bitmap image [2508.17275]. Preprocessing then re-oriented all 3-D volumes to RAS, resampled voxel spacing to \(1 \text{ mm}\times1 \text{ mm}\times1 \text{ mm}\), clipped HU values to \([-175,+250]\), linearly mapped intensities to \([0,1]\), and extracted the L3 slice for downstream processing [2508.17275].

In longitudinal ultrasound, the localization problem is different. Skeletal muscle cross-sectional area measurement over time consists of three steps: retrieving the same anatomical slice at each visit, segmenting the muscle boundary, and converting that segmentation into an area measurement [2306.04739]. The first step is especially critical because manual ultrasound measurements are time-consuming and suffer large inter- and intra-operator variability due to inconsistent probe placement, angulation, and subjective tracing of the muscle border [2306.04739].

## 3. Computational measurement pipelines

Ultrasound SMA tracking in ICU patients has been formulated as a retrieval-plus-segmentation problem. A self-supervised contrastive model inspired by SimCLR augments each frame \(X\) twice by random crop and flip to form a positive pair \((X_i,X_j)\). A shared encoder \(f(\cdot)\), consisting of four \(3\times3\) convolution layers with ReLU and two max-pooling layers, maps each view to a representation \(h=f(X)\), and a two-layer MLP projection head \(g(\cdot)\) with 512 units per layer yields \(z=g(h)\). The contrastive loss is
$$
\mathcal{L}_{contra}
= -\log\frac{\exp(\operatorname{sim}(z_{i},z_{j})/\tau)}
{\sum_{k=1}^{2N}\mathbbm{1}_{[k\neq i]}\exp(\operatorname{sim}(z_{i},z_{k})/\tau)},
$$
where \(\operatorname{sim}(\cdot,\cdot)\) is cosine similarity and \(\tau\) is a temperature parameter. After 500 epochs of self-supervised pretraining, the encoder is frozen and its outputs feed a classifier with four dense layers of \(2048\rightarrow1024\rightarrow512\rightarrow2\) units and softmax to predict whether two frames correspond to the same anatomical view [2306.04739]. The retrieved frame is then segmented by a standard U-Net trained on 1,000 manually annotated rectus-femoris cross-sections using Adam optimization, \(\ell_2\) regularization, and data augmentation [2306.04739].

The CT oncology pipeline in SMAART-AI is fully automated from image retrieval through prognostic modeling. After DICOM retrieval and L3 identification, HU values are windowed to \([-29,150]\) to isolate muscle density, converted to PNG, and normalized to zero mean and unit variance. A 2D nnU-Net receives one \(512\times512\) normalized PNG slice and outputs a pixel-wise probability map and binary mask using a threshold of 0.5. The architecture is a 2D U-Net with 7 encoder and 7 decoder stages, skip connections, instance normalization, and LeakyReLU activations; two variants are described, PlainConvUNet and ResidualEncoderUNet, the latter using deeper residual blocks in the encoder. Training uses the nnU-Net default composite Dice plus cross-entropy loss and stochastic gradient descent with initial learning rate \(1\times10^{-2}\), under 5-fold cross-validation, with average pixel probabilities across folds used at inference [2503.16556].

DAFS Express applies a related MRI segmentation strategy to single L3 slices. The software imports DICOM data, applies a proprietary intensity-normalization routine to standardize fat-water contrast, and uses a U-Net-style CNN to classify each pixel as SKM, VAT, or SAT. For SKM, the network output is thresholded to a binary mask, small isolated islands under 50 pixels are removed, and holes within the SKM mask are filled. The final mask is converted into a tag file, the cross-sectional area is computed, mean signal intensity within the SKM mask is recorded for potential fat-infiltration analyses, and a one-page PDF report is generated per slice [2409.06942].

A further CT workflow couples transfer learning for sarcopenia classification with self-supervised segmentation for quantitative SMA estimation. The classification arm replaces the final layer of Densenet-121, InceptionV3, InceptionResNetV2, and ResNet-50 backbones pretrained on RadImageNet with a two-node SoftMax head. The segmentation arm, termed SMIT, uses a Vision Transformer in a teacher-student self-distillation framework with masked image modeling, pretrained on more than 3,500 unlabeled CT scans and then fine-tuned on 79 labeled CT volumes with masks. This student network outputs a pixel-wise skeletal-muscle mask \(M_{pred}\) [2508.17275].

## 4. Quantification formulas, evaluation metrics, and reported performance

Across modalities, SMA is computed from a binary muscle mask and pixel geometry. In the ultrasound formulation,
$$
A \;=\; \sum_{i=1}^{N} p_i \,\times s^2,
$$
where \(N\) is the total number of pixels, \(p_i\in\{0,1\}\) indicates whether pixel \(i\) belongs to muscle, and \(s\) is the known pixel spacing; no further coordinate transform is needed beyond pixel-to-millimeter conversion [2306.04739]. The MRI formulation is equivalent:
$$
A_{SKM} \;=\; N_{pixels}^{SKM}\;\times\;A_{pixel},
$$
with \(A_{pixel}\) obtained from DICOM metadata [2409.06942].

Reported quantitative results span retrieval, segmentation overlap, area error, and reliability metrics:

| Setting | Reported results | Study |
|---|---|---|
| ICU ultrasound longitudinal exams | AUC \(=73.52\%\); precision \(=67.2\%\); recall \(=68.3\%\); F1 \(=67.7\%\); mean relative absolute CSA error \(=5.7\%\pm0.24\%\) | [2306.04739] |
| CT cancer cachexia pipeline | Mean Dice \(=97.80\%\pm0.93\%\); median absolute error \(=2.48\%\); uncertainty correlations: variance \(r\approx0.83\), entropy \(r\approx0.76\), CV \(r\approx0.73\) | [2503.16556] |
| MRI L3 body-composition analysis | Mean SKM area automated \(=132.51\ \text{cm}^2\); manual \(=132.36\ \text{cm}^2\); bias \(=0.15\ \text{cm}^2\); Bland-Altman limits \([ -3.73, +4.03]\ \text{cm}^2\); Jaccard \(=99.03\%\); Dice \(=99.51\%\); ICC \(=0.998\); runtime \(=18\) s per DICOM | [2409.06942] |
| CT sarcopenia segmentation | Mean Dice \(=0.93\pm0.05\); mean SMA error \(=3.0\pm7.8\) percentage-points; average fold errors from \(-3.2\%\) to \(10.0\%\); after cut-off thresholding, held-out sarcopenia detection accuracy \(=100\%\) | [2508.17275] |

The overlap metrics are standard. The Dice similarity coefficient is written as
$$
\mathrm{Dice} \;=\;\frac{2\,|A\cap B|}{|A|+|B|},
$$
where \(A\) is the predicted mask and \(B\) is the manual mask [2503.16556]. The MRI study additionally reports the Jaccard coefficient, Intraclass Correlation Coefficient under a two-way mixed absolute-agreement model, and Bland-Altman bias and limits of agreement, thereby distinguishing pixel overlap from agreement in area measurements [2409.06942]. This distinction matters because reliable SMA quantification is not exhausted by boundary overlap alone.

## 5. Clinical roles in critical illness, sarcopenia, and cancer cachexia

In critical care, skeletal muscle atrophy is common in ICU patients who spend long periods in bed, and muscle mass must be recovered through physiotherapy before discharge. Ultrasound is frequently used to assess recovery by measuring muscle size over time, but manual measurements suffer from large variability. Automated retrieval of the same cross-section at each visit, combined with automated U-Net segmentation, reduces random and systematic errors in longitudinal SMA tracking and enables more precise monitoring of muscle atrophy and recovery [2306.04739]. The same study states that objective, reproducible quantification may guide the timing and intensity of physiotherapy, help predict functional outcomes, and facilitate large-scale studies of muscle wasting interventions [2306.04739].

In oncology, progressive loss of SMA is presented as a hallmark of cancer cachexia, a multifactorial syndrome characterized by involuntary weight and muscle loss. Early detection of declining SMA can guide nutritional and exercise interventions, improve treatment tolerance, and is strongly prognostic of survival [2503.16556]. SMAART-AI extends SMA estimation into a multimodal cachexia-prediction workflow by combining SMA, SMI, BMI, and clinical variables such as age, sex, race, ethnicity, and cancer stage in a multi-layer perceptron for binary cachexia classification at diagnosis. The described architecture uses three hidden linear layers with LeakyReLU and dropout followed by a sigmoid output; reported evaluation includes \(79\%\) precision for predicting cachexia, \(70\%\) overall accuracy, and an F1-score of \(76.9\%\) [2503.16556].

For sarcopenia detection in CT, one study reports that quantitative assessment of SMA is more precise and informative than qualitative assessment models for detecting sarcopenia [2508.17275]. After segmentation, SMA can be thresholded by published cut-offs of male \(<144\ \text{cm}^2\) and female \(<92\ \text{cm}^2\), yielding \(100\%\) sarcopenia detection accuracy on held-out test sets in that dataset [2508.17275]. This frames SMA not merely as a descriptive measurement but as a quantitative intermediate phenotype for diagnostic stratification.

## 6. Reliability, variability, and recurring methodological issues

The principal threats to reliable SMA estimation differ by modality. In ultrasound, variability arises from inconsistent probe placement, angulation, and subjective tracing of the muscle border; the retrieval stage is therefore indispensable, not optional [2306.04739]. In MRI, documented variability sources include poor image quality such as motion blur and stitching artifacts, unusual anatomy including hip prostheses and severe scoliosis, varying voxel sizes, and rater inconsistency during manual correction, although the reported ICC of 0.998 indicates that rater-induced disagreement was largely mitigated in that study [2409.06942].

Uncertainty quantification has become an explicit mechanism for reliability control in automated CT SMA pipelines. SMAART-AI computes variance, entropy, and coefficient of variation from an ensemble of models, and flags cases above dataset-specific uncertainty thresholds for expert correction. The reported Pearson correlations between these uncertainty measures and absolute SMA error—approximately \(0.83\) for variance, \(0.76\) for entropy, and \(0.73\) for coefficient of variation—indicate that elevated uncertainty reliably flags likely high-error cases [2503.16556]. This directly addresses the misconception that high mean segmentation performance eliminates the need for human review.

Another recurring methodological issue is conflation of slice selection with segmentation accuracy. The workflows summarized here make slice or view localization a separate stage: TotalSegmentator-based L3 identification in CT, expert L3 selection in MRI and in one CT study, and contrastive retrieval of corresponding ultrasound views over time [2503.16556] [2409.06942] [2508.17275] [2306.04739]. A plausible implication is that reproducible SMA estimation depends jointly on anatomical standardization, acquisition consistency, and segmentation quality, rather than on the segmentation network alone.

Source: https://www.emergentmind.com/topics/skeletal-muscle-area-sma