[Proj] Standard projection for Mediterranean basin

[strebe]

Jan & colleagues:

Regarding the question of optimal standard parallels for Albers, I found my notes from 8 years ago. To summarize, there is no good rule, so perhaps the exact calculation I describe below has more importance than just academic.

Regards,
— daan Strebe

________________________


The problem: given two latitudinal extents phi[N] and phi[S], find the optimal standard parallels for the Albers projection. Optimal is defined to mean that the maximum angular deformation between the two standard parallels equals the maximum angular deformation at the two latitudinal extents. In this way one is assured that the maximum deviation from true remains minimal throughout the region of interest.

In a breathtakingly brief treatment, Melluish touches on the issue by describing how to select standard parallels such that the "error" on the extreme parallels is equal in magnitude and opposite in sign to the error of the central parallel. The distortion measurement is equivalent to mine, but the central parallel (the average of the north and south extreme parallels) generally does not carry the greatest distortion, so his method does not yield an optimal map. Hinks discusses the problem precisely, even noting that the central parallel is not the parallel of greatest distortion... but his treatment applies to the equidistant conic, not the equal-area, contrary to Snyder's comment in USGS P.P 1395.

We define a family of Albers projections to be all Albers projections characterized by a common 'central' latitude at which angular deformation reaches a peak between the two standard parallels, or sinks to zero in the case of one standard parallel. I have placed 'central' between quote marks because it is not any simple average of the standard parallels. We can demonstrate that the ratio of the angular deformation of any two points on the map remains constant for every member of the family. This is the beautiful and useful characteristic of the family. The property arises out of a happy mathematical coincidence.

First I declare the member of the family with only one standard parallel to be the canonical form. Because of the family's property of invariant scale ratios, we can first reduce the problem of choosing an optimal Albers projection to the single standard parallel case, which is much more tractable than trying to deal with two free parameters. This gives us a projection with zero distortion along the canonical parallel and identical distortion at the northern and southern extremities. Once we have found the canonical projection, we can choose the member of its family that produces the same magnitude of distortion at the canonical parallel as at the extreme parallels. We are guaranteed that the projection that fulfills this condition is one from the family of the canonical projection.

Let phi[S] denote the southern parallel of interest; phi[N] the northern parallel of interest; phi[0] the canonical parallel; phi[1] a standard parallel; and phi[2] the other standard parallel.

phi[S] and phi[N] are given. From them we must determine phi[0], phi[1], and phi[2], such that the magnitude of the angular deformation at phi[0] is the same as the magnitude of the angular deformation at phi[N] and phi[S]. Furthermore, there must be no greater angular deformation registered anywhere between the phi[S] and phi[N]. Given those conditions:

1) sin (phi[0]) = A +/- sqrt (A^2 - 1)

where

2) A = (sin (phi[N]) * cos^2 (phi[S]) - sin (phi[S]) * cos^2 (phi[N])) / (cos^2 (phi[S] - cos^2 (phi[N]))

That establishes the canonical parallel of the canonical form, and therefore the family, wherein lies the solution.

If we arbitrarily select some phi[1], one may ask how to arrive at the correct phi[2] within the family, given phi[0].

3) sin (phi[2]) = [sin (phi[1]) - 2 sin (phi[N]) + sin (phi[1]) sin^2 (phi[N])] / [2 sin (phi[N]) * sin (phi[1]) - sin^2 (phi[N]) - 1]

However, we do not yet know the proper phi[1]. With much manipulation, we discover the iterative equation:

4) cos^2 (phi[1]) * sqrt (1 + sin^2 (phi[0]) - 2 sin (phi[0]) * sin (phi[S])) - cos (phi[S]) (1 + sin^2 (phi[0] - 2 sin (phi[0]) sin (phi[1]) = 0

We solve for phi[1] with (4) and then use (3) to compute phi[2]. Now have the optimal Albers for any circumstance.

The literature fixates quite a bit on rules of thumb for choosing the standard parallels, with some texts recommending inward 1/5th the meridional separation of the extreme parallels. Other texts recommend 1/7th, and yet others recommend 1/6th. My own calculations based on the above criteria show that the correct values range from 5% to 25% for reasonable conics, and the symmetry between north and south is not necessarily very good. Hence, no rule of thumb is very good. The higher the latitude of the northern parallel, the greater the asymmetry. The asymmetry remains remarkable even if both latitudes are high. For instance, extreme latitudes at 50° and 70° result in standard parallels placed optimally at roughly 10% below phi[N] and 18% above phi[S].

Deetz's and Adams's recommendation of 29°30' and 45°30' for the US are pretty good, of course. Given maximum extents of 24°34' and 49°25', the optimal standard parallels are 28°42' and 46°19'. If we just can't bring ourselves to treat Lake of the Woods equitably then the upper standard parallel drops to 45°56'. The lower parallel only drops by 5'. The former map yield a maximum east-west scale distortion of 1.19%; the latter yields 1.15% (excluding Lake of the Woods, of course). Deetz's and Adams's recommendation yields 0.98% in the heartland; 1.25% at the northern border; 1.45% at the northern tip of the Lake of the Woods boundary; and 1.37% at Key West.

________________________


On Feb 15, 2010, at 1:05:19 PM, "Jan Hartmann" <j.l.h.hartmann at uva.nl> wrote:
From:   "Jan Hartmann" <j.l.h.hartmann at uva.nl>
Subject:    Re: [Proj] Standard projection for Mediterranean basin
Date:   February 15, 2010 1:05:19 PM PST
To: "PROJ.4 and general Projections Discussions" <proj at lists.maptools.org>
Cc: strebe <strebe at aol.com>
What's the 1/6 rule, Daan, and how much does it deviate from the optimal formula? 

BTW Nice Java applet to study projection deformations here:

http://www.uff.br/mapprojections/Albers_en.html

Pity they didn't make it possible to alter projection parameters like standard parallels.

Jan 

On 15-2-2010 21:04, strebe wrote:

I do not understand Dr. Mugnier's response. Using the authalic latitude does not yield an equal-area projection if you project using a conformal projection. I cannot think of any use for mixing authalic latitude with conformal projection.

The purpose of the authalic latitude is to provide an equal-area projection of the ellipsoid onto the sphere. From there you may project to the plane using any spherical equal-area projection, yielding an equal-area projection of the ellipsoid onto the plane. While this technique will not preserve other properties of the spherical projection (such as standard parallels or constant scale along some meridian or somesuch), at least it drastically increases the available projections for ellipsoidal use.

I have, by the way, formulæ somewhere that I derived for computing the optimal standard parallels for Albers given the north/south extents of the map. I will post them if I happen to find them. The 1/6th rule would be fine for any real use, of course. The optimal is more of academic interest.

Regards,
— daan Strebe


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20100215/af79532e/attachment-0001.htm