Skip to content

Commit eb2350b

Browse files
rossbarjarrodmillman
authored andcommitted
Brief accounting of eht-imaging and pycbc libraries.
Added a few sentences on how components of the scientific python ecosystem are used by pycbc and eht-imaging for the analysis of gravitational wave data and black hole imaging, respectively. Focused mostly on numpy, scipy, and matplotlib.
1 parent 39b84e0 commit eb2350b

File tree

2 files changed

+39
-0
lines changed

2 files changed

+39
-0
lines changed

body.tex

+27
Original file line numberDiff line numberDiff line change
@@ -165,6 +165,33 @@ \section*{Scientific Python ecosystem}
165165
% we should have an example of how these tools may all be used in an
166166
% example analysis. maybe one of the ones mentioned in the abstract?
167167

168+
The \tt{eht-imaging} library developed by the Event Horizon Telescope
169+
collabortion relies on many components of the scientific Python ecosystem.
170+
NumPy arrays were used to store and manipulate numerical data at every step
171+
in the processing chain: from raw data through calibration and image
172+
reconstruction.
173+
Scipy was used to provide tools used in general image processing such as
174+
filtering and image alignment, while scikit\_image provided higher-level
175+
functionality such as edge filters and Hough transforms.
176+
Other components of Scipy are used throughout the library, such as the
177+
\tt{scipy.optimize} package for handling general optimization tasks.
178+
Matplotlib is used for visualizing data throughout the analysis pipeline,
179+
including the generation of the final image of the black hole.
180+
% N.B. - most of this comes from introspecting the Image class defined in
181+
% ehtim/image.py.
182+
183+
Tools for analyzing data from gravitational wave observatories such as LIGO
184+
and Virgo are provided by the \tt{pycbc} package, which was used in the first
185+
detection of gravitationaly waves \cite{abbott2016observation} and in
186+
on-going analysis of data from LIGO and Virgo.
187+
\tt{pycbc} make extensive use of the scientific Python ecosystem.
188+
The time-series data from the interferometers are stored in NumPy arrays,
189+
while the \tt{scipy.signal} is used to construct and apply filters to the
190+
data.
191+
Matplotlib is used to visualize data at all points along the analysis chain
192+
including raw data from the instrumentation and the time-frequency
193+
visualization of the ``chirp'' from the binary black hole merger.
194+
168195
Exposing array programming primitives, as well as the surrounding ecosystem of
169196
tools, in Python---an interpreted language---creates an ideal environment for
170197
interactive, exploratory data analysis where users may iteratively inspect,

references.bib

+12
Original file line numberDiff line numberDiff line change
@@ -2213,6 +2213,18 @@ @article{eht-imaging
22132213
publisher={IOP Publishing}
22142214
}
22152215

2216+
@article{abbott2016observation,
2217+
title={Observation of gravitational waves from a binary black hole merger},
2218+
author={Abbott, Benjamin P and Abbott, Richard and Abbott, TD and Abernathy, MR and Acernese, Fausto and Ackley, Kendall and Adams, Carl and Adams, Thomas and Addesso, Paolo and Adhikari, RX and others},
2219+
journal={Physical review letters},
2220+
volume={116},
2221+
number={6},
2222+
pages={061102},
2223+
year={2016},
2224+
publisher={APS}
2225+
}
2226+
2227+
22162228
@incollection{pytorch,
22172229
title = {PyTorch: An Imperative Style, High-Performance Deep Learning Library},
22182230
author = {Paszke, Adam and Gross, Sam and Massa, Francisco and Lerer, Adam and Bradbury, James and Chanan, Gregory and Killeen, Trevor and Lin, Zeming and Gimelshein, Natalia and Antiga, Luca and Desmaison, Alban and Kopf, Andreas and Yang, Edward and DeVito, Zachary and Raison, Martin and Tejani, Alykhan and Chilamkurthy, Sasank and Steiner, Benoit and Fang, Lu and Bai, Junjie and Chintala, Soumith},

0 commit comments

Comments
 (0)