Lambert Conformal Conical Projection

• Lambert conformal conical projection is a cartographic method that preserves local angles by leveraging two standard parallels to ensure scale accuracy.
• Its methodology employs closed-form forward and inverse formulas that optimize the exponent n and constant F to minimize bi-Lipschitz distortion.
• Empirical analysis, such as in mapping the Russian Empire, shows this projection achieves less than 0.87% maximum scale error, outperforming competing conical projections.

The Lambert conformal conical projection is a cartographic method for mapping a subset of the sphere—often representing a designated terrestrial region such as a country—onto the Euclidean plane. Its principal aim is to minimize specific geometric distortions, notably by providing local conformality, i.e., precise preservation of angles at every point. This projection is characterized by the utilization of two standard parallels, selected so that the projection surface (a cone) is tangent to the sphere along these parallels, resulting in optimal scale distribution within the targeted latitudinal belt. Optimized for both analytical tractability and minimal bi-Lipschitz distortion, the Lambert conformal conical projection is extensively validated in its application to the mapping of wide latitudinal zones, including historical instances such as Euler’s general map of the Russian Empire (Miyachi et al., 13 Jan 2026).

1. Mathematical Formulation

The forward and inverse projection formulas derive from explicit analytical expressions. Denote by $\phi$ the geographic latitude, $\lambda$ the longitude, $R$ the mean Earth radius, $\lambda_{0}$ the central meridian, and $\phi_{1}<\phi_{2}$ the selected standard parallels. The exponent $n$ is defined as:

$$n = \frac{\ln\cos\phi_{1} - \ln\cos\phi_{2}}{\ln\tan(\frac{\pi}{4} + \frac{\phi_{2}}{2}) - \ln\tan(\frac{\pi}{4} + \frac{\phi_{1}}{2})}$$

A scale factor constant $F$ is then introduced:

$$F = \frac{\cos\phi_{1}\,[\tan(\frac{\pi}{4}+\frac{\phi_{1}}{2})]^n}{n}$$

For a point $(\phi,\lambda)$, the radial (distance from cone apex) and angular coordinates on the developed cone are:

$$\rho(\phi) = F \left[\tan\left(\frac{\pi}{4} + \frac{\phi}{2}\right)\right]^{-n}$$

$\theta = n(\lambda - \lambda_{0})$

The corresponding planar coordinates are:

$$x = \rho(\phi)\,\sin\theta, \quad y = \rho(\phi)\,\cos\theta$$

Inverse formulas recover $(\phi,\lambda)$ from $(x,y)$:

$\rho = \sqrt{x^2 + y^2},\ \theta = \arctan2(x, y)$

$$\phi = 2\,\arctan\left[\left(F/\rho\right)^{1/n}\right] - \frac{\pi}{2},\quad \lambda = \lambda_{0} + \frac{\theta}{n}$$

All relevant quantities admit closed-form representations, supporting direct computational implementation and analytical optimization (Miyachi et al., 13 Jan 2026).

2. Role of Standard Parallels and Central Meridian

The two standard parallels $(\phi_{1}, \phi_{2})$ are chosen to ensure the cone intersects the sphere precisely along these lines of latitude. By construction, the scale factor along these parallels is unity. The central meridian $\lambda_{0}$ determines the origin for angular measurements on the cone—typically the axis of symmetry for the mapped region. The parameter $n$ governs the convergence of meridians in the conical model, with $n = \sin\alpha$ in the traditional apex-angle parameterization. The constant $F$ ensures the radial scale matches that of the sphere along the standard parallels (Miyachi et al., 13 Jan 2026).

3. Distortion Analysis: Tissot Indicatrix and Scale Variation

This projection is exactly conformal: infinitesimal circles on the sphere are mapped to circles on the plane cone, yielding zero angular distortion. The scale factor is isotropic at each point. The local scale $k(\phi)$, identical in both meridional and parallel directions, is given by:

$$k(\phi) = \frac{nF\,[\tan(\frac{\pi}{4} + \frac{\phi}{2})]^{-n}}{R\cos\phi}$$

In Tissot’s terminology, $a(\phi) = b(\phi) = k(\phi)$, and the angular shear $\theta_{\rm max}=0$. Area distortion is $k(\phi)^2$. The scale equals one precisely at the standard parallels; elsewhere, maximum deviation from true scale is minimized by optimizing $n$ (Miyachi et al., 13 Jan 2026).

4. Methodology for Distortion Optimization and Comparative Metrics

Evaluation proceeds by regarding each projection as a homeomorphism from the spherical annulus $\phi_{1} \leq \phi \leq \phi_{2}$ onto a planar annulus. The primary measure of distortion is the bi-Lipschitz constant:

$$\delta = \ln(\sup_{p} L(p)) - \ln(\inf_{p} L(p)) = \ln(\max k / \min k), \quad \text{(for conformal projections)}$$

where $L(p)$ is the local Lipschitz constant at $p$. Projection parameters (notably, $n$) are then optimized to minimize $\delta$, providing the tightest possible scale uniformity over the region. Classically, this approach is equivalent to setting the standard parallels as scale-true and adjusting $n$ to confine interior scale deviations (Miyachi et al., 13 Jan 2026).

5. Empirical Results for Russian Empire Mapping

A comprehensive numerical analysis was conducted on the latitudinal span of Euler’s general map of the Russian Empire, with $\phi_{1}=47^\circ 30'$, $\phi_{2}=67^\circ 30'$. The optimal exponent obtained is:

$$n_{0} = \sin \alpha_{0},\quad \alpha_{0} \approx 55^\circ 14',\quad n_{0} \approx 0.8215294$$

$\delta_{\rm Lambert} \approx 0.0086263354$

The maximum scale error is thus $\exp(\delta_{\rm Lambert}) \approx 1.00866$, corresponding to a deviation of less than $0.87\%$ from true scale anywhere within the belt. The numerical comparison is summarized below:

Projection	$\delta$
Central (gnomonic on cone)	0.0171839
Delisle–Euler (true on $\phi_{1},\phi_{2}$)	0.00862621
Delisle equidistant cone	0.00921812
Orthogonal cone	0.00866925
Teichmüller extremal cone	0.0115244
Lambert conformal cone	0.0086263354

Graphs of the bi-Lipschitz constant $\sigma(\phi) = \max\{k(\phi), 1/k(\phi)\}$ show the Lambert conformal cone to be uniformly below or coincident with the best competing schemes, demonstrating minimal distortion (Miyachi et al., 13 Jan 2026).

6. Advantages Over Competing Projections

The Lambert conformal conical projection possesses several distinguishing features for cartographic representation over wide latitudinal belts:

• Local conformality: Absence of angular distortion, a property not available in Delisle–Euler or equidistant conical maps.
• Uniformly minimized scale-variation: Bi-Lipschitz measure $\delta$ is minimized through choice of $n$, ensuring both meridional and parallel scale-factors remain within $0.87\%$ of unity.
• Analytical tractability: Closed-form solutions for all relevant functions facilitate direct calculus-based optimization.
• Superior empirical performance: The distortion metric $\delta$ for the Lambert scheme is lower than, or matching, every other evaluated classical conical projection—including the Delisle–Euler map—while also guaranteeing exact preservation of local angles.

For the mapping of the Russian Empire, this combination yields the optimal trade-off: strict conformality and nearly uniform scale, substantiated by detailed numerical comparison and analytic foundations (Miyachi et al., 13 Jan 2026).