Install Insaseis on MacOS (py3)

conda裡面的pysiede2沒有包qt binary,要執行instaseis的GUI會造成.ui無法編譯。
若用pip去裝pyside2,生活變得更簡單。 cela vous simplifiera la vie

If we want to launch instaseis GUI, install Pyside2 from pip instead of conda.

python3 environment
conda create --name py37 python=3.7
conda activate py37

pip install instaseis
pip install PySide2
pip install pyqtgraph
conda install -c anaconda basemap
python -m instaseis.gui

UPDATE 18/11/2020:
On macOS Big Sur, add the following command to make qt work again.
export QT_MAC_WANTS_LAYER=1

Phase picking model training video

This video shows a convolutional neural network transferring seismic wave into P phase probability density function during the training process. The predict label looks like the original trace at first, drop down rapidly and slowly matches the training label.

Time-lapse video for delay-time kernel

Author: Li-Yu Kan (1) & Li Zhao (2)
(1) Department of Earth Sciences, National Central University, Taoyuan, Taiwan
(2) School of Earth and Space Sciences, Peking University, Beijing, China

A movie showing the Vp and Vs delay-time kernel. Made with mc kernel.

Original:

Source depth Station Distance Link
500 km 105 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep500_105deg.mp4
100 km 150 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep100_150deg.mp4
20 km 60 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep20_60deg.mp4

High-pass 100s filtered (zero phase)

Source depth Station Distance Link
500 km 105 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep500_105deg_hp.mp4
500 km
(Radial)
130 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep500_130deg_hpR.mp4
100 km 150 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep100_150deg_hp.mp4
100 km 180 deg https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep100_180deg_hp.mp4
20 km 105 deg (not finished)
https://sgylab.gl.ntu.edu.tw/kernel_video/PS_dep20_102deg_hp.mp4

Draw the outline of the text in GMT

While we choose a bright color for filling the text, a black outline is usually used to emphasize it (or it is unnoticeable). For plotting text on the canvas, we would immediately come up with pstext. But when we want to depict the outline of the text, there is no intuitive way or independent arguments for this. The -G and -W of pstext are actually for the background textbox rather than the text itself. Here, we provide two methods to achieve this: append =<pen> to the font specification in pstext or use psxy -Sl to plot the outline of the text.

閱讀全文

Plot SAC with pssac in the cartesian time axis

photo_2019-03-20 18.19.24

After the release of GMT5.3.0, pssac has been included into the GMT package. PSSAC operates a more convenient way than pswiggle or psxy to plot with SAC file (a common file format used in seismology community). We no longer have to convert SAC files into XY files and plot them by sorting them manually. PSSAC provides several types of regular profile like sorting by azimuth, epicentral distance, or user-defined profile through -E arguments. However, there is no direct argument in -E to plot such a figure with a cartesian time of X-axis (or Y if given -Q). In this article, we provide the trick for plotting SAC file in an axis of time.

閱讀全文