FAILURE MAP
← Case archive

FA-69911 / Map projection transforms / Open access

Slippy map tile index from coordinates: column rounding · case 01

Points in the eastern half of a tile are assigned to the next column.

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

ROOT CAUSE

The fractional column is rounded to nearest instead of floored.

VERIFIED REPAIR

At the column rounding step restore `math.floor(xf)`, leaving the rest of the model unchanged.

Unsuccessful approach: ceil(xf) - 1 differs from floor(xf) on exact tile edges, pushing edge points one column west.

Case contract

Input [lon, lat, z] with lon in [-180, 180], any latitude, z >= 0. With n = 2**z tiles per axis, column = floor((lon+180)/360*n) and row = floor((1 - asinh(tan(lat))/pi)/2*n) (row 0 at the north edge); both are clamped to [0, n-1] so lon=180 falls in the last column and latitudes beyond the mercator limit fall in the edge rows. Return [z, column, row].

Why this case matters

Tile servers and caches key storage by these indices; an off-by-one or flipped axis fetches the wrong image.

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, z = x
    n = 2 ** z
    phi = math.radians(max(-89.9, min(89.9, lat)))
    xf = (lon + 180.0) / 360.0 * n
    yf = (1.0 - math.asinh(math.tan(phi)) / math.pi) / 2.0 * n
    xt = min(max(int(round(xf)), 0), n - 1)
    yt = min(max(math.floor(yf), 0), n - 1)
    return [z, xt, yt]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('boundary #0', [180.0, 0.0, 3], [3, 7, 4]), ('boundary #1', [-180.0, 10.0, 4], [4, 0, 7]), ('boundary #2', [10.0, 89.0, 5], [5, 16, 0]), ('boundary #3', [10.0, -89.0, 5], [5, 16, 31]), ('control #4', [13.4, 52.5, 10], [10, 550, 335]), ('control #5', [-74.0, 40.7, 12], [12, 1206, 1540]), ('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1])], [('boundary #3', [10.0, -89.0, 5], [5, 16, 31]), ('control #4', [13.4, 52.5, 10], [10, 550, 335]), ('control #5', [-74.0, 40.7, 12], [12, 1206, 1540]), ('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('boundary #23', [-67.5, 20.0, 4], [4, 5, 7])], [('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('control #9', [-122.4, 37.8, 9], [9, 81, 197]), ('control #10', [100.5, 13.7, 6], [6, 49, 29]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0])], [('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('control #10', [100.5, 13.7, 6], [6, 49, 29]), ('control #11', [-58.4, -34.6, 11], [11, 691, 1234]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('control #13', [37.6, 55.75, 8], [8, 154, 80]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0]), ('boundary #23', [-67.5, 20.0, 4], [4, 5, 7])], [('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #9', [-122.4, 37.8, 9], [9, 81, 197]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('control #15', [18.4, -33.9, 13], [13, 4514, 4916]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0]), ('control #17', [139.7, 35.7, 14], [14, 14549, 6450]), ('boundary #18', [180.0, 45.0, 6], [6, 63, 23]), ('boundary #19', [180.0, -60.0, 10], [10, 1023, 726])]]
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 fixtureActualExpectedOutcome
boundary #0[3, 7, 4][3, 7, 4]Passed
boundary #1[4, 0, 7][4, 0, 7]Passed
boundary #2[5, 17, 0][5, 16, 0]Failed
boundary #3[5, 17, 31][5, 16, 31]Failed
control #4[10, 550, 335][10, 550, 335]Passed
control #5[12, 1206, 1540][12, 1206, 1540]Passed
control #6[7, 118, 76][7, 117, 76]Failed
control #7[1, 1, 1][1, 1, 1]Passed

SHA-256 / c3492ca6abd3b95373ecc87d97923a0f09d89aae0b3d2d06a1ff5da99f621ae6

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, z = x
    n = 2 ** z
    phi = math.radians(max(-89.9, min(89.9, lat)))
    xf = (lon + 180.0) / 360.0 * n
    yf = (1.0 - math.asinh(math.tan(phi)) / math.pi) / 2.0 * n
    xt = min(max(math.ceil(xf) - 1, 0), n - 1)
    yt = min(max(math.floor(yf), 0), n - 1)
    return [z, xt, yt]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('boundary #0', [180.0, 0.0, 3], [3, 7, 4]), ('boundary #1', [-180.0, 10.0, 4], [4, 0, 7]), ('boundary #2', [10.0, 89.0, 5], [5, 16, 0]), ('boundary #3', [10.0, -89.0, 5], [5, 16, 31]), ('control #4', [13.4, 52.5, 10], [10, 550, 335]), ('control #5', [-74.0, 40.7, 12], [12, 1206, 1540]), ('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1])], [('boundary #3', [10.0, -89.0, 5], [5, 16, 31]), ('control #4', [13.4, 52.5, 10], [10, 550, 335]), ('control #5', [-74.0, 40.7, 12], [12, 1206, 1540]), ('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('boundary #23', [-67.5, 20.0, 4], [4, 5, 7])], [('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('control #9', [-122.4, 37.8, 9], [9, 81, 197]), ('control #10', [100.5, 13.7, 6], [6, 49, 29]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0])], [('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('control #10', [100.5, 13.7, 6], [6, 49, 29]), ('control #11', [-58.4, -34.6, 11], [11, 691, 1234]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('control #13', [37.6, 55.75, 8], [8, 154, 80]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0]), ('boundary #23', [-67.5, 20.0, 4], [4, 5, 7])], [('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #9', [-122.4, 37.8, 9], [9, 81, 197]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('control #15', [18.4, -33.9, 13], [13, 4514, 4916]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0]), ('control #17', [139.7, 35.7, 14], [14, 14549, 6450]), ('boundary #18', [180.0, 45.0, 6], [6, 63, 23]), ('boundary #19', [180.0, -60.0, 10], [10, 1023, 726])]]
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 fixtureActualExpectedOutcome
boundary #0[3, 7, 4][3, 7, 4]Passed
boundary #1[4, 0, 7][4, 0, 7]Passed
boundary #2[5, 16, 0][5, 16, 0]Passed
boundary #3[5, 16, 31][5, 16, 31]Passed
control #4[10, 550, 335][10, 550, 335]Passed
control #5[12, 1206, 1540][12, 1206, 1540]Passed
control #6[7, 117, 76][7, 117, 76]Passed
control #7[1, 0, 1][1, 1, 1]Failed

SHA-256 / 9324ac15a8071c1d7bab5686f5e41d366b965d9385399291d7097ce471489535

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, z = x
    n = 2 ** z
    phi = math.radians(max(-89.9, min(89.9, lat)))
    xf = (lon + 180.0) / 360.0 * n
    yf = (1.0 - math.asinh(math.tan(phi)) / math.pi) / 2.0 * n
    xt = min(max(math.floor(xf), 0), n - 1)
    yt = min(max(math.floor(yf), 0), n - 1)
    return [z, xt, yt]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('boundary #0', [180.0, 0.0, 3], [3, 7, 4]), ('boundary #1', [-180.0, 10.0, 4], [4, 0, 7]), ('boundary #2', [10.0, 89.0, 5], [5, 16, 0]), ('boundary #3', [10.0, -89.0, 5], [5, 16, 31]), ('control #4', [13.4, 52.5, 10], [10, 550, 335]), ('control #5', [-74.0, 40.7, 12], [12, 1206, 1540]), ('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1])], [('boundary #3', [10.0, -89.0, 5], [5, 16, 31]), ('control #4', [13.4, 52.5, 10], [10, 550, 335]), ('control #5', [-74.0, 40.7, 12], [12, 1206, 1540]), ('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('boundary #23', [-67.5, 20.0, 4], [4, 5, 7])], [('control #6', [151.2, -33.9, 7], [7, 117, 76]), ('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('control #9', [-122.4, 37.8, 9], [9, 81, 197]), ('control #10', [100.5, 13.7, 6], [6, 49, 29]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0])], [('control #8', [2.35, 48.85, 15], [15, 16597, 11274]), ('control #10', [100.5, 13.7, 6], [6, 49, 29]), ('control #11', [-58.4, -34.6, 11], [11, 691, 1234]), ('boundary #12', [179.99, -85.0, 2], [2, 3, 3]), ('control #13', [37.6, 55.75, 8], [8, 154, 80]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0]), ('boundary #23', [-67.5, 20.0, 4], [4, 5, 7])], [('control #7', [0.0, 0.0, 1], [1, 1, 1]), ('control #9', [-122.4, 37.8, 9], [9, 81, 197]), ('control #14', [-0.12, 51.5, 16], [16, 32746, 21794]), ('control #15', [18.4, -33.9, 13], [13, 4514, 4916]), ('boundary #16', [0.0, 0.0, 0], [0, 0, 0]), ('control #17', [139.7, 35.7, 14], [14, 14549, 6450]), ('boundary #18', [180.0, 45.0, 6], [6, 63, 23]), ('boundary #19', [180.0, -60.0, 10], [10, 1023, 726])]]
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 fixtureActualExpectedOutcome
boundary #0[3, 7, 4][3, 7, 4]Passed
boundary #1[4, 0, 7][4, 0, 7]Passed
boundary #2[5, 16, 0][5, 16, 0]Passed
boundary #3[5, 16, 31][5, 16, 31]Passed
control #4[10, 550, 335][10, 550, 335]Passed
control #5[12, 1206, 1540][12, 1206, 1540]Passed
control #6[7, 117, 76][7, 117, 76]Passed
control #7[1, 1, 1][1, 1, 1]Passed

SHA-256 / f62230f3bd9d4efd3b4de4373f6c678b16435b4281eacbfc4412ba7732746c52

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:15.817904+00:00.

Case digest / 1cfb4cb1f90e0de5b2b8a97a260c3a3baab8e024ed1dde0d4e7163267d3486a7