Misplaced Pages

Mollweide projection: Difference between revisions

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Browse history interactively← Previous editContent deleted Content addedVisualWikitext
Revision as of 03:42, 19 September 2013 editStrebe (talk | contribs)Extended confirmed users8,417 edits Edits resulted in unparsable formulæ; please try again.← Previous edit Latest revision as of 16:45, 8 December 2024 edit undoStrebe (talk | contribs)Extended confirmed users8,417 edits It’s not about equivalence; it’s that it appears in records both ways. Undid revision 1261825089 by Haush2909 (talk)Tag: Undo 
(85 intermediate revisions by 54 users not shown)
Line 1: Line 1:
{{Short description|Pseudocylindrical equal-area map projection}}
]
] of deformation]]
The '''Mollweide projection''' is a ] generally used for global maps of the world (or sky). Also known as the '''Babinet projection''', '''homalographic projection''', '''homolographic projection''', and '''elliptical projection'''. As its more explicit name '''Mollweide equal area projection''' indicates, it sacrifices accuracy of angle and shape in favor of accurate proportions in area. It is used primarily where accurate representation of area takes precedence over shape, for instance small maps depicting global distributions.


]
The projection was first published by mathematician and astronomer ] (1774 &ndash; 1825) of ] in 1805. It was popularized by ] in 1857, giving it the name '''homalographic projection'''. The variation '''homolographic''' arose from frequent nineteenth century usage in star atlases.<ref>''Flattening the Earth: Two Thousand Years of Map Projections'', John P. Snyder, 1993, pp.&nbsp;112&ndash;113, ISBN 0-226-76747-7.</ref>
] of deformation]]


The '''Mollweide projection''' is an ], ] generally used for maps of the world or ]. It is also known as the '''Babinet projection''', '''homalographic projection''', '''homolographic projection''', and '''elliptical projection'''. The projection trades accuracy of angle and shape for accuracy of proportions in area, and as such is used where that property is needed, such as maps depicting global distributions.
] image (2012) of the ].<ref name="Space-20121221">{{cite web |last=Gannon |first=Megan |title=New 'Baby Picture' of Universe Unveiled |url=http://www.space.com/19027-universe-baby-picture-wmap.html|date=December 21, 2012 |publisher=] |accessdate=December 21, 2012 }}</ref><ref name="arXiv-20121220">{{cite journal |last=Bennett |first=C.L. |last2=Larson |first2=L.|last3=Weiland |first3=J.L. |last4=Jarosk |first4= N. |last5=Hinshaw |first5=N. |last6=Odegard|first6=N. |last7=Smith |first7=K.M. |last8=Hill |first8=R.S. |last9=Gold |first9=B.|last10=Halpern |first10=M. |last11=Komatsu |first11=E. |last12=Nolta |first12=M.R.|last13=Page |first13=L. |last114=Spergel |first14=D.N. |last15=Wollack |first15=E.|last16=Dunkley |first16=J. |last17=Kogut |first17=A. |last18=Limon |first18=M. |last19=Meyer|first19=S.S. |last20=Tucker |first20=G.S. |last21=Wright |first21=E.L. |title=Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results|url=http://arxiv.org/abs/1212.5225 |arxiv=1212.5225 |date=December 20, 2012|accessdate=December 22, 2012 }}</ref> Projected using the Mollweide projection.]]

The projection was first published by mathematician and astronomer ] (1774–1825) of ] in 1805. It was reinvented and popularized in 1857 by ], who gave it the name '''homalographic projection'''. The variation '''homolographic''' arose from frequent nineteenth-century usage in star atlases.<ref>''Flattening the Earth: Two Thousand Years of Map Projections'', John P. Snyder, 1993, pp.&nbsp;112&ndash;113, {{ISBN|0-226-76747-7}}.</ref>

] image (2012) of the ].<ref name="Space-20121221">{{cite web |last=Gannon |first=Megan |title=New 'Baby Picture' of Universe Unveiled |url=http://www.space.com/19027-universe-baby-picture-wmap.html|date=December 21, 2012 |publisher=] |access-date=December 21, 2012 }}</ref><ref name="arXiv-20121220">{{cite journal |last1=Bennett |first1=C.L. |last2=Larson |first2=L.|last3=Weiland |first3=J.L. |last4=Jarosk |first4= N. |last5=Hinshaw |first5=N. |last6=Odegard|first6=N. |last7=Smith |first7=K.M. |last8=Hill |first8=R.S. |last9=Gold |first9=B.|last10=Halpern |first10=M. |last11=Komatsu |first11=E. |last12=Nolta |first12=M.R.|last13=Page |first13=L. |last14=Spergel |first14=D.N. |last15=Wollack |first15=E.|last16=Dunkley |first16=J. |last17=Kogut |first17=A. |last18=Limon |first18=M. |last19=Meyer|first19=S.S. |last20=Tucker |first20=G.S. |last21=Wright |first21=E.L. |title=Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results|journal=]|volume=208|issue=2|page=20|year=2013|arxiv=1212.5225|bibcode=2013ApJS..208...20B|doi=10.1088/0067-0049/208/2/20|s2cid=119271232 }}</ref> Projected using the Mollweide projection.]]


] levels measured by the ]. Projected using the Mollweide projection.]] ] levels measured by the ]. Projected using the Mollweide projection.]]
The projection is:


==Properties==
:<math>x = \frac{2 \sqrt 2}{\pi} \lambda \cos\left(\theta \right),</math>
The Mollweide is a ] projection in which the ] is represented as a straight horizontal line perpendicular to a central ] that is one-half the equator's length. The other ]s compress near the poles, while the other meridians are equally spaced at the equator. The meridians at 90 degrees east and west form a perfect circle, and the whole earth is depicted in a proportional 2:1 ellipse. The proportion of the area of the ellipse between any given parallel and the equator is the same as the proportion of the area on the globe between that parallel and the equator, but at the expense of shape distortion, which is significant at the perimeter of the ellipse, although not as severe as in the ].


Shape distortion may be diminished by using an ''interrupted'' version. A ''sinusoidal interrupted'' Mollweide projection discards the central meridian in favor of alternating half-meridians which terminate at right angles to the equator. This has the effect of dividing the globe into lobes. In contrast, a ''parallel interrupted'' Mollweide projection uses multiple disjoint central meridians, giving the effect of multiple ellipses joined at the equator. More rarely, the projection can be drawn obliquely to shift the areas of distortion to the oceans, allowing the continents to remain truer to form.
:<math>y = \sqrt 2 \sin\left(\theta \right),\,</math>


The Mollweide, or its properties, has inspired the creation of several other projections, including the ], ] and the ].<ref>, ] Professional Paper 1395, John P. Snyder, 1987, pp.&nbsp;249&ndash;252</ref>
where <math>\theta\,</math> is an auxiliary angle defined by


==Mathematical formulation==
:<math>2 \theta + \sin(2 \theta) = \pi \sin(\varphi)\qquad (1)</math>
The projection transforms from latitude and longitude to map coordinates ''x'' and ''y'' via the following equations:<ref name="MathWorldMollweide">{{MathWorld|title=Mollweide Projection|urlname=MollweideProjection|author=Weisstein, Eric W.}}</ref>


:<math>\begin{align} x &= R \frac{2 \sqrt 2}{\pi} \left( \lambda - \lambda_{0} \right) \cos \theta, \\ y &= R \sqrt 2 \sin \theta ,\end{align}</math>
and ''λ'' is the longitude from the central meridian, and ''φ'' is the latitude.


where ''θ'' is an auxiliary angle defined by
Equation (1) may be solved with rapid convergence (but slow near the poles) using ] iteration:


:<math> \theta_0 = \varphi,\,</math> :<math>2\theta + \sin 2\theta = \pi \sin \varphi \qquad (1)</math>
:<math> \theta_{n+1} = \theta_n - \frac{(2\theta_n + \sin(2\theta_n) - \pi \sin(\varphi))}{2 + 2\cos(2\theta_n)}.\,</math>


and ''λ'' is the longitude, ''λ''{{sub|0}} is the central meridian, ''φ'' is the latitude, and ''R'' is the radius of the globe to be projected. The map has area 4{{pi}}''R''{{sup|2}}, conforming to the surface area of the generating globe. The ''x''-coordinate has a range of , and the ''y''-coordinate has a range of .
If ''φ'' = ±π/2, then also θ = ±π/2. In that case the iteration should be bypassed; otherwise, ] may result.


Equation (1) may be solved with rapid convergence (but slow near the poles) using ] iteration:<ref name="MathWorldMollweide"/>
==Properties==
The Mollweide is a ] projection in which the ] is represented as a straight horizontal line perpendicular to a central ] one-half its length. The other ]s compress near the poles, while the other meridians are equally spaced at the equator. The meridians at 90 degrees east and west form a perfect circle, and the whole earth is depicted in a proportional 2:1 ellipse. The proportion of the area of the ellipse between any given parallel and the equator is the same as the proportion of the area on the globe between that parallel and the equator, but at the expense of shape distortion, which is significant at the perimeter of the ellipse, although not as severe as in the ].


:<math>\begin{align} \theta_0 &= \varphi, \\ \theta_{n+1} &= \theta_n - \frac{2 \theta_n + \sin 2\theta_n - \pi \sin \varphi}{2 + 2 \cos 2\theta_n}.\end{align}</math><ref group=note>The formula in the text helps the reader confirm that the formula is correct. For ] the denominator should be changed, starting with the double angle identity.
Shape distortion may be diminished by using an ''interrupted'' version. A ''sinusoidal interrupted'' Mollweide projection discards the central meridian in favor of alternating half-meridians which terminate at right angles to the equator. This has the effect of dividing the globe into lobes shape. In contrast, a ''parallel interrupted'' Mollweide projection uses multiple disjoint central meridians, giving the effect of multiple ellipses joined at the equator. More rarely, the project can be drawn obliquely to shift the areas of distortion to the oceans, allowing the continents to remain truer to form.


:<math>\begin{align} \cos 2\theta_n &= 2 \cos^2 \theta_n - 1, \\ 1 + \cos 2\theta_n &= 2 \cos^2 \theta_n, \\ 2 + 2 \cos 2\theta_n &= 4 \cos^2 \theta_n, \\ \theta_{n+1} &= \theta_n - \frac{2 \theta_n + \sin 2\theta_n - \pi \sin \varphi}{4 \cos^2 \theta_n}.\end{align}</math>
The Mollweide, or its properties, has inspired the creation of several other projections, including the ], ] and the ].<ref>, ] Professional Paper 1395, John P. Snyder, 1987, pp.&nbsp;249&ndash;252</ref>

In numerical computation, the original denominator could result in zero for ''θ'' near ±{{sfrac|{{pi}}|2}} (catastrophic cancellation). This substitution is true for all angles and avoids the problem near ''θ'' = ±{{sfrac|{{pi}}|2}} without making it a special case.</ref>


If ''φ'' = ±{{sfrac|{{pi}}|2}}, then also ''θ'' = ±{{sfrac|{{pi}}|2}}. In that case the iteration should be bypassed; otherwise, ] may result.

There exists a ] inverse transformation:<ref name="MathWorldMollweide"/>

:<math>\begin{align} \varphi &= \arcsin \frac{2 \theta + \sin 2\theta}{\pi}, \\ \lambda &= \lambda_{0} + \frac{\pi x}{2 R \sqrt{2} \cos \theta}, \end{align}</math>

where ''θ'' can be found by the relation

:<math>\theta = \arcsin \frac{y}{R \sqrt{2}}. \,</math>

The inverse transformations allow one to find the latitude and longitude corresponding to the map coordinates ''x'' and ''y''.

==Alterations==
]]]


==See also== ==See also==

* ]
* ] * ]
* ] * ]
* ] family

==Notes==
{{reflist|group=note}}


==References== ==References==
Line 42: Line 68:


==External links== ==External links==
{{Commonscat|Mollweide projection}} {{Commons category|Mollweide projection}}
* *
* *


{{Map Projections}} {{Map projections}}


]
] ]

Latest revision as of 16:45, 8 December 2024

Pseudocylindrical equal-area map projection
Mollweide projection of the world
The Mollweide projection with Tissot's indicatrix of deformation

The Mollweide projection is an equal-area, pseudocylindrical map projection generally used for maps of the world or celestial sphere. It is also known as the Babinet projection, homalographic projection, homolographic projection, and elliptical projection. The projection trades accuracy of angle and shape for accuracy of proportions in area, and as such is used where that property is needed, such as maps depicting global distributions.

The projection was first published by mathematician and astronomer Karl (or Carl) Brandan Mollweide (1774–1825) of Leipzig in 1805. It was reinvented and popularized in 1857 by Jacques Babinet, who gave it the name homalographic projection. The variation homolographic arose from frequent nineteenth-century usage in star atlases.

Nine-year WMAP image (2012) of the cosmic microwave background radiation. Projected using the Mollweide projection.
Sea-surface freon levels measured by the Global Ocean Data Analysis Project. Projected using the Mollweide projection.

Properties

The Mollweide is a pseudocylindrical projection in which the equator is represented as a straight horizontal line perpendicular to a central meridian that is one-half the equator's length. The other parallels compress near the poles, while the other meridians are equally spaced at the equator. The meridians at 90 degrees east and west form a perfect circle, and the whole earth is depicted in a proportional 2:1 ellipse. The proportion of the area of the ellipse between any given parallel and the equator is the same as the proportion of the area on the globe between that parallel and the equator, but at the expense of shape distortion, which is significant at the perimeter of the ellipse, although not as severe as in the sinusoidal projection.

Shape distortion may be diminished by using an interrupted version. A sinusoidal interrupted Mollweide projection discards the central meridian in favor of alternating half-meridians which terminate at right angles to the equator. This has the effect of dividing the globe into lobes. In contrast, a parallel interrupted Mollweide projection uses multiple disjoint central meridians, giving the effect of multiple ellipses joined at the equator. More rarely, the projection can be drawn obliquely to shift the areas of distortion to the oceans, allowing the continents to remain truer to form.

The Mollweide, or its properties, has inspired the creation of several other projections, including the Goode's homolosine, van der Grinten and the Boggs eumorphic.

Mathematical formulation

The projection transforms from latitude and longitude to map coordinates x and y via the following equations:

x = R 2 2 π ( λ λ 0 ) cos θ , y = R 2 sin θ , {\displaystyle {\begin{aligned}x&=R{\frac {2{\sqrt {2}}}{\pi }}\left(\lambda -\lambda _{0}\right)\cos \theta ,\\y&=R{\sqrt {2}}\sin \theta ,\end{aligned}}}

where θ is an auxiliary angle defined by

2 θ + sin 2 θ = π sin φ ( 1 ) {\displaystyle 2\theta +\sin 2\theta =\pi \sin \varphi \qquad (1)}

and λ is the longitude, λ0 is the central meridian, φ is the latitude, and R is the radius of the globe to be projected. The map has area 4πR, conforming to the surface area of the generating globe. The x-coordinate has a range of , and the y-coordinate has a range of .

Equation (1) may be solved with rapid convergence (but slow near the poles) using Newton–Raphson iteration:

θ 0 = φ , θ n + 1 = θ n 2 θ n + sin 2 θ n π sin φ 2 + 2 cos 2 θ n . {\displaystyle {\begin{aligned}\theta _{0}&=\varphi ,\\\theta _{n+1}&=\theta _{n}-{\frac {2\theta _{n}+\sin 2\theta _{n}-\pi \sin \varphi }{2+2\cos 2\theta _{n}}}.\end{aligned}}}


If φ = ±⁠π/2⁠, then also θ = ±⁠π/2⁠. In that case the iteration should be bypassed; otherwise, division by zero may result.

There exists a closed-form inverse transformation:

φ = arcsin 2 θ + sin 2 θ π , λ = λ 0 + π x 2 R 2 cos θ , {\displaystyle {\begin{aligned}\varphi &=\arcsin {\frac {2\theta +\sin 2\theta }{\pi }},\\\lambda &=\lambda _{0}+{\frac {\pi x}{2R{\sqrt {2}}\cos \theta }},\end{aligned}}}

where θ can be found by the relation

θ = arcsin y R 2 . {\displaystyle \theta =\arcsin {\frac {y}{R{\sqrt {2}}}}.\,}

The inverse transformations allow one to find the latitude and longitude corresponding to the map coordinates x and y.

Alterations

Allen K. Philbrick (1953) Sinu-Mollweide uninterrupted projection, with Tissot indicatrices

See also

Notes

  1. The formula in the text helps the reader confirm that the formula is correct. For numerical computation the denominator should be changed, starting with the double angle identity.
    cos 2 θ n = 2 cos 2 θ n 1 , 1 + cos 2 θ n = 2 cos 2 θ n , 2 + 2 cos 2 θ n = 4 cos 2 θ n , θ n + 1 = θ n 2 θ n + sin 2 θ n π sin φ 4 cos 2 θ n . {\displaystyle {\begin{aligned}\cos 2\theta _{n}&=2\cos ^{2}\theta _{n}-1,\\1+\cos 2\theta _{n}&=2\cos ^{2}\theta _{n},\\2+2\cos 2\theta _{n}&=4\cos ^{2}\theta _{n},\\\theta _{n+1}&=\theta _{n}-{\frac {2\theta _{n}+\sin 2\theta _{n}-\pi \sin \varphi }{4\cos ^{2}\theta _{n}}}.\end{aligned}}}
    In numerical computation, the original denominator could result in zero for θ near ±⁠π/2⁠ (catastrophic cancellation). This substitution is true for all angles and avoids the problem near θ = ±⁠π/2⁠ without making it a special case.

References

  1. Flattening the Earth: Two Thousand Years of Map Projections, John P. Snyder, 1993, pp. 112–113, ISBN 0-226-76747-7.
  2. Gannon, Megan (December 21, 2012). "New 'Baby Picture' of Universe Unveiled". Space.com. Retrieved December 21, 2012.
  3. Bennett, C.L.; Larson, L.; Weiland, J.L.; Jarosk, N.; Hinshaw, N.; Odegard, N.; Smith, K.M.; Hill, R.S.; Gold, B.; Halpern, M.; Komatsu, E.; Nolta, M.R.; Page, L.; Spergel, D.N.; Wollack, E.; Dunkley, J.; Kogut, A.; Limon, M.; Meyer, S.S.; Tucker, G.S.; Wright, E.L. (2013). "Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results". The Astrophysical Journal Supplement Series. 208 (2): 20. arXiv:1212.5225. Bibcode:2013ApJS..208...20B. doi:10.1088/0067-0049/208/2/20. S2CID 119271232.
  4. Map Projections – A Working Manual, USGS Professional Paper 1395, John P. Snyder, 1987, pp. 249–252
  5. ^ Weisstein, Eric W. "Mollweide Projection". MathWorld.

External links

Map projection
By surface
Cylindrical
Mercator-conformal
Equal-area
Pseudocylindrical
Equal-area
Conical
Pseudoconical
Azimuthal
(planar)
General perspective
Pseudoazimuthal
By metric
Conformal
Equal-area
Bonne
Bottomley
Cylindrical
Tobler hyperelliptical
Equidistant in
some aspect
Gnomonic
Loxodromic
Retroazimuthal
(Mecca or Qibla)
By construction
Compromise
Hybrid
Perspective
Planar
Polyhedral
See also
Category: