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
Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib
Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...
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 you can plot the shear-wave splitting measurements from splitting database using PyGMT.
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...
Is the signal denoising using a wavelet effective? We will implement the wavelet-based signal denoising using the MATLAB
Follow the instructions to make your python script executable from anywhere in Linux system.
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.
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. ...
Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib
A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.
This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...
This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...
We learn how to make the three-dimensional map using both GMT and PyGMT
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.
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...
In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...
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 Python app for system CPU and RAM usage monitoring in real time. The app is build mainly using the PyQt5, pyqtgraph and psutil.
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 is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...
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...
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
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...
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...
Numerical tests to understand the degrees of freedom in time series
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 ...
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...
This post demonstrate how to use Python to set up clip topographic map based on coastlines.
Generating GMT style shaded relief map in Python
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 ...
We read the station information from a file to plot the 1 arc-minute relief map
Plotting 1 arc-minute topographic map in python
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 ...
A basic to advanced guide to making interactive plots in Bokeh.
Generators don’t hold the entire result in memory. It yields one result at a time.
In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.
Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib
How you can plot the shear-wave splitting measurements from splitting database using PyGMT.
Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.
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...
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. ...
Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB
In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...
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...
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 ...
Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...
Simple earthquake location problem and its solution using Geiger’s method
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 ...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
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.
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Python code to automatically plot the record section for the highest magnitude earthquake in the given time range
Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
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...
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.
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Python code to automatically plot the record section for the highest magnitude earthquake in the given time range
Introduction to the concepts of tomography with equations and codes. Introduction to the concepts of overdetermined, underdetermined and mix-determined probl...
We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...
Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom
Using Randomization to test the disprove the null hypothesis
Generating GMT style shaded relief map in Python
We read the station information from a file to plot the 1 arc-minute relief map
How to fit data to non-linear model
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 ...
Plotting publishable quality maps using generic mapping tools (GMT)
Generators don’t hold the entire result in memory. It yields one result at a time.
Plotting publishable quality maps using generic mapping tools (GMT)
We learn how to make the three-dimensional map using both GMT and PyGMT
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...
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...
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...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.
This post demonstrate how to use Python to set up clip topographic map based on coastlines.
Generating GMT style shaded relief map in Python
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 ...
We read the station information from a file to plot the 1 arc-minute relief map
Plotting 1 arc-minute topographic map in python
Plotting publishable quality maps using generic mapping tools (GMT)
A basic to advanced guide to making interactive plots in Bokeh.
Plotting publishable quality maps using generic mapping tools (GMT)
Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...
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...
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 Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem
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 is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...
Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...
Factor Analysis is an exploratory data analysis method used to search influential underlying factors or latent variables from a set of observed variables
We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...
Using Randomization to test the disprove the null hypothesis
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...
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. ...
Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram
Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB
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...
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 ...
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 ...
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Python code to automatically plot the record section for the highest magnitude earthquake in the given time range
We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib
This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...
This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...
We learn how to make the three-dimensional map using both GMT and PyGMT
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...
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
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...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
This post demonstrate how to use Python to set up clip topographic map based on coastlines.
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...
In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.
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. ...
We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB
Transfer learning using the pre-trained deep learning networks from MATLAB can be easily implemented to achieve fast and impressive results
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...
We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection
We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...
Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom
Using Randomization to test the disprove the null hypothesis
How to fit data to non-linear model
How to fit data to non-linear model
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 is essential in most fields of science, including geophysics, economics, etc. Most of the geophysical data comes in a time-series format...
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 ...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
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.
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Python code to automatically plot the record section for the highest magnitude earthquake in the given time range
Plotting publishable quality maps using generic mapping tools (GMT)
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. ...
We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB
In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...
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 ...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
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.
Python code to automatically plot the record section for the highest magnitude earthquake in the given time range
Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib
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...
We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.
Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB
In geophysics, it is important to understand and identify the complex and unknown relationships between two time-series. Cross-correlation is an established ...
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...
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...
Numerical tests to understand the degrees of freedom in time series
Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...
We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...
Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom
Using Randomization to test the disprove the null hypothesis
Plot the counties of Taiwan using shapefile data on a map created by PyGMT
How you can plot the shear-wave splitting measurements from splitting database using PyGMT.
Uses pandas to read the html page and extract the html data into pandas dataframe
Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.
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...
We learn how to make the three-dimensional map using both GMT and PyGMT
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
We learn how to make the three-dimensional map using both GMT and PyGMT
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...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
Generating GMT style shaded relief map in Python
Plotting publishable quality maps using generic mapping tools (GMT)
Plotting publishable quality maps using generic mapping tools (GMT)
We learn how to make the three-dimensional map using both GMT and PyGMT
This post demonstrate how to use Python to set up clip topographic map based on coastlines.
Generating GMT style shaded relief map in Python
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 ...
We read the station information from a file to plot the 1 arc-minute relief map
Plotting 1 arc-minute topographic map in python
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...
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...
Numerical tests to understand the degrees of freedom in time series
Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...
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 ...
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 Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. 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 is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...
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...
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 Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Shortcut code for quickly logging temperature in Apple health app
Quick action for mac to easily love, dislike, rate songs in apple music app
Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib
Read the seismic traces from the miniseed files and compute the cross-correlation and spectrogram
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are 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...
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...
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...
We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection
Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
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...
Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...
Simple earthquake location problem and its solution using Geiger’s method
Plot the counties of Taiwan using shapefile data on a map created by PyGMT
Uses pandas to read the html page and extract the html data into pandas dataframe
Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.
Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
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...
Plotting publishable quality maps using generic mapping tools (GMT)
Plotting publishable quality maps using generic mapping tools (GMT)
Numerical tests to understand the degrees of freedom in time series
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...
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 ...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.
Plotting 1 arc-minute topographic map in python
Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Quick action for mac to easily love, dislike, rate songs in apple music app
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...
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...
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 ...
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...
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...
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 ...
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...
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...
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 ...
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...
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...
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 ...
An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem
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 is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...
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...
An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem
Least-squares method is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...
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. ...
We will learn how to convert a mseed data file into mat format and then read and analyze it using MATLAB
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 ...
We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib
This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...
This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...
Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...
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...
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...
Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...
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...
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...
Plot the counties of Taiwan using shapefile data on a map created by PyGMT
Uses pandas to read the html page and extract the html data into pandas dataframe
Read the earthquake data from a csv file and overlay on a topographic map using PyGMT.
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 ...
Generators don’t hold the entire result in memory. It yields one result at a time.
This post demonstrate how to use Python to set up clip topographic map based on coastlines.
Plotting publishable quality maps using generic mapping tools (GMT)
Numerical tests to understand the degrees of freedom in time series
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 ...
Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc
Quick action for mac to easily love, dislike, rate songs in apple music app
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...
Python code to automatically plot the record section for the highest magnitude earthquake in the given time range
This post is aimed to resolve the issues regarding the conflicts of using obspy and basemap libraries together.
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are 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...
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.
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...
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.
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...
A simple tutorial on how to plot high resolution topographic map using GMT tools in Python
NetCDF file format has been designed for storing multidimensional scientific data such as temperature, rainfall, humidity, etc. In this post, we will see how...
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...
NetCDF file format has been designed for storing multidimensional scientific data such as temperature, rainfall, humidity, etc. In this post, we will see how...
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...
A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.
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.
A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.
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.
A PyQt5 application for retrieving and visualizing sound waveforms in real time. Codes included.
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.
We learn how to make the three-dimensional map using both GMT and PyGMT
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...
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...
An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!
We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.
An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. 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...
An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. 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
This tutorial gives a brief description of scientific computing using Pandas by introducing Series, DataFrame, Pandas common operations, methods, conditional...
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...
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib
We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.
Runge-Kutta methods are most popular method to solve ordinary differential equations (ODEs) with a better approximation than the Euler method. We compare the...
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...
In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.
In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.
In this tutorial post, I give a quick demo of how to install Python (using anaconda) and then getting started with writing simple scripts.
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 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 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 don’t hold the entire result in memory. It yields one result at a time.
A basic to advanced guide to making interactive plots in Bokeh.
A basic to advanced guide to making interactive plots in Bokeh.
A basic to advanced guide to making interactive plots in Bokeh.
Plotting 1 arc-minute topographic map in python
Plotting 1 arc-minute topographic map in python
How to fit data to non-linear model
How to fit data to non-linear model
How to fit data to non-linear model
How to fit data to non-linear model
How to fit data to non-linear model
We read the station information from a file to plot the 1 arc-minute relief map
We read the station information from a file to plot the 1 arc-minute relief map
We read the station information from a file to plot the 1 arc-minute relief map
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 ...
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 ...
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 ...
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 ...
This post demonstrate how to use Python to set up clip topographic map based on coastlines.
Using Randomization to test the disprove the null hypothesis
Using Randomization to test the disprove the null hypothesis
Visualize the statistics of the data using MATLAB: mean, median, std, interquartile range, skewness, kurtosis, t-statistic, degrees of freedom
We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...
We test for the correlation coefficients or the covariance between two sets of random numbers selected from normal distribution using the Monte Carlo simulat...
Introduction to the concepts of tomography with equations and codes. Introduction to the concepts of overdetermined, underdetermined and mix-determined probl...
Introduction to the concepts of tomography with equations and codes. Introduction to the concepts of overdetermined, underdetermined and mix-determined probl...
Quick action for mac to easily love, dislike, rate songs in apple music app
Quick action for mac to easily love, dislike, rate songs in apple music app
Quick action for mac to easily love, dislike, rate songs in apple music app
Quick action for mac to easily love, dislike, rate songs in apple music app
Quick action for mac to easily love, dislike, rate songs in apple music app
Shortcut code for quickly logging temperature in Apple health app
Shortcut code for quickly logging temperature in Apple health app
Shortcut code for quickly logging temperature in Apple health app
Shortcut code for quickly logging temperature in Apple health app
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Mac can be easily automated by the help of several tools such as automator, quick actions, applescripts
Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc
Some handy tweaks for mac like relocating default screenshot location, renaming batch files etc
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Short demonstration of the ppsd class defined in Obspy using 3 days of data for station PB-B075
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are included.
Short demonstration of how to plot the track or trajectory of a hurricane on a map. Codes are 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.
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 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 ...
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 ...
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 ...
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...
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...
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...
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...
NetCDF file format has been designed for storing multidimensional scientific data such as temperature, rainfall, humidity, etc. In this post, we will see how...
We demonstrate how to model a simple wave, obtain its frequencies, apply Hilbert transform, and perform edge detection
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 ...
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 is essential in most fields of science, including geophysics, economics, etc. Most of the geophysical data comes in a time-series format...
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 is an exploratory data analysis method used to search influential underlying factors or latent variables from a set of observed variables
Simple earthquake location problem and its solution using Geiger’s method
Seismic resolution and fidelity are the two important measures of the quality of the seismic record and the seismic images. Seismic resolution quantifies the...
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...
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...
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...
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...
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...
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...
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...
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...
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
Codes for plotting advanced 2D plots using matplotlib library in Python. Includes simple 2D plot, error bars, bar graphs, histograms, multiple plots, etc
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...
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...
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 is a popular approach in geophysical inversion to estimate the parameters of a postulated Earth model from given observations. This meth...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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...
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...
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...
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...
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...
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...
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...
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...
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...
An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem
An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem
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.
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.
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.
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.
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 or generic mapping tools have become synonymous with plotting maps in Earth, Ocean, and Planetary sciences. It can be used for processing data, generatin...
An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!
An introduction to the wavelet analysis for a real geophysical data set. I compared the analysis to the Fourier analysis. Codes included!
We learn how to make the three-dimensional map using both GMT and PyGMT
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...
This tutorial gives a brief description of scientific computing using numpy by introducing arrays, methods, attributes, random numbers, indexing, broadcastin...
We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib
We learn how to read huge csv file containing time series data by breaking it into chunks and then visualizing it with matplotlib
In this introduction to the concepts of Pytorch data structures, we will learn about how to create and reshape tensors using Pytorch and compare it with the ...
In this introduction to the concepts of Pytorch data structures, we will learn about how to create and reshape tensors using Pytorch and compare it with the ...
Transfer learning using the pre-trained deep learning networks from MATLAB can be easily implemented to achieve fast and impressive results
Transfer learning using the pre-trained deep learning networks from MATLAB can be easily implemented to achieve fast and impressive results
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
We will learn the basics of Fourier analysis and implement it to remove noise from the synthetic and real signals
Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated u...
We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.
We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.
We will learn the basic concepts of wavelet tranform and multi-resolution analysis starting from the Fourier Transform, and Gabor Transform.
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.
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.
Follow the instructions to make your python script executable from anywhere in Linux system.
Is the signal denoising using a wavelet effective? We will implement the wavelet-based signal denoising using the MATLAB
Is the signal denoising using a wavelet effective? We will implement the wavelet-based signal denoising using the MATLAB
Is the signal denoising using a wavelet effective? We will implement the wavelet-based signal denoising using the MATLAB
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...
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...
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...
Uses pandas to read the html page and extract the html data into pandas dataframe
Uses pandas to read the html page and extract the html data into pandas dataframe
How can you link unsplash images directly to your blog without hosting it locally.
How can you link unsplash images directly to your blog without hosting it locally.
How can you link unsplash images directly to your blog without hosting it locally.
How you can plot the shear-wave splitting measurements from splitting database using PyGMT.
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...
Plot the counties of Taiwan using shapefile data on a map created by PyGMT
Plot the counties of Taiwan using shapefile data on a map created by PyGMT
Librosa can efficiently compute the spectrogram for large time series data in seconds. We will use that to plot the spectrogram using matplotlib