关于如何利用obspy求出穿透点的经纬度信息

1. 介绍

在Taup中可以利用taup_pierce求出不连续界面穿透点的经纬度信息，在obspy中可以利用TauPyModel中的get_pierce_points_geo获取穿透点的经纬度信息，在此记录一下。注意其返回值是Arrival类型数据，可以从Arrival中得到获取到穿透点的震相的名称、走时、深度，穿透点的经纬度等信息。

2. 实现

利用get_pierce_points_geo获取穿透点信息后在通过Arrival索引并利用pandas获取经纬度信息。

其他信息可以参照obspy官网手册。

    import pandas as pd
    from obspy.taup import TauPyModel
    model = TauPyModel(model='prem')

    evdp=20;evlo=0;evla=15      # Event information
    stlo=40;stla=40             # Receiver information
    # Get information about pierce point 
    info = model.get_pierce_points_geo(source_depth_in_km=evdp, source_latitude_in_deg=evla, source_longitude_in_deg=evlo, receiver_latitude_in_deg=stla,\
                     receiver_longitude_in_deg=stlo, phase_list=['PcP'])
    # Data type
    print('Data type is :\n %s\n' % (info[0].pierce.dtype))
    # Data of pierce point 
    data_info = pd.DataFrame(info[0].pierce)
    print('Data information is : \n %s' % data_info)

    # Get information of (lon,lat) on Core Mantle Boundary (PcP Phase)
    depth_max_value = data_info.depth.max();depth_max_idnex = data_info.depth.idxmax()
    lon = data_info.lon[depth_max_idnex]; lat = data_info.lat[depth_max_idnex]
    print('\nPierce point information of PcP on CMB is : (%f,%f)' % (lon,lat))

3. 输出

Data type is :
 [('p', '<f8'), ('time', '<f8'), ('dist', '<f8'), ('depth', '<f8'), ('lat', '<f8'), ('lon', '<f8')]

Data information is : 
              p        time      dist   depth        lat        lon
0   191.402313    0.000000  0.000000    20.0  15.000000   0.000000
1   191.402313    0.661101  0.000145    24.4  15.005736   0.006235
2   191.402313   25.732053  0.008111   220.0  15.320268   0.348925
3   191.402313   47.180361  0.016635   400.0  15.656238   0.716749
4   191.402313   76.335193  0.032359   670.0  16.274339   1.398427
5   191.402313  291.410855  0.373603  2891.0  28.953970  17.579481
6   191.402313  506.486516  0.714848   670.0  39.208814  37.775856
7   191.402313  535.641348  0.730571   400.0  39.591747  38.831133
8   191.402313  557.089656  0.739095   220.0  39.795326  39.408068
9   191.402313  582.160608  0.747061    24.4  39.982988  39.950332
10  191.402313  582.821709  0.747207    20.0  39.986385  39.960240
11  191.402313  583.572935  0.747371    15.0  39.990237  39.971484
12  191.402313  586.199420  0.747789     0.0  40.000000  40.000000

Pierce point information of PcP on CMB is : (17.579481,28.953970)