US Speckle Models — Part 3: m as a Tissue Map

Two posts ago I argued that ultrasound speckle isn't random noise — it's the interference pattern of many sub-resolution scatterers, and its statistics carry information about the tissue that produced it. Rayleigh distribution was the classical model: mathematically elegant, physically motivated, but built on assumptions that break in real tissue.

Nakagami generalized Rayleigh. Same physical intuition, one extra parameter, and the freedom to adapt to whatever scatterer density the tissue actually has.

Last post I laid out the Nakagami PDF and what each term does. This post is where it stops being a math exercise and starts being clinically useful.

From distribution to map

The Nakagami PDF is:

$$p(x;\, m,\, \omega) \;=\; \frac{2 m^m}{\Gamma(m)\, \omega^m}\; x^{2m-1}\, e^{-m x^2 / \omega}, \qquad x \geq 0$$

Two parameters carry the tissue information:

• m — shape parameter. Low m spreads the distribution out (few scatterers, high uncertainty in the envelope amplitude). High m sharpens the peak (many scatterers, predictable behavior).
• ω — scale parameter, equal to
  $$E[x^2]$$
  . Sets the local mean intensity, normalizes for depth attenuation and gain.

At m = 1, the whole expression collapses to Rayleigh. Nakagami includes Rayleigh as a single point in a continuum, then extends both directions — sub-Rayleigh (fewer scatterers, edges and boundaries) and super-Rayleigh (denser, more organized tissue).

Here's the pivot that makes it clinically useful: m can be estimated locally. Slide a small window across the image, compute the method-of-moments estimator inside it:

$$\hat{m} \;=\; \frac{\big(\mathbb{E}[X^2]\big)^2}{\mathrm{Var}(X^2)}, \qquad \hat{\omega} \;=\; \mathbb{E}[X^2]$$

The output isn't a denoised image. It's a map — one m value per window, spatially varying, tracking whatever tissue characteristic changes as you move across the scan.

Canonical envelope-domain values

Published literature converges on a rough three-regime interpretation, all measured on raw envelope amplitude before any display processing:

• m ≈ 0.5 — lesion boundaries, cyst walls, tissue interfaces. Few scatterers per resolution cell, high variance in the envelope.
• m ≈ 1.0 — fully developed speckle. What Rayleigh assumed everywhere. Homogeneous soft tissue at sufficient depth.
• m > 1.5 — dense parenchyma, fibrous tissue, structured collagen. Many scatterers, predictable statistics.

Same image, three regimes, one model. That's the promise.

From the thesis: real B-mode, real fit

Here's what happens when you actually run this on a real hospital scan.

Three 64×64 windows on a B-mode frame from my PhD data:

• A — inside a dense, near-uniform lesion. High m: histogram is narrow, PDF nearly Gaussian-looking.
• B — a homogeneous parenchyma region. Mid-range m, classic developed-speckle shape.
• C — a heterogeneous zone with varying intensity. Low m, right-skewed distribution.

The relative ordering matches the canonical interpretation: dense → high m, heterogeneous → low m. But the absolute m values are much higher than the 0.5–2 range the literature reports. The reason is what the image is: this is not envelope amplitude. This is the display frame — log-compressed, gain-adjusted, gamma-corrected, 8-bit quantized.

Log compression squashes the dynamic range. A distribution that looked broad and skewed in the envelope domain becomes narrower and more symmetric after compression. Method-of-moments picks that up as a higher m estimate. It's not wrong — it's just measuring m on a different variable than the textbook does.

For most practical work, the relative m ordering carries the tissue information. Absolute values shift with the specific compression curve of the scanner, but the ranking of regions doesn't.

I'll come back to why this matters at the end.

For contrast: what the textbook world looks like

If I simulate the same three regimes in the envelope domain — Nakagami-sampled directly, no compression, on a synthetic sector scan geometry — this is what I get:

Two things worth noting from the contrast:

First, the m values here (0.5, 1.0, 1.7) sit exactly where the canonical literature puts them. That's because this is synthetic envelope amplitude, drawn directly from Nakagami distributions with those m values. No log compression, no display processing in between. The values are the ground truth, not an estimate contaminated by a display pipeline.

Second, the PDF shapes are visibly cleaner — sharper mode separation, less overlap between regimes. Real clinical data doesn't look like this. Once you go through display processing, distributions compress toward each other, m estimates shift, and the model has to accommodate what the scanner actually outputs, not what a textbook assumes.

Where a lot of QUS quietly breaks

This is where a lot of published quantitative ultrasound work quietly breaks when it hits real hospital data.

Most method development happens on simulated envelope data or on lab-grade research scanners that expose raw RF or beamformed envelope. That's the world where the canonical m values live. It's a clean world, and the math works there.

Real clinics don't run research scanners. They run commercial systems that output DICOM images with log compression, dynamic range mapping, and a proprietary post-processing pipeline you don't have access to. If you take a method validated on envelope-domain simulation and drop it on a hospital scanner's video output, three things happen:

• Absolute m values shift out of the calibration range the method assumes
• Distributions overlap more than the training data, so region classification degrades
• The relationship between m and underlying tissue reflectivity is no longer direct — there's a monotonic-but-nonlinear display curve sitting in between

The math is right. The domain assumption isn't.

You can still do useful QUS on log-compressed data. But you have to design for it — pick estimators that are compression-tolerant, work with relative rather than absolute m, and calibrate on the specific scanner's output curve. That's a different research program than "assume envelope, apply Nakagami."

Next: the other thing Nakagami by itself misses

There's a second problem, independent of log compression. Nakagami models the amplitude distribution at each pixel. It says nothing about the relationship between pixels. And in ultrasound, that relationship is never zero: the PSF spreads every scatterer's signal across a neighborhood, and neighboring pixels are always weighted combinations of overlapping scatterer sets.

Standard local m estimation treats those pixels as if they were independent samples. They aren't. Next post: why the iid assumption breaks, what it costs you, and what you have to fold into the observation model to fix it.

---

References

• My paper on Bayesian restoration for B-mode ultrasound, accounting for correlated speckle noise: Ultrasonic Imaging, 2019
• ORCID (all publications): 0000-0003-2633-9371

---

I take on selective remote consulting on ultrasound image quality, quantitative ultrasound, and restoration problems — particularly the messy real-clinical-data variety this post is about. Reach out via LinkedIn if that's the kind of problem you're working on.