gw开放科学数据存档的python接口
gwopensci的Python项目详细描述
gwopensci包提供了查询 从ligo和 处女座引力波观测站。
要安装:
pip install gwopensci
搜索数据集
搜索可用数据集(截至2018年3月14日正确):
>>>fromgwopensciimportdatasets>>>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:
>>>fromgwopensci.locateimportget_event_urls>>>get_event_urls('GW150914')['https://losc.ligo.org//s/events/GW150914/H-H1_LOSC_4_V2-1126259446-32.hdf5','https://losc.ligo.org//s/events/GW150914/L-L1_LOSC_4_V2-1126259446-32.hdf5','https://losc.ligo.org//s/events/GW150914/H-H1_LOSC_4_V2-1126257414-4096.hdf5','https://losc.ligo.org//s/events/GW150914/L-L1_LOSC_4_V2-1126257414-4096.hdf5']
您可以使用关键字参数向下选择URL:
>>>get_event_urls('GW150914',detector='L1',duration=32)['https://losc.ligo.org//s/events/GW150914/L-L1_LOSC_4_V2-1126259446-32.hdf5']
按GPS间隔定位数据URL
您可以根据GPS时间间隔搜索远程数据URL,如下所示 如下:
>>>fromgwopensci.locateimportget_urls>>>get_urls('L1',968650000,968660000)['https://losc.ligo.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968646656-4096.hdf5','https://losc.ligo.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968650752-4096.hdf5','https://losc.ligo.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968654848-4096.hdf5','https://losc.ligo.org/archive/data/S6/967835648/L-L1_LOSC_4_V1-968658944-4096.hdf5']
此函数的参数如下
- detector:相关引力波的前缀 干涉仪,'H1'用于ligo hanford,或'L1'用于 利戈利文斯顿,
- start:感兴趣区间的GPS开始时间
- end:感兴趣区间的GPS结束时间
默认情况下,此方法将返回4的hdf5文件的路径 kHz采样率数据,这些可以指定为关键字参数。为了 完整信息,运行
>>>help(get_urls)
查询时间线段
您还可以基于标记名搜索时间线段,以及 GPS时间间隔如下:
>>>fromgwopensci.timelineimportget_segments>>>get_segments('H1_DATA',1126051217,1126151217)[(1126073529,1126114861),(1126121462,1126123267),(1126123553,1126126832),(1126139205,1126139266),(1126149058,1126151217)]
输出是一个list的(start, end)2元组,每个元组 表示半开放时间间隔。
有关可用标志的文档,例如O1 科学运行,请参见the O1 data release page(data quality)。