# Posts by Tag

## python

#### Understanding the metropolis-hastings algorithm: a step-by-step guide with python examples

The Metropolis-Hastings algorithm is a cornerstone of Markov Chain Monte Carlo (MCMC) methods, enabling us to generate samples from complex probability distr...

#### Exploring the multi-taper method: a robust approach for time-frequency analysis

Explore the Multi-Taper Method’s unique ability to refine spectral estimates in seismology. Leveraging multiple orthogonal tapers, this approach minimizes va...

#### What is pink noise and how is it different from white noise?

Pink noise, also known as 1/f noise or flicker noise, is a type of random signal that has equal energy per octave. It is called pink because it is analogous ...

#### Object oriented programming (oop) with examples in python

This post explains the concepts of object-oriented programming (OOP) and its key features, including encapsulation, abstraction, inheritance, and polymorphis...

#### Using mpi4py for parallel computing in python on supercomputers

The article discusses the mpi4py module, which is a Python wrapper for MPI, used for parallel computing in supercomputing environments. It provides higher-le...

#### How to plot the boundaries of the states on the basemap of the usa

We will plot the boundaries of the states of the USA on a basemap figure

#### Read yaml input file in bash, c/c++ and python

We will see how to read a YAML file in Bash, C/C++ and Python.

#### What are linked lists in data structures?

We look into a quick overview of the idea of linked list data structure with some examples.

#### Numerical optimization based on the l-bfgs method

We will inspect the L-BFGS optimization method using one example and compare its performance with the gradient descent method

#### Create beamer slides using python

We will use the Python package to create beamer presentation and append existing figures to each slide

#### How to deal with irregularly sampled time series data

While analyzing time series data, we often come across data that is non-uniformly sampled, i.e., they have non-equidistant time-steps. Infact, most of the re...

#### The easiest way to publish a python package on pypi using poetry

We will learn how to use the poetry package to create, manage, build and publish Python package on PyPI

#### Encrypt your data and software with python

Learn to encrypt and decrypt any files, data or software with python

#### Perform distributed computing easily using ray in python

We will introduce the concepts of distributed computing and then use the open-source Python library Ray to write scalable codes that can work on distributed ...

#### Convolution of waveforms using python compared to matlab

We will see compare the convolution functions in Python (Numpy) with the conv function in MATLAB. If you have tried them both then you would know that its no...

#### Convert any text to lifelike speech using amazon polly

We will see how we can use amazon web services, specifically amazon polly to convert any text into a speech

#### Maximum likelihood estimation for the regression parameters

We will learn the basics of the maximum likelihood method, and then apply it on a regression problem. We will also compare it with the least-squares estimati...

#### Efficiently compute spectrogram for large dataset in python

Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib

#### Numerically solving initial value problems using the runge-kutta method

Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...

#### Numerically solving initial value problems using the euler method

The simplest algorithm to solve a system of differential equations is the Euler method. We understand the Euler method by looking into a simple heat transfer...

#### How to plot shear-wave splitting measurements using pygmt

How you can plot the shear-wave splitting measurements from splitting database using PyGMT.

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

#### How to host a http server on your linux computer to share files

Is the signal denoising using a wavelet effective? We will implement the wavelet-based signal denoising using the MATLAB

#### Make your python script executable from anywhere in linux

Follow the instructions to make your python script executable from anywhere in Linux system.

#### How to securely send data from local client to remote server using paramiko (codes included)

Paramiko module can be used in Python to securely send data from the local client to the remote server. It is analogous to the SSH and SCP in Linux.

#### Concatenating daily seismic traces into one miniseed file (codes included)

I concatenate the daily seismic traces for 15 days into one miniseed file for further analysis. Then I obtained the spectrogram of the 15 days seismic data. ...

#### Computing cross-correlation and spectrogram of two seismic traces (codes included)

Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram

#### Principal component analysis to decompose signals and reduce dimensionality (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### How to analyze a huge data file with pandas (codes included)

We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib

#### Pyqt5 application for visualizing ambient sound in real-time (codes included)

A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.

#### How to start using pandas immediately for earth data analysis (codes included)

This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...

#### Introduction to numpy for scientific computing

This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### The right way to loop in python (codes included)

What is the fastest and most efficient way to loop in Python. We found that the numpy is fastest and python builtins are the most memory efficient.

#### Easily integrate custom functions in matlab with python (codes included)

How can we use the MATLAB functions in Python? MATLAB implementation are usually reliable as it is developed by the professionals. But the advantages of usin...

#### The easy way to compute and visualize the time & frequency domain correlation (codes included)

In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...

#### Monte carlo methods and earthquake location problem (codes included)

The common geophysical problems most often have multimodal objective function with many possible minima. In this post, we will look into the Monte Carlo meth...

#### A simple system monitor app in python (codes included)

less than 1 minute read      DESKTOPAPPS

A simple Python app for system CPU and RAM usage monitoring in real time. The app is build mainly using the PyQt5, pyqtgraph and psutil.

#### Build a flask web application: sea level rise monitoring (codes included)

This post gives a quick introduction on how to build a web application using Flask and deploy on Heroku server. Then, I share my codes for building advanced ...

#### Least-squares method in geosciences (codes included)

Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

#### Hypothesis test for the significance of linear trend using the monte carlo simulations (codes included)

We pose a null hypothesis and enquire that given that the null hypothesis is true, how likely is the observed pattern of results? This likelihood is known as...

#### How to avoid common mistakes in analyzing correlations of two time-series

Most often data analyst consider correlation between two time-series as a causation effect. Two time-series are correlated that does not imply that one cause...

#### Estimation of the degrees of freedom for time series (codes included)

Numerical tests to understand the degrees of freedom in time series

#### Exploratory factor analysis (codes included)

Factor Analysis is an exploratory data analysis method used to search influential underlying factors or latent variables from a set of observed variables

Obspy is an open-source Python framework developed for the processing of seismological data. In this post, I will introduce how to use Obspy along with some ...

#### Time series analysis in python: filtering or smoothing data (codes included)

In this post, we will see how we can use Python to low-pass filter the 10 year long daily fluctuations of GPS time series. We need to use the “Scipy” package...

#### Topographic map clipped by coastlines in python (codes included)

This post demonstrate how to use Python to set up clip topographic map based on coastlines.

#### Topographic map with shading in python (codes included)

Generating GMT style shaded relief map in Python

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

#### Plotting 1 arc-minute global relief map in python (codes included)

Plotting 1 arc-minute topographic map in python

#### Estimation of the degrees of freedom of time series in python (codes included)

Two time series with predominant linear trends (very low DOF) can have a very high correlation coefficient, which can hardly be construed as an evidence for ...

#### Interactive data visualization in python using bokeh

A basic to advanced guide to making interactive plots in Bokeh.

#### Why we should use generators? [python]

Generators don’t hold the entire result in memory. It yields one result at a time.

#### Introduction to python for beginners

In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.

## numericalMethodsScientificComputation

#### Understanding the metropolis-hastings algorithm: a step-by-step guide with python examples

The Metropolis-Hastings algorithm is a cornerstone of Markov Chain Monte Carlo (MCMC) methods, enabling us to generate samples from complex probability distr...

#### Numerical optimization based on the l-bfgs method

We will inspect the L-BFGS optimization method using one example and compare its performance with the gradient descent method

#### How to deal with irregularly sampled time series data

While analyzing time series data, we often come across data that is non-uniformly sampled, i.e., they have non-equidistant time-steps. Infact, most of the re...

#### Convolution of waveforms using python compared to matlab

We will see compare the convolution functions in Python (Numpy) with the conv function in MATLAB. If you have tried them both then you would know that its no...

#### Maximum likelihood estimation for the regression parameters

We will learn the basics of the maximum likelihood method, and then apply it on a regression problem. We will also compare it with the least-squares estimati...

#### Solving boundary value problems using the shooting method

The boundary value problems require information at the present time and a future time. We will see how we can use shooting method to solve problems where we ...

#### Numerically solving initial value problems using the runge-kutta method

Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...

#### Numerically solving initial value problems using the euler method

The simplest algorithm to solve a system of differential equations is the Euler method. We understand the Euler method by looking into a simple heat transfer...

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

#### Principal component analysis to decompose signals and reduce dimensionality (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### Genetic algorithm: a highly robust inversion scheme for geophysical applications (codes included)

An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem

#### Monte carlo methods and earthquake location problem (codes included)

The common geophysical problems most often have multimodal objective function with many possible minima. In this post, we will look into the Monte Carlo meth...

#### Least-squares method in geosciences (codes included)

Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...

#### Numerical tests for seismic resolution (codes included)

Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...

#### Exploratory factor analysis (codes included)

Factor Analysis is an exploratory data analysis method used to search influential underlying factors or latent variables from a set of observed variables

#### Monte carlo simulation to test for the correlation between two dataset in matlab (codes included)

We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...

#### Hypothesis test in matlab (codes included)

Using Randomization to test the disprove the null hypothesis

## obspyTutorial

#### Efficiently compute spectrogram for large dataset in python

Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib

#### How to plot shear-wave splitting measurements using pygmt

How you can plot the shear-wave splitting measurements from splitting database using PyGMT.

#### How to automatically enquire the availability of seismic data using obspy

In this post, we will see how can we retrive the available seismic waveforms information for a given network, station, channel and client in a given period o...

#### Concatenating daily seismic traces into one miniseed file (codes included)

I concatenate the daily seismic traces for 15 days into one miniseed file for further analysis. Then I obtained the spectrogram of the 15 days seismic data. ...

#### Computing cross-correlation and spectrogram of two seismic traces (codes included)

Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram

#### Signal denoising using fourier analysis in python (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### Analyzing miniseed seismic data in matlab (codes included)

We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB

#### The easy way to compute and visualize the time & frequency domain correlation (codes included)

In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Write ascii data to mseed file using obspy (codes included)

In this post, I will read a ASCII file whose first few lines contains the header information and then the three-component data. I will read using the pandas ...

#### Numerical tests for seismic resolution (codes included)

Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...

#### Locating earthquakes using geiger’s method (codes included)

Simple earthquake location problem and its solution using Geiger’s method

#### Working with obspy and basemap (codes included)

This post is aimed to resolve the issues regarding the conflicts of using obspy and basemap libraries together.

Obspy is an open-source Python framework developed for the processing of seismological data. In this post, I will introduce how to use Obspy along with some ...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

#### Automatically plotting record section for an earthquake in the given time range in python (codes included)

Python code to automatically plot the record section for the highest magnitude earthquake in the given time range

## geophysics

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

#### Automatically plotting record section for an earthquake in the given time range in python (codes included)

Python code to automatically plot the record section for the highest magnitude earthquake in the given time range

#### Numerical tests on travel time tomography [matlab]

Introduction to the concepts of tomography with equations and codes. Introduction to the concepts of overdetermined, underdetermined and mix-determined probl...

#### Monte carlo simulation to test for the correlation between two dataset in matlab (codes included)

We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...

#### Easy statistical analysis using the tools of matlab (codes included)

Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom

#### Hypothesis test in matlab (codes included)

Using Randomization to test the disprove the null hypothesis

#### Topographic map with shading in python (codes included)

Generating GMT style shaded relief map in Python

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

#### Estimation of the degrees of freedom of time series in python (codes included)

Two time series with predominant linear trends (very low DOF) can have a very high correlation coefficient, which can hardly be construed as an evidence for ...

#### Gmt advanced tutorial, part ii (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

#### Why we should use generators? [python]

Generators don’t hold the entire result in memory. It yields one result at a time.

#### Gmt advanced tutorial (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

## geospatial data visualization

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### High-quality maps using the modern interface to the generic mapping tools (codes included)

GMT or generic mapping tools have become synonymous with plotting maps in Earth, Ocean, and Planetary sciences. It can be used for processing data, generatin...

#### Gmt tutorial for beginners (codes included)

The Generic Mapping Tools is widely used across Earth and Planetary and other fields of studies to process data and generate high-quality illustrations. This...

#### Reading netcdf4 data in python (codes included)

In Earth Sciences, we often deal with multidimensional data structures such as climate data, GPS data. It ‘s hard to save such data in text files as it would...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

#### Topographic map clipped by coastlines in python (codes included)

This post demonstrate how to use Python to set up clip topographic map based on coastlines.

#### Topographic map with shading in python (codes included)

Generating GMT style shaded relief map in Python

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

#### Plotting 1 arc-minute global relief map in python (codes included)

Plotting 1 arc-minute topographic map in python

#### Gmt advanced tutorial, part ii (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

#### Interactive data visualization in python using bokeh

A basic to advanced guide to making interactive plots in Bokeh.

#### Gmt advanced tutorial (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

## obspy

#### How to automatically enquire the availability of seismic data using obspy

In this post, we will see how can we retrive the available seismic waveforms information for a given network, station, channel and client in a given period o...

#### Concatenating daily seismic traces into one miniseed file (codes included)

I concatenate the daily seismic traces for 15 days into one miniseed file for further analysis. Then I obtained the spectrogram of the 15 days seismic data. ...

#### Computing cross-correlation and spectrogram of two seismic traces (codes included)

Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram

#### Computing cross-correlation and spectrogram of two seismic traces (codes included)

Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram

#### Signal denoising using fourier analysis in python (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### Analyzing miniseed seismic data in matlab (codes included)

We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Write ascii data to mseed file using obspy (codes included)

In this post, I will read a ASCII file whose first few lines contains the header information and then the three-component data. I will read using the pandas ...

#### Working with obspy and basemap (codes included)

This post is aimed to resolve the issues regarding the conflicts of using obspy and basemap libraries together.

Obspy is an open-source Python framework developed for the processing of seismological data. In this post, I will introduce how to use Obspy along with some ...

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

#### Automatically plotting record section for an earthquake in the given time range in python (codes included)

Python code to automatically plot the record section for the highest magnitude earthquake in the given time range

## introEarthDataAnalysis

#### How to analyze a huge data file with pandas (codes included)

We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib

#### How to start using pandas immediately for earth data analysis (codes included)

This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...

#### Introduction to numpy for scientific computing

This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### High-quality maps using the modern interface to the generic mapping tools (codes included)

GMT or generic mapping tools have become synonymous with plotting maps in Earth, Ocean, and Planetary sciences. It can be used for processing data, generatin...

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

#### Time series analysis in python: filtering or smoothing data (codes included)

In this post, we will see how we can use Python to low-pass filter the 10 year long daily fluctuations of GPS time series. We need to use the “Scipy” package...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Topographic map clipped by coastlines in python (codes included)

This post demonstrate how to use Python to set up clip topographic map based on coastlines.

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

#### Introduction to python for beginners

In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.

## matlab

#### Convolution of waveforms using python compared to matlab

We will see compare the convolution functions in Python (Numpy) with the conv function in MATLAB. If you have tried them both then you would know that its no...

#### Concatenating daily seismic traces into one miniseed file (codes included)

I concatenate the daily seismic traces for 15 days into one miniseed file for further analysis. Then I obtained the spectrogram of the 15 days seismic data. ...

#### Analyzing miniseed seismic data in matlab (codes included)

We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB

#### Transfer learning applied on the unsplash data using alexnet pretrained network (codes included)

Transfer learning using the pre-trained deep learning networks from MATLAB can be easily implemented to achieve fast and impressive results

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

#### Simple wave modeling and hilbert transform in matlab (codes included)

We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection

#### Monte carlo simulation to test for the correlation between two dataset in matlab (codes included)

We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...

#### Easy statistical analysis using the tools of matlab (codes included)

Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom

#### Hypothesis test in matlab (codes included)

Using Randomization to test the disprove the null hypothesis

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

## pygmt

#### How to plot great circle path through your region using pygmt

In this article, we will learn how to visualize the great circle paths that traverse a designated region of interest using PyGMT. Understanding these paths i...

#### A quick overview on geospatial data visualization using pygmt

We will see how to plot a topographic map, overlay earthquake data on topographic maps, plot focal mechanism solutions and plot tomographic results on a geog...

#### How to plot earthquakes data on a three-dimensional topographic map

Read the earthquake data from a data file and overlay on a three-dimensional topographic map using PyGMT.

#### How to overlay shapefile data on pygmt maps

We learn how to plot selected shapefile data using geopandas on top of PyGMT maps

#### How to plot shear-wave splitting measurements using pygmt

How you can plot the shear-wave splitting measurements from splitting database using PyGMT.

#### How to extract station data from the web page

Uses pandas to read the html page and extract the html data into pandas dataframe

#### How to plot the earthquakes data on a topographic map

Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.

#### How to automatically enquire the availability of seismic data using obspy

In this post, we will see how can we retrive the available seismic waveforms information for a given network, station, channel and client in a given period o...

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

## visualization

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Time series analysis in geophysics (codes included)

Time-series analysis is essential in most fields of science, including geophysics, economics, etc. Most of the geophysical data comes in a time-series format...

#### Gmt tutorial for beginners (codes included)

The Generic Mapping Tools is widely used across Earth and Planetary and other fields of studies to process data and generate high-quality illustrations. This...

Obspy is an open-source Python framework developed for the processing of seismological data. In this post, I will introduce how to use Obspy along with some ...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

#### Automatically plotting record section for an earthquake in the given time range in python (codes included)

Python code to automatically plot the record section for the highest magnitude earthquake in the given time range

#### Gmt advanced tutorial, part ii (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

## seismology

#### How to plot shear-wave splitting measurements using pygmt

How you can plot the shear-wave splitting measurements from splitting database using PyGMT.

#### Concatenating daily seismic traces into one miniseed file (codes included)

I concatenate the daily seismic traces for 15 days into one miniseed file for further analysis. Then I obtained the spectrogram of the 15 days seismic data. ...

#### Analyzing miniseed seismic data in matlab (codes included)

We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB

#### The easy way to compute and visualize the time & frequency domain correlation (codes included)

In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...

#### Locating earthquakes using geiger’s method (codes included)

Simple earthquake location problem and its solution using Geiger’s method

Obspy is an open-source Python framework developed for the processing of seismological data. In this post, I will introduce how to use Obspy along with some ...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

#### Automatically plotting record section for an earthquake in the given time range in python (codes included)

Python code to automatically plot the record section for the highest magnitude earthquake in the given time range

## time-frequency analysis

#### Exploring the multi-taper method: a robust approach for time-frequency analysis

Explore the Multi-Taper Method’s unique ability to refine spectral estimates in seismology. Leveraging multiple orthogonal tapers, this approach minimizes va...

#### Efficiently compute spectrogram for large dataset in python

Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib

#### How effective is the signal denoising using the matlab based wavelet analysis

We see how to download seismic waveforms, convert them into mat format from mini-seed and then perform denoising using wavelet analysis. We first performed w...

#### Towards multi-resolution analysis with wavelet transform

We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.

#### Computing cross-correlation and spectrogram of two seismic traces (codes included)

Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram

#### Principal component analysis to decompose signals and reduce dimensionality (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### Signal denoising using fourier analysis in python (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### Analyzing miniseed seismic data in matlab (codes included)

We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB

#### The easy way to compute and visualize the time & frequency domain correlation (codes included)

In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...

## statistics

#### Understanding the metropolis-hastings algorithm: a step-by-step guide with python examples

The Metropolis-Hastings algorithm is a cornerstone of Markov Chain Monte Carlo (MCMC) methods, enabling us to generate samples from complex probability distr...

#### Hypothesis test for the significance of linear trend using the monte carlo simulations (codes included)

We pose a null hypothesis and enquire that given that the null hypothesis is true, how likely is the observed pattern of results? This likelihood is known as...

#### How to avoid common mistakes in analyzing correlations of two time-series

Most often data analyst consider correlation between two time-series as a causation effect. Two time-series are correlated that does not imply that one cause...

#### Estimation of the degrees of freedom for time series (codes included)

Numerical tests to understand the degrees of freedom in time series

#### Numerical tests for seismic resolution (codes included)

Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...

#### Monte carlo simulation to test for the correlation between two dataset in matlab (codes included)

We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...

#### Easy statistical analysis using the tools of matlab (codes included)

Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom

#### Hypothesis test in matlab (codes included)

Using Randomization to test the disprove the null hypothesis

## geospatial data

#### How to plot great circle path through your region using pygmt

In this article, we will learn how to visualize the great circle paths that traverse a designated region of interest using PyGMT. Understanding these paths i...

#### How to plot the boundaries of the states on the basemap of the usa

We will plot the boundaries of the states of the USA on a basemap figure

#### A quick overview on geospatial data visualization using pygmt

We will see how to plot a topographic map, overlay earthquake data on topographic maps, plot focal mechanism solutions and plot tomographic results on a geog...

#### How to plot earthquakes data on a three-dimensional topographic map

Read the earthquake data from a data file and overlay on a three-dimensional topographic map using PyGMT.

#### How to overlay shapefile data on pygmt maps

We learn how to plot selected shapefile data using geopandas on top of PyGMT maps

#### How to extract station data from the web page

Uses pandas to read the html page and extract the html data into pandas dataframe

#### How to plot the earthquakes data on a topographic map

Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

## GMT

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### High-quality maps using the modern interface to the generic mapping tools (codes included)

GMT or generic mapping tools have become synonymous with plotting maps in Earth, Ocean, and Planetary sciences. It can be used for processing data, generatin...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Topographic map with shading in python (codes included)

Generating GMT style shaded relief map in Python

#### Gmt advanced tutorial, part ii (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

#### Gmt advanced tutorial (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

## geospatial data visualization python

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### Topographic map clipped by coastlines in python (codes included)

This post demonstrate how to use Python to set up clip topographic map based on coastlines.

#### Topographic map with shading in python (codes included)

Generating GMT style shaded relief map in Python

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

#### Plotting 1 arc-minute global relief map in python (codes included)

Plotting 1 arc-minute topographic map in python

## parallel computing

#### Mastering large data processing with mpi4py in python

Explore the power of mpi4py for simplifying data distribution in parallel computing with its efficient broadcast functionality, which seamlessly sends data f...

#### Using mpi4py for parallel computing in python on supercomputers

The article discusses the mpi4py module, which is a Python wrapper for MPI, used for parallel computing in supercomputing environments. It provides higher-le...

Using multiple threads in C for concurrent process flow

#### Perform distributed computing easily using ray in python

We will introduce the concepts of distributed computing and then use the open-source Python library Ray to write scalable codes that can work on distributed ...

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

## signal processing

#### Exploring the multi-taper method: a robust approach for time-frequency analysis

Explore the Multi-Taper Method’s unique ability to refine spectral estimates in seismology. Leveraging multiple orthogonal tapers, this approach minimizes va...

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

#### How effective is the signal denoising using the matlab based wavelet analysis

We see how to download seismic waveforms, convert them into mat format from mini-seed and then perform denoising using wavelet analysis. We first performed w...

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

#### Principal component analysis to decompose signals and reduce dimensionality (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

#### Signal denoising using fourier analysis in python (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

## techniques

#### Hypothesis test for the significance of linear trend using the monte carlo simulations (codes included)

We pose a null hypothesis and enquire that given that the null hypothesis is true, how likely is the observed pattern of results? This likelihood is known as...

#### How to avoid common mistakes in analyzing correlations of two time-series

Most often data analyst consider correlation between two time-series as a causation effect. Two time-series are correlated that does not imply that one cause...

#### Estimation of the degrees of freedom for time series (codes included)

Numerical tests to understand the degrees of freedom in time series

#### Numerical tests for seismic resolution (codes included)

Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

## apple

#### How to install jekyll on apple m1 macbook

If you have recently bought a M1 mac and have been doing blogging using Jekyll, then you must have experienced that installing Jekyll on the M1 architecture ...

#### Some handy mac tweaks [macos]

Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

#### Iphone shortcut to quickly log body temperature [ios]

Shortcut code for quickly logging temperature in Apple health app

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## matplotlib

#### How to plot the boundaries of the states on the basemap of the usa

We will plot the boundaries of the states of the USA on a basemap figure

#### Efficiently compute spectrogram for large dataset in python

Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib

#### Computing cross-correlation and spectrogram of two seismic traces (codes included)

Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

## datascience

#### Mastering large data processing with mpi4py in python

Explore the power of mpi4py for simplifying data distribution in parallel computing with its efficient broadcast functionality, which seamlessly sends data f...

#### Using mpi4py for parallel computing in python on supercomputers

The article discusses the mpi4py module, which is a Python wrapper for MPI, used for parallel computing in supercomputing environments. It provides higher-le...

#### Perform distributed computing easily using ray in python

We will introduce the concepts of distributed computing and then use the open-source Python library Ray to write scalable codes that can work on distributed ...

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

## time series analysis

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

#### Wavelet analysis applied to the real dataset in a quick and easy way (codes included)

An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!

#### Genetic algorithm: a highly robust inversion scheme for geophysical applications (codes included)

An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem

#### Least-squares method in geosciences (codes included)

Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...

## plotting map

#### How to plot the boundaries of the states on the basemap of the usa

We will plot the boundaries of the states of the USA on a basemap figure

#### How to plot earthquakes data on a three-dimensional topographic map

Read the earthquake data from a data file and overlay on a three-dimensional topographic map using PyGMT.

#### How to overlay shapefile data on pygmt maps

We learn how to plot selected shapefile data using geopandas on top of PyGMT maps

#### How to extract station data from the web page

Uses pandas to read the html page and extract the html data into pandas dataframe

#### How to plot the earthquakes data on a topographic map

Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.

## time-series

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

#### Hypothesis test for the significance of linear trend using the monte carlo simulations (codes included)

We pose a null hypothesis and enquire that given that the null hypothesis is true, how likely is the observed pattern of results? This likelihood is known as...

#### How to avoid common mistakes in analyzing correlations of two time-series

Most often data analyst consider correlation between two time-series as a causation effect. Two time-series are correlated that does not imply that one cause...

#### Simple wave modeling and hilbert transform in matlab (codes included)

We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection

## MATLAB

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

#### Easily integrate custom functions in matlab with python (codes included)

How can we use the MATLAB functions in Python? MATLAB implementation are usually reliable as it is developed by the professionals. But the advantages of usin...

#### Numerical tests for seismic resolution (codes included)

Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...

#### Locating earthquakes using geiger’s method (codes included)

Simple earthquake location problem and its solution using Geiger’s method

## optimization method

#### Maximum likelihood estimation for the regression parameters

We will learn the basics of the maximum likelihood method, and then apply it on a regression problem. We will also compare it with the least-squares estimati...

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

#### Genetic algorithm: a highly robust inversion scheme for geophysical applications (codes included)

An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem

#### Least-squares method in geosciences (codes included)

Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...

## iterative methods

#### Solving boundary value problems using the shooting method

The boundary value problems require information at the present time and a future time. We will see how we can use shooting method to solve problems where we ...

#### Numerically solving initial value problems using the runge-kutta method

Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...

#### Numerically solving initial value problems using the euler method

The simplest algorithm to solve a system of differential equations is the Euler method. We understand the Euler method by looking into a simple heat transfer...

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

## gmt mapping

#### Gmt tutorial for beginners (codes included)

The Generic Mapping Tools is widely used across Earth and Planetary and other fields of studies to process data and generate high-quality illustrations. This...

#### Gmt advanced tutorial, part ii (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

#### Gmt advanced tutorial (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

## data science

#### What are linked lists in data structures?

We look into a quick overview of the idea of linked list data structure with some examples.

#### Estimation of the degrees of freedom of time series in python (codes included)

Two time series with predominant linear trends (very low DOF) can have a very high correlation coefficient, which can hardly be construed as an evidence for ...

#### Why we should use generators? [python]

Generators don’t hold the entire result in memory. It yields one result at a time.

## time series

#### Estimation of the degrees of freedom for time series (codes included)

Numerical tests to understand the degrees of freedom in time series

#### Time series analysis in geophysics (codes included)

Time-series analysis is essential in most fields of science, including geophysics, economics, etc. Most of the geophysical data comes in a time-series format...

#### Estimation of the degrees of freedom of time series in python (codes included)

Two time series with predominant linear trends (very low DOF) can have a very high correlation coefficient, which can hardly be construed as an evidence for ...

## topography

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

#### Plotting 1 arc-minute global relief map in python (codes included)

Plotting 1 arc-minute topographic map in python

## mac

#### Some handy mac tweaks [macos]

Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## basemap

#### How to plot the boundaries of the states on the basemap of the usa

We will plot the boundaries of the states of the USA on a basemap figure

#### Working with obspy and basemap (codes included)

This post is aimed to resolve the issues regarding the conflicts of using obspy and basemap libraries together.

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

## MS word

#### Writing and formatting a scientific manuscript in microsoft word

less than 1 minute read      UTILITIES

If you are ready to use the Microsoft Word as your favourite tool for writing your awesome scientific thoughts and ideas into a manuscript, then I would like...

#### How to insert equation numbering in ms word

It is essential to insert equation numbers in your thesis and/or any scientific paper. In this post, I will show you some of the easiest ways to insert equat...

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## microsoft word

#### Writing and formatting a scientific manuscript in microsoft word

less than 1 minute read      UTILITIES

If you are ready to use the Microsoft Word as your favourite tool for writing your awesome scientific thoughts and ideas into a manuscript, then I would like...

#### How to insert equation numbering in ms word

It is essential to insert equation numbers in your thesis and/or any scientific paper. In this post, I will show you some of the easiest ways to insert equat...

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## writing manuscript

#### Writing and formatting a scientific manuscript in microsoft word

less than 1 minute read      UTILITIES

If you are ready to use the Microsoft Word as your favourite tool for writing your awesome scientific thoughts and ideas into a manuscript, then I would like...

#### How to insert equation numbering in ms word

It is essential to insert equation numbers in your thesis and/or any scientific paper. In this post, I will show you some of the easiest ways to insert equat...

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## equations

#### Writing and formatting a scientific manuscript in microsoft word

less than 1 minute read      UTILITIES

If you are ready to use the Microsoft Word as your favourite tool for writing your awesome scientific thoughts and ideas into a manuscript, then I would like...

#### How to insert equation numbering in ms word

It is essential to insert equation numbers in your thesis and/or any scientific paper. In this post, I will show you some of the easiest ways to insert equat...

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## earthquake location problem

#### Genetic algorithm: a highly robust inversion scheme for geophysical applications (codes included)

An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem

#### Monte carlo methods and earthquake location problem (codes included)

The common geophysical problems most often have multimodal objective function with many possible minima. In this post, we will look into the Monte Carlo meth...

#### Least-squares method in geosciences (codes included)

Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...

## mseed

#### Concatenating daily seismic traces into one miniseed file (codes included)

I concatenate the daily seismic traces for 15 days into one miniseed file for further analysis. Then I obtained the spectrogram of the 15 days seismic data. ...

#### Analyzing miniseed seismic data in matlab (codes included)

We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB

#### Write ascii data to mseed file using obspy (codes included)

In this post, I will read a ASCII file whose first few lines contains the header information and then the three-component data. I will read using the pandas ...

## fourier transform

#### Exploring the multi-taper method: a robust approach for time-frequency analysis

Explore the Multi-Taper Method’s unique ability to refine spectral estimates in seismology. Leveraging multiple orthogonal tapers, this approach minimizes va...

#### Towards multi-resolution analysis with wavelet transform

We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.

#### Wavelet analysis applied to the real dataset in a quick and easy way (codes included)

An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!

## scientific computing

#### How to analyze a huge data file with pandas (codes included)

We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib

#### How to start using pandas immediately for earth data analysis (codes included)

This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...

#### Introduction to numpy for scientific computing

This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...

## spectrogram

#### Exploring the multi-taper method: a robust approach for time-frequency analysis

Explore the Multi-Taper Method’s unique ability to refine spectral estimates in seismology. Leveraging multiple orthogonal tapers, this approach minimizes va...

#### Efficiently compute spectrogram for large dataset in python

Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib

#### Towards multi-resolution analysis with wavelet transform

We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.

## initial value problem

#### Numerically solving initial value problems using the runge-kutta method

Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...

#### Numerically solving initial value problems using the euler method

The simplest algorithm to solve a system of differential equations is the Euler method. We understand the Euler method by looking into a simple heat transfer...

#### How to implement the iterative newton–raphson method to find roots of a function in python

The Newton–Raphson method (commonly known as Newton’s method) is developed for finding roots of a given function or polynomial iteratively. We show two examp...

## c

#### How to load dynamic libraries in c/c++?

We will see how to load shared libraries in C/C++. We will write a library to convert km to degrees and vice-versa. Then we create a utility program to conve...

#### What are linked lists in data structures?

We look into a quick overview of the idea of linked list data structure with some examples.

#### How to create a simple makefile in linux

We learn how to write a Makefile to automate the compilation of our source code. We will use one example from Fortran.

## clone

#### How to set up ssh-keys for github

You will learn how to generate and set up an SSH key for github so that you don’t need to always type your username and password when accessing github repo.

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## maps

#### Topographic map clipped by coastlines in python (codes included)

This post demonstrate how to use Python to set up clip topographic map based on coastlines.

#### Gmt advanced tutorial, part ii (codes included)

Plotting publishable quality maps using generic mapping tools (GMT)

## degrees of freedom

#### Estimation of the degrees of freedom for time series (codes included)

Numerical tests to understand the degrees of freedom in time series

#### Estimation of the degrees of freedom of time series in python (codes included)

Two time series with predominant linear trends (very low DOF) can have a very high correlation coefficient, which can hardly be construed as an evidence for ...

## quick action

#### Some handy mac tweaks [macos]

Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## macbook

#### How to install jekyll on apple m1 macbook

If you have recently bought a M1 mac and have been doing blogging using Jekyll, then you must have experienced that installing Jekyll on the M1 architecture ...

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## record section

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Automatically plotting record section for an earthquake in the given time range in python (codes included)

Python code to automatically plot the record section for the highest magnitude earthquake in the given time range

## geospatial data visualization python python

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

## epicentral distance

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

## waveforms

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

## Generic Mapping Tools

#### Gmt tutorial for beginners (codes included)

The Generic Mapping Tools is widely used across Earth and Planetary and other fields of studies to process data and generate high-quality illustrations. This...

#### Pygmt: high-resolution topographic map in python (codes included)

A simple tutorial on how to plot high resolution topographic map using GMT tools in Python

## netcdf

#### Writing netcdf4 data using python (codes included)

NetCDF file format has been designed for storing multidimensional scientific data such as temperature, rainfall, humidity, etc. In this post, we will see how...

#### Reading netcdf4 data in python (codes included)

In Earth Sciences, we often deal with multidimensional data structures such as climate data, GPS data. It ‘s hard to save such data in text files as it would...

## data analysis

#### Writing netcdf4 data using python (codes included)

NetCDF file format has been designed for storing multidimensional scientific data such as temperature, rainfall, humidity, etc. In this post, we will see how...

#### Time series analysis in python: filtering or smoothing data (codes included)

In this post, we will see how we can use Python to low-pass filter the 10 year long daily fluctuations of GPS time series. We need to use the “Scipy” package...

## spectral analysis

#### Exploring the multi-taper method: a robust approach for time-frequency analysis

Explore the Multi-Taper Method’s unique ability to refine spectral estimates in seismology. Leveraging multiple orthogonal tapers, this approach minimizes va...

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

Using multiple threads in C for concurrent process flow

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

## least-squares method

#### Maximum likelihood estimation for the regression parameters

We will learn the basics of the maximum likelihood method, and then apply it on a regression problem. We will also compare it with the least-squares estimati...

#### Least-squares method in geosciences (codes included)

Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...

## desktopApp

#### Pyqt5 application for visualizing ambient sound in real-time (codes included)

A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.

#### A simple system monitor app in python (codes included)

less than 1 minute read      DESKTOPAPPS

A simple Python app for system CPU and RAM usage monitoring in real time. The app is build mainly using the PyQt5, pyqtgraph and psutil.

## PyQt5

#### Pyqt5 application for visualizing ambient sound in real-time (codes included)

A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.

#### A simple system monitor app in python (codes included)

less than 1 minute read      DESKTOPAPPS

A simple Python app for system CPU and RAM usage monitoring in real time. The app is build mainly using the PyQt5, pyqtgraph and psutil.

## pyqtgraph

#### Pyqt5 application for visualizing ambient sound in real-time (codes included)

A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.

#### A simple system monitor app in python (codes included)

less than 1 minute read      DESKTOPAPPS

A simple Python app for system CPU and RAM usage monitoring in real time. The app is build mainly using the PyQt5, pyqtgraph and psutil.

## monte carlo methods

#### Understanding the metropolis-hastings algorithm: a step-by-step guide with python examples

The Metropolis-Hastings algorithm is a cornerstone of Markov Chain Monte Carlo (MCMC) methods, enabling us to generate samples from complex probability distr...

#### Monte carlo methods and earthquake location problem (codes included)

The common geophysical problems most often have multimodal objective function with many possible minima. In this post, we will look into the Monte Carlo meth...

## GMT6

#### Three-dimensional perspective map of taiwan using gmt and pygmt (codes included)

We learn how to make the three-dimensional map using both GMT and PyGMT

#### High-quality maps using the modern interface to the generic mapping tools (codes included)

GMT or generic mapping tools have become synonymous with plotting maps in Earth, Ocean, and Planetary sciences. It can be used for processing data, generatin...

## wavelets

#### How effective is the signal denoising using the matlab based wavelet analysis

We see how to download seismic waveforms, convert them into mat format from mini-seed and then perform denoising using wavelet analysis. We first performed w...

#### Wavelet analysis applied to the real dataset in a quick and easy way (codes included)

An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!

## indian monsoon rainfall data

#### Empirical orthogonal function analysis to inspect the spatial coherency in the geospatial data (codes included)

Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...

#### Wavelet analysis applied to the real dataset in a quick and easy way (codes included)

An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!

## pandas dataframe

#### How to analyze a huge data file with pandas (codes included)

We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib

#### How to start using pandas immediately for earth data analysis (codes included)

This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...

## denoising

#### How effective is the signal denoising using the matlab based wavelet analysis

We see how to download seismic waveforms, convert them into mat format from mini-seed and then perform denoising using wavelet analysis. We first performed w...

#### Signal denoising using fourier analysis in python (codes included)

We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals

## Runge Kutta method

#### Numerically solving initial value problems using the runge-kutta method

Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...

#### Numerically solving initial value problems using the euler method

The simplest algorithm to solve a system of differential equations is the Euler method. We understand the Euler method by looking into a simple heat transfer...

## compile

#### Read and write text files in modern fortran

We will read and write text and numeric data into a file using modern fortran.

#### How to create a simple makefile in linux

We learn how to write a Makefile to automate the compilation of our source code. We will use one example from Fortran.

## programs

#### Encrypt your data and software with python

Learn to encrypt and decrypt any files, data or software with python

#### How to create a simple makefile in linux

We learn how to write a Makefile to automate the compilation of our source code. We will use one example from Fortran.

## fortran

#### Read and write text files in modern fortran

We will read and write text and numeric data into a file using modern fortran.

#### How to create a simple makefile in linux

We learn how to write a Makefile to automate the compilation of our source code. We will use one example from Fortran.

## fortran 90

#### Read and write text files in modern fortran

We will read and write text and numeric data into a file using modern fortran.

#### How to create a simple makefile in linux

We learn how to write a Makefile to automate the compilation of our source code. We will use one example from Fortran.

## gmt

#### How to plot great circle path through your region using pygmt

In this article, we will learn how to visualize the great circle paths that traverse a designated region of interest using PyGMT. Understanding these paths i...

#### A quick overview on geospatial data visualization using pygmt

We will see how to plot a topographic map, overlay earthquake data on topographic maps, plot focal mechanism solutions and plot tomographic results on a geog...

## tomographic map

#### How to plot great circle path through your region using pygmt

In this article, we will learn how to visualize the great circle paths that traverse a designated region of interest using PyGMT. Understanding these paths i...

#### A quick overview on geospatial data visualization using pygmt

We will see how to plot a topographic map, overlay earthquake data on topographic maps, plot focal mechanism solutions and plot tomographic results on a geog...

## focal mechanism

#### How to plot great circle path through your region using pygmt

In this article, we will learn how to visualize the great circle paths that traverse a designated region of interest using PyGMT. Understanding these paths i...

#### A quick overview on geospatial data visualization using pygmt

We will see how to plot a topographic map, overlay earthquake data on topographic maps, plot focal mechanism solutions and plot tomographic results on a geog...

## data structures

#### How to load dynamic libraries in c/c++?

We will see how to load shared libraries in C/C++. We will write a library to convert km to degrees and vice-versa. Then we create a utility program to conve...

#### What are linked lists in data structures?

We look into a quick overview of the idea of linked list data structure with some examples.

## cpp

#### How to load dynamic libraries in c/c++?

We will see how to load shared libraries in C/C++. We will write a library to convert km to degrees and vice-versa. Then we create a utility program to conve...

#### What are linked lists in data structures?

We look into a quick overview of the idea of linked list data structure with some examples.

## C++

#### Read yaml input file in bash, c/c++ and python

We will see how to read a YAML file in Bash, C/C++ and Python.

Using multiple threads in C for concurrent process flow

## supercomputing

#### Mastering large data processing with mpi4py in python

Explore the power of mpi4py for simplifying data distribution in parallel computing with its efficient broadcast functionality, which seamlessly sends data f...

#### Using mpi4py for parallel computing in python on supercomputers

The article discusses the mpi4py module, which is a Python wrapper for MPI, used for parallel computing in supercomputing environments. It provides higher-le...

## high performance computing

#### Mastering large data processing with mpi4py in python

Explore the power of mpi4py for simplifying data distribution in parallel computing with its efficient broadcast functionality, which seamlessly sends data f...

#### Using mpi4py for parallel computing in python on supercomputers

The article discusses the mpi4py module, which is a Python wrapper for MPI, used for parallel computing in supercomputing environments. It provides higher-le...

## MEMS accelerometers

#### Integrating mems accelerometers and broadband seismometers: challenges and opportunities in modern seismology

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

#### Understanding mems accelerometers

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

#### Integrating mems accelerometers and broadband seismometers: challenges and opportunities in modern seismology

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

#### Understanding mems accelerometers

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

## seismic monitoring

#### Integrating mems accelerometers and broadband seismometers: challenges and opportunities in modern seismology

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

#### Understanding mems accelerometers

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

## low-frequency sensing

#### Integrating mems accelerometers and broadband seismometers: challenges and opportunities in modern seismology

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

#### Understanding mems accelerometers

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

## earthquake detection

#### Integrating mems accelerometers and broadband seismometers: challenges and opportunities in modern seismology

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

#### Understanding mems accelerometers

While MEMS accelerometers offer advantages in cost, size, and deployment flexibility compared to traditional broadband seismometers, they face limitations in...

## python3

#### Introduction to python for beginners

In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.

## anaconda

#### Introduction to python for beginners

In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.

## miniconda

#### Introduction to python for beginners

In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.

## git

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## github

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## merge

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## branch

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## github desktop

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## tutorial

#### Using github for team collaboration

Tutorial on how to use Git and GitHub for team collaboration on a project. Content includes installing, setting up, creating a repository, making commits, un...

## generators

#### Why we should use generators? [python]

Generators don’t hold the entire result in memory. It yields one result at a time.

## bokeh

#### Interactive data visualization in python using bokeh

A basic to advanced guide to making interactive plots in Bokeh.

## streaming data

#### Interactive data visualization in python using bokeh

A basic to advanced guide to making interactive plots in Bokeh.

## web data visualization

#### Interactive data visualization in python using bokeh

A basic to advanced guide to making interactive plots in Bokeh.

## noaa

#### Plotting 1 arc-minute global relief map in python (codes included)

Plotting 1 arc-minute topographic map in python

## global relief

#### Plotting 1 arc-minute global relief map in python (codes included)

Plotting 1 arc-minute topographic map in python

## non-linear

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

## curve-fitting

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

## optimization

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

## lsqcurvefit

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

## non-linear regression

#### Non-linear curve fitting to a model with multiple observational variables in matlab (codes included)

How to fit data to non-linear model

## station map

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

## topographic

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

## etopo1

#### How to plot topographic station map in python (codes included)

We read the station information from a file to plot the 1 arc-minute relief map

## mapping

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

## kriging

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

## pykrige

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

## Ordinary Kriging

#### Plotting the geospatial data clipped by coastlines in python (codes included)

In geosciences, we most frequently have to make geospatial plots, but the available data is unevenly distributed and irregular. We like to show the data, in ...

## clipped

#### Topographic map clipped by coastlines in python (codes included)

This post demonstrate how to use Python to set up clip topographic map based on coastlines.

## hypothesis testing

#### Hypothesis test in matlab (codes included)

Using Randomization to test the disprove the null hypothesis

## null hypothesis

#### Hypothesis test in matlab (codes included)

Using Randomization to test the disprove the null hypothesis

## histograms

#### Easy statistical analysis using the tools of matlab (codes included)

Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom

## monte-carlo simulations

#### Monte carlo simulation to test for the correlation between two dataset in matlab (codes included)

We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...

## correlation

#### Monte carlo simulation to test for the correlation between two dataset in matlab (codes included)

We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...

## tomography

#### Numerical tests on travel time tomography [matlab]

Introduction to the concepts of tomography with equations and codes. Introduction to the concepts of overdetermined, underdetermined and mix-determined probl...

## seismic tomography

#### Numerical tests on travel time tomography [matlab]

Introduction to the concepts of tomography with equations and codes. Introduction to the concepts of overdetermined, underdetermined and mix-determined probl...

## apple music

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## Rating song

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## love song

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## dislike song

#### Mac quick action to rate the songs in apple music app [macos]

Quick action for mac to easily love, dislike, rate songs in apple music app

## ios

#### Iphone shortcut to quickly log body temperature [ios]

Shortcut code for quickly logging temperature in Apple health app

## log data

#### Iphone shortcut to quickly log body temperature [ios]

Shortcut code for quickly logging temperature in Apple health app

## iphone

#### Iphone shortcut to quickly log body temperature [ios]

Shortcut code for quickly logging temperature in Apple health app

## health

#### Iphone shortcut to quickly log body temperature [ios]

Shortcut code for quickly logging temperature in Apple health app

## automation

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

## automator

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

## services

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

## quick actions

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

## applescripts

#### Introduction to automating mac [macos]

Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts

## screenshot

#### Some handy mac tweaks [macos]

Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc

## batch renaming

#### Some handy mac tweaks [macos]

Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc

## psd

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

## probabilistic power spectral density

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

## temporal plot

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

## spectrogram ppsd

#### Visualizing power spectral density using obspy in python (codes included)

Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075

## hurricane track

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

## typhoon track

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

## cyclone track

#### Plotting track and trajectory of tropical cyclones on a topographic map in python (codes included)

Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.

## P-arrival

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

## S-arrival

#### Plotting seismograms with increasing epicentral distance using python (codes included)

Short demostration of how to plot the distance vs seismic waveforms and mark the P and S arrival times using the IASP91 earth model. Codes are included.

## Ulysses

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## Papers

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## Mendeley

#### Writing scientific papers in ulysses

Ulysses is a natural, freestyle way of writing. If you got any idea, just write it down, worry about the format and other things when you’re done. Don’t let ...

## hierarchical data analysis

#### Reading netcdf4 data in python (codes included)

In Earth Sciences, we often deal with multidimensional data structures such as climate data, GPS data. It ‘s hard to save such data in text files as it would...

## xarray

#### Reading netcdf4 data in python (codes included)

In Earth Sciences, we often deal with multidimensional data structures such as climate data, GPS data. It ‘s hard to save such data in text files as it would...

## lowpass filter

#### Time series analysis in python: filtering or smoothing data (codes included)

In this post, we will see how we can use Python to low-pass filter the 10 year long daily fluctuations of GPS time series. We need to use the “Scipy” package...

## smoothing

#### Time series analysis in python: filtering or smoothing data (codes included)

In this post, we will see how we can use Python to low-pass filter the 10 year long daily fluctuations of GPS time series. We need to use the “Scipy” package...

## writing data

#### Writing netcdf4 data using python (codes included)

NetCDF file format has been designed for storing multidimensional scientific data such as temperature, rainfall, humidity, etc. In this post, we will see how...

## wave modeling

#### Simple wave modeling and hilbert transform in matlab (codes included)

We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection

## hilbert transform

#### Simple wave modeling and hilbert transform in matlab (codes included)

We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection

Obspy is an open-source Python framework developed for the processing of seismological data. In this post, I will introduce how to use Obspy along with some ...

## analysis

#### Time series analysis in geophysics (codes included)

Time-series analysis is essential in most fields of science, including geophysics, economics, etc. Most of the geophysical data comes in a time-series format...

#### Time series analysis in geophysics (codes included)

Time-series analysis is essential in most fields of science, including geophysics, economics, etc. Most of the geophysical data comes in a time-series format...

## factor analysis

#### Exploratory factor analysis (codes included)

Factor Analysis is an exploratory data analysis method used to search influential underlying factors or latent variables from a set of observed variables

## factor analysis python example

#### Exploratory factor analysis (codes included)

Factor Analysis is an exploratory data analysis method used to search influential underlying factors or latent variables from a set of observed variables

## earthquake location

#### Locating earthquakes using geiger’s method (codes included)

Simple earthquake location problem and its solution using Geiger’s method

## seismic resolution

#### Numerical tests for seismic resolution (codes included)

Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...

## cross-correlations

#### How to avoid common mistakes in analyzing correlations of two time-series

Most often data analyst consider correlation between two time-series as a causation effect. Two time-series are correlated that does not imply that one cause...

## monte-carlo

#### Hypothesis test for the significance of linear trend using the monte carlo simulations (codes included)

We pose a null hypothesis and enquire that given that the null hypothesis is true, how likely is the observed pattern of results? This likelihood is known as...

## Fourier Transform

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

## Welchs method

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

## Kaiser window

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

## plomb

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

## spectral analysis matlab

#### Time-frequency analysis in matlab (codes included)

A signal has one or more frequency components in it and can be viewed from two different standpoints: time-domain and frequency domain. In general, signals a...

## 2D plots

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

## bar plot

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

## log scale

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

## multiple plots

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

## histogram

#### Advanced 2d plots with matplotlib in python (codes included)

Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc

## multiprocessing

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

## joblib

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

## concurrent futures

#### Speed-up your codes by parallel computing in python (codes included)

Parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often only use single process to do most of...

## ascii

#### Write ascii data to mseed file using obspy (codes included)

In this post, I will read a ASCII file whose first few lines contains the header information and then the three-component data. I will read using the pandas ...

## text

#### Write ascii data to mseed file using obspy (codes included)

In this post, I will read a ASCII file whose first few lines contains the header information and then the three-component data. I will read using the pandas ...

#### Build a flask web application: sea level rise monitoring (codes included)

This post gives a quick introduction on how to build a web application using Flask and deploy on Heroku server. Then, I share my codes for building advanced ...

## webapp

#### Build a flask web application: sea level rise monitoring (codes included)

This post gives a quick introduction on how to build a web application using Flask and deploy on Heroku server. Then, I share my codes for building advanced ...

## deploy

#### Build a flask web application: sea level rise monitoring (codes included)

This post gives a quick introduction on how to build a web application using Flask and deploy on Heroku server. Then, I share my codes for building advanced ...

## heroku

#### Build a flask web application: sea level rise monitoring (codes included)

This post gives a quick introduction on how to build a web application using Flask and deploy on Heroku server. Then, I share my codes for building advanced ...

## slrm

#### Build a flask web application: sea level rise monitoring (codes included)

This post gives a quick introduction on how to build a web application using Flask and deploy on Heroku server. Then, I share my codes for building advanced ...

## seismograms

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...

## SAC

#### Plotting a record section using obspy (codes included)

Quickly plot record section of a stream using Obspy. I will introduce you how to make a stream from a set of SAC data, plot the record section and store it a...