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))) )