Theory, strict formula derivation and algorithm development for the computation of a geodesic polygon area
Edward Nowak,Joanna Nowak Da Costa
DOI: https://doi.org/10.1007/s00190-022-01606-z
IF: 4.809
2022-03-24
Journal of Geodesy
Abstract:Abstract The paper presents a coherent, original theory and an effective algorithm to calculate the surface area of a geodesic polygon, i.e. the area on an ellipsoid of revolution limited by geodesics, connecting pairs of subsequent vertices of known geographical coordinates. Using the cylindrical coordinate system and introducing an equatorial geodesic triangle, the novel strict recursive formulas were derived for the longitude on the ellipsoid and the distance along the geodesic and the area of a geodesic polygon. The area is determined by constant coefficients, not the values of finite series. The orthodromic properties of geodesics, which led to the detection of bifurcation lines composed of the double solutions of the inverse geodetic problem were analysed. The novel algorithm presented, integrating the problem with geodesic polygon area calculation, is based on the newly proposed parameterisation with two alternative working unknowns and a new formula to obtain the initial value of the geodesic C parameter. It allows solutions for all cases, including near-antipodal or near-vertices, and offers two solutions if bifurcation occurs. The algorithm presented was tested and verified using numerical examples from known publications. The numerical tests proved the accuracy of 1m 2 of the computed geodesic area using standard IEEE 754 Double Floating Point PC arithmetic, for geodesic lengths up to ten thousand km.
geochemistry & geophysics,remote sensing