depthdat = config.get('FileIO', 'depthdat')
stackfile = config.get('FileIO', 'stackfile')
Velmod = config.get('FileIO', 'Velmod')
stalist = config.get('FileIO', 'stalist')
domperiod = float(config.get('para', 'domperiod'))
Profile_width = float(config.get('para', 'Profile_width'))
bin_radius = float(config.get('para', 'bin_radius'))
Stack_range = np.arange(350, 751)
azi = distaz.distaz(lat1, lon1, lat2, lon2).baz
dis = distaz.distaz(lat1, lon1, lat2, lon2).delta
Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width)
Profile_lat = []
Profile_lon = []
for Profile_loca in Profile_range:
(lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi,
distaz.km2deg(Profile_loca))
Profile_lat = np.append(Profile_lat, [lat_loca], axis=1)
Profile_lon = np.append(Profile_lon, [lon_loca], axis=1)
# ----- Read depth .mat file -----#
depthmat = sio.loadmat(depthdat)
RFdepth = depthmat['YN_RFdepth']
# find stations beside the profile
stalat_all = RFdepth[0, 0::]['stalat']
stalon_all = RFdepth[0, 0::]['stalon']
stalat_full = []
stalon_full = []
staidx = []
for i in range(stalon_all.shape[0]):
stalat_full = np.append(stalat_full, stalat_all[i][0])
lat2 = float(line.split('/')[2])
lon2 = float(line.split('/')[3])
depthdat = config.get('FileIO', 'depthdat')
stackfile = config.get('FileIO', 'stackfile')
Velmod = config.get('FileIO', 'Velmod')
stalist = config.get('FileIO', 'stalist')
domperiod = float(config.get('para', 'domperiod'))
Profile_width = float(config.get('para', 'Profile_width'))
Stack_range = np.arange(1, 101)
azi = distaz.distaz(lat1, lon1, lat2, lon2).baz
dis = distaz.distaz(lat1, lon1, lat2, lon2).delta
Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width)
Profile_lat = []
Profile_lon = []
for Profile_loca in Profile_range:
(lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi, distaz.km2deg(Profile_loca))
Profile_lat = np.append(Profile_lat, [lat_loca], axis=1)
Profile_lon = np.append(Profile_lon, [lon_loca], axis=1)
# read depth data
depthmat = sio.loadmat(depthdat)
RFdepth = depthmat['YN_RFdepth']
# find stations beside the profile
stalat_all = RFdepth[0, 0::]['stalat']
stalon_all = RFdepth[0, 0::]['stalon']
stalat_full = []
stalon_full = []
staidx = []
for i in range(stalon_all.shape[0]):
stalat_full = np.append(stalat_full, stalat_all[i][0])
lon2 = float(line.split('/')[3])
depthdat = config.get('FileIO', 'depthdat')
stackfile = config.get('FileIO', 'stackfile')
Velmod = config.get('FileIO', 'Velmod')
stalist = config.get('FileIO', 'stalist')
domperiod = float(config.get('para', 'domperiod'))
Profile_width = float(config.get('para', 'Profile_width'))
bin_radius = float(config.get('para', 'bin_radius'))
Stack_range = np.arange(350, 751)
azi = distaz.distaz(lat1, lon1, lat2, lon2).baz
dis = distaz.distaz(lat1, lon1, lat2, lon2).delta
Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width)
Profile_lat = []
Profile_lon = []
for Profile_loca in Profile_range:
(lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi, distaz.km2deg(Profile_loca))
Profile_lat = np.append(Profile_lat, [lat_loca], axis=1)
Profile_lon = np.append(Profile_lon, [lon_loca], axis=1)
# ----- Read depth .mat file -----#
depthmat = sio.loadmat(depthdat)
RFdepth = depthmat['YN_RFdepth']
# find stations beside the profile
stalat_all = RFdepth[0, 0::]['stalat']
stalon_all = RFdepth[0, 0::]['stalon']
stalat_full = []
stalon_full = []
staidx = []
for i in range(stalon_all.shape[0]):
stalat_full = np.append(stalat_full, stalat_all[i][0])
lon2 = float(line.split('/')[3])
depthdat = config.get('FileIO', 'depthdat')
stackfile = config.get('FileIO', 'stackfile')
Velmod = config.get('FileIO', 'Velmod')
stalist = config.get('FileIO', 'stalist')
domperiod = float(config.get('para', 'domperiod'))
Profile_width = float(config.get('para', 'Profile_width'))
Stack_range = np.arange(1, 101)
azi = distaz.distaz(lat1, lon1, lat2, lon2).baz
dis = distaz.distaz(lat1, lon1, lat2, lon2).delta
Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width)
Profile_lat = []
Profile_lon = []
for Profile_loca in Profile_range:
(lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi,
distaz.km2deg(Profile_loca))
Profile_lat = np.append(Profile_lat, [lat_loca], axis=1)
Profile_lon = np.append(Profile_lon, [lon_loca], axis=1)
# read depth data
depthmat = sio.loadmat(depthdat)
RFdepth = depthmat['YN_RFdepth']
# find stations beside the profile
stalat_all = RFdepth[0, 0::]['stalat']
stalon_all = RFdepth[0, 0::]['stalon']
stalat_full = []
stalon_full = []
staidx = []
for i in range(stalon_all.shape[0]):
stalat_full = np.append(stalat_full, stalat_all[i][0])
