---
title: 'EMC²: Motion-Compensated Ultrasound Imaging'
url: https://www.emergentmind.com/topics/estimated-motion-compensated-compounding-emc2
type: topic
---

# EMC²: Motion-Compensated Ultrasound Imaging

Searching arXiv for the specified EMC² papers and closely related context.
Searching arXiv for the specified EMC² papers and closely related context.
Estimated Motion-Compensated Compounding (EMC²) is an ultrasound image formation strategy in which partially formed images acquired across multiple transmit events are first used to estimate inter-acquisition motion, then spatially realigned, and finally combined by coherent summation. In the literature summarized here, the term appears in two closely related but architecturally distinct settings: as an extension to Recursive Aperture Decoded Ultrasound Imaging (READI) for Fast Orthogonal Row-Column Electronic Scanning (FORCES), and as a real-time framework for ultrafast echocardiography that combines motion-corrected coherent compounding, harmonic imaging, and angular-coherence weighting [2509.08781] [2312.01628]. In both settings, EMC² addresses the same central problem: motion between constituent acquisitions degrades compounding gain, tissue speckle, contrast, and sharpness unless the sub-images are aligned before summation.

## 1. Terminological scope and problem setting

EMC² is motivated by the motion sensitivity of coherent compounding. In FORCES, motion sensitivity arises from ensemble size and aperture encoding; in ultrafast echocardiography, it arises because successive steered diverging-wave transmissions are acquired at different angles while tissue continues to move. The common remedy is to estimate displacement between sub-images or between successive pulse-echo data, warp the data into a common reference geometry, and preserve complex phase so that coherent summation remains valid.

A concise comparison of the two documented formulations is useful.

| Context | Motion-estimation basis | Final compounded image |
|---|---|---|
| READI with FORCES | Block-matching with normalized cross-correlation on low-resolution READI sub-images | \(F_{\mathrm{EMC}^2}(r)=\sum_{s=1}^S \hat{R}_s(r)\) |
| Ultrafast echocardiography | Cross-correlation of complex baseband signals or Doppler-based ensemble autocorrelation | \(y_{\mathrm{EMC}^2}(r)=\mathrm{ACF}(r)\,y_{\mathrm{comp}}(r)\) |

The literature therefore does not present EMC² as a single immutable algorithm. Rather, it designates a class of motion-aware coherent-compounding procedures whose concrete implementation depends on the acquisition scheme. This suggests that the acronym names a pipeline role—estimated motion followed by motion-compensated compounding—more than a unique estimator or beamformer.

## 2. READI-based EMC² in FORCES imaging

In the READI formulation, EMC² is built on a decomposition of an \(N\)-line FORCES acquisition \(G(t)\in\mathbb{R}^N\) with \(N=S\cdot Q\) into \(S\) groups of \(Q\) Hadamard-encoded receive signals. For the \(s\)-th block \(g'_s(t)\in\mathbb{R}^Q\), partial decoding is performed as
\[
d'_s(t)=H_Q^{-1} g'_s(t),
\]
where \(H_Q\in\{\pm1\}^{Q\times Q}\) is the Sylvester Hadamard matrix of order \(Q\). Each partially decoded group is then beamformed into \(S\) sub-images by a modified delay-and-sum operator
\[
D'_l\{d'_s,r\}=\sum_{v=1}^Q d'_{s,v}\!\bigl(\tau_{v+Q(l-1)}(r)\bigr),
\]
with \(\tau_i(r)=|d_i(r)|/c\). These sub-images are combined using the \(s\)-th row of \(H_S^{-1}\),
\[
R_s(r)=\sum_{l=1}^S \hat{h}^{(S)}_{l,s}\,D'_l\{d'_s,r\},
\]
so that, by linearity, the full FORCES image is recovered as
\[
F\{G,r\}=\sum_{s=1}^S R_s(r).
\]

EMC² modifies this workflow by retaining each \(R_s(r)\) in complex analytic (RF) form. That choice is structurally important: later coherent summation is intended to preserve phase, so motion correction is applied before envelope detection. The pipeline proceeds by selecting one low-resolution image, typically \(R_1\), as a reference; estimating motion from each remaining \(R_s\) relative to that reference; interpolating the resulting sparse displacement estimates to a dense field; warping each low-resolution image; and summing the aligned complex images coherently [2509.08781].

This architecture exploits a specific property of READI: the low-resolution images are less susceptible to motion than the full FORCES reconstruction, yet they sum to form the complete FORCES image. The reported purpose of EMC² is therefore not merely post hoc registration, but restoration of a motion-corrupted encoded-aperture reconstruction by operating on intermediate sub-images that remain physically tied to the original acquisition.

## 3. Motion estimation, outlier rejection, and warping in READI EMC²

For READI-based EMC², motion estimation is performed by block matching with normalized cross-correlation (NCC). At each motion-grid location \(r_k\) in the reference image, a template patch
\[
T(u)=R_{\mathrm{ref}}(r_k+u), \qquad u\in\Omega_{\mathrm{ref}}
\]
is compared against candidate patches within a target-image search window
\[
\Omega_{\mathrm{search}}=\{r_k+p:\; p\in[-W,\ldots,+W]\times[-H,\ldots,+H]\}.
\]
The NCC score at offset \(p\) is
\[
\mathrm{NCC}(r_k;p)=
\frac{
\sum_{u\in\Omega_{\mathrm{ref}}}[R_{\mathrm{ref}}(r_k+u)-\mu_{\mathrm{ref}}]\,[R_{\mathrm{target}}(r_k+u+p)-\mu_{\mathrm{target}}]
}{
\sqrt{
\sum |R_{\mathrm{ref}}(r_k+u)-\mu_{\mathrm{ref}}|^2 \cdot
\sum |R_{\mathrm{target}}(r_k+u+p)-\mu_{\mathrm{target}}|^2
}
},
\]
and the integer-pixel displacement is chosen as
\[
p^*=\arg\max_{p\in\Omega_{\mathrm{search}}}\mathrm{NCC}(r_k;p).
\]

Sub-pixel refinement is obtained by fitting a \(2\times 2\)-D paraboloid to the \(5\times 5\) neighborhood around the integer-pixel peak and taking the maximum of that fitted surface. Spurious or unreliable vectors are rejected if any of three criteria hold: peak \(\mathrm{NCC}(r_k;p^*)<T_{\mathrm{abs}}\), \(\mathrm{NCC}(r_k;p^*)/\mathrm{NCC}(r_k;0)<T_{\mathrm{rel}}\), or the smallest eigenvalue of the paraboloid Hessian at \(p^*\) is below \(\kappa_{\min}\). The paper reports typical parameter ranges tuned per experiment: motion grid spacing \(4\)–\(16\) px, template size \(16\)–\(64\) px, search extension \(+8\)–\(64\) px, \(T_{\mathrm{abs}}=0.01\)–\(0.5\), \(T_{\mathrm{rel}}=1.01\)–\(1.20\), and \(\kappa_{\min}=0.005\)–\(0.1\) [2509.08781].

After sparse displacements \(u_s(r_k)\) are estimated, they are interpolated to full image resolution \(u_s(r)\). Each low-resolution image is then warped according to
\[
\hat{R}_s(r)=R_s(T_s(r)),\qquad T_s(r)=r+u_s(r),
\]
with complex-valued bilinear interpolation. If \(y=T_s(r)\) lies between integer samples, one writes \(y=\lfloor y\rfloor+\delta\) and evaluates
\[
\hat{R}_s(r)=\sum_{i=0}^1\sum_{j=0}^1
w_{i,j}(\delta_x,\delta_z)\,
R_s(\lfloor y_x\rfloor+i,\lfloor y_z\rfloor+j),
\]
using the usual bilinear weights. The final motion-compensated reconstruction is the phase-coherent sum
\[
F_{\mathrm{EMC}^2}(r)=\sum_{s=1}^S R_s(r+u_s(r))=\sum_{s=1}^S \hat{R}_s(r),
\]
followed optionally by envelope detection \(|F_{\mathrm{EMC}^2}(r)|\) and post-processing such as coherence-factor weighting or log compression.

## 4. EMC² in ultrafast echocardiography

A separate EMC² formulation was introduced for ultrafast echocardiography as a unified framework combining coherent compounding, second-harmonic imaging, and angular coherence for simultaneous high-quality B-mode and tissue Doppler. Here the acquisition consists of \(N\) steered diverging-wave transmissions emitted at successive angles. Motion between the \(n\)-th and \((n+1)\)-th acquisitions misaligns the echo fields and degrades compounding gain, so EMC² estimates inter-pulse motion either by cross-correlation of complex baseband signals or, more efficiently, by Doppler-based ensemble autocorrelation [2312.01628].

The Doppler-based estimator uses slow-time autocorrelations over ascending and descending halves of a triangular angular sweep,
\[
R_1(r)=\sum_{m=1}^{N/2-1} s_m(r)s_{m+1}^*(r),
\qquad
R_2(r)=\sum_{m=N/2}^{N-1} s_m(r)s_{m+1}^*(r).
\]
From the resulting phase \(\Phi(r)\), the radial Doppler velocity is implemented in pseudocode as
\[
V_{\mathrm{dop}}=\frac{c}{4\pi f_0 T_{\mathrm{prp}}}\,\Phi,
\]
and the spatial delay per transmit is
\[
\hat d_n(r)=\Bigl(n-\frac{N+1}{2}\Bigr)\hat V(r)T_{\mathrm{PRP}}.
\]
After beamforming each diverging wave into a complex IQ field \(y_n(r)\), motion-compensated coherent compounding is
\[
y_{\mathrm{comp}}(r)=\frac{1}{N}\sum_{n=1}^N y_n\!\bigl(r+\hat d_n(r)\bigr).
\]
In practice, the realignment is implemented by \(1\)D linear interpolation along the radial coordinate followed by phase correction \(\exp\{j\,(n-\tfrac{N+1}{2})\,\Phi(r)\}\).

The same framework incorporates second-harmonic imaging via pulse inversion. The receive signal is modeled as
\[
p(r,t)\approx p_1 e^{j(\omega t+kr)}+p_2 e^{j(2\omega t+2kr)},
\]
so that transmitting alternating positive and negative pulses and summing the echoes yields
\[
s^+(r,t)+s^-(r,t)\propto 2p_2 e^{j(2\omega t+2kr)},
\]
after which a band-pass filter around \(2\omega\) isolates the second harmonic. Both fundamental and second-harmonic channels undergo the same motion estimation and compounding steps.

Residual misalignment and clutter are further handled through a lag-1 angular coherence factor,
\[
\mathrm{ACF}(r)=
\frac{\sum_{n=1}^{N-1} y_{n,\mathrm{moco}}(r)\,y_{n+1,\mathrm{moco}}^*(r)}
{\sum_{n=1}^{N} |y_{n,\mathrm{moco}}(r)|^2},
\]
which weights the motion-compensed IQ image to produce
\[
y_{\mathrm{EMC}^2}(r)=\mathrm{ACF}(r)\,y_{\mathrm{comp}}(r).
\]
For B-mode display, \(|y_{\mathrm{EMC}^2}(r)|\) is envelope-detected and logarithmically compressed. In this formulation, EMC² therefore denotes not only motion-compensated compounding but a compound pipeline in which motion correction, harmonic clutter suppression, and angular-coherence weighting are explicitly coupled.

## 5. Quantitative performance across reported implementations

The two EMC² instantiations report quantitative gains under different motion regimes and acquisition models. In the READI/FORCES setting, the reported benchmarks emphasize recovery from probe motion, performance of low-resolution READI sub-images relative to sparse STA baselines, and speckle preservation in flow. In the ultrafast echocardiography setting, the reported benchmarks emphasize gCNR under controlled in vitro motion, Doppler accuracy, and in vivo performance during low- and high-velocity phases of the cardiac cycle.

| Setting | Condition | Reported result |
|---|---|---|
| READI/FORCES | Lateral probe motion at \(\sim 15\) cm/s | Standard FORCES gCNR fell by \(30\)–\(40\%\); EMC² restored gCNR to within \(2\)–\(5\%\) of the static reference |
| READI/FORCES | Same \(Q\) transmit count as uFORCES | A single READI sub-image achieved \(10\)–\(20\%\) higher generalized CNR than uFORCES sparse-Hadamard |
| READI/FORCES | \(Q=8\), \(S=16\), \(8\) MHz, PRF \(=4\) kHz | EMC² recovered coherent blood-speckle flow at \(42\) cm/s |
| Ultrafast echocardiography | In vitro spinning disk | Standard compounding gCNR fell from \(\approx 0.8\) at \(0\) cm/s to \(<0.1\) at \(15\) cm/s |
| Ultrafast echocardiography | In vitro EMC² | gCNR \(\gtrsim 0.8\) up to \(24\) cm/s in fundamental and \(17\) cm/s in second harmonic |
| Ultrafast echocardiography | Doppler estimation | NRMSE \(<12\%\) for \(V\leq 21\) cm/s in fundamental and \(V\leq 17\) cm/s in second harmonic |
| Ultrafast echocardiography | In vivo ventricular filling | Fundamental gCNR \(0.47 \to 0.79\); second-harmonic gCNR \(0.53 \to 0.87\) |
| Ultrafast echocardiography | In vivo diastasis | Fundamental gCNR \(0.67 \to 0.83\); second-harmonic gCNR \(0.80 \to 0.92\) |

In the READI paper, EMC² is also described qualitatively as fully recovering images corrupted by probe motion and restoring tissue speckle and sharpness to an image of a beating heart. In the echocardiography paper, the unified framework is reported to increase gCNR from \(0.47\) to \(0.87\) during ventricular filling when compared against coherent compounding, and the conclusion states that real-time implementation at \(32\)-angle compounding and \(140\) Hz yields high-contrast B-mode images and reliable tissue Doppler quantification even during fast myocardial motion up to \(17\) cm/s [2509.08781] [2312.01628].

## 6. Computational profile, methodological implications, and recurring misconceptions

The computational cost of READI-based EMC² can be decomposed into partial decoding, Hilbert transforms, beamforming, NCC block matching, and warping with coherent compounding. The paper specifies: \(S\) matrix multiplies of size \(Q\times Q\) using cublasSgemm, with complexity \(O(S\cdot Q^3)\); Hilbert transforms costing \(S\cdot M\cdot (N_{\mathrm{samp}}\log N_{\mathrm{samp}})\) per channel; beamforming over \(S\cdot Q\) subsets with \(O(S\cdot Q\cdot N_{\mathrm{pix}})\) delay-and-sum operations; NCC block matching over \((S-1)\) targets times \((N_{\mathrm{grid}}\times \mathrm{search\_area}\times \mathrm{template\_area})\); and warping plus compounding with \(O(S\cdot N_{\mathrm{pix}})\) interpolations and additions. On a modern NVIDIA GPU using CUDA, cuBLAS, cuFFT, and NPP, real-time performance greater than \(20\) Hz was reported for \(128\times128\) TOBE arrays and \(128\) transmit events [2509.08781].

In ultrafast echocardiography, the real-time pipeline is summarized at the algorithmic level rather than through asymptotic complexity terms. The implemented parameters are a \(75^\circ\) sector, \(N=32\) diverging waves, angular step \(\Delta\theta=1.25^\circ\), transmit PRF \(=4700\) Hz, frame rate \(\mathrm{PRF}/N\simeq147\) Hz, demodulation to IQ at \(f_0=2.98\) MHz for the fundamental or \(3.78\) MHz for the second harmonic, and grid \(\Delta r=\lambda\), \(\Delta\theta=0.4^\circ\). The processing sequence consists of beamforming each diverging wave into IQ lines, computing the autocorrelations \(R_1\) and \(R_2\), deriving \(\Phi(r)\), estimating \(V_{\mathrm{dop}}\), warping and phase compensating each angle, summing the motion-corrected fields, computing the angular coherence factor, and applying envelope detection with log compression [2312.01628].

Several recurring misconceptions can be addressed directly from these formulations. First, EMC² is not synonymous with a particular motion estimator: one implementation uses NCC block matching on low-resolution sub-images, whereas another uses ensemble autocorrelation and Doppler phase. Second, EMC² is not always limited to geometric realignment alone: in ultrafast echocardiography, the final reported EMC² image includes angular-coherence weighting, and the framework supports both fundamental and second-harmonic imaging. Third, motion compensation is not applied after irreversible magnitude formation in the READI setting; the low-resolution images are retained in complex analytic form specifically so that coherent compounding can preserve phase. A plausible implication is that EMC² is best understood as a phase-aware alignment-and-summation paradigm whose performance depends on the reliability of motion estimation in the relevant signal domain.

Source: https://www.emergentmind.com/topics/estimated-motion-compensated-compounding-emc2