def test_distaz(self):
"""
Tests distance and azimuth calculation between two points on a sphere.
"""
client = Client()
# normal request
result = client.distaz(stalat=1.1, stalon=1.2, evtlat=3.2, evtlon=1.4)
self.assertAlmostEquals(result["distance"], 2.09554)
self.assertAlmostEquals(result["backazimuth"], 5.46946)
self.assertAlmostEquals(result["azimuth"], 185.47692)
# w/o kwargs
result = client.distaz(1.1, 1.2, 3.2, 1.4)
self.assertAlmostEquals(result["distance"], 2.09554)
self.assertAlmostEquals(result["backazimuth"], 5.46946)
self.assertAlmostEquals(result["azimuth"], 185.47692)
# missing parameters
self.assertRaises(Exception, client.distaz, stalat=1.1)
self.assertRaises(Exception, client.distaz, 1.1)
self.assertRaises(Exception, client.distaz, stalat=1.1, stalon=1.2)
self.assertRaises(Exception, client.distaz, 1.1, 1.2)
def test_distaz(self):
"""
Tests distance and azimuth calculation between two points on a sphere.
"""
client = Client()
# normal request
result = client.distaz(stalat=1.1, stalon=1.2, evtlat=3.2, evtlon=1.4)
self.assertAlmostEquals(result['distance'], 2.09554)
self.assertAlmostEquals(result['backazimuth'], 5.46946)
self.assertAlmostEquals(result['azimuth'], 185.47692)
# w/o kwargs
result = client.distaz(1.1, 1.2, 3.2, 1.4)
self.assertAlmostEquals(result['distance'], 2.09554)
self.assertAlmostEquals(result['backazimuth'], 5.46946)
self.assertAlmostEquals(result['azimuth'], 185.47692)
# missing parameters
self.assertRaises(Exception, client.distaz, stalat=1.1)
self.assertRaises(Exception, client.distaz, 1.1)
self.assertRaises(Exception, client.distaz, stalat=1.1, stalon=1.2)
self.assertRaises(Exception, client.distaz, 1.1, 1.2)
def az2(sfile):
azDic = {}
client = Client() #Creando objeto cliente
evtlat,evtlon = evcor(sfile) #Obteniendo las coordenadas del evento del sfile en decimal
infile = open(sfile,'rU')
for line in infile:
if line[9:11] == 'ES':
#print 'Dentro del if'
station_name = line[1:6].split(' ')[0]
Baz = client.distaz(stalat=stacor('corstations.txt')[station_name][0], stalon=stacor('corstations.txt')[station_name][1], evtlat=evtlat, evtlon=evtlon)
azDic[station_name] = Baz['backazimuth']
return azDic
#!/usr/bin/python from obspy.iris import Client #from obspy.fdsn import Client client = Client() result = client.distaz(stalat=1.1, stalon=1.2, evtlat=3.2, evtlon=1.4) print(result['distance']) print(result['backazimuth']) print(result['azimuth'])
