Title Abstract Introduction Related Work Methodology Algorithm 1 Inference Experiments Results Conclusions
BMVC 2026 · Submission #131

Composed Historical Image
Retrieval by Modeling
Temporal Representations

Can a neural embedding space preserve temporal structure while remaining effective for both date estimation and object retrieval — simultaneously?

Method
TDIR
Dataset
DEW · 7.2M detections
Best Date P@10
86.3% (ViT)
Zero-label mode
Image + Image
Explore the paper
Abstract

Orthogonal Subspaces for Date and Content

While time evolves linearly, neural embeddings are chaotic and multi-dimensional. TDIR asks: can we decompose a single embedding into independent temporal and categorical subspaces — so the year of one image can be surgically transplanted onto another, label-free?

Disentangle the space — drag the slider or press play

Entangled TDIR
Car Jacket Poster point size encodes the year (small=1930s, large=1990s)

We propose Temporally Decomposable Image Representations (TDIR) — the first formal treatment of temporal decomposability in visual embeddings. We prove conditions under which decomposition holds, that orthogonality between subspaces emerges from joint optimisation without being explicitly imposed, and characterise the error when conditions are only partially met.

Beyond geometry, TDIR enables transitive operations: the year of one image can be extracted and injected into another with no label supervision — enabling the first zero-label composed historical image retrieval.

Chapter 1 · Introduction

When an Archivist Queries by Composition

Around 5–10% of archival material is photographic. Unlike text, photographs encode temporal information implicitly — in the design of a car, the cut of a jacket, the typography of a poster. Archivists infer dates from objects; TDIR makes this machine-readable.

Probative Value — objects as temporal vessels

Not all objects carry date information equally

In archival science, the probative value of an image is its capacity to uniquely illustrate a phenomenon. Objects evolve at different temporal rhythms. A tire or a pair of shorts are strong temporal signals; a building or a wheel are ambiguous. TDIR learns which objects are reliable "time capsules."

Temporal reliability — Average Year MAE by object (lower = more reliable)
The Three Query Modes

How should an ideal archive search work?

TDIR supports three retrieval modes of increasing expressiveness. Select each to see the geometric operation on the embedding space.

Mode A · q = µc + µy
Label + Label
Mode B · q = f(xc) + µy
Image + Label
Mode C · q = f(xc) + ry
Image + Image
Car col.Jacket col.Poster col. faded glyph = older year, vivid = newer gold ring = retrieved (inside radius)
Chapter 3 · Methodology

Building the Mathematical Foundation

TDIR rests on two definitions, two propositions, a lemma, a theorem, and its corollary — each with a geometric illustration inline. No appendix needed.

§3.1 · Definition 1

A representation is temporally independent from object classes when:

p(µc, c, µy, y) = p(µc, c) · p(µy, y)

Category centroids carry no date information; date centroids carry no category information.

Independence, illustrated — pick a category, then drag the year through it
§3.1 · Definition 2

The year subspace is temporally transitive if for any triple (A, B, C):

∆yA→B + ∆yB→C = ∆yA→C   where  ∆yA→B = µyB − µyA

Necessary for zero-label inference: a year residual from any image acts as a consistent temporal operator when transplanted.

Drag the year proxies — watch the displacement identity hold
Drag µyB — the displacement ΔA→C = ΔA→B + ΔB→C updates live.
§3.2 · Proxy Sets

Three learned proxy sets

The ideal decomposition v = f(x) = µc + µy requires three proxy sets, each with a distinct role at training and inference time.

𝒫c — CATEGORY
µc1 … µcC
One per object category. Integrates over all years → year-agnostic centroid. Operates on unnormalised embeddings.
𝒫y — YEAR
µy1 … µyY
One per year, shared across categories. Operates on ℓ₂-normalised residuals — temporal info purely in direction.
𝒦y — DISPLACEMENT
ky1 … kyY
Learnable temporal offsets for the swapping trick. Used only at training. Discarded at inference.
§3.2 · Proposition 1

Under Definition 1:

I(µc, µy; c, y) = I(µc; c) + I(µy; y)

The total information in the joint embedding decomposes into two independent sums — no cross-information between categorical and temporal components.

The proof is a two-line cancellation. Step through it — read the diagram as an information map: the total area is everything the proxies tell us about (category, year). Watch the cross-terms vanish, leaving two disjoint pools.

START
Joint info
CHAIN ×1
Peel off µᶜ
DEF 1
Kill cross-terms
RESULT
Disjoint sum

So v = f(x) is temporally decomposable: knowing both proxies captures everything v encodes about (c, y), with no shared (overlapping) information left between the temporal and categorical parts.

§3.2 · Lemma 1

Under Definition 1, the category and temporal centroids are orthogonal:

Cov(µc, µy) = 0    Var(µc + µy) = Var(µc) + Var(µy)

Independence → zero covariance → Pythagorean variance decomposition.

The two subspaces are perpendicular — drag the noise slider to see the effect of residual entanglement

⟨µc, µy⟩ = 0.00

As I(ε;c,y) grows, µc picks up a projection onto µy — but Theorem 1 shows joint optimisation drives this projection back to zero.

§3.2 · Proposition 2

In real data, residual noise ε entangles year and category. But:

I(v; c,y) ≤ I(µc; c) + I(µy; y) + δ   where  δ ≜ I(ε; c,y)

δ is a data-dependent constant — it does not depend on the proxy vectors µc, µy. The orthogonality constraint ⟨µc, µy⟩ = 0 can be satisfied without affecting this bound.

The error floor δ ignores how you orient the proxies — drag both knobs

An embedding splits as v = µc + µy + ε. You are free to rotate the orthogonal (µc, µy) frame to satisfy ⟨µc,µy⟩ = 0 — but that only re-shuffles information between the two proxies. The residual δ = I(ε; c,y) is fixed by the data alone, so the proxies can be learned freely without touching the bound.

PROXY ORIENTATION θ
DATA ENTANGLEMENT I(ε;c,y)
δ ≤ I(ε; c) + I(ε; y)

When ε is disentangled: δ → 0
⟹ I(v; c,y) ≈ I(µc;c) + I(µy;y)
Circularity resolved by Theorem 1: Proposition 2 assumes orthogonality to guarantee free learning — but orthogonality hasn't been enforced yet. Theorem 1 shows it emerges implicitly from the training strategy.
§3.4 · Theorem 1 (Emergent TDIR)

Part I (exact). If I(ε;c,y) = 0:

minθ ℒc + ℒy ⟹ p(µc, µy, c, y) = p(µc, c) · p(µy, y)

Part II (asymptotic). If I(ε;c,y) ≠ 0:

minθ ℒc + ℒy ⟹ ⟨µc, µy⟩ → O(I(ε;c,y) / ‖µy‖²) ∀c,y

Orthogonality is an emergent property of the joint optimisation — it does not need to be imposed explicitly.

Interactive geometric proof (Part II) — drag the category proxy µc

Decompose µc = µ̃c + αµyA + βµyB, where (α,β) is its projection onto the year subspace. Drag µc: watch the two year residuals rA, rB drift away from the year proxies they are supposed to match, and the year loss climb. Only at α=β=0 do the residuals land exactly on µyA, µyB — that is the unique minimum.

LEMMA · YEAR PROXIES: click to force them parallel and watch the unique minimum collapse
GEOMETRY — drag µᶜ
YEAR-LOSS ENERGY Q(α,β)

The right panel is the actual loss surface Q(α,β) = 2α²‖µyA‖² + 2β²‖µyB‖² + 4αβ⟨µyA,µyB⟩ = 2‖α µyA + β µyB‖². It is a positive-definite bowl with a single minimum at α=β=0 — but only while the year proxies are linearly independent (the Lemma). Force them parallel and the bowl flattens into a valley: the minimum is no longer unique. Double-click the left panel to reset.

§3.4 · Corollary 1
Projµy(µc) → 0    ∀ c, y

The shared displacement vectors ky cannot simultaneously compensate projections α ≠ α' across two categories — since ε is not a control variable, the unique globally consistent solution is α = α' = 0 for all c, c'.

The swapping trick forces residuals from different categories to be displaced by the same ky. This makes the optimizer implicitly require all category proxies to be orthogonal to all year proxies — not just category by category. It's the key to scalability.

What is α? Write each category centroid as µc = µ̃c + α·µyA + β·µyB. The coefficient α is simply how far µc leans into the year axis — the length of its shadow on µy. Orthogonality means α = 0 (no shadow). Below, α and α' are drawn as the coloured offsets that pull each residual off the year proxy.

Why one shared ky forces it globally — adjust each category's projection

Two categories (blue, red) must both have their year-A residual mapped onto the same target µyB by the same displacement kyB. If their projections α, α' differ, a single kyB cannot hit the target for both. The only globally consistent solution is α = α' = 0.

CATEGORY c (blue) · α
CATEGORY c' (red) · α'

Set both sliders to 0 to see both residuals land exactly on µyB with the shared displacement — the consistent optimum.

Chapter 3 · Algorithm 1

The TDIR Training Loop

Click each step — or press Play — to step through Algorithm 1. The canvas shows the geometric operation on the embedding space at each point. Color-coded to match the pseudocode.

■ Forward / Backward ■ Losses ■ Residuals ■ Normalise
Init
Forward
ℒc
Center
Swap
Norm
ℒy direct
ℒy trans
Backward
Chapter 3 · Inference

Three Modes of Composed Retrieval

At inference, 𝒦y is discarded. The three retrieval modes follow from f(x) ≈ µc + µy. Each is a different geometric operation on the learned embedding space.

q = µc + µy
Label + Label
q = f(xc) + µy
Image + Label
q = f(xc) + ry
Image + Image
Car col.Jacket col.Poster col. faded glyph = older year, vivid = newer gold ring = retrieved (inside radius)
Chapter 4 · Experimental Setup

Dataset, Architecture, Protocol

7.2M
Object Detections
586
Object Categories
50
Top Categories (focus)
70yr
1930–1999
Dataset Construction

Object-Centric DEW

Starting from DEW (1M natural-scene photographs, 1930–1999, one year per image), we apply OWLv2 with a 30% threshold — yielding 7.2M detections across 586 classes with year-level annotation. Top-50 categories used for evaluation (4.4M detections).

DEW object detections by decade — click a bar
Evaluation Protocol

Two precision axes — object and date

For each inference mode, retrieve top-K=10 images for every query. A date retrieval is correct if within the same lustrum (5-year window). Random baseline for date: 7.14% (70yr / 14 bins).

Object Precision@K:
  #{xᵢ : cᵢ = c} / K

Date Precision@K:
  #{xᵢ : yᵢ ≈ y} / K
  (y ≈ yᵢ iff same 5-year lustrum)
Random date baseline: 7.14%

Best Label+Label: Date 86.3% (ViT)
Best Image+Image: Obj 72.5% (ConvNeXt)
Chapter 5 · Results

Evaluation: Backbones, Ablations, Calibration

Table 2 · Backbone Comparison (K=10)

All TDIR architectures dominate CLIP on date precision

CLIP is competitive on objects but fails temporally — date composability requires TDIR's explicit geometric structure.
BackboneLabel+LabelImage+LabelImage+Image
Obj P@KDate P@KObj P@KDate P@KObj P@KDate P@K
CLIP (arithmetic).152.115.387.076.389.102
CLIP (prompting).415.105.372.079.520.081
VGG + TDIR.397.712.388.579.669.234
ResNet + TDIR.489.803.444.665.721.306
ConvNeXt + TDIR.443.824.416.664.725.393
ViT + TDIR.413.863.355.708.646.429
Table 3 · Ablation Study (ConvNeXt, K=10)

All three components are necessary

The paper ablates exactly one component at a time, so this is a single choice: pick what to remove and the bars compare it against the full model (faded bars). Notably, removing the swapping trick hurts Label+Label too — confirming Corollary 1 (implicit orthogonality regularisation).

Figure 4 · Calibration Analysis

Category preserved; date degrades with temporal gap

Object precision remains stable across source years, target years, and gaps — the category subspace is insulated from temporal injection. Date precision follows a U-shaped curve over target year and degrades monotonically with gap.

Chapter 6 · Conclusions

What TDIR Establishes

THEORETICAL
First formal treatment of temporal decomposability

Propositions 1 & 2 bound the decomposition error. Theorem 1 proves orthogonality emerges from joint optimisation. Corollary 1 shows the swapping trick enforces this cross-category.

EMPIRICAL
First zero-label composed historical retrieval

DEW extended with 7.2M detections across 50 categories. Image+Image mode transplants the year of one image onto the content of another — no labels required.

LIMITATIONS
Shared displacement; additive operator

A single temporal displacement across all categories may not reflect uneven historical evolution. Additive transplantation loses fidelity at large temporal distances.

FUTURE WORK
Non-linear operators; category-specific subspaces

Replace additive operator with a non-linear mechanism. Explore category-specific temporal subspaces. Extend to zero-shot evaluation over unseen categories.

The geometry of a well-structured embedding is not merely a tool for retrieval — it is a formal language for historical reasoning. Every µc + µy is a sentence: "Show me this object from this era."