Acerca de cómo usar obspy para encontrar la información de latitud y longitud del punto de penetración

1. Introducción

En Taup, puede usar taup_pierce para obtener la información de longitud y latitud del punto de penetración de la interfaz discontinua. En obspy, puede usar get_pierce_points_geo en TauPyModel para obtener la información de longitud y latitud del punto de penetración. Regístrelo aquí. Tenga en cuenta que el valor de retorno son datos de tipo de llegada. Puede obtener el nombre, el tiempo de viaje, la profundidad, la longitud y la latitud del punto de penetración y otra información de la llegada.

2. Implementación

Utilice get_pierce_points_geo para obtener la información del punto de penetración, luego utilice el índice de llegada y utilice pandas para obtener la información de latitud y longitud.

Para obtener más información, consulte el manual del sitio web oficial de 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. Salida

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)