Color Tracer Mapping: Methods & Applications

• Color Tracer Mapping is a set of methods that transform color signals into structured representations for identifying coherent sets, inferring hidden variables, and supporting analysis.
• Techniques span spectral graph clustering for unsteady flows, physics-constrained reconstruction in cardiac imaging, and similarity-preserving encodings in visualization.
• Applications include remote sensing, LiDAR mapping, and image editing where color-derived data is mapped into domains with enhanced inferential or visual capabilities.

Color Tracer Mapping denotes a family of technically distinct mapping problems in which color, colored tracers, or color-derived signals are transformed into more structured representations in order to identify coherence, reconstruct hidden variables, encode multivariate relations, or support analysis and visualization. In the literature, the phrase is associated with reference-trajectory coherent-set recovery in unsteady flows, physics-constrained reconstruction of intraventricular blood-flow vectors from color Doppler, similarity-preserving multivariate color encodings, mappings between color images and vector fields, lithologic mapping from multispectral color ratios, color-assisted 3D mapping, and RGB-conditioned image editing (Schlueter-Kuck et al., 2017, Vixège et al., 2021, Cheng et al., 2016). Taken together, these usages indicate that the term does not denote a single standardized algorithm; rather, it names a recurring strategy in which color or tracer information is treated as a structured signal rather than as a purely decorative display variable.

1. Scope and conceptual variants

Across the cited work, Color Tracer Mapping appears in at least four recurrent forms. In dynamical systems, it refers to trajectory-centered identification of the set of particles that move coherently with a chosen reference tracer, using spectral graph methods rather than Euclidean proximity in physical space (Schlueter-Kuck et al., 2017). In medical flow imaging, it denotes reconstruction of a full planar velocity field from the one-component scalar measurements returned by color Doppler, under explicit physical constraints (Vixège et al., 2021). In visualization, it refers to data-driven mappings from multivariate samples, vector fields, or compositions into perceptually organized color spaces, so that similarity, direction, or mixture can be read directly from color (Cheng et al., 2016, Waters et al., 2020, Lu et al., 2023). In image analysis and mapping systems, it can denote coordinate-to-color neural encodings, forward and backward mappings between color images and vector fields, trace-transform descriptors for image-domain identification, colorized LiDAR map construction, and RGB-to-embedding control modules for diffusion editing (Bricman et al., 2018, Snarskii et al., 2022, Olaizola et al., 2012, Lu et al., 24 Feb 2025, Yang et al., 17 Sep 2025).

A common misconception is that these works all solve the same task. They do not. Some infer hidden kinematics from sparse tracers, some convert scalar or categorical measurements into vector or compositional fields, and others use color as an auxiliary constraint for registration or editing. The unifying principle is narrower: each method maps color-associated information into a domain in which structure is more separable, more constrained, or more interpretable.

Another misconception is that color is merely a visualization layer added after the analysis. Several of the cited methods treat color as part of the inference problem itself. In intraventricular vector flow mapping, the color Doppler field is the observation from which the missing velocity component must be recovered (Vixège et al., 2021). In CAR-LOAM, color differences enter the scan-to-map optimization as robust weights (Lu et al., 24 Feb 2025). In multivariate visualization, perceptual color distance and color-name consistency are explicit optimization terms rather than aesthetic afterthoughts (Waters et al., 2020, Lu et al., 2023).

2. Reference-centered coherent-set identification in unsteady flows

In the trajectory-centered formulation based on Coherent Structure Coloring (CSC), the input is a set of Lagrangian particle trajectories, potentially sparse, and the output is not a full segmentation of the flow but the coherent set associated with one selected reference tracer (Schlueter-Kuck et al., 2017). The method begins by constructing a graph whose weighted adjacency matrix encodes pairwise trajectory dissimilarity. The graph Laplacian is

$L=D-A,$

and the generalized eigenvalue problem

$LX=\lambda DX$

produces eigenvectors $X_1,X_2,\dots$, with the largest eigenvalue/eigenvector pair $X_1,\lambda_1$ defining the original CSC field. Each entry of $X_1$ is the CSC value assigned to a particle at the final time, and it highlights the strongest trajectory-to-trajectory dissimilarities in the flow (Schlueter-Kuck et al., 2017).

The crucial extension is that the method does not stop at $X_1$. Additional eigenvectors are used because different coherent sets can be indistinguishable in $X_1$ alone. Since $A$ is real and symmetric, the eigenvectors are orthogonal, and the dominant generalized eigenvectors define a spectral embedding in which kinematically similar trajectories remain close even when they are far apart in physical space. For a chosen reference trajectory $\mathrm{ref}$, the simplest one-dimensional distance is

$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$

and in a $LX=\lambda DX$0-dimensional eigenspace the distance becomes

$LX=\lambda DX$1

The paper also proposes a weight choice that emphasizes larger eigenvalues, although in the examples the unweighted metric $LX=\lambda DX$2 works well (Schlueter-Kuck et al., 2017).

The number of eigenvectors is chosen by identifying a critical eigenspace dimension from the plateau of

$LX=\lambda DX$3

This is analogous to an elbow criterion in a scree plot. Once a suitable dimension has been selected, the method applies agglomerative hierarchical clustering to separate the trajectories into two groups: those inside the coherent structure containing the reference particle, and those outside. The cluster containing the reference trajectory just before the last merge is taken as the coherent set associated with that trajectory (Schlueter-Kuck et al., 2017).

The method was demonstrated on the quadruple gyre and the Bickley jet. In the quadruple gyre, a reference particle in the upper-left vortex core yielded a compact coherent region using about 6 dimensions, and the recovered set matched the gyre core seen in dense-particle FTLE calculations. In the Bickley jet, the method identified both a vortex core and the elongated jet itself as coherent structures, including the nonconvex jet that spans much of the domain. A key limitation was also made explicit: when a background-flow particle is chosen as the reference, there may be no clear eigenspace plateau, and clustering can absorb nearly all particles into a single set. The paper therefore frames the method as reference-trajectory based rather than as a guarantee that every trajectory belongs to a meaningful coherent set (Schlueter-Kuck et al., 2017).

The same work uses the shape of the $LX=\lambda DX$4 curve as a diagnostic of relative coherence. A clear elbow suggests a genuine coherent set, whereas the absence of a plateau suggests incoherent background flow. The authors also note that true coherence should be commutative: if particle $LX=\lambda DX$5 is coherent with $LX=\lambda DX$6, then $LX=\lambda DX$7 should be coherent with $LX=\lambda DX$8. The background-flow example fails this intuition, reinforcing that it is not a true coherent set (Schlueter-Kuck et al., 2017). This directly addresses a common misunderstanding of trajectory-centered mapping: proximity in spectral tracer space is not synonymous with universal set membership.

3. Physics-constrained reconstruction from color Doppler

In cardiac flow imaging, Color Tracer Mapping appears as physics-constrained intraventricular vector flow mapping (iVFM), where the objective is to reconstruct the missing tangential velocity component from 2-D fan-shaped color Doppler measurements (Vixège et al., 2021). Color Doppler is a one-component velocimetric technique. If the actual velocity is $LX=\lambda DX$9, the measured signal is

$X_1,X_2,\dots$0

so the inverse problem is to recover $X_1,X_2,\dots$1 and thereby the full planar field $X_1,X_2,\dots$2 (Vixège et al., 2021).

The reconstruction is formulated in polar coordinates centered on the probe. The method seeks $X_1,X_2,\dots$3 that minimize the mismatch to Doppler data,

$X_1,X_2,\dots$4

subject to two equality constraints. The first is planar mass conservation,

$X_1,X_2,\dots$5

and the second is the free-slip boundary condition,

$X_1,X_2,\dots$6

The first enforces a divergence-free 2-D velocity field, and the second imposes no normal penetration through the endocardium while allowing tangential slip, which the paper argues is appropriate because Doppler resolution is too coarse to resolve the boundary layer (Vixège et al., 2021).

The constrained optimization is solved with Lagrange multipliers. The discretization is performed on the polar Doppler grid before scan conversion, with constant radial and angular steps $X_1,X_2,\dots$7, and differential operators approximated by second-order central finite differences using three-point stencils. The resulting stationarity conditions produce a real, sparse, symmetric linear system. After removing null rows and columns due to incomplete ROI coverage, the system becomes full-rank and positive-definite, and is solved by Cholesky decomposition in MATLAB. Because the Doppler data and wall velocities are noisy, the method adds a smoothing penalty,

$X_1,X_2,\dots$8

leading to

$X_1,X_2,\dots$9

The paper emphasizes that this version uses one regularization parameter $X_1,\lambda_1$0, selected automatically by the L-curve method and then reused for the remainder of the cardiac cycle after being determined once at the end of early filling (Vixège et al., 2021).

Validation used a patient-specific CFD model with geometry from CT-based cardiac anatomy, an ALE method for moving walls, a simulated apical three-chamber Doppler sector, 50 scanlines and 160 samples per scanline, radial spacing about 0.61 mm and angular spacing $X_1,\lambda_1$1, 100 color Doppler frames across a cardiac cycle, and Gaussian noise added to simulate SNRs from 10 to 50 dB. The iVFM-derived velocities agreed well with CFD: for the radial component $X_1,\lambda_1$2, the coefficients of determination were $X_1,\lambda_1$3 with normalized RMSE $X_1,\lambda_1$4–$X_1,\lambda_1$5; for the angular component $X_1,\lambda_1$6, the coefficients of determination were $X_1,\lambda_1$7 with normalized RMSE $X_1,\lambda_1$8–$X_1,\lambda_1$9. The method also preserved large-scale flow organization: mean vorticity agreed with CFD at $X_1$0 and mean difference $X_1$1, while the stream function agreed at $X_1$2 (Vixège et al., 2021).

In vivo, the method revealed intraventricular flow structures that are hardly visible in raw color Doppler. In a normal patient, it clearly showed the large vortex forming during early filling, and the vortex persisted into diastasis. The paper suggests that this vortex behavior could become a biomarker of diastolic function and concludes that the algorithm is ready for pilot clinical studies (Vixège et al., 2021). The broader significance is methodological: color is the measurement domain, but the clinically useful output is a reconstructed vector field constrained by fluid mechanics.

4. Similarity-preserving color encodings and perceptual uniformity

A distinct branch of Color Tracer Mapping concerns the design of color encodings that preserve relationships among variables, samples, classes, or vector directions. In the multivariate case, a data-driven method first gauges the similarity of the attributes, arranges them according to the periphery of a convex 2D color space such as HSL, and then assigns each multivariate sample a color via generalized barycentric coordinate (GBC) interpolation (Cheng et al., 2016). The method fixes lightness to a constant, typically

$X_1$3

and places the variables on the periphery of the HSL slice so that similar attributes are nearby and dissimilar attributes are separated. For a point $X_1$4 inside the convex polygon with vertices $X_1$5, the generalized barycentric weights satisfy

$X_1$6

and

$X_1$7

The same weights define the color interpolation,

$X_1$8

The paper also describes three feature-extraction or warping modes based on an elliptical distribution region identified by PCA: color preserving enhancement, contrast enhancement coloring, and comparison compression coloring (Cheng et al., 2016).

Perceptual colorimetry generalizes this concern from multivariate composition to scalar fields, vector fields, and three-component composition fields. A central claim is that color is not merely decorative in scientific figures but a data-encoding channel, and that mappings should therefore be constructed in perceptually uniform color spaces such as CIELAB, CAM02-UCS, or CAM16-UCS rather than raw sRGB or HSV (Waters et al., 2020). In CIELAB, perceptual distance is measured by

$X_1$9

and perceptual uniformity is expressed through a constant derivative $X_1$0. The paper notes that a typical perceptually detectable difference is about $X_1$1. For 2D vectors, PAPUC maps vector angle to hue angle in the $X_1$2 plane and vector magnitude to lightness $X_1$3, using an inverted cone in CIELAB with axis along $X_1$4, tip at $X_1$5, base at $X_1$6, and radius 40. For three-component composition fields, CMPUC maps total intensity to $X_1$7 and relative composition to an equilateral triangle in the $X_1$8 plane, with the ideal triangle chosen as large as possible while remaining inside the sRGB gamut: $X_1$9, circumradius $X_1$0, and hue angle $X_1$1 (Waters et al., 2020).

Interactive context-preserving highlighting extends these ideas to dynamic categorical visualization. The method precomputes two coordinated palettes for the same classes: a salient palette $X_1$2 used for highlighted points and a faint palette $X_1$3 used for de-emphasized points (Lu et al., 2023). The stated design requirements are: highlight selected data points as much as possible; maintain visual discrimination among all classes; and preserve color consistency during dynamic emphasis changes. Corresponding colors across the two palettes are constrained to preserve hue and saturation,

$X_1$4

while differing mainly in lightness. The faint palette is further constrained by

$X_1$5

with $X_1$6 after experimentation. Optimization uses a customized simulated annealing algorithm. Evaluation consisted of two crowd-sourced experiments on Amazon Mechanical Turk with 150 participants total, plus two case studies. The paper reports that the method significantly outperformed most benchmark approaches in the highlighting task and supported lower error and response time in the matching task than Tableau Highlighter variants (Lu et al., 2023).

This literature directly counters the idea that “good” color mapping is a subjective matter. In all three cases, similarity structure, perceptual distance, or semantic consistency enters the mathematical construction of the map itself (Cheng et al., 2016, Waters et al., 2020, Lu et al., 2023).

5. Image-to-structure mappings: implicit functions, vector fields, and domain descriptors

Another major usage of Color Tracer Mapping converts images into compact functions, vector fields, or discriminative descriptors. CocoNet treats a single image as a continuous coordinate-to-color function. Rather than learning from many labeled images, it learns a function

$X_1$7

that maps a six-dimensional coordinate representation

$X_1$8

to the normalized RGB triplet $X_1$9 (Bricman et al., 2018). The network is a fully connected feedforward neural network with 6 input neurons, tanh hidden activations, 3 sigmoid output neurons, and task-dependent depth and width. It is trained with MSE and Adam, with learning rates between $A$0 and $A$1. On CIFAR-10 denoising with Gaussian noise variance 10, CocoNet achieved 31.25 PSNR and 0.9503 SSIM; with variance 20, it achieved 27.22 PSNR and 0.8925 SSIM. On Set5 $A$2 upsampling, it outperformed bicubic interpolation on some images in PSNR and SSIM, and in completion experiments it reconstructed global shape structure well but struggled with fine texture details (Bricman et al., 2018). The method is explicit about its main limitation: one network must be trained individually for each input image.

A different image-to-structure construction is the fiber bundle color space mapping between a color image and a 2D vector field (Snarskii et al., 2022). The paper defines a color vector

$A$3

with

$A$4

so that brightness corresponds to

$A$5

Any color vector is decomposed into a spectral component and an achromatic component directed along the achromatic axis $A$6. The forward mapping orthogonally projects the spectral part onto a 2D plane; the achromatic component has zero projection and is preserved separately as a scalar field. The backward mapping reconstructs a color image from a 2D vector field under the simplifying assumption $A$7 and the constraint $A$8, which ensures a unique solution. Sequential application of forward and backward mapping defines an adaptive image filter that removes the achromatic component completely, or partially if the inverse mapping is modified. The mapped vector field then supports divergence, curl, gradients, and related vector-analysis tools; the paper illustrates this with a still-life image and with color rendering of a 2D magnetic vector field scattered by three corrosion pits in a pipeline defect scenario (Snarskii et al., 2022).

For image-domain identification, DITEC uses the trace transform as a global descriptor for color image context categorization (Olaizola et al., 2012). The method applies the trace transform to each color channel, computes a global 2D DCT of the transformed signal, and summarizes groups of DCT coefficients by the mean $A$9 and kurtosis

$\mathrm{ref}$0

or, for a discrete set,

$\mathrm{ref}$1

The first descriptor pair for the DC term is replaced by the mean and variance of the original image in HSV space. The paper gives a numerical example in which $\mathrm{ref}$2 coefficients are reduced to $\mathrm{ref}$3 using mean and kurtosis, and then to $\mathrm{ref}$4 after channel-wise pruning. Feature ranking uses the square of the SVM weight, and the best SVM-ranking results are reported as about $\mathrm{ref}$5 higher than PCA. On the reported experiments, precision was typically $\mathrm{ref}$6–$\mathrm{ref}$7 on the Corel dataset and above $\mathrm{ref}$8 on the satellite dataset (Olaizola et al., 2012). This is a mapping from color images to low-dimensional global descriptors rather than to displayed colors.

These three cases share a structural move: the image is re-expressed in a domain where continuity, vector calculus, or statistical discrimination becomes tractable. This suggests that “tracer” in Color Tracer Mapping can refer not only to physical particles but also to the propagation of color information through a representation pipeline.

6. Remote sensing, 3D mapping, and controllable color transformation

In planetary remote sensing, color parameters are used as compositional tracers. Dawn Framing Camera data on Vesta were used to compute three 1-$\mathrm{ref}$9m pyroxene band parameters:

$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$0

$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$1

$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$2

BC measures curvature around the 1-$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$3m absorption, BT measures the slope or tilt across the feature, and BS is a proxy for band depth (Thangjam et al., 2013). Laboratory calibration used 239 HED meteorite spectra from RELAB: 41 howardites, 157 eucrites, and 41 diogenites. The resulting parameter ranges were $d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$4–$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$5 for eucrite BC and $d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$6–$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$7 for diogenite BC; $d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$8–$d_1(i)=|CSC(i)-CSC(\mathrm{ref})|,$9 for eucrite BT and $LX=\lambda DX$00–$LX=\lambda DX$01 for diogenite BT; and $LX=\lambda DX$02–$LX=\lambda DX$03 for eucrite BS and $LX=\lambda DX$04–$LX=\lambda DX$05 for diogenite BS. The paper concludes that BC and BT are the most robust discriminants, whereas BS is less reliable because it is strongly affected by grain size and opaques. The resulting global FC maps showed eucrite-rich regions mainly in the equatorial and northern hemisphere, diogenite-rich regions concentrated in the southern hemisphere, and widespread howardite-rich mixed regolith. The average color spectrum of Vesta was found to be essentially identical to that of howardite-rich regions, supporting an extensive mixing of surface regolith due to impact gardening (Thangjam et al., 2013).

In robotics, CAR-LOAM uses camera color as an auxiliary signal for LiDAR odometry and mapping. LiDAR points are projected into the camera image and assigned RGB values, yielding colored features of the form

$LX=\lambda DX$06

Edge and planar features are matched to a global map, and the pose is estimated by minimizing a weighted sum of geometric residuals. Positional correspondence outliers are handled by the Welsch function

$LX=\lambda DX$07

with $LX=\lambda DX$08, while color correspondence outliers are suppressed by Gaussian weights built from the CIEDE2000 perceptual color difference,

$LX=\lambda DX$09

with $LX=\lambda DX$10 (Lu et al., 24 Feb 2025). The overall optimization is performed on $LX=\lambda DX$11 with Gauss-Newton. In a campus outdoor field test, the ablation study reported ATE reduced to $LX=\lambda DX$12 with both color weighting and Welsch loss, compared to $LX=\lambda DX$13 without both, corresponding to about $LX=\lambda DX$14 ATE reduction and about $LX=\lambda DX$15 RTE reduction. In building reconstruction, the paper reported 23% of nearest neighbors below 0.1 mm, 60% below 0.5 mm, and 72% below 1 mm. The method also notes an important limitation: it may degrade in featureless environments or in scenes with similar or missing color (Lu et al., 24 Feb 2025).

In diffusion-based image editing, color mapping becomes an explicit control mechanism. A color mapping module is trained to predict a CLIP text embedding from an RGB value so that a diffusion model can perform continuous color editing without relying on ambiguous prompt wording or uncalibrated interpolation coefficients (Yang et al., 17 Sep 2025). The method first interpolates between two color-related prompt embeddings,

$LX=\lambda DX$16

uses InstructPix2Pix to generate edited images, and samples an RGB value $LX=\lambda DX$17 from the edited region. The original embedding $LX=\lambda DX$18 of shape $LX=\lambda DX$19 is compressed by PCA to $LX=\lambda DX$20, with PCA dimension 15, and an MLP learns

$LX=\lambda DX$21

under the objective

$LX=\lambda DX$22

The predicted embedding is then reconstructed by inverse PCA and fed back into the diffusion model. Training uses 30 embedding–RGB pairs per image, Adam, learning rate 0.001, 500 epochs, and a single NVIDIA RTX 4090 GPU. A SAM-derived binary mask restricts the edit region. The paper reports smoother transitions, a more linear relationship in RGB distribution plots, and improved semantic consistency relative to interpolation-based control, while also noting dependencies on sampled training pairs, single-pixel RGB supervision, PCA compression, and segmentation quality (Yang et al., 17 Sep 2025).

These remote-sensing, robotic, and generative cases illustrate a broad but precise pattern: color-derived quantities can serve as tracers of composition, correspondence reliability, or intended appearance. A plausible implication is that Color Tracer Mapping is best understood not as a discipline-specific method, but as a recurrent computational motif in which color-associated measurements are mapped into a latent domain where inference, registration, classification, or control becomes better posed.