FA-70069 / Map projection transforms / Member archive

Equirectangular projection with standard parallel: inverse origin latitude · case 04

Inverse latitudes are offset when the grid origin is not on the equator.

Member previewVariant 4 · 3 implementations · 8 checks per implementation

Case contract

Input [mode, a, b, lon0, lat0, lat1] on a sphere R = 6371007. mode "fwd": a=lon, b=lat; dlon = lon-lon0 wrapped into [-180, 180); x = R*rad(dlon)*cos(rad(lat1)), y = R*rad(lat-lat0). mode "inv": a=x, b=y; lon = lon0 + deg(x/(R*cos(rad(lat1)))) wrapped into [-180, 180), lat = lat0 + deg(y/R). Forward results round to 2 decimals, inverse results to 6.

Why this case matters

Plate carree style grids with a true-scale parallel are common for regional rasters and simple web previews.

One recorded failure

Sample boundary fixture

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

Boundary fixture	Actual	Expected	Outcome
control #10	`[-111.774572, -62.486373]`	`[-111.774572, -42.486373]`	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.