<html><body name="Mail Message Editor"><div><br></div><div>I wote:</div><div><br></div><div><div>> yielding two solutions, depending on the sign chosen for the radical in B. Both are valid solutions.</div><div><br></div><div>Which is not correct. Only one of the two solutions will produce a legal value of φ₁, the other being outside the range of ±90°. Hence, choose the one within range. There is probably some analytic way of deciding which it will be, but I haven't tried to work it out.</div><div><br></div></div><div>Regards,</div><div>-- daan Strebe</div><div><br></div><br>On Oct 20, 2008, at 4:59:34 PM, strebe <strebe@aol.com> wrote:<br><blockquote style="padding-left: 5px; margin-left: 5px; border-left-width: 2px; border-left-style: solid; border-left-color: blue; color: blue; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><div><br class="Apple-interchange-newline">Mikæl:</div><div><br></div><div>The usual spherical formulæ for computing a point at a known distance and azimuth are:<br></div><div><br></div><div><span class="Apple-style-span" style="font-family: Monaco; ">φ₂ = arcsin (sin φ₁ cos c + cos φ₁ sin c cos Az)<br></span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">λ₂ = λ₁ + atan2 [sin c sin Az, (cos φ₁ cos c - sin φ₁ sin c cos Az)]</span></div><div><br></div><div>where the start point is subscripted with 1, the end point is subscripted with 2, and the arc distance is c. We know φ₂, λ₂, Az, and c. Therefore these reduce to</div><div><br></div><div><div><span class="Apple-style-span" style="font-family: Monaco; ">(1) n₄ = n₁ sin φ₁ + n₂ cos φ₁</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">(2) λ₂ = λ₁ + atan2 [n₃, (n₁ cos φ₁ - n₂ sin φ₁)]</span></div><div><br></div><div>where we have made substitutions for symbolic clarity:</div><div><br></div><div><span class="Apple-style-span" style="font-family: Monaco; ">n₁ = cos c</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">n₂ = sin c cos Az</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">n₃ = sin c sin Az</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">n₄ = sin φ₂</span><br></div><div><br></div><div>Since equation 1 <span class="Apple-style-span" style="font-family: Monaco; "><span class="Apple-style-span" style="font-family: Helvetica; "><span class="Apple-style-span" style="background-color: transparent; ">ha</span></span><span class="Apple-style-span" style="font-family: Helvetica; ">s only φ₁ unknown, we solve with trigonometric and algebraic manipulations:</span></span></div></div><div><br></div><div><span class="Apple-style-span" style="font-family: Monaco; ">tan φ₁ = A / B</span></div><div><br></div><div>with φ₁ resolvable using the usual atan2 function in order to preserve quadrants. Here</div><div><br></div><div><span class="Apple-style-span" style="font-family: Monaco; "> n₄ - n₂ B</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">A = _________</span></div><div><span class="Apple-style-span" style="font-family: Monaco; "> n₁</span></div><div><span class="Apple-style-span" style="font-family: Monaco; "><br></span></div><div><span class="Apple-style-span" style="font-family: Monaco; "> n₂ n₄ ± √(n₂² n₁² + n₁⁴ - n₁² n₄²)</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">B = ----------------------------------</span></div><div><span class="Apple-style-span" style="font-family: Monaco; "> n₁² + n₂²</span></div><div><br></div><div><br></div><div>yielding two solutions, depending on the sign chosen for the radical in B. Both are valid solutions.</div><div><br></div><div>Thus armed with φ₁, we can easily compute<span class="Apple-style-span" style="background-color: transparent; "> </span><span class="Apple-style-span" style="font-family: Monaco; "><span class="Apple-style-span" style="font-family: Helvetica; "><span class="Apple-style-span" style="background-color: transparent; ">λ₁ </span></span><span class="Apple-style-span" style="font-family: Helvetica; ">as</span></span></div><div><br></div><div><div><span class="Apple-style-span" style="font-family: Monaco; ">λ₁ = λ₂ - atan2 [n₃, (n₁ cos φ₁ - n₂ sin φ₁)]</span></div><div><span class="Apple-style-span" style="font-family: Monaco; "><br></span></div></div><div>The ellipsoidal case is much harder. Your iterative process will work, of course, though it will be terribly inefficient and likely comes with problematic regions. In principle it seems like the problem you are trying to solve is intermediate between the forward and the inverse problems. The reason it is not as simple as the forward problem is that, in the forward case, you know both the starting point and the azimuth. This wholly characterizes which geodesic you are dealing with, since a geodesic is completely characterized by its constant and any point along it:</div><div><span class="Apple-style-span" style="font-family: Monaco; "><br></span></div><div><span class="Apple-style-span" style="font-family: Monaco; "> cos φ sin Az</span></div><div><span class="Apple-style-span" style="font-family: Monaco; ">K = ________________</span></div><div><span class="Apple-style-span" style="font-family: Monaco; "> √(1 - e² sin² φ)</span></div><div><br></div><div>In your case, you cannot compute the geodesic's constant because you don't know both the latitude and azimuth from any known point; instead you know the azimuth from an unknown point. On the other hand, you already know the distance, the calculation of which is what makes the inverse problem so difficult.</div><div><br></div><div>Warning: Do NOT use Pearson's analysis of geodesics (Frederick Pearson II 1984: "Map Projection Methods", pp. 81–82, and 1990: "Map Projections: Theory and Applications", pp. 80–81). He makes an error early on in the derivation, making the final, easily computed integral nothing but a fiction. Sadly, Donald Fenna has repeated this error in his 2007 volume, "A Compendium of Map Projections, with Derivations", pp. 394–395. </div><div><br></div><div>Regards,</div><div>-- daan Strebe</div><div><br></div></span><br class="Apple-interchange-newline"></blockquote><br><div><br></div><div class="aol_ad_footer" id="u037837880DB048AF887CB1FE63291D4D"></div></body></html>