gwosc 客户端API

《code>gwosc 软件包提供了从GEO、LIGO和Virgo引力波观测站托管在 https://gwosc.org 的开放数据发布进行查询的接口。

发布状态

开发状态

安装

要安装

conda install -c conda-forge gwosc

或

pip install gwosc

搜索数据集

要搜索可用的数据集（截至2018年3月14日正确）

>>> from gwosc import datasets
>>> datasets.find_datasets()
['GW150914', 'GW151226', 'GW170104', 'GW170608', 'GW170814', 'GW170817', 'LVT151012', 'O1', 'S5', 'S6']
>>> datasets.find_datasets(detector='V1')
['GW170814', 'GW170817']
>>> datasets.find_datasets(type='run')
['O1', 'S5', 'S6']

要查询事件数据集的GPS时间（或反之亦然）

>>> datasets.event_gps('GW170817')
1187008882.43
>>> datasets.event_at_gps(1187008882)
'GW170817'

对于观测运行数据集也有类似的查询可用

>>> datasets.run_segment('O1')
(1126051217, 1137254417)
>>> datasets.run_at_gps(1135136350)  # event_gps('GW151226')
'O1'

通过事件名称定位数据URL

您可以根据事件名称搜索远程数据URL

>>> from gwosc.locate import get_event_urls
>>> get_event_urls('GW150914')
['https://gwosc.org/eventapi/json/GWTC-1-confident/GW150914/v3/H-H1_GWOSC_4KHZ_R1-1126259447-32.hdf5', 'https://gwosc.org/eventapi/json/GWTC-1-confident/GW150914/v3/H-H1_GWOSC_4KHZ_R1-1126257415-4096.hdf5', 'https://gwosc.org/eventapi/json/GWTC-1-confident/GW150914/v3/L-L1_GWOSC_4KHZ_R1-1126259447-32.hdf5', 'https://gwosc.org/eventapi/json/GWTC-1-confident/GW150914/v3/L-L1_GWOSC_4KHZ_R1-1126257415-4096.hdf5']

您可以使用关键字参数来下载数据URL

>>> get_event_urls('GW150914', detector='L1', duration=32)
['https://gwosc.org/eventapi/json/GWTC-1-confident/GW150914/v3/L-L1_GWOSC_4KHZ_R1-1126259447-32.hdf5']

通过GPS间隔定位数据URL

您可以根据GPS时间间隔搜索远程数据URL，如下所示

>>> from gwosc.locate import get_urls
>>> get_urls('L1', 968650000, 968660000)
['https://gwosc.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968646656-4096.hdf5', 'https://gwosc.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968650752-4096.hdf5', 'https://gwosc.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968654848-4096.hdf5', 'https://gwosc.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968658944-4096.hdf5']

此函数的参数如下

• detector：相关引力波干涉仪的前缀，可以是LIGO-Hanford的'H1'，或者LIGO Livingston的'L1'，
• start：感兴趣时间间隔的GPS开始时间
• end：感兴趣时间间隔的GPS结束时间

默认情况下，此方法将返回4 kHz采样率数据的HDF5文件路径，这些可以指定为关键字参数。获取完整信息，请运行

>>> help(get_urls)

查询时间线段

您还可以根据标志名称和GPS时间间隔搜索时间线段，如下所示

>>> from gwosc.timeline import get_segments
>>> get_segments('H1_DATA', 1126051217, 1126151217)
[(1126073529, 1126114861), (1126121462, 1126123267), (1126123553, 1126126832), (1126139205, 1126139266), (1126149058, 1126151217)]

输出是一个包含(start, end)二元组的list，每个二元组代表一个半开时间间隔。

有关可用标志的文档，例如O1科学运行，请参阅O1数据发布页面（数据质量）。