Next Article in Journal
Holocene Vegetation Dynamics, Landscape Change and Human Impact in Western Ireland as Revealed by Multidisciplinary, Palaeoecological Investigations of Peat Deposits and Bog-Pine in Lowland Connemara
Next Article in Special Issue
The Spatiotemporal Characteristics and Dynamic Changes of Tidal Flats in Florida from 1984 to 2020
Previous Article in Journal
The Importance of Connected and Interspersed Urban Green and Blue Space for Biodiversity: A Case Study in Cork City, Ireland
Previous Article in Special Issue
The Utility of Sentinel-2 MSI Data to Estimate Wetland Vegetation Leaf Area Index in Natural and Rehabilitated Wetlands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Linear Scale Factors in Map Projections of an Ellipsoid

Faculty of Geodesy, University of Zagreb, Kačićeva 26, 10000 Zagreb, Croatia
Submission received: 24 August 2021 / Revised: 22 October 2021 / Accepted: 28 October 2021 / Published: 9 November 2021
(This article belongs to the Special Issue Early Career Scientists’ (ECS) Contributions to Geographies)

Abstract

:
The main problem in cartography is that it is not possible to map/project/transform a spherical or ellipsoidal surface into a plane without distortions. The distortions of areas, angles, and/or distances are immanent to all maps. It is known that scale changes from point to point, and at certain points, the scale usually depends on the direction. The local linear scale factor c is one of the most important indicators of distortion distribution in the theory of map projections. It is not possible to find out the values of the local linear scale factor c in directions of coordinate axes x and y immediately from the definition of c. To solve this problem, in this paper, we derive new formulae for the calculation of c for a rotational ellipsoid. In addition, we derive the formula for computing c in any direction defined by dx and dy. We also considered the position and magnitude of the extreme values of c and derived new formulae for a rotational ellipsoid.

1. Introduction

A map is a result of mapping data, usually from the Earth, a celestial body, or an imagined world to a plane of representation on a piece of paper or on a digital display such as a computer monitor. Usually, maps are created by transforming data to a spherical or ellipsoidal surface and then to a plane. The mapping from a curved surface into a plane is known as map projection and can take a variety of forms [1].
Since no map projection maintains the correct scale throughout, it is important to determine the extent to which it varies on a map. On a world map, qualitative distortion is evident to an eye familiar with maps, after noting the extent to which landmasses are improperly sized or out of shape, and the extent to which meridians and parallels do not intersect at right angles or are not spaced uniformly along a given meridian or given parallel. On maps of countries or even of continents, distortion may not be evident to the eye, but it becomes apparent upon careful measurement and analysis [2].
All map projections involve the distortion of areas, angles, and/or distances. The types of distortion can be controlled to preserve specific characteristics, but map projections must distort other characteristics of the represented object. The main problem in cartography is that it is not possible to map a spherical or ellipsoidal surface into a plane without distortions. Euler first proved as early as 1772 that a sphere cannot be mapped into a plane with zero distortion [3,4,5].
The principal (linear) scale PS is the ratio of the length in the plane of projection and its origin on the surface (sphere, ellipsoid) to be projected/mapped. Notice that the scale of a map is not the ratio of a distance on the map to the corresponding distance on the ground [6].
The PS is usually indicated on maps because it determines the general degree of reduction in the length on the map. On most maps, it is usually simply called ‘scale’ and is known as the map scale.
The scale changes from point to point, and usually depends on direction. This is the local scale. The local linear scale factor c is the ratio of the differential of the curve arc in the plane of projection and the differential of the corresponding curve arc on an ellipsoid or spherical surface (see details in Section 2).
Tissot’s indicatrix is often used to illustrate the variation of point scale across a map. Each function can be locally approximated by a linear function. Each cartographic projection can be locally approximated by affine mapping. This was noticed by N. Tissot in the second half of the 19th century [7]. By mapping a circle of infinitesimal radius, an ellipse is obtained, which is called the distortion ellipse or Tissot’s indicatrix. The directions of this ellipse are called the main directions, and along them, the local linear scale is both the largest and the smallest.
A single indicatrix describes the distortion at a single point. Because distortion varies across a map, generally, Tissot’s indicatrices are placed across a map to illustrate the spatial change in distortion. A common scheme places them at each intersection of the displayed meridians and parallels (Figure 1). These schematics are important in the study of map projections, both to illustrate distortion and to provide the basis for the calculations that represent the magnitude of distortion precisely at each point.
Tissot’s indicatrices illustrate linear, angular, and areal distortions of maps, defined as follows:
  • A map distorts distances (linear distortion) wherever the quotient between the lengths of an infinitesimally short line as projected onto the projection surface and as it originally is on the Earth model deviates from 1.
  • A map distorts angles wherever the angles measured on the model of the Earth are not conserved in the projection. This is expressed by an ellipse of distortion, which is not a circle.
  • A map distorts areas wherever areas measured in the model of the Earth are not conserved in the projection. This is expressed by ellipses of distortion whose areas vary across the map.
In conformal maps, where each point preserves angles projected from the geometric model, the Tissot’s indicatrices are all circles whose size varies by location. In equal-area projections, where area proportions between objects are conserved, the Tissot’s indicatrices all have the same area, though their shapes and orientations vary with location. In arbitrary projections, both the area and shape vary across the map.
Since it is not possible to read a local linear scale factor in the direction of a coordinate axis immediately from the definition of a local linear scale factor, this paper considers the derivation of new formulae that enable local linear scale factors in the direction of coordinate x and y axes to be calculated for mapping a rotational ellipsoid. The formula for computing the local linear scale factor in any direction defined by dx and dy is also derived. Furthermore, the position and magnitude of the extreme values of the local linear scale factor are considered and new formulas are derived. The paper is a generalization to the ellipsoid of a previous paper developed for a sphere [8].

2. Ellipsoid, Map Projection, and Local Linear Scale Factor

This chapter does not contain novelties in the theory of cartographic projections, but it is necessary for an understanding of what follows. The geodetic parameterization of a rotational ellipsoid with semiaxes a and b and the center located in the origin of the coordinate system is a mapping defined by the following formulae [2]:
x = N cos φ cos λ ,   y = N cos φ sin λ ,   z = N ( 1 e 2 ) sin φ
where:
e 2 = a 2 b 2 a 2
N = a 1 e 2 sin 2 φ ,
φ [ π 2 , π 2 ] , λ [ π , π ] . In this case, φ is the latitude, and λ is the longitude. It is not difficult to derive that the first fundamental form of this mapping reads as:
d s 2 = M 2 d φ 2 + N 2 cos 2 φ d λ 2
where:
M = a ( 1 e 2 ) ( 1 e 2 sin 2 φ ) 3
A map projection is mapping given by the functions:
x = x ( φ , λ ) ,   y = y ( φ , λ ) ,
where the geodetic coordinates are φ [ π 2 , π 2 ] , λ [ π , π ] , as usual, and x and y are the coordinates of a point in a rectangular (mathematical, right-oriented) coordinate system in a plane. The first fundamental form of such a mapping is [9]:
d s 2 = E d φ 2 + 2 F d φ d λ + G d λ 2 ,
with the coefficients:
E = ( x φ ) 2 + ( y φ ) 2 ,   F = x φ x λ + y φ y λ ,   G = ( x λ ) 2 + ( y λ ) 2 .
Figure 2 on the left-hand side shows a differential quadrangel ABCD on the ellipsoid, while on the right-hand side is the image A′B′C′D′ of that quadrangle in the plane of projection. The following labels were used:
ds: differential arc on the ellipsoid surface;
ds′: image of ds in the plane of projection;
α: the angle between the differential arc on the ellipsoid surface and a meridian;
α′: image of α;
θ: the angle between the images of a meridian and a parallel.
The local linear scale factor c for mapping (3) of a rotational ellipsoid is usually defined in the theory of map projections by using the following relation:
c 2 = d s 2 d s 2 = E d φ 2 + 2 F d φ d λ + G d λ 2 M 2 d φ 2 + N 2 cos 2 φ d λ 2
which can also be written as follows in [2,3]:
c 2 ( α ) = E M 2 cos 2 α + F M N cos φ sin 2 α + G N 2 cos 2 φ sin 2 α
where:
tan α = N cos φ d λ M d φ
The poles are singular points of geographic parameterization (1), and therefore, expression (9) and all subsequent ones should be interpreted in the poles as limiting cases when φ π 2 or φ π 2 .
If α = 0 or, more generally, α = z π , z Z , where Z denotes the set of all integers, then the local linear scale factor c along a meridian ( d λ = 0 ) is:
h = c ( d λ = 0 ) = E M
and if α = π 2 or, more generally, α = π 2 + z π , z Z , then the local linear scale factor c along a parallel ( d φ = 0 ) is given by:
k = c ( d φ = 0 ) = G N cos φ .

3. Local Linear Scale Factors in the Directions of Coordinate Axes

It is not possible to read a local linear scale factor in the direction of a coordinate axis immediately from the definition of local linear scale factor (9). The same is true for Equation (10), where α denotes the azimuth, i.e., the angle between a meridian and any direction in a point in question. To be able to obtain a local linear scale factor in a direction defined by dx and dy, we need to modify Equation (9) or (10) in the appropriate way.
Let us start with the general Equations (6) of a map projection. Then, we can write the following:
d x = x φ d φ + x λ d λ ,   d y = y φ d φ + y λ d λ .
From (14), we have:
d φ = 1 H ( y λ d x x λ d y ) ,   d λ = 1 H ( y φ d x x φ d y ) ,
where:
H = E G F 2 = | y φ x λ y λ x φ |
and
H > 0 .
If we suppose that:
d y = 0
then:
d φ = 1 H y λ d x ,   d λ = 1 H y φ d x ,
and by substituting (19) in (9), we obtain the local linear scale factor in the direction of the x-axis as follows:
c ( d y = 0 ) = H M 2 ( y λ ) 2 + N 2 cos 2 φ ( y φ ) 2
If we suppose that:
d x = 0
then:
d φ = 1 H x λ d y ,   d λ = 1 H x φ d y ,
and by substituting (19) in (9), we obtain the local linear scale factor in the direction of the y-axis as follows:
c ( d x = 0 ) = H M 2 ( x λ ) 2 + N 2 cos 2 φ ( x φ ) 2
From the general theory of map projections, it is known that the local linear scale factor can be visualized as an ellipse. It is the ellipse of distortion or the Tissot’s indicatrix. Figure 3 represents a general case of Tissot’s indicatrix with local linear scale factors along a meridian c ( d λ = 0 ) , along a parallel c ( d φ = 0 ) , in the direction of the x-axis c ( d y = 0 ) , in the direction of the y-axis c ( d x = 0 ) , and the extremal values c m i n and c m a x .

4. Local Linear Scale Factor in a Given Direction

Let us suppose that we need a local linear scale factor in a given direction. If the direction is defined by d φ and d λ , we can use Equation (9), and our problem will be solved. If the direction is defined by dx and dy, we can use the following procedure. Let us denote that ψ is known as a meridian convergence:
tan ψ = d y d x .
Then:
cos ψ d y = sin ψ d x
and by using (15), (9) can be transformed into:
c 2 = H 2 ( d x 2 + d y 2 ) M 2 ( y λ d x x λ d y ) 2 + N 2 cos 2 φ ( y φ d x x φ d y ) 2 = H 2 a 1 cos 2 ψ + a 2 sin ψ cos ψ + a 3 sin 2 ψ 2 ,
where:
a 1 = M 2 ( y λ ) 2 + N 2 cos 2 φ ( y φ ) 2 , a 2 = 2 ( M 2 x λ y λ + N 2 cos 2 φ x φ y φ ) , a 3 = M 2 ( x λ ) 2 + N 2 cos 2 φ ( x φ ) 2 .
It follows that a local linear scale factor c in the direction ψ defined by (24) can be calculated by the formula:
c = H a 1 cos 2 ψ + a 2 sin ψ cos ψ + a 3 sin 2 ψ ,
where the coefficients a1, a2, and a3 are given by (27) and H by (16).
In a special case, when dy = 0, then ψ = 0 , and (28) reduces to (29):
c ( d y = 0 ) = H a 1 ,
which is equivalent to (20).
If dx = 0, then ψ = π 2 , and (28) takes the form of (30):
c ( d x = 0 ) = H a 3 ,
which is equivalent to (23).
The extremal values of c = c ( ψ ) given by (28) can be obtained in the usual way:
d c d ψ = 0
which gives:
tan 2 ψ = a 2 a 1 a 3
and by substituting (27) in (26), the extremal values of c are as follows:
c 1 , 2 = H 2 a 1 + a 3 ( a 1 a 3 ) 2 + a 2 2 = a 1 + a 3 ± ( a 1 a 3 ) 2 + a 2 2 2 M N cos φ = = M 2 G + E N 2 cos 2 φ ± ( M 2 G + E N 2 cos 2 φ ) 2 4 H 2 M 2 N 2 cos 2 φ 2 M N cos φ = = h 2 + k 2 ± ( h 2 k 2 ) 2 + 4 h 2 k 2 cos 2 β 2 .
In (33), c 1 , 2 = h and k are the local linear scale factors along a meridian and a parallel, respectively (see (12) and (13)).
The angle β (Figure 1) is the angle between the images of a meridian and a parallel at one point (see [10]):
sin 2 β = H 2 E G ,   cos 2 β = F 2 E G .
The extremal values of the local linear scale factor (33) are the semiaxes of the Tissot indicatrix or the ellipse of distortion. In order to be able to visualize an ellipse of distortions, in addition to its size, the direction of its axes is also required. The formula is known from trigonometry:
tan 2 ψ = 2 tan ψ 1 tan 2 ψ .
If we express tan ψ from that formula, we obtain the following quadratic equation:
tan 2 ψ + 2 tan 2 ψ tan ψ 1 = 0 .
From (36), we see that tan ψ 1 tan ψ 2 = 1 holds true, which means that the directions ψ 1 and ψ 2 will be orthogonal.
Relation (33) is not sufficient to determine the direction of the axis of the ellipse because it gives the direction of both axes since tan 2 ψ = tan ( π 2 ψ ) = tan 2 ( π 2 ψ ) . So, we will solve Equation (36). The solutions of this quadratic equation are as follows:
tan ψ 1 , 2 = cot 2 ψ ± cot 2 2 ψ + 1 ,
i.e.,
tan ψ 1 , 2 = a 3 a 1 ± ( a 3 a 1 ) 2 + a 2 2 a 2 .
The question of which of the directions ψ 1 and ψ 2 corresponds to the major and which to the minor axis of the ellipse remains open. The larger semiaxis is calculated as:
c 1 = H 2 a 1 + a 3 ( a 1 a 3 ) 2 + a 2 2 = a 1 + a 3 + ( a 1 a 3 ) 2 + a 2 2 2 M N cos φ
and the corresponding direction as:
tan ψ 1 = a 3 a 1 + ( a 3 a 1 ) 2 + a 2 2 a 2 .
The smaller semiaxis is as follows:
c 2 = H 2 a 1 + a 3 + ( a 1 a 3 ) 2 + a 2 2 = a 1 + a 3 ( a 1 a 3 ) 2 + a 2 2 2 M N cos φ
with the corresponding direction:
tan ψ 2 = a 3 a 1 ( a 3 a 1 ) 2 + a 2 2 a 2   or   ψ 2 = ψ 1 + π 2 .
We can now write Equation (28) in the following form:
c ( ψ ) = c 1 c 2 c 1 2 cos 2 ( ψ ψ 1 ) + c 2 2 sin 2 ( ψ ψ 1 ) ,
which is the equation of an ellipse in the polar coordinate system. In parametric form in the coordinate system x, y will be:
x = c ( ψ ) sin ( ψ ψ 1 ) ,   y = c ( ψ ) cos ( ψ ψ 1 )
and in canonical form:
x 2 c 1 2 + y 2 c 2 2 = 1 .

5. Examples

All derived formulas in the previous sections will be applied to several map projections to illustrate their functionality and validity.

5.1. Local Linear Scale Factors in the Mercator Projection

The equations for the normal aspect conformal cylindrical or Mercator projection of an ellipsoid are as follows (according to [2], but in different notation):
x = n ( λ λ 0 ) ,   y = n [ tanh 1 ( sin φ ) e tanh 1 ( e   sin φ ) ]
where λ 0 [ π , π ] represents the longitude of the central meridian, n > 0 , and e is defined by (2). From (46) we can obtain partial derivatives:
x λ = n ,   x φ = 0 ,   y λ = 0 ,   y φ = n M N cos φ
and then:
H = M n 2 N cos φ ,   a 1 = a 3 = M 2 n 2 ,   a 2 = 0
and by using Formulae (12), (13), (29), and (30):
c ( φ ) = c ( d λ = 0 ) = c ( d φ = 0 ) = c ( d x = 0 ) = c ( d y = 0 ) = n N cos φ
as expected, because the Mercator projection is conformal.
Tissot’s indicatrices are circles (Figure 4) whose radii do not depend on the longitude λ and increase with latitude φ or with the ordinate y = y ( φ ) . The standard parallel corresponds to the latitude for which N cos φ = n . For example, if n = a , then the equator is the standard parallel.

5.2. Local Linear Scale Factors in the Transverse Mercator Projection

The equations for the transverse aspect conformal cylindrical or transverse Mercator projection of a rotational ellipsoid are [12] as follows:
x = l N cos φ + l 3 6 N cos 3 φ ( 1 tan 2 φ ) +
y = S φ + l 2 2 N sin φ cos φ + ,
where l = λ λ 0 , λ 0 represents the longitude of the central meridian, and S φ is the length of the meridian arc from the equator to the point in question. From (49), we can obtain the following partial derivatives:
x λ = N cos φ ( 1 + l 2 2 cos 2 φ ) ,   y λ = l N sin φ cos φ .
The transverse Mercator projection is a conformal projection, and the Cauchy-Riemann conditions for mapping the ellipsoid read as follows:
x φ = M N cos φ y λ ,   y φ = M N cos φ x λ
and then:
H = M N cos φ ( 1 + l 2 cos 2 φ + l 4 4 cos 2 2 φ ) M N cos φ ( 1 + l 2 2 cos 2 φ ) 2 ,
a 1 = a 3 = M 2 N 2 cos 2 φ ( 1 + l 2 cos 2 φ + l 4 4 cos 2 2 φ ) = M 2 N 2 cos 2 φ ( 1 + l 2 2 cos 2 φ ) 2 , a 2 = 0
and then by using Formulae (12), (13), (29), and (30):
c ( φ , λ ) = c ( d λ = 0 ) = c ( d φ = 0 ) = c ( d x = 0 ) = c ( d y = 0 ) 1 + l 2 2 cos 2 φ
as expected because the transverse Mercator projection is also conformal. Tissot’s indicatrices are circles (Figure 5) whose radii increase with the longitude λ . From (52), we see that the local linear scale factor will be equal to 1 for l = 0 . Therefore, the central meridian of the mapping area is the standard meridian.

5.3. Local Linear Scale Factors in the Web-Mercator Projection

The equations for the Web-Mercator projection of an ellipsoid are [14] as follows:
x = a λ ,   y = a tanh 1 ( sin φ ) ,
where φ and λ are the ellipsoidal latitude and longitude, respectively. From (53), we can obtain partial derivatives:
x λ = a ,   x φ = 0 ,   y λ = 0 ,   y φ = a cos φ
and then:
E = a 2 cos 2 φ ,   F = 0 ,   G = a 2 ,   H = a 2 cos φ
a 1 = N 2 a 2 ,   a 2 = 0 ,   a 3 = M 2 a 2 .
By using Formulae (12), (13), (29), and (30), we have the following:
h = h ( φ ) = c ( d λ = 0 ) = E M = a M cos φ = ( 1 e 2 sin 2 φ ) 3 ( 1 e 2 ) cos φ
k = k ( φ ) = c ( d φ = 0 ) = G N cos φ = a N cos φ = 1 e 2 sin 2 φ cos φ
c ( d x = 0 ) = a M cos φ = h
c ( d y = 0 ) = a N cos φ = k
Furthermore:
h k = 1 e 2 sin 2 φ 1 e 2
and
( h k ) m a x = 1 1 e 2 1 + e 2
which is achieved for φ = 0 . The Web-Mercator projection is not conformal. Tissot’s indicatrices are ellipses very close to circles. Their semiaxes do not depend on the longitude, and the largest ratio of these semiaxes is at the equator and is approximately 1 + e 2 = 1.0067 . There is no standard parallel, i.e., a parallel along which h = k = 1 would be valid [14]. The distribution of distortions is very similar to that of the Mercator projection (Figure 4). The difference cannot be seen on a small-scale map.

5.4. Local Linear Scale Factors in the Albers Equal-Area Conic Projection

The equations for the normal aspect Albers or equal-area conic projection of an ellipsoid are the following [2]:
x = ρ sin θ ,   y = ρ 0 ρ cos θ
where φ and λ are the ellipsoidal latitude and longitude, respectively, calculated as follows:
ρ = a n C n q ( φ ) ,   θ = n ( λ λ 0 )
q ( φ ) = ( 1 e 2 ) ( sin φ 1 e 2 sin 2 φ 1 2 e ln 1 e sin φ 1 + e sin φ )
ρ 0 = a n C n q ( φ 0 )
m ( φ ) = cos φ 1 e 2 sin 2 φ
C = m 2 ( φ 1 ) + n q ( φ 1 ) , n = m 2 ( φ 1 ) m 2 ( φ 2 ) q ( φ 2 ) q ( φ 1 )
where φ 0 and λ 0 are the latitude and longitude, respectively, for the origin of the rectangular coordinates, and φ 1 and φ 2 are the latitudes of standard parallels.
From (62) and (63), we can obtain the following partial derivatives:
x λ = n ρ cos n ( λ λ 0 ) ,   x φ = M r n ρ sin n ( λ λ 0 ) ,
y λ = n ρ sin n ( λ λ 0 ) ,   y φ = M r n ρ cos n ( λ λ 0 ) ,
and then:
a 1 = M 2 n 2 ρ 2 [ n 4 ρ 4 sin 2 n ( λ λ 0 ) + r 4 cos 2 n ( λ λ 0 ) ]
a 2 = 2 M 2 n 2 ρ 2 ( n 4 ρ 4 r 4 ) sin n ( λ λ 0 ) cos n ( λ λ 0 )
a 3 = M 2 n 2 ρ 2 [ n 4 ρ 4 cos 2 n ( λ λ 0 ) + r 4 sin 2 n ( λ λ 0 ) ] ,
and then, by using Formulae (12), (13), (29), and (30), we can obtain:
h = c ( d λ = 0 ) = 1 k = m C n q ,
k = c ( d φ = 0 ) = n ρ r = C n q m
c ( d y = 0 ) = n r ρ n 4 ρ 4 sin 2 n ( λ λ 0 ) + r 4 cos 2 n ( λ λ 0 ) ,
c ( d x = 0 ) = n r ρ n 4 ρ 4 cos 2 n ( λ λ 0 ) + r 4 sin 2 n ( λ λ 0 ) .
According to (39) and (41), the larger and smaller semiaxes of the Tissot’s indicatrices are:
c 1 = c 1 ( φ ) = n ρ r = a C n q ( φ ) N cos φ = k ,   c 2 = c 2 ( φ ) = r n ρ = N cos φ a C n q ( φ ) = h .
We can see that the local linear scale factors (65)–(66) are different and depend on the direction. The distribution of distortions is depicted in Figure 6, although the choice of standard parallels for the map in that Figure is not in accordance with the theory [10,15].

6. Instead of Conclusions

All map projections involve distortion of areas, angles, and/or distances. The types of distortion can be controlled to preserve specific characteristics, but map projections must distort other characteristics of the represented object. The main problem in cartography is that it is not possible to map/project/transform a spherical or ellipsoidal surface into a plane with zero distortions.
Noting distortions is not enough. Cartographers need a way to measure the amount of distortion in a given projection. Some specific computational procedure must be followed. Such a procedure, typically described by a mathematical formula, is called the distortion measure. Therefore, in order to evaluate any kind of distortions, we need some kind of measure. There is no perfect choice because every map projection uniquely alters some aspect of space during the transformation process. Detailed information about the type, amount, and distribution of distortion is essential for choosing the best projection for a particular map or data set [17]. The distortion inherent in projections can be measured and symbolized much like any other map variable [18]. Map projections and the distortions they carry need to be better understood by spatial data developers, distributors, and users. Map distortions should be included along with map data as confidence layers, and easily accessible distortion displays should be available to help in the selection of map projections.
It is well known that scale changes from point to point, and at certain points, usually depends on direction. This is the local scale. The local linear scale factor c is the ratio of the differential of the curve arc in the plane of projection and the differential of the corresponding curve arc on an ellipsoid or spherical surface. The local linear scale factor c is one of the most important indicators of distortion distribution in the theory of map projections.
Knowing the distribution of distortion is important when choosing a map projection [19]. In any software that serves to select a map projection, as well as in any GIS, distortion display options should be built in. Formulae for determining the local linear scale factor in a direction of a coordinate axis are especially important in working with raster data [8,20,21,22].
It is not possible to read the local linear scale factor in a direction of a coordinate axis immediately from the definition (9). The same is true for Equation (10), where α is the angle between the meridian and any direction in a point in question. In the paper, we derive new formulae that enable calculation of a local linear scale factor in the direction of coordinate axes x and y for a rotational ellipsoid. Moreover, we derive the formula for computing the local linear scale factor in any direction defined by dx and dy. The position and magnitude of the extreme values of the local linear scale factor were also considered, and new formulae were derived for a rotational ellipsoid.

Funding

This research received no external funding.

Data Availability Statement

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Lapaine, M. Map projection article on Wikipedia. Adv. Cartogr. GIScience ICA 2019, 1, 1–8. [Google Scholar] [CrossRef] [Green Version]
  2. Snyder, J.P. Map Projections–A Working Manual: U.S. Geological Survey Professional Paper 1395; United States Government Printing Office: Washington, DC, USA, 1987. [Google Scholar]
  3. Euler, L. De repraesentatione superficiei sphaericae super plano. Acta Acad. Sci. Imp. Petropolitanae 1778, 1777, 107–132. [Google Scholar]
  4. Biernacki, F. Teoria Odwzorowań Powierzchni dla Geodetów i Kartografów; Główny Urząd Pomiarów Kraju, Prace Geodezyjnego Instytutu Naukowo-Badawczego: Warsaw, Poland, 1949; Translated into English as Theory of Representation of Surfaces for Surveyors and Cartographers: U.S. Dept. of Commerce, 1965. [Google Scholar]
  5. Lapaine, M. Kartografske projekcije i njihove deformacije (in Croatian). In Proceedings of the Peti hrvatski kongres o katastru, Zagreb, Croatia, 8–9 May 2014; pp. 15–32, ISBN 978-953-97081-9-9. [Google Scholar]
  6. Lapaine, M. Mjerilo i istinito mjerilo?/Scale and True Scale? Kartogr. Geoinformacije/Cartogr. Geoinf. 2021, 20, 110–111. [Google Scholar]
  7. Tissot, N.A. Mémoire sur la Représentation des Surfaces et les Projections des Cartes Géographiques; Gauthier Villars: Paris, France, 1881. [Google Scholar]
  8. Lapaine, M.; Frančula, N.; Tutek, Ž. Local linear scale factors in map projections in the direction of coordinate axes. Geo-Spat. Inf. Sci. 2021, 1–8. [Google Scholar] [CrossRef]
  9. Canters, F.; Decleir, H. The World in Perspective, a Directory of World Map Projections; John Wiley & Sons: Chichester, UK, 1989. [Google Scholar]
  10. Bugayevskiy, L.M.; Snyder, J.P. Map Projections: A Reference Manual; Taylor and Francis: London, UK, 1995. [Google Scholar]
  11. Tissot Mercator. Available online: https://commons.wikimedia.org/wiki/File:Tissot_mercator.png (accessed on 8 August 2021).
  12. Thomas, P.D. Conformal Projections in Geodesy and Cartography; U.S. Coast and Geodetic Survey Spec. Pub. No. 251; United States Government Printing Office: Washington, DC, USA, 1952.
  13. Transverse Mercator with Tissot’s Indicatrices of Distortion. Available online: https://commons.wikimedia.org/wiki/File:Transverse_Mercator_with_Tissot%27s_Indicatrices_of_Distortion.svg (accessed on 8 August 2021).
  14. Lapaine, M.; Frančula, N. Web mercator projection—One of cylindrical projections of an ellipsoid to a plane/web-mercatorova projekcija–jedna od cilindričnih projekcija. Kartogr. Geoinformacije/Cartogr. Geoinf. 2021, 20, 30–47. [Google Scholar]
  15. Каврайский, В.В. Избранные труды. Т. 2., Вып. 2. Математическая Картoграфия. Кoнические и Цилиндрические Прoекции, их Применение; Управления начальника Гидрoграфическoй службы ВМФ: Мoсква, Russia, 1959. [Google Scholar]
  16. The World on an Albers Projection of the Earth’s Sphere, with 10° Graticule and Tissot’s Indicatrices Overlaid. Standard Parallels are at 45° N and 15° N. Available online: https://commons.wikimedia.org/wiki/File:Albers_with_Tissot%27s_Indicatrices_of_Distortion.svg (accessed on 8 August 2021).
  17. Laskowski, P. Part 1: Distortion-spectrum fundamentals—A new tool for analyzing and visualizing map distortions. Cartographica 1997, 34, 3–18. [Google Scholar] [CrossRef] [Green Version]
  18. Mulcahy, K.A.; Clarke, K.C. Symbolization of map projection distortion: A review. Cartogr. Geogr. Inf. Sci. 2001, 28, 167–181. [Google Scholar] [CrossRef]
  19. Lapaine, M.; Usery, E.L. (Eds.) Choosing a Map Projection; Lecture Notes in Geoinformation and Cartography; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  20. Seong, J.C.; Mulcahy, K.A.; Usery, E.L. The sinusoidal projection: A new importance in relation to global image data. Prof. Geogr. 2002, 54, 218–225. [Google Scholar] [CrossRef]
  21. Seong, J.C.; Usery, E.L. Assessing raster representation accuracy using a scale factor model. Photogramm. Eng. Remote. Sens. 2001, 67, 1185–1191. [Google Scholar]
  22. Usery, E.L.; Seong, J.C. All equal-area map projections are created equal, but some are more equal than others. Cartogr. Geogr. Inf. Sci. 2001, 28, 183–194. [Google Scholar] [CrossRef]
Figure 1. The Hammer projection with Tissot’s indicatrices.
Figure 1. The Hammer projection with Tissot’s indicatrices.
Geographies 01 00014 g001
Figure 2. An ellipsoidal differential quadrangle (left) and its image in the plane of projection (right).
Figure 2. An ellipsoidal differential quadrangle (left) and its image in the plane of projection (right).
Geographies 01 00014 g002
Figure 3. General case of Tissot’s indicatrix showing local linear scale factors along a meridian h = c ( d λ = 0 ) , along a parallel k = c ( d φ = 0 ) , in the direction of the x-axis c ( d y = 0 ) , in the direction of the y-axis c ( d x = 0 ) , and extremal values c m i n and c m a x . The angle β is the angle between the images of a meridian and a parallel at a point under consideration.
Figure 3. General case of Tissot’s indicatrix showing local linear scale factors along a meridian h = c ( d λ = 0 ) , along a parallel k = c ( d φ = 0 ) , in the direction of the x-axis c ( d y = 0 ) , in the direction of the y-axis c ( d x = 0 ) , and extremal values c m i n and c m a x . The angle β is the angle between the images of a meridian and a parallel at a point under consideration.
Geographies 01 00014 g003
Figure 4. Mercator projection map of the earth’s sphere with Tissot’s indicatrices. Source [11].
Figure 4. Mercator projection map of the earth’s sphere with Tissot’s indicatrices. Source [11].
Geographies 01 00014 g004
Figure 5. The world on a transverse Mercator projection of the Earth’s sphere, with 10° graticule and Tissot’s indicatrices overlaid. Source [13].
Figure 5. The world on a transverse Mercator projection of the Earth’s sphere, with 10° graticule and Tissot’s indicatrices overlaid. Source [13].
Geographies 01 00014 g005
Figure 6. The world on an Albers projection of the Earth’s sphere, with 10° graticule and Tissot’s indicatrices overlaid. Standard parallels are at 45° N and 15° N. Source [16].
Figure 6. The world on an Albers projection of the Earth’s sphere, with 10° graticule and Tissot’s indicatrices overlaid. Standard parallels are at 45° N and 15° N. Source [16].
Geographies 01 00014 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lapaine, M. Local Linear Scale Factors in Map Projections of an Ellipsoid. Geographies 2021, 1, 238-250. https://0-doi-org.brum.beds.ac.uk/10.3390/geographies1030014

AMA Style

Lapaine M. Local Linear Scale Factors in Map Projections of an Ellipsoid. Geographies. 2021; 1(3):238-250. https://0-doi-org.brum.beds.ac.uk/10.3390/geographies1030014

Chicago/Turabian Style

Lapaine, Miljenko. 2021. "Local Linear Scale Factors in Map Projections of an Ellipsoid" Geographies 1, no. 3: 238-250. https://0-doi-org.brum.beds.ac.uk/10.3390/geographies1030014

Article Metrics

Back to TopTop