重力波調查員

鄉民都來分析重力波觀測資料

蘇羿豪@PyCon Taiwan 2020
創用CC 姓名標示 4.0

前情提要

黑蛋被牆上銀幕播放的「重力波調查員」招募廣告所吸引…

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(p1)

  • 事件名稱: GW150914
  • 事件描述: 35.6
    Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →
    和30.6
    Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →
    相撞
  • 事件照片連結
  • 事件發生時的GPS時間: _____
  • 事件查詢系統連結

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(p2)

現場鑑識人員: 用GWpy工具包可取得蒐證資料





from gwpy.timeseries import TimeSeries

data = TimeSeries.fetch_open_data(interferometer, GPS_start_time, GPS_end_time)

鑑識人員的蒐證報告連結

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(p3)

現場鑑識人員: 用GWpy工具包可萃取資料







from gwpy.signal import filter_design

bp = filter_design.bandpass(low_f, high_f, sample_rate)
notche = filter_design.notch(frequency, sample_rate)
zpk = filter_design.concatenate_zpks(bp, notche)
filtered_data = data.filter(zpk)

鑑識人員的資料萃取報告連結

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(p4)

現場鑑識人員: 用GWpy工具包可轉換資料

qspecgram = data.q_transform()

鑑識人員的資料轉換報告連結

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(p5)

現場鑑識人員: 用PyCBC工具包可從資料搜索事件



from pycbc.waveform import get_td_waveform
from pycbc.filter import matched_filter

鑑識人員的事件搜索報告連結

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(p6)

GW150914是由兩個質量分別為35.6

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
和30.6
Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
的黑洞合併釋出重力波的事件，重力波源距離地球\(1.3 * 10 ^ {22}\)公里…[載入更多]

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(A0)

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(A1)

如何用GWpy取得重力波時序資料並畫圖?

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(A2)

如何濾除資料中的雜訊？

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(A3)

如何將時序資料轉換成時頻資料？

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
GW150914案發現場調查報告(A4)

如何從資料中搜尋緻密星體合併事件?

黑蛋連續劇，明年想聽哪一齣？

  1. 以天文為主題的程式教育
  2. 與天文相關Python套件
  3. 基於Python web framework的天文網站

謝謝，請幫我

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
也歡迎來中央大學天文所Python使用者討論會的FB社團交流

Select a repo