本文整理汇总了Python中gwpy.segments.DataQualityFlag.write方法的典型用法代码示例。如果您正苦于以下问题:Python DataQualityFlag.write方法的具体用法?Python DataQualityFlag.write怎么用?Python DataQualityFlag.write使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类gwpy.segments.DataQualityFlag
的用法示例。
在下文中一共展示了DataQualityFlag.write方法的3个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: test_write_hdf5
# 需要导入模块: from gwpy.segments import DataQualityFlag [as 别名]
# 或者: from gwpy.segments.DataQualityFlag import write [as 别名]
def test_write_hdf5(self, delete=True):
flag = DataQualityFlag(FLAG1, active=ACTIVE, known=KNOWN)
hdfout = self.tmpfile % 'hdf'
try:
flag.write(hdfout)
except ImportError as e:
self.skipTest(str(e))
else:
if delete:
os.remove(hdfout)
return hdfout
示例2: argwhere
# 需要导入模块: from gwpy.segments import DataQualityFlag [as 别名]
# 或者: from gwpy.segments.DataQualityFlag import write [as 别名]
idxhighscat = argwhere(scatf2>=thresh)
highscatf2 = scatf2[idxhighscat]
highscattimes = times[idxhighscat]
highscattimesgps = highscattimes+start_time
# save text file with values above threshold [GPS f2 index]
outdata = hstack((highscattimesgps,highscatf2,idxhighscat))
savetxt('%s-ALL-TIMES-SCATTER-GT%dHZ-%d-%d.txt' % (ifo,thresh,start_time,dur),outdata,fmt='%f %f %i')
# save segments XML file with segments (based off code from Duncan Macleod)
from math import (floor, ceil)
from gwpy.segments import (Segment, DataQualityFlag)
flag = '%s:DCH-SCATTERED_LIGHT_GT%dHZ:1' % (ifo,thresh)
flag = DataQualityFlag(flag)
segs = []
append = segs.append
for gps in highscattimesgps:
if len(segs) and gps in segs[-1]:
continue
seg = Segment(floor(gps), ceil(gps))
append(seg)
flag.active = segs
flag.known = [Segment(start_time, end_time)]
flag.coalesce()
flag.write('%s-%s_%d-%d-%d.xml.gz' % (flag.ifo, flag.tag.replace('-', '_'), flag.version,start_time, dur))
#EOF
示例3: range
# 需要导入模块: from gwpy.segments import DataQualityFlag [as 别名]
# 或者: from gwpy.segments.DataQualityFlag import write [as 别名]
else:
print filename + " does not exist. Looking for the segment file in next time increment."
break
f.close()
#construct flag and filename
flag_name = 'H1:UPVh-RND:1' #NEEDS TO BE CHANGED
name = 'segments_UPVh_RND.xml' #NEEDS TO BE CHANGED
try: knownsegments = numpy.loadtxt('total_UPVh_segs.txt')
except:
print 'No total_UPVh_segs.txt file in current working directory. It should have been produced from last loop.'
#knownsegments = numpy.loadtxt(total_UPVh_segs.txt'
known_start = [knownsegments[i,0] for i in range(len(knownsegments))]
known_end = [knownsegments[i,1] for i in range(len(knownsegments))]
# read the data
data = numpy.loadtxt('total_UPVh_trigs.txt', dtype=float)
# get an array for the start_time and end_time of each segment
start_time = [data[i,0] for i in range(len(data))]
end_time = [data[i,1] for i in range(len(data))]
# create a data quality flag object
#zip will truncate the start and end time. is this OK?
flag = DataQualityFlag(flag_name, active=zip(start_time, end_time), known=zip(known_start, known_end))
# write flag
flag.write(name)