Colour Science & Perception
what the eye settles, and what it cannot
Concept

Quadrature — where it appears

Computing an integral by sampling it, with an error whose order follows from the geometry rather than from a fit to the data. Stating the order lets a computed result be checked against a refinement rather than believed, which is how every integral on this site is tested.

Named by 22 essays across 6 fields — each of them below, with the objects they name alongside it.

The form-factor matrix of the box. F_ij is the fraction of everything leaving face i that arrives at face j. The diagonal is zero because a flat face sees none of itself; each row sums to exactly 1 because the cavity is closed; and A_i F_ij = A_j F_ji, which is reciprocity and is checked to 10⁻⁹. For a cube the opposite face takes 19.98% and each of the four adjacent faces 20.00%, and the near-equality of those two numbers is a coincidence of the cube rather than a rule.

A corner is not a wall

A form factor is the fraction of everything leaving one surface that arrives at another, and it is the only place geometry enters the colour of a room. It is also the one number here with a published closed form to check against — and the check turned out to converge at two different rates for two cases that look identical.

scene · Scene What a press and a display can reach at L* = 50. A slice through both solids at lightness 50, with the press dashed and the display solid. The boundaries cross: the press reaches past sRGB in 15 of 48 directions and falls inside it in 31. That is the shape of every conversion between them — colours are lost in one direction and gained in the other, and the picture cannot show the gained ones, because it is being displayed on the gamut that cannot reach them. Those are hatched.

Neither gamut contains the other

Converting a picture to CMYK is described everywhere as a reduction, as though the press were a smaller version of the screen. Measured as solids, the press reaches eleven per cent of its volume outside sRGB while sRGB reaches fifty-four per cent of its volume outside the press — and in the cyans the press wins by seventy per cent of chroma.

applied · Delivery Every census row under five constructions of the same test set. A slope chart with 5 columns — lattice, coarse, fine, uniform, natural — and one line per change of light in the census, each line joining that row's mean residual under each construction. Four of the five columns describe the same region of surfaces walked at different densities or against different measures; the last is the clamped, realistic family, which is not linear in its parameters and is therefore answering a slightly different question. The levels move: between the coarse and fine lattices every row shifts by seven to nine per cent, in the same direction, which is a common-mode factor no published residual here has ever carried. The order almost survives. Inside the region exactly one pair crosses, and it is the pair the standard error had already flagged; under the clamped set two more cross, including one the error separates by nearly nine standard errors. The crossing lines are drawn heavy.

A lattice is a quadrature rule

Walking a set of test surfaces more finely does not converge on a better answer, because refining a lattice under a constraint changes which corners of the region get sampled and not only how densely. The lattice used here turns out to be a two per cent biased estimate of the integral it stands for.

difference · Metric Five fields, by how much light arrives from each elevation. The radiance arriving at a surface from each direction in one vertical plane, for five ways of lighting it. The vertical axis is logarithmic, spanning the three decades between a sun and the sky around it. The number beside each name is the share of the light that would have to be moved to make the field uniform: zero for the overcast sky, 0.93 for a lamp on a stand. A uniform field is the condition under which a reading is the sample's own reflectance, and the only place it exists is inside an instrument.

A room is not a sphere

An integrating sphere reads a surface's own reflectance exactly, and the exactness is a theorem rather than good engineering — under a hemisphere of constant radiance, reciprocity makes the reading the sample's directional-hemispherical reflectance whatever the surface is. Every room fails that condition, and a viewing booth and a window get the sign of the error wrong in opposite directions.

scene · Scene Six functions of wavelength, and the six different places they stop. Every table this collection integrates against, drawn over the range the body that published it defined it on. The scale is logarithmic so that the ultraviolet and the near infrared both fit. The bottom row is the range used here before the infrared band was added, and it is the intersection of the two rows that matter for an eye looking at a reflector — which is the right answer only while everything in the integral is being multiplied together. The daylight basis runs 80 nanometres further down than that intersection, and it was published that way because the ultraviolet in daylight is what makes a brightened sheet of paper glow. The analytic row is drawn to the edge of the plot because it has no edge: Planck's law is a formula and is exact at every wavelength, which is why illuminant A needs no table at all.

The grid is a range, not an index

Three of the four departures in this round restore an argument the model dropped. The fourth does not — the 380-to-780-nanometre grid is the range of the one argument the model kept, chosen in the first weeks and never revisited. It costs 6.70 ΔE₀₀ on a coated printing paper, it is the cheapest of the four to fix, and it is the one still unfixed.

light · Light The six arguments a surface's response has, and the one this model keeps. A surface's response to light is a function of six arguments: the wavelength, direction and place the light arrives with, and the wavelength, direction and place it leaves with. The model every colour here is computed from keeps one number per wavelength, which means it takes the diagonal of the first pair, integrates the second away, and assumes the third pair equal. Each departure drawn here restores one of them. The fourth departure is not on the diagram: the wavelength grid is the range of the index that was kept rather than an index that was dropped, which is why it is the cheapest of the four to fix and was still not fixed.

The audit that changed the object

Three rounds audited this collection's numbers, its sets and its conventions, and each one ended by naming the same thing it could not reach — that a surface is a reflectance. This round reached it, found that four departures share one algebraic form, had two of its predictions refused, and discovered that its own hemisphere quadrature had been passing a convergence check by luck.

limits · Limits a fluorescent tube, mercury lines on a phosphor bed, with a 5-nanometre grid marked on it. The light drawn at a fifth of a nanometre, with the 5-nanometre tabulation points marked beneath. 38 of the 53 points carry more than a twentieth of the peak. What a summation over those points computes is not an approximation to the area under this curve; on a spectrum with features narrower than the spacing it is a different quantity, and the difference depends on where the points fall rather than on how many there are.

The index is a choice too

Every colour in this collection is a sum over eighty-one numbers running from 380 to 780 nanometres in steps of five. That index is not a property of the eye, the light or the sample — it is a tabulation, and it holds three separable decisions that behave completely differently from one another.

light · Light What an instrument's slit width does to a tabulated colour. The horizontal axis is the full width of a triangular slit, from zero — perfect point sampling — to twenty nanometres; the vertical is the distance from the true colour, logarithmic. For a smooth light the lines are flat: a slit narrower than any feature changes nothing. For a line spectrum they fall off a cliff at the left. Point-sampling a mercury line at five nanometres costs 1.26 ΔE₀₀ and integrating the same spectrum through a five-nanometre slit costs 0.014. A spectrometer does not sample a spectrum; it integrates one, and the blur everybody would remove if they could is what makes a five-nanometre table safe.

The slit is what makes it legal

Point-sampling a mercury line at five nanometres costs a colour difference of one unit, and a laser projector thirty-seven. Integrating the same spectrum through the five-nanometre slit every spectrometer already has costs 0.014 and 0.19. The blur anybody would remove if they could is what makes a coarse table honest.

light · Light What three fill-in rules cost when a five-nanometre table is read at one. A five-nanometre table read at one nanometre by three rules — hold the value, straight lines, a cubic through four points — each compared with the same tenth-nanometre reference. The dashed rule is the answer obtained by summing the table as it stands, at 0.00000 ΔE₀₀. Two of the three interpolations are worse than not interpolating. That is not a paradox: the summation's error is already small because the normaliser cancels most of it, and an interpolator introduces a shape the original curve did not have, which the cancellation cannot touch. A finer grid is not more resolution when the table is not finer.

A finer reading of a coarser table

Interpolating a five-nanometre spectrum to one nanometre helps a daylight calculation by a factor of five and harms a three-emitter LED by a factor of a hundred and twenty thousand. Both are the same operation on the same table, and which one happens is decided by a property of the light nobody records.

light · Light The two tabulation choices over forty-two surfaces, under a tungsten lamp at 2856 K. Each column is one choice, measured over a family of forty-two analytic reflectances rather than on a single example: an absorption band of stated centre, width and depth. The four marks are the smallest, the median, the ninety-fifth percentile and the largest cost in ΔE₀₀, logarithmically. Under a smooth light the range is worth 6.3 times the step at the median, so a collection wanting one repair should widen its range rather than refine its step — and under a fluorescent tube the ranking reverses outright.

A neutral has no grid

A perfectly flat reflectance computes to exactly the same colour on every wavelength grid, through every slit, at every origin, and for every observer — not nearly, but to the last bit of a floating-point number. The condition is an identity rather than a limit, and what makes it one is the white point.

difference · Metric How much the answer moves when the 5-nanometre grid is slid through one cell. Each bar is the spread of one light's colour across five grid origins, all at the same 5-nanometre step, in ΔE₀₀. A smooth light barely moves, and what movement it has is the end cells rather than the sampling. The fluorescent tube moves by 3.18 units and the laser projector by 35.0, because their emission lines are narrower than the step and whether a sample lands on one is a coincidence of arithmetic. This is the measurement that separates a quadrature error from an aliasing error, and no average over origins can substitute for it.

Where the grid starts

Holding the step at five nanometres and sliding the grid's origin through one cell moves a fluorescent tube's computed colour by 3.18 ΔE₀₀ and a laser projector's by 35.0. Refining the step does not fix it and averaging over origins hides it. It is the one tabulation fault with no smooth error to cancel against.

light · Light The two tabulation choices over forty-two surfaces, under a 6500 K thermal radiator. Each column is one choice, measured over a family of forty-two analytic reflectances rather than on a single example: an absorption band of stated centre, width and depth. The four marks are the smallest, the median, the ninety-fifth percentile and the largest cost in ΔE₀₀, logarithmically. Under a smooth light the range is worth 9.1 times the step at the median, so a collection wanting one repair should widen its range rather than refine its step — and under a fluorescent tube the ranking reverses outright.

Which end to buy

Refining a five-nanometre grid to one buys a daylight calculation 0.05 ΔE₀₀ and widening its range buys 0.54. Under a fluorescent tube the same two purchases are worth 0.83 and 0.0001. The ranking reverses completely, and what decides it is one length compared against one other length.

light · Light What the normaliser cancels, per light. Two bars per light, logarithmic. The upper is the colour error a 5-nanometre sum makes when the white it is divided by is computed finely; the lower is the same sum divided by the white computed on the same coarse grid, which is what every colorimetric calculation actually does. The ratio is between 1.3 and 4.1. The grid appears twice in a tristimulus value and the two errors are the same error, so most of it divides out — which is why five nanometres has been good enough for a century without anybody having to be careful about it.

The normaliser carries the error too

A five-nanometre sum gets a red pigment's tristimulus value wrong by two hundredths of a per cent and its colour wrong by six hundredths of a unit. Those two numbers are not the same size because the grid appears twice in a colour — once in the sample and once in the white — and the two errors are largely the same error.

difference · Metric The rectangle sum against the trapezoid sum, under a 6500 K thermal radiator. Colorimetry's summation is the rectangle rule at the tabulated points. The trapezoid rule differs from it by exactly one thing — half a cell at each end of the range — and the gap between these two lines is therefore that term and nothing else. At five nanometres it is a factor of 14.3, which means the number everybody calls a sampling error is mostly a truncation error wearing the step's clothes. On a light whose lines are narrower than the step the two rules agree to three decimal places, because there the error really is the sampling.

The endpoint term has a name

The five-nanometre error on a smooth light falls linearly with the step, which is not what a sampling error does. It is the half-cell at each end of a truncated range, it is first order where the sampling is second, and halving two weights removes fourteen fifteenths of it for nothing.

light · Light What a tabulation step costs, by light, on paper. The horizontal axis is the tabulation step in nanometres, from one to twenty; the vertical is how far the resulting colour is from the same integral taken at a tenth of a nanometre over the same range, in ΔE₀₀, on a logarithmic scale. Each line is one light. The three with no feature narrower than the step fall smoothly and stay below a tenth of a unit at five nanometres, which is the grid used throughout. The fluorescent tube and the laser projector do not fall at all: their lines are narrower than any step drawn here, so the answer depends on where the samples land rather than on how many there are. The sample is held at paper throughout.

A grid is not a resolution

Ten essays into this collection there is one sentence about wavelength sampling, and it is that five nanometres is enough. Ten measurements later there are three decisions, three mechanisms, three repairs and two rankings, and the word resolution names none of them.

light · Light The directional solver reduces to the radiosity solver exactly. A solver with a new unknown in it is worth nothing until it reproduces the one it replaces. Setting every wall's bidirectional distribution to ρ/π collapses all thirty ordered-pair radiances onto their patch's radiosity divided by π, and the answer agrees with this collection's existing radiosity solution to 9.8e-16 relative — the floating-point floor. That is the check that makes every other number in this family a statement about lobes rather than about a new piece of arithmetic, and it is the reason the reduction is drawn rather than mentioned.

Thirty unknowns instead of six

A directional transport solver is worth nothing until it reproduces the one it replaces. Setting every wall's bidirectional distribution to ρ over π collapses thirty ordered-pair radiances onto six radiosities and reproduces this collection's existing answer to 9.8 × 10⁻¹⁶ relative — which is the only reason anything else it says can be believed.

scene · Scene Where a point-sampled patch stops resolving a lobe. The horizontal axis is the wall's roughness, logarithmic; the vertical is how far the answer moves when the cone quadrature is refined from eight directions to sixteen, also logarithmic. A patch in a cube subtends a cone of angular radius 25.8° at the face opposite, and a microfacet lobe of roughness a is about a radians wide, so a lobe below about 0.15 is narrower than the quadrature that samples it. The movement at 0.05 is 21.2 ΔE₀₀ and at 0.3 it is 0.033. This is where the method stops, not where the paint does: matt, eggshell and satin finishes are inside it and a high-gloss varnish is not.

Where a patch stops being a point

A patch in a cube subtends a cone of 25.8° at the face opposite. A microfacet lobe of roughness 0.05 is about three degrees wide. Point-sampling the second inside the first returned a chroma of thirty-four thousand and a negative lightness, and the boundary between working and not working is measured rather than declared.

limits · Limits The two tabulation choices over forty-two surfaces, under a 6500 K thermal radiator. Each column is one choice, measured over a family of forty-two analytic reflectances rather than on a single example: an absorption band of stated centre, width and depth. The four marks are the smallest, the median, the ninety-fifth percentile and the largest cost in ΔE₀₀, logarithmically. Under a smooth light the range is worth 9.1 times the step at the median, so a collection wanting one repair should widen its range rather than refine its step — and under a fluorescent tube the ranking reverses outright.

The grid under the census

The adaptation census is computed on eighty-one wavelengths from 380 to 780 nanometres. Under the daylight and blackbody sources it uses, the range is worth about half a colour difference on ordinary surfaces and the step about a twentieth — so the census carries a tabulation term as well as an observer one, and they are not the same size.

brain · Appearance A 2-nanometre notch, and what two five-nanometre grids record of it. The reflectance of a sample with a 2-nanometre notch at 552.3 nm, drawn finely from 535 to 570 nm. The dark ticks are the samples of a five-nanometre grid starting at 380 nm and the pale ticks those of the same grid started half a step later. The true minimum is 0.06; the first grid's deepest sample reads 0.70 and the second's 0.08. The sample has put a line into a calculation whose light has none.

The line can be in the sample

The rule for which tabulation defect to fix compares the light's narrowest feature with the grid's step, and it named its own failure case — a sample with structure narrower than the step. Given one, a two-nanometre notch under a thermal source with no feature at all costs 2.1 colour differences at five nanometres and moves 2.0 when the grid slides, which is what a fluorescent tube costs on a smooth pigment. Under that tube a notch twelve nanometres wide, more than twice the step, moves 8.1 when it sits on the mercury line.

light · Light What a five-nanometre grid costs a steep-sided notch, against its width, under a 6500 K source. The cost of a five-nanometre grid starting at 380 nm, against a reference at two hundredths of a nanometre, for a sample with a flat-bottomed notch at 552.3 nm under a 6500 K thermal radiator, against the notch's width from 2 to 30 nm, for edges rising in 0.4, 2.2, 6.6 nm. With the steepest edges the cost is 0.02 at 10 nm, 1.99 at 12.5 and 0.03 at 20: it rises and falls with the step as its period and does not die away as the notch widens. With the softest edges it stays under 0.10 at every width.

The cost of a steep notch repeats every step

A Gaussian notch is safe on a five-nanometre grid once it is a couple of steps wide. A flat-bottomed notch with steep sides never is. Its cost on the grid rises and falls with its width, with the step as its period — 0.02 colour differences at ten nanometres, 1.99 at twelve and a half, 0.03 at twenty — and it does not die away as the notch widens. What sets its size is how fast the edges rise, and an interference filter's edges rise in under a nanometre.

light · Light Six departures, and three ways of adding them up. The six audited observer departures on 120 smooth reflectances, across the dial. Summing them assumes they all point the same way and is an overestimate everywhere; taking the largest alone assumes only one matters and is an underestimate everywhere. Quadrature — the usual way of combining contributions taken to be independent — assumes they are mutually perpendicular, and the measured combination crosses it at a degree of 0.8618. Below that the departures are on balance pointing together and quadrature is too small; above it they are on balance pointing apart and quadrature is too large. It is exactly right in one room.

Quadrature is exact in one room

An observer allowance is built by adding the departures in quadrature, which assumes they are mutually perpendicular. Over fifteen pairs their angles run from 18 degrees to 179 and hardly any are perpendicular. Measured against the real combination, quadrature is too small in a cinema by eight per cent and too large under the sky by fourteen, crossing at an adapting luminance of 22 candelas a square metre — and on individual surfaces it is out by a third in both directions in every room.

difference · Metric The glossiest finish the solver can report, against what it costs to report it. Six quadratures, each with the roughness at which its answer stops being stable, on logarithmic axes. The line is a fit and its slope is -0.350: the reachable roughness falls as the cost to the power of about a third, so reaching a finish twice as glossy costs about 7 times the work. The solver used here sits at 108 directions and reports down to a roughness of 0.145, which is where its own note put the boundary by inspection.

The boundary belongs to the quadrature

The directional solver stops at a roughness of about 0.15, and below that its answers are not imprecise but unphysical. The boundary is where the lobe stops being resolved by the sampling, so it belongs to the discretisation rather than to the room — and moving it is a purchase. Measured across six quadratures the reachable roughness falls as the cost to the power of a third, so a finish twice as glossy costs seven times the work and a polished varnish costs two hundred and thirty-six times.

scene · Scene

Named alongside it

The objects these essays reach for when they reach for this one.

Wavelength gridMeasurement errorAuditAliasingConvergenceModelling assumptionSpecificationStructural choiceTest setIntegrationReflectanceSampling

All concepts