ReadLIGO

此模块提供读取LIGO数据文件的工具。数据及相关文档可以从GWOSC网站下载。

安装

要安装，请下载仓库并使用pip进行安装

git clone https://git.ligo.org/gwosc/readligo.git; cd readligo
pip install .

文档

以下展示了某些可能的用例。

示例 0

从单个文件加载所有数据

strain, time, dq = rl.loaddata("ligo_data/H-H1_LOSC_4_V1-842653696-4096.hdf5", "H1")

某些GWF文件需要参数来命名应变、DQ和硬件注入通道

strain, time, dq = rl.loaddata(
    "H-H1_LOSC_16_V1-1127415808-4096.gwf",
    "H1",
    strain_chan="H1:GWOSC-16KHZ_R1_STRAIN", 
    dq_chan="H1:GWOSC-16KHZ_R1_DQMASK",
    inj_chan="H1:GWOSC-16KHZ_R1_INJMASK"
)

示例 1

此默认配置假定所需的LIGO数据文件在当前工作目录或子目录中可用。LIGO数据在输入GPS时间之间被加载到STRAIN中。META是包含gps开始、gps停止和样本时间的字典。DQ是包含数据质量标志的字典。

segList = getsegs(842657792, 842658792, 'H1')
for (start, stop) in segList:
  strain, meta, dq = getstrain(start, stop, 'H1')
  # -- Analysis code here
  ...

示例 2

在示例2中，H1_segs.txt是从GWOSC网站使用时间线应用下载的段列表。它可以像示例1中的segList一样使用。

segList = SegmentList("H1_segs.txt")

示例 3

在此示例中，第一个命令搜索指定目录及其子目录中的LIGO数据文件。然后使用这些数据文件列表构建段列表并加载所需数据。

filelist = FileList(directory="/home/ligodata")
segList = getsegs(842657792, 842658792, "H1", filelist=filelist)
for start, stop in segList:
  strain, meta, dq = getstrain(start, stop, "H1", filelist=filelist)
  # -- Analysis code here

段列表

您可以使用时间线查询表单从GWOSC网站下载段列表，或者直接从数据文件构建。

使用SegmentList读取从GWOSC网站的时间线应用程序下载的段列表。

from readligo import SegmentList
seglist = SegmentList("H1_segs.txt")

或者，使用getsegs直接从LIGO数据文件构建段列表。

from readligo import getsegs
seglist = getsegs(842657792, 842658792, 'H1', flag='DATA', filelist=None)