A wall's twist can be tested against the error the adjustment reports
Worth reading first: Another picture of the same sweep · The track and the scene together.
The bundle keeps the needles turned, if the wall has three rows followed a survey of a facade’s flatness from six cameras in two groups, sixty degrees either side of the wall’s normal. With the cameras’ poses known, the arrangement measured each mark’s departure from the wall’s plane four and a half times better than six cameras on a narrow arc. With the poses found from the same pictures in one bundle adjustment, it still did, provided the wall carried its marks in three rows. With ten marks in two rows, the adjustment’s uncertainty collapsed onto one shape: a twist, the top row’s left end forward and right end back and the bottom row the reverse, carrying seventy-one per cent of all the flatness error. A surveyor whose reconstructed wall came out as a saddle should suspect the reconstruction before the wall.
The essay ended on what the surveyor actually has. An adjustment does not only return the wall; it returns the covariance of everything it fitted, and so an uncertainty for the twist itself. The question was whether that uncertainty is worth believing: how often a fitted twist clears twice its own reported error when the wall is flat and when it is truly twisted, and whether the third row of marks that shrinks the twist also lets the adjustment say when to believe it.
It does both, and the report is honest at every layout. What changes with the layout is not the report’s truthfulness but its size.
The twist, estimated and reported
The rig is the earlier essay’s: a wall six metres wide with marks from half a metre to two and a half metres up, eight metres from two groups of three cameras centred sixty degrees either side of its normal, every mark found in every picture and read with half a pixel of error, every camera’s pose found with the points. A twisted wall is a saddle: each mark pushed out of the wall in proportion to its place across the marked area times its place up it, so the four corners move by a stated amount, two forward and two back, and the middle does not move.
The twist is estimated by least squares along that saddle’s pattern of departures, weighted by the departures’ full covariance from the adjustment, and reported with the standard deviation the same covariance implies. Both the estimate and its reported error come out of numbers the adjustment already computes; nothing is added to the survey but a line of arithmetic.
One property of the estimate matters before any number is read. A reconstruction from pictures alone can be placed, turned and scaled freely — seven numbers no picture can name — and every coordinate the adjustment reports moves with that choice. The twist does not. It is read from each mark’s departure from the plane fitted through all the marks, and a translation, a turn or a change of scale of the whole reconstruction moves every mark along a plane the fit absorbs. So the fitted twist and its reported error are the same whichever seven numbers the adjustment held, and a surveyor comparing two reports from two programs that fix the frame differently is comparing like with like. That is not true of most numbers an adjustment prints, and it is why the twist is a quantity a test can be built on.
The bars are the reported error, doubled, around each fitted twist, and the dashed line is the truth. For ten marks in two rows the adjustment states its twist to ±5.05 millimetres — ±10.1 doubled — and the twelve reconstructions of a wall truly twisted five millimetres scatter from −3.1 to 12.4 millimetres around it. One of the twelve clears zero by twice its error. For fifteen marks in three rows the stated error is ±2.31 millimetres, the reconstructions sit between 0.1 and 7.2, and seven clear zero.
The slider moves the truth. At zero, neither layout should call a twist, and over these twelve draws the two-row wall never does while the three-row wall does twice, which is the size of chance that twelve draws allow. At twenty millimetres both see it every time. The interesting range is between, and it is set by the reported error: five millimetres is one stated error on the two-row wall and more than two on the three-row one.
The report is honest at every layout
A reported error is worth believing only if a test against it has the false-alarm rate it promises. A test at twice the reported error promises about one false alarm in twenty.
On a flat wall every layout calls the twist real four or five times in a hundred, as a test at twice its error should. That is the half of the question with the cleanest answer: the adjustment’s own number is honest about the twist, at ten marks as at thirty-nine, in two rows as in three. It is honest because the adjustment’s covariance is the covariance of a linear estimate under the error model it was given, and the error model here is right — the marks really are read with independent errors of half a pixel. When it is not right the test misleads, and a later section takes that apart.
What the layout decides is the power. A real twist of five millimetres is called real 17 times in a hundred from ten marks in two rows, 30 from twenty-six marks in two rows, 56 from fifteen in three rows and 91 from thirty-nine in three. To be called real nine times in ten a twist must be about sixteen millimetres on the ten-mark wall, twelve on twenty-six marks in two rows, eight on fifteen in three, and five on thirty-nine.
All four curves are one curve read at different scales. The fitted twist is the true twist plus a normal scatter of exactly the stated error, so the chance of clearing twice that error depends on the true twist only through its ratio to the stated error: about one in six at one error, a half at two, nineteen in twenty at a little under four. A layout’s whole contribution is the stated error it produces, and a surveyor can read the power of the test off the adjustment’s report before deciding what a saddle of a given size would mean.
The earlier essay guessed that a real five-millimetre twist would stand at one stated error from the two-row wall and at three from the three-row one. The two-row guess was exact: 5.05 millimetres stated, one error. The three-row wall of fifteen marks states 2.31, so five millimetres is 2.2 errors; thirty-nine marks state 1.50, so it is 3.3. The third row is what lets the survey see a twist of this size at all, and adding marks along the two rows a two-row wall already has buys much less — twenty-six marks in two rows are worse than fifteen in three.
What the groups’ opening buys
The twist is the two groups’ relative turn about the wall’s vertical showing through the marks: if one group’s estimated pose is turned a little against the other’s, every point they triangulate together is pushed forward on one side and back on the other, more at the top than the bottom of the wall’s marked area. The groups’ separation and the marks’ spread both pin that turn, as the spread a point gets found the spread of directions, not the count of cameras, pinning a single point.
Opening the groups helps every layout and helps steadily: from fifteen degrees to seventy-five, the two-row wall’s stated error falls by a factor of three and the three-row wall’s by more than five. The third row helps at every opening, and its advantage grows as the groups open, since a mark off the two rows’ lines is exactly what a turn of one group against the other moves differently from the rest. Split the track and the needles turn found the groups’ opening paying for each point’s depth; the twist is the bundle’s version of the same purchase, made with the marks’ layout as well as the cameras’.
So the earlier essay’s advice about the wall — three rows, a few dozen marks — is also advice about the adjustment’s honesty in the useful sense. A survey whose stated twist error is a millimetre and a half can tell a five-millimetre twist from a flat wall; one whose stated error is five millimetres can only report that it does not know.
A test against an assumed error tests the assumption
Every figure so far has given the adjustment the right reading error. A surveyor rarely knows it. Marks are read to half a pixel on a good day with sharp targets; on a soft wall with natural features they may be read to two.
The reported error is computed from the reading error the surveyor supplies. Marks read twice as badly as assumed make every fitted twist twice as noisy and leave the reported error where it was, so a flat wall clears twice its stated error 32 times in a hundred; read four times as badly, 61 times. The layout does not help, because the failure is in the assumption and not in the geometry: the three-row wall’s false alarms climb exactly as the two-row wall’s do. A wall read better than assumed calls fewer twists than it should, real ones included.
The repair is already in the adjustment. Its residuals — how far each reading sits from where the fitted cameras and points put it — measure the reading error directly, and with sixty-one redundant observations on the ten-mark wall and a hundred and six on the fifteen-mark one, they measure it well. Rescale the reported error by the reading error the residuals imply and the false-alarm rate returns to about one in twenty at every true error. Where the adjustment stops found the reading error setting the floor of every residual; here that floor is the instrument the twist test needs.
This is the same lesson an uncertainty is quoted from something drew about the gauge: a stated uncertainty means what its inputs mean. A twist error computed from an assumed half pixel is a statement about half a pixel. One computed from the residuals is a statement about the survey.
A mismatch that fakes a twist shows in its own residual
The earlier essay warned that a single mismatched mark carries far more weight among ten than among thirty-nine. A mark matched to the wrong feature in one picture is not a small error, and an adjustment that absorbs it could spread it into the weakest shape it has, which is the twist.
A wrong match is not a small error made the point about a single point: a mismatch moves it somewhere else rather than slightly off its own place. Inside an adjustment the same mismatch is shared out. Part of it stays in the bad reading’s residual, part moves the point that reading belongs to, and part moves the cameras that saw it — and a camera moved is a group turned, which is a twist.
It does spread. One reading ten pixels wrong makes a flat two-row wall look twisted 24 times in a hundred; twenty pixels wrong, 33 times, and on the fifteen-mark three-row wall 53 times, because a smaller stated error is easier for a bias to clear. A mismatch is a twist-maker, and a surveyor who took the test at face value would be fooled about once in three by a wrong match of that size.
But it is never fooled silently. The wrong reading’s own residual, divided by its expected size, exceeds three every time once the mistake reaches five pixels, on both layouts, and 85 to 93 times in a hundred at three pixels. A twist called and the bad reading not caught happens 1.5 per cent of the time at three pixels on the two-row wall and never at five. The reason is the readings’ leverage: no reading in either layout has more than about three-quarters of its own error absorbed by the adjustment, so a mismatch big enough to push the twist past twice its error leaves about a quarter of itself or more in the residual, where an ordinary test for outliers sees it.
The order of operations follows. Before reading the twist, test the residuals; remove or re-read any reading whose residual is more than three times its expected size; then read the twist against the error the remaining residuals imply. Done in that order, the test keeps its promise of one false alarm in twenty.
What the adjustment can say about its own saddle
Put together, a bundle adjustment can be asked whether the twist it found is real, and its answer is trustworthy on three conditions: the reading error is measured from the residuals rather than assumed, gross mismatches are removed first, and the layout gives the test enough power to matter. On those conditions a flat wall is called twisted one time in twenty and a real twist is called real as often as its size over the stated error allows — five millimetres seventeen times in a hundred from ten marks in two rows, ninety-one from thirty-nine in three.
The earlier essay’s warning about saddles stands, sharpened. A two-row wall that comes out as a five-millimetre saddle is, on these numbers, a wall the adjustment could not have distinguished from a flat one most of the time, and its own report says so: one stated error. A three-row wall that comes out the same way is a twist the adjustment believes at better than two errors, and on a well-marked wall at better than three. The difference between those two statements is a row of marks.
A survey is trusted at its own accuracy found that an adjustment does best when its inputs are weighted by their true accuracy and is not much hurt by being told they are worse. The twist test is the mirror case: a test of the outputs is trustworthy when its error is the true error, and a surveyor who must guess should take the residuals’ word over their own.
The linear world the numbers live in
The errors are small. Every estimate here is the linearised adjustment, exact for errors small beside the geometry, which half a pixel on a six-metre wall at eight metres is. The mismatch figures leave that world by design, but they test only whether a gross error is visible, which the linear residual answers correctly.
The twist is a saddle. A wall that is bowed, bulged or racked in some other pattern projects onto the saddle only partly, and the test reads only that part. A surveyor looking for a different shape should estimate along that shape with the same arithmetic; the reported error will be different, and on the two-row wall it will usually be smaller, since the saddle is the adjustment’s loosest shape.
Every mark is found in every picture. A real matcher loses oblique marks, and a mark seen by one group only ties nothing between the groups. The twist error then grows, and the test’s power falls with it, but its honesty does not: the covariance of a survey with missing marks is still the right covariance for that survey.
The camera calibration is known. Every camera shares one focal length, held fixed. An adjustment that also fits the focal length has a further loose direction, scale against distance, which is a plane and is taken out by the flatness reading, but its coupling to the groups’ relative turn was not measured here.
Still open: whether a third group of cameras does what a third row does
The twist is loose because two groups are tied to each other only through the marks both see, and the third row of marks tightens the tie. A third group of cameras would tighten it differently: a few pictures taken square to the wall see every mark the two oblique groups see, and so tie each to the middle rather than to each other.
The measurement that settles what that is worth adds one, two or three cameras facing the wall squarely to the two groups at sixty degrees, on the ten-mark two-row wall, and asks how far the twist’s reported error falls compared with adding a third row of marks instead — and whether three square-on pictures, which cost a photographer a few seconds, buy the five-millimetre twist the three-row wall sees, so that a plain wall with only a string course and a cornice to mark can still be surveyed for a saddle.
Shares its objects with
Essays that name at least two of the same things, and that neither author linked.
- A length fixes the scale where it lies, and the worst camera nowhere — both name bundle adjustment, covariance, gauge freedom
- A loop's far side is a length — both name bundle adjustment, covariance, gauge freedom
- A mismatch on its own line needs a third eye — both name least squares, outlier, residual
- A scale bar is worth its ends, not its tape — both name bundle adjustment, covariance, gauge freedom
- A stepped hand passes for a tiring one on a short strip — both name degrees of freedom, least squares, residual
- A stepped hand's walk is in the creeps it was given — both name degrees of freedom, least squares, residual
Named objects
A flat tag is an object no other essay names yet.
bundle adjustmentCovariancedegrees of freedomgauge freedomleast squaresOutlierResidual