FA-70207 / Map projection transforms / Member archive

Azimuthal equidistant range and bearing: antipode exclusion · case 02

Points at or beside the antipode are plotted at absurd distances instead of being rejected.

Member previewVariant 2 · 3 implementations · 8 checks per implementation

Case contract

Input [lon, lat, lon0, lat0]; sphere R = 6371000 centred on (lon0, lat0). cos c is clamped to [-1, 1], c = acos. Points with c > pi - 1e-6 (the antipode neighbourhood) return None. k = 1 at the centre (c < 1e-12) else c/sin c. x = R k cos(p) sin(dl); y = R k (cos(p0) sin(p) - sin(p0) cos(p) cos(dl)). Return [x rounded to 1 decimal, y rounded to 1 decimal, hypot(x, y)/1000 km rounded to 3, azimuth clockwise from north in [0, 360) rounded to 4].

Why this case matters

Radio coverage, flight-range rings and seismology plots use this projection to read range and bearing directly.

One recorded failure

Sample boundary fixture

This sample comes from the broken implementation of a controlled reproducer.

Boundary fixture	Actual	Expected	Outcome
boundary #24	`[0.0, 20028776.7, 20028.777, 0.0]`	`null`	Failed

MEMBER ARCHIVE

The complete case is available to members.

This record includes three runnable implementations, regression fixtures, execution results, and source hashes.

Member access is invitation-based. Sign in with your invited account to inspect the sources.