FA-70246 / Map projection transforms / Open access

Mollweide forward with Newton solve: easting auxiliary angle · case 01

The outline is not an ellipse; high latitudes are too wide.

Verified by executionVariant 1 · 8 checks per implementationDownload source bundle ↓JSON ↗

ROOT CAUSE

The easting uses cos(latitude) instead of cos of the auxiliary angle.

VERIFIED REPAIR

At the easting auxiliary angle step restore `dl * math.cos(theta)`, leaving the rest of the model unchanged.

Unsuccessful approach: The patched term uses half the auxiliary angle.

Case contract

Input [lon, lat, lon0]; sphere R = 6371000, dlon wrapped to [-180, 180). Solve 2t + sin 2t = pi sin(phi) by Newton from t = phi, at most 50 steps, stopping when the residual is below 1e-12, with derivative 2 + 2 cos 2t. x = R (2 sqrt 2 / pi) dlon cos t, y = R sqrt 2 sin t, rounded to 2 decimals.

Why this case matters

Global thematic maps use Mollweide for equal area; an under-converged auxiliary angle skews high latitudes.

1 / The failure

Exit 1

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    lon, lat, lon0 = x
    R = 6371000.0
    phi = math.radians(lat)
    dl = math.radians(((lon - lon0 + 180.0) % 360.0) - 180.0)
    theta = phi
    for _ in range(50):
        f = 2 * theta + math.sin(2 * theta) - math.pi * math.sin(phi)
        if abs(f) < 1e-12:
            break
        theta -= f / (2 + 2 * math.cos(2 * theta))
    px = R * 2 * math.sqrt(2) / math.pi * dl * math.cos(phi)
    py = R * math.sqrt(2) * math.sin(theta)
    return [round(px, 2), round(py, 2)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('control #0', [101.361, 11.005, 160.0], [-5803523.67, 1355983.62]), ('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #7', [114.805, 18.464, -90.0], [-15037610.8, 2265254.5]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #10', [-134.836, -65.77, 160.0], [3720190.56, -7401319.17]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #13', [-49.492, 24.945, 0.0], [-4663467.36, 3043363.54]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #16', [120.053, 58.23, -90.0], [-10039242.4, 6698546.76]), ('control #17', [164.893, -9.572, 10.0], [15372850.89, -1180095.58]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])]]
for label, args, expected in fixtures[N-1]:
    check(label, solve(args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)

Boundary fixture	Actual	Expected	Outcome
control #0	`[-5762432.57, 1355983.62]`	`[-5803523.67, 1355983.62]`	Failed
control #1	`[6011366.97, 1966970.83]`	`[6102963.12, 1966970.83]`	Failed
control #2	`[-6219644.27, -6684729.63]`	`[-7888806.65, -6684729.63]`	Failed
control #3	`[-3942442.81, -6673720.01]`	`[-4994731.15, -6673720.01]`	Failed
control #4	`[2245159.06, 8144275.14]`	`[3649809.48, 8144275.14]`	Failed
control #5	`[8223217.54, 6828509.25]`	`[10593963.16, 6828509.25]`	Failed
boundary #20	`[0.0, 9009954.61]`	`[0.0, 9009954.61]`	Passed
control #22	`[0.0, 0.0]`	`[0.0, 0.0]`	Passed

SHA-256 / 83323f8b499d1cb32b58752c3ffa69ab7184203ba91da8cf6f1a7b51c2463006

2 / The unsuccessful fix

Exit 1

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    lon, lat, lon0 = x
    R = 6371000.0
    phi = math.radians(lat)
    dl = math.radians(((lon - lon0 + 180.0) % 360.0) - 180.0)
    theta = phi
    for _ in range(50):
        f = 2 * theta + math.sin(2 * theta) - math.pi * math.sin(phi)
        if abs(f) < 1e-12:
            break
        theta -= f / (2 + 2 * math.cos(2 * theta))
    px = R * 2 * math.sqrt(2) / math.pi * dl * math.sin(math.pi / 2 - theta / 2)
    py = R * math.sqrt(2) * math.sin(theta)
    return [round(px, 2), round(py, 2)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('control #0', [101.361, 11.005, 160.0], [-5803523.67, 1355983.62]), ('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #7', [114.805, 18.464, -90.0], [-15037610.8, 2265254.5]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #10', [-134.836, -65.77, 160.0], [3720190.56, -7401319.17]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #13', [-49.492, 24.945, 0.0], [-4663467.36, 3043363.54]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #16', [120.053, 58.23, -90.0], [-10039242.4, 6698546.76]), ('control #17', [164.893, -9.572, 10.0], [15372850.89, -1180095.58]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])]]
for label, args, expected in fixtures[N-1]:
    check(label, solve(args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)

Boundary fixture	Actual	Expected	Outcome
control #0	`[-5853646.45, 1355983.62]`	`[-5803523.67, 1355983.62]`	Failed
control #1	`[6215983.5, 1966970.83]`	`[6102963.12, 1966970.83]`	Failed
control #2	`[-10753029.74, -6684729.63]`	`[-7888806.65, -6684729.63]`	Failed
control #3	`[-6797256.68, -6673720.01]`	`[-4994731.15, -6673720.01]`	Failed
control #4	`[7209953.41, 8144275.14]`	`[3649809.48, 8144275.14]`	Failed
control #5	`[14760239.6, 6828509.25]`	`[10593963.16, 6828509.25]`	Failed
boundary #20	`[0.0, 9009954.61]`	`[0.0, 9009954.61]`	Passed
control #22	`[0.0, 0.0]`	`[0.0, 0.0]`	Passed

SHA-256 / 02b2f4becc2761095ee368c0534f146a9ceafe2be10e1d0ec8d2f8c20f1aa6df

3 / The verified repair

Exit 0

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    lon, lat, lon0 = x
    R = 6371000.0
    phi = math.radians(lat)
    dl = math.radians(((lon - lon0 + 180.0) % 360.0) - 180.0)
    theta = phi
    for _ in range(50):
        f = 2 * theta + math.sin(2 * theta) - math.pi * math.sin(phi)
        if abs(f) < 1e-12:
            break
        theta -= f / (2 + 2 * math.cos(2 * theta))
    px = R * 2 * math.sqrt(2) / math.pi * dl * math.cos(theta)
    py = R * math.sqrt(2) * math.sin(theta)
    return [round(px, 2), round(py, 2)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('control #0', [101.361, 11.005, 160.0], [-5803523.67, 1355983.62]), ('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #7', [114.805, 18.464, -90.0], [-15037610.8, 2265254.5]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #10', [-134.836, -65.77, 160.0], [3720190.56, -7401319.17]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #13', [-49.492, 24.945, 0.0], [-4663467.36, 3043363.54]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #16', [120.053, 58.23, -90.0], [-10039242.4, 6698546.76]), ('control #17', [164.893, -9.572, 10.0], [15372850.89, -1180095.58]), ('boundary #20', [0.0, 90.0, 0.0], [0.0, 9009954.61]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])]]
for label, args, expected in fixtures[N-1]:
    check(label, solve(args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)

Boundary fixture	Actual	Expected	Outcome
control #0	`[-5803523.67, 1355983.62]`	`[-5803523.67, 1355983.62]`	Passed
control #1	`[6102963.12, 1966970.83]`	`[6102963.12, 1966970.83]`	Passed
control #2	`[-7888806.65, -6684729.63]`	`[-7888806.65, -6684729.63]`	Passed
control #3	`[-4994731.15, -6673720.01]`	`[-4994731.15, -6673720.01]`	Passed
control #4	`[3649809.48, 8144275.14]`	`[3649809.48, 8144275.14]`	Passed
control #5	`[10593963.16, 6828509.25]`	`[10593963.16, 6828509.25]`	Passed
boundary #20	`[0.0, 9009954.61]`	`[0.0, 9009954.61]`	Passed
control #22	`[0.0, 0.0]`	`[0.0, 0.0]`	Passed

SHA-256 / bbd817eccdf9b0441c0fb033c6950436a8b7653c4bc1a6442b5804df72192d10

Verification & scope

Stipulated deterministic toy contract on a bounded input domain; results are rounded as stated and no conformance with any published standard or library is claimed. This reproducer isolates one failure mechanism. Results cover the supplied fixtures. Variants within a family share a test contract and should remain grouped when constructing evaluation splits. Related mechanisms with a shared evaluation_group must also remain together; these controlled models are not independent production incidents.

Observations recorded using Python 3.12.14 at 2026-09-29T14:48:19.073333+00:00.

Case digest / 600e715ca35f9dc537f267aad8018f8493bd7c09537019673bc47b467af83705