# 重力波調查員 ## 鄉民都來分析重力波觀測資料 蘇羿豪@PyCon Taiwan 2020 [創用CC 姓名標示 4.0](https://creativecommons.org/licenses/by/4.0/deed.zh_TW) <iframe frameborder="0" src="https://button.like.co/in/embed/astrobackhacker/button?referrer=https://hackmd.io/@astrobackhacker/rkZ4n-0MD"></iframe> --- 前情提要 [黑蛋被牆上銀幕播放的「重力波調查員」招募廣告所吸引](https://matters.news/@astrobackhacker/%E4%BA%94%E5%80%8B%E5%8E%BB%E8%99%95-%E8%AE%93%E4%BD%A0%E5%B0%8B%E8%A6%93%E5%A4%A9%E6%96%87%E7%9B%B8%E9%97%9Cpython%E5%A5%97%E4%BB%B6-%E5%9B%9B-ads%E8%97%8F%E7%B6%93%E9%96%A3-bafyreiagagbykv5zdschaj3s4wy4cja5gu34ohcujgamrf7sndwjcsm2yi)... --- ### :book: GW150914案發現場調查報告(p1) * 事件名稱: GW150914 * 事件描述: 35.6:sun_with_face: 和30.6:sun_with_face:相撞 * [事件照片連結](https://journals.aps.org/prl/article/10.1103/PhysRevLett.116.061102/figures/1/medium) * 事件發生時的GPS時間: _____ * [事件查詢系統連結](https://www.gw-openscience.org/eventapi/) --- ### :book: GW150914案發現場調查報告(p2) 現場鑑識人員: 用[GWpy](https://gwpy.github.io/docs/stable/)工具包可取得蒐證資料 ```python= from gwpy.timeseries import TimeSeries data = TimeSeries.fetch_open_data(interferometer, GPS_start_time, GPS_end_time) ``` [鑑識人員的蒐證報告連結](https://colab.research.google.com/github/gw-odw/odw-2020/blob/master/Day_1/Tuto%201.2%20Open%20Data%20access%20with%20GWpy.ipynb) --- ### :book: GW150914案發現場調查報告(p3) 現場鑑識人員: 用[GWpy](https://gwpy.github.io/docs/stable/)工具包可萃取資料 ```python= 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) ``` [鑑識人員的資料萃取報告連結](https://gwpy.github.io/docs/stable/examples/signal/gw150914.html) --- ### :book: GW150914案發現場調查報告(p4) 現場鑑識人員: 用[GWpy](https://gwpy.github.io/docs/stable/)工具包可轉換資料 ```python= qspecgram = data.q_transform() ``` [鑑識人員的資料轉換報告連結](https://gwpy.github.io/docs/stable/examples/timeseries/qscan.html) --- ### :book: GW150914案發現場調查報告(p5) 現場鑑識人員: 用[PyCBC](https://pycbc.org/)工具包可從資料搜索事件 ```python= from pycbc.waveform import get_td_waveform from pycbc.filter import matched_filter ``` [鑑識人員的事件搜索報告連結](https://colab.research.google.com/github/gwastro/pycbc-tutorials/blob/master/tutorial/3_WaveformMatchedFilter.ipynb) --- ### :book: GW150914案發現場調查報告(p6) GW150914是由兩個質量分別為35.6:sun_with_face:和30.6:sun_with_face:的黑洞合併釋出重力波的事件，重力波源距離地球$1.3 * 10 ^ {22}$公里...[[載入更多]](https://doi.org/10.1103/PhysRevLett.116.061102) --- ### :book: GW150914案發現場調查報告(A0) * 如何查詢已被證實的重力波事件? * [GWOSC - Event Portal](https://www.gw-openscience.org/eventapi/) * 有哪些用來分析重力波資料的軟體? * [GWOSC - Software](https://www.gw-openscience.org/software/) * 有哪些重力波分析教材? * [GWOSC - Tutorials](https://www.gw-openscience.org/tutorials/) * [GWOSC - GW Open Data Workshops ](https://www.gw-openscience.org/workshops/) --- ### :book: GW150914案發現場調查報告(A1) 如何用GWpy取得重力波時序資料並畫圖? * [GWpy Docs - Remote data access](https://gwpy.github.io/docs/stable/timeseries/remote-access.html#open-data-releases) * [GWpy Docs - Plotting time-domain data](https://gwpy.github.io/docs/stable/timeseries/plot.html) * [GWpy Docs - Plotting public LIGO data](https://gwpy.github.io/docs/stable/examples/timeseries/public.html) * [GWpy Docs - Filtering a TimeSeries to detect gravitational waves](https://gwpy.github.io/docs/stable/examples/signal/gw150914.html) * [Gravitational-Wave Open Data Workshop Tutorials](https://github.com/gw-odw/odw-2020) --- ### :book: GW150914案發現場調查報告(A2) 如何濾除資料中的雜訊？ * [GWpy Docs - Filtering a TimeSeries to detect gravitational waves](https://gwpy.github.io/docs/stable/examples/signal/gw150914.html) * [GWpy Docs - Signal processing](https://gwpy.github.io/docs/stable/signal/index.html) * [Gravitational-Wave Open Data Workshop Tutorials](https://github.com/gw-odw/odw-2020) --- ### :book: GW150914案發現場調查報告(A3) 如何將時序資料轉換成時頻資料？ * [GWpy Docs - Generate the Q-transform of a TimeSeries](https://gwpy.github.io/docs/stable/examples/timeseries/qscan.html) * [GWpy Docs - Compute the raw Q-transform of a TimeSeries](https://gwpy.github.io/docs/stable/examples/signal/qscan.html) * [Gravitational-Wave Open Data Workshop Tutorials](https://github.com/gw-odw/odw-2020) --- ### :book: GW150914案發現場調查報告(A4) 如何從資料中搜尋緻密星體合併事件? * [GWpy Docs - Calculating the SNR associated with a given astrophysical signal model](https://gwpy.github.io/docs/stable/examples/timeseries/pycbc-snr.html) * [PyCBC Tutorials - Generating Waveforms and Matched Filtering](https://colab.research.google.com/github/gwastro/pycbc-tutorials/blob/master/tutorial/3_WaveformMatchedFilter.ipynb) * [Gravitational-Wave Open Data Workshop Tutorials](https://github.com/gw-odw/odw-2020) --- 黑蛋連續劇，明年想聽哪一齣？ 1. 以天文為主題的程式教育 2. 與天文相關Python套件 3. 基於Python web framework的天文網站 --- 謝謝，請幫我:clap: 也歡迎來[中央大學天文所Python使用者討論會的FB社團](https://www.facebook.com/groups/1022708484514663/)交流 <iframe frameborder="1" src="https://button.like.co/in/embed/astrobackhacker/button?referrer=https://hackmd.io/@astrobackhacker/rkZ4n-0MD"></iframe>