def test_should_compute_equator( self ): """somewhere near -5d 50m for 4/18/2012""" p1= LatLon( lat=Angle2.fromstring("0.0N"), lon=Angle2.fromstring( "0.0E" ) ) #print( "p1 = {0}".format( p1.dms ) ) d= declination( p1, date=datetime.date(2012,4,18) ) #print( "declination= {0}".format(d.dm) ) self.assertAlmostEqual( -5, d.dm[0] ) self.assertTrue( 45 <= d.dm[1] <= 55 )
def test_should_compute( self ): p1= LatLon( lat=Angle2.fromstring("51 07 32N"), lon=Angle2.fromstring("001 20 17E") ) #print( "p1 = {0}".format( p1.dms ) ) bearing= Angle( "116 38 10" ) #print( "theta = {0}".format( bearing.dms ) ) p2= destination( p1, 40.23, bearing, R=KM ) self.assertEqual( ("50 57 48.1N", "001 51 08.8E"), p2.dms ) self.assertEqual( ("50.963N", "001.852E"), p2.d )
def test_should_compute( self ): p1= LatLon( lat=Angle2.fromstring("50 21 50N"), lon=Angle2.fromstring("004 09 25W") ) p2= LatLon( lat=Angle2.fromstring("42 21 04N"), lon=Angle2.fromstring("071 02 27W") ) #print( "p1 = {0}".format( p1.dms ) ) #print( "p2 = {0}".format( p2.dms ) ) d, brg = range_bearing( p1, p2, R=KM ) self.assertEqual( 5196, int(d) ) self.assertEqual( (260, 7, 38), tuple(map(int,map(round,brg.dms))) ) d, brg = range_bearing( p1, p2, R=NM ) self.assertEqual( 2805, int(d) ) self.assertEqual( (260, 7, 38), tuple(map(int,map(round,brg.dms))) )