# 重力波調查員
## 鄉民都來分析重力波觀測資料

蘇羿豪@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/)工具包可取得蒐證資料

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/)工具包可萃取資料

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/)工具包可轉換資料

qspecgram = data.q_transform()

[鑑識人員的資料轉換報告連結](https://gwpy.github.io/docs/stable/examples/timeseries/qscan.html)

---

### :book: GW150914案發現場調查報告(p5)

現場鑑識人員: 用[PyCBC](https://pycbc.org/)工具包可從資料搜索事件

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>