This article was automatically translated from the original Turkish version.
Map projection is the process of transforming a curved reference surface such as a sphere or ellipsoid defined for the Earth’s surface onto a plane through a geometric transformation. Mathematically, it is defined by a pair of functions that establish the relationship between the reference surface coordinate system and the plane coordinate system. This transformation is known as the “forward transformation”; the transition from plane coordinates back to reference surface coordinates is called the “inverse transformation”.
Projecting the Earth’s surface onto a plane offers advantages in terms of map usage, perception, and storage. However, since the reference surfaces used for the Earth are closed surfaces, their transfer to a plane is not a simple scaling problem.
Since it is geometrically impossible to unfold a curved surface such as a sphere or ellipsoid onto a plane, changes in area, angle, and length are unavoidable during the transformation process. These changes are termed “distortion” or “deformation”.
It is impossible to preserve all three fundamental cartographic properties—area, angle, and length—simultaneously. Therefore, map projections have been developed to preserve one of these properties according to the intended use of the map.
Distortions are visualized and analyzed using a method known as Tissot’s Indicatrix (or distortion ellipse). It is based on the projection of an infinitesimally small circle on the original surface as an ellipse on the projection plane. The major axis (a) and minor axis (b) of this ellipse represent the maximum and minimum length distortions at that point.
The fundamental classifications of projections according to the properties they preserve are as follows:
Map projections are classified according to the projection surface used, the orientation of that surface, its contact condition with the reference surface, and the property preserved.
A projection can be made directly from the reference surface onto a plane, or it can be performed via auxiliary surfaces that can be unfolded into a plane.

Planar Projection (Generated by Artificial Intelligence)

Cylindrical Projection (Generated by Artificial Intelligence)

Conic Projection (Generated by Artificial Intelligence)
This refers to the position of the projection surface (plane, cylinder, or cone axis) relative to the reference surface.
The choice of an appropriate projection depends on factors such as the map’s intended use, target audience, location, shape, size of the area being mapped, and the map scale. The primary goal is to select the method that produces the smallest possible distortion for the region of interest.
Some widely accepted selection rules are as follows:
By Region Location: Azimuthal (planar) projections are preferred for polar regions, conic projections for mid-latitudes, and cylindrical projections for equatorial regions.
By Region Shape (Young’s Rule): If the area to be projected is approximately circular, a planar (azimuthal) projection is recommended. If the area is elongated (rectangular strip), a cylindrical or conic projection aligned with the longer dimension is preferred.
By Map Purpose: Conformal (angle-preserving) projections are generally preferred for large-scale topographic maps and medium-scale physical maps. For small-scale thematic maps covering large areas (e.g., population distribution), equal-area projections are recommended.
By Region Size:
Theoretically, an infinite number of map projections can be defined, and advances in information technology and Geographic Information Systems (GIS) have created a need for a standardized system of definition and coding.
The most widely used coding system in this field is EPSG. EPSG (European Petroleum Survey Group) has published a dataset containing geodetic parameters, coordinate reference systems (CRS), units of measurement, and ellipsoid definitions. These codes are known as SRID (Spatial Reference System Identifier).
Some commonly used EPSG codes are:
Map projections are grounded in mathematical principles. The complexity of projections necessitates the use of numerical analysis methods in modern cartography.
To calculate projection distortions (such as Tissot’s Indicatrix elements), the partial derivatives of the projection equations (∂x/∂φ, ∂y/∂λ, etc.) must be known. Analytically computing these derivatives can be complex. Numerical differentiation methods allow approximate values of these partial derivatives to be calculated directly from the function itself, without analytical derivation—for example, using <span class="katex"><span class="katex-html" aria-hidden="true"><span class="base"><span class="strut" style="height:1.0019em;vertical-align:-0.25em;"></span><span class="mord"><span class="mord mathnormal" style="margin-right:0.10764em;">f</span><span class="msupsub"><span class="vlist-t"><span class="vlist-r"><span class="vlist" style="height:0.7519em;"><span style="top:-3.063em;margin-right:0.05em;"><span class="pstrut" style="height:2.7em;"></span><span class="sizing reset-size6 size3 mtight"><span class="mord mtight"><span class="mord mtight">′</span></span></span></span></span></span></span></span></span><span class="mopen">(</span><span class="mord mathnormal">x</span><span class="mclose">)</span><span class="mspace" style="margin-right:0.2778em;"></span><span class="mrel">≈</span><span class="mspace" style="margin-right:0.2778em;"></span></span><span class="base"><span class="strut" style="height:1.355em;vertical-align:-0.345em;"></span><span class="mord"><span class="mopen nulldelimiter"></span><span class="mfrac"><span class="vlist-t vlist-t2"><span class="vlist-r"><span class="vlist" style="height:1.01em;"><span style="top:-2.655em;"><span class="pstrut" style="height:3em;"></span><span class="sizing reset-size6 size3 mtight"><span class="mord mtight"><span class="mord mtight">2</span><span class="mord mathnormal mtight">t</span></span></span></span><span style="top:-3.23em;"><span class="pstrut" style="height:3em;"></span><span class="frac-line" style="border-bottom-width:0.04em;"></span></span><span style="top:-3.485em;"><span class="pstrut" style="height:3em;"></span><span class="sizing reset-size6 size3 mtight"><span class="mord mtight"><span class="mord mathnormal mtight" style="margin-right:0.10764em;">f</span><span class="mopen mtight">(</span><span class="mord mathnormal mtight">x</span><span class="mbin mtight">+</span><span class="mord mathnormal mtight">t</span><span class="mclose mtight">)</span><span class="mbin mtight">−</span><span class="mord mathnormal mtight" style="margin-right:0.10764em;">f</span><span class="mopen mtight">(</span><span class="mord mathnormal mtight">x</span><span class="mbin mtight">−</span><span class="mord mathnormal mtight">t</span><span class="mclose mtight">)</span></span></span></span></span><span class="vlist-s"></span></span><span class="vlist-r"><span class="vlist" style="height:0.345em;"><span></span></span></span></span></span><span class="mclose nulldelimiter"></span></span></span></span></span>. A small parameter value such as t=10−6 has been found to yield satisfactory results for map projection equations.
Some pseudo-cylindrical projections (e.g., Mollweide projection) do not define a direct relationship between latitude (φ) and the plane coordinate (y), but instead use an auxiliary variable (θ) defined by an implicit function. For example, for Mollweide, the equation <span class="katex"><span class="katex-html" aria-hidden="true"><span class="base"><span class="strut" style="height:0.7778em;vertical-align:-0.0833em;"></span><span class="mord">2</span><span class="mord mathnormal" style="margin-right:0.02778em;">θ</span><span class="mspace" style="margin-right:0.2222em;"></span><span class="mbin">+</span><span class="mspace" style="margin-right:0.2222em;"></span></span><span class="base"><span class="strut" style="height:1em;vertical-align:-0.25em;"></span><span class="mop">sin</span><span class="mopen">(</span><span class="mord">2</span><span class="mord mathnormal" style="margin-right:0.02778em;">θ</span><span class="mclose">)</span><span class="mspace" style="margin-right:0.2778em;"></span><span class="mrel">=</span><span class="mspace" style="margin-right:0.2778em;"></span></span><span class="base"><span class="strut" style="height:1em;vertical-align:-0.25em;"></span><span class="mord mathnormal" style="margin-right:0.03588em;">π</span><span class="mspace" style="margin-right:0.1667em;"></span><span class="mop">sin</span><span class="mopen">(</span><span class="mord mathnormal">φ</span><span class="mclose">)</span></span></span></span> holds. Finding the value of θ for a given φ requires iterative numerical methods such as the Newton-Raphson method.
For some projections (e.g., Winkel Tripel), the inverse transformation equations that convert plane coordinates (x,y) back to geographic coordinates (φ,λ) form nonlinear equation systems that cannot be solved analytically. In such cases, the inverse transformation can only be performed using numerical methods such as Newton-Raphson. This approach, combined with numerical differentiation, enables the calculation of inverse projections using only the forward projection equations.
As a theoretical approach, cartograms differ from traditional map projections. Their purpose is not to preserve geographic area (in square meters) but to “adjust” the representation of a phenomenon (population, income, etc.) within that area. The goal of these transformations is to scale the area of a region on the map proportionally to a value assigned to it (e.g., population or retail sales volume). Mathematically, this can be expressed as a transformation where Tissot’s area distortion measure (S) is set equal to the given density distribution (D): <span class="katex"><span class="katex-html" aria-hidden="true"><span class="base"><span class="strut" style="height:0.6833em;"></span><span class="mord mathnormal" style="margin-right:0.02778em;">D</span><span class="mspace" style="margin-right:0.2778em;"></span><span class="mrel">=</span><span class="mspace" style="margin-right:0.2778em;"></span></span><span class="base"><span class="strut" style="height:0.6833em;"></span><span class="mord mathnormal" style="margin-right:0.05764em;">S</span></span></span></span>.
A normal-position, conformal cylindrical projection introduced by Gerardus Mercator in 1569. Its key feature is the accurate representation of loxodromes—curves that cross meridians at constant angles—as straight lines, preserving angles. This property made it useful for compass navigation. However, it causes extreme area distortion near the poles and is unsuitable for representing the entire Earth.
The Web Mercator variant (commonly EPSG:3857) is widely used in online mapping services such as Google Maps. Its conformal nature ensures that square or rectangular shapes like buildings remain undistorted at high zoom levels and that geographic grids intersect at right angles, making it favorable for computer graphics. Nevertheless, this adaptation has faced criticism for severe area distortion at global scales.
This is the application of the Mercator projection in transverse position using an ellipsoidal reference surface. It is a conformal projection used for large-scale mapping. The projection is constructed on a cylinder tangent to the ellipsoid along a chosen central meridian. To avoid increasing distortion away from the central meridian, the projection is applied in narrow “zones”.
A common problematic practice is the direct use of geographic coordinates (latitude, longitude) as if they were planar (x,y) coordinates in a map projection. This system, commonly known by the EPSG:4326 (WGS84) code, is primarily intended for data storage.
When these data are displayed on a map, the result is equivalent to an equidistant cylindrical projection (also called Plate Carrée). This projection is only suitable for regions near the equator. When used in mid-latitude regions such as Türkiye, it causes severe compression in the north-south direction, resulting in an inaccurate representation. Additionally, when horizontal units are in degrees and vertical units (if elevation data is present) are in meters, software such as QGIS produces incorrect results in three-dimensional operations like hillshading or hydrological analysis. For such analyses, data must be transformed into a metric CRS.
In Türkiye, standard projection systems defined by legal regulations and guidelines are used for large and medium-scale map production. Both standards are based on the conformal transverse cylindrical Gauss-Krüger projection.
TM (Transverse Mercator) - 3° Zones:
UTM (Universal Transverse Mercator) - 6° Zones:
Core Issue: Distortion
Classification of Projections
By Projection Surface (Auxiliary Surfaces)
By Surface Orientation
By Contact Condition (Points of Contact)
By Mathematical Structure
Projection Selection
Standardization and Coordinate Reference Systems (CRS)
Theoretical Approaches and Numerical Analysis
Numerical Differentiation
Iterative Solutions and Auxiliary Variables
Numerical Inverse Projection
Cartograms (Density-Preserving Projections)
Common Projections and Application Errors
Mercator Projection
Gauss-Krüger (Transverse Mercator) Projection
Direct Use of Geographic Coordinates (EPSG:4326)
Applications in Türkiye