def testGeodCalc(self): self.subtitle(karney, 'GeodCalc') # <https://GeographicLib.sourceforge.io/scripts/geod-calc.html#area> p = [ LatLon_(*t) for t in ((-63.1, -58), (-72.9, -74), (-71.9, -102), (-74.9, -102), (-74.3, -131), (-77.5, -163), (-77.4, 163), (-71.7, 172), (-65.9, 140), (-65.7, 113), (-66.6, 88), (-66.9, 59), (-69.8, 25), (-70.0, -4), (-71.0, -14), (-77.3, -33), (-77.9, -46), (-74.7, -61)) ] self.test('area', ellipsoidalKarney.areaOf(p), 13662703680020.125, fmt='%.0f') self.test('perimeter', ellipsoidalKarney.perimeterOf(p, closed=True), 16831067.89279071, fmt='%.6f')
def areaOf(points, datum=Datums.WGS84, wrap=True): # PYCHOK no cover '''DEPRECATED, use function C{ellipsoidalKarney.areaOf}. ''' from pygeodesy.ellipsoidalKarney import areaOf return areaOf(points, datum=datum, wrap=wrap)