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.
Posts by Tag
Tags in Chronological Order
- python 52
- obspyTutorial 18
- numericalMethodsScientificComputation 17
- geophysics 16
- geospatial data visualization 14
- obspy 12
- introEarthDataAnalysis 11
- matlab 11
- visualization 9
- seismology 9
- pygmt 9
- time-frequency analysis 8
- statistics 7
- GMT 6
- geospatial data visualization python 6
- geospatial data 6
- techniques 5
- apple 5
- time series analysis 5
- signal processing 5
- matplotlib 4
- time-series 4
- MATLAB 4
- optimization method 4
- iterative methods 4
- plotting map 4
- gmt mapping 3
- data science 3
- time series 3
- topography 3
- mac 3
- MS word 3
- microsoft word 3
- writing manuscript 3
- equations 3
- parallel computing 3
- earthquake location problem 3
- mseed 3
- scientific computing 3
- initial value problem 3
- c 3
- maps 2
- degrees of freedom 2
- quick action 2
- macbook 2
- record section 2
- basemap 2
- geospatial data visualization python python 2
- epicentral distance 2
- waveforms 2
- Generic Mapping Tools 2
- netcdf 2
- data analysis 2
- threading 2
- least-squares method 2
- desktopApp 2
- PyQt5 2
- pyqtgraph 2
- GMT6 2
- wavelets 2
- fourier transform 2
- indian monsoon rainfall data 2
- pandas dataframe 2
- denoising 2
- spectrogram 2
- Runge Kutta method 2
- compile 2
- programs 2
- fortran 2
- fortran 90 2
- data structures 2
- cpp 2
- C++ 2
- python3 1
- anaconda 1
- miniconda 1
- git 1
- github 1
- clone 1
- merge 1
- branch 1
- github desktop 1
- tutorial 1
- generators 1
- bokeh 1
- streaming data 1
- web data visualization 1
- noaa 1
- global relief 1
- non-linear 1
- curve-fitting 1
- optimization 1
- lsqcurvefit 1
- non-linear regression 1
- station map 1
- topographic 1
- etopo1 1
- mapping 1
- kriging 1
- pykrige 1
- Ordinary Kriging 1
- clipped 1
- hypothesis testing 1
- null hypothesis 1
- histograms 1
- monte-carlo simulations 1
- correlation 1
- tomography 1
- seismic tomography 1
- apple music 1
- Rating song 1
- love song 1
- dislike song 1
- ios 1
- log data 1
- iphone 1
- health 1
- automation 1
- automator 1
- services 1
- quick actions 1
- applescripts 1
- screenshot 1
- batch renaming 1
- psd 1
- probabilistic power spectral density 1
- temporal plot 1
- spectrogram ppsd 1
- hurricane track 1
- typhoon track 1
- cyclone track 1
- P-arrival 1
- S-arrival 1
- Ulysses 1
- Papers 1
- Mendeley 1
- hierarchical data analysis 1
- xarray 1
- lowpass filter 1
- smoothing 1
- writing data 1
- wave modeling 1
- hilbert transform 1
- downloading-data 1
- analysis 1
- loading 1
- factor analysis 1
- factor analysis python example 1
- earthquake location 1
- seismic resolution 1
- cross-correlations 1
- monte-carlo 1
- spectral analysis 1
- Fourier Transform 1
- Welchs method 1
- Kaiser window 1
- plomb 1
- spectral analysis matlab 1
- 2D plots 1
- bar plot 1
- log scale 1
- multiple plots 1
- histogram 1
- multiprocessing 1
- joblib 1
- concurrent futures 1
- ascii 1
- text 1
- flask 1
- webapp 1
- deploy 1
- heroku 1
- slrm 1
- seismograms 1
- SAC 1
- bandpass filter 1
- instaseis 1
- global optimization method 1
- eigen functions 1
- matlab engine api 1
- matlab compiler 1
- matlab runtime 1
- Genetic Algorithm 1
- coffee cup size problem 1
- looping 1
- while loop 1
- for loop 1
- speedup 1
- subplots in GMT 1
- GMT basics 1
- Wavelet transform 1
- Sea surface temperature data 1
- three-dimensional map 1
- sktime 1
- forecasting 1
- NaiveForecaster 1
- KNeighborsRegressor 1
- ExponentialSmoothing 1
- ARIMA 1
- arrays 1
- methods and attributes 1
- big data 1
- read_csv 1
- PyTorch 1
- Tensors 1
- Alexnet 1
- transfer learning 1
- fft 1
- filter 1
- empirical orthogonal function analysis 1
- space-time field 1
- signals processing 1
- wavelet transform 1
- Gabor transform 1
- scp in python 1
- ssh in python 1
- executable 1
- http.server. webserver 1
- hosting blogs 1
- hosting server locally 1
- sharing files 1
- inversion method 1
- get stations 1
- plotting stations map 1
- pandas 1
- web extract 1
- unsplash 1
- unsplash api 1
- jekyll blog 1
- SKS splitting 1
- plotting polarization directions 1
- Euler method 1
- shapefiles 1
- geopandas 1
- Taiwan 1
- North America 1
- librosa 1
- Shooting method 1
- boundary value problem 1
- regression analysis 1
- linux 1
- arm processor 1
- jekyll 1
- M1 1
- install 1
- amazon 1
- text to speech 1
- ai 1
- aws polly 1
- audio 1
- read 1
- write 1
- convolution 1
- dask 1
- data 1
- earthquake catalog 1
- Ray 1
- distributed computing 1
- data encryption 1
- decryption 1
- security 1
- 3D 1
- publish package 1
- pip 1
- poetry 1
- irregular data 1
- Lomb–Scargle periodogram 1
- interpolate 1
- tabular data 1
- earth model 1
- seismic earth model 1
- read text file in fortran 1
- Fortran 90 1
- Fortran 95 1
- SEMUCB 1
- gmt 1
- tomographic map 1
- focal mechanism 1
- topographic map 1
- plot earthquakes 1
- latex 1
- slides 1
- beamer 1
- pylatex 1
- L-BFGS 1
- gradient descent method 1
- multithreading in C 1
- concurrency 1
- bash 1
- shell 1
- yml 1
- yaml 1
- loading libraries 1
- shared library 1
- dynamic library 1
- getopt_long 1
Posts with Each Tags
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 ...
Using dask python library to read a huge global earthquake catalog file
We will see how we can read a large data file with earthquake catalog in Python using the Dask library
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)
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
Getting started with obspy: downloading waveform data (codes included)
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.
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.
Getting started with obspy: downloading waveform data (codes included)
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
numericalMethodsScientificComputation
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
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.
Getting started with obspy: downloading waveform data (codes included)
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
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...
Getting started with obspy: downloading waveform data (codes included)
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
Getting started with obspy: downloading waveform data (codes included)
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
pygmt
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
time-frequency analysis
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
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
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
geospatial data
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...
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
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...
signal processing
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
matplotlib
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.
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...
plotting map
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.
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
MS word
Writing and formatting a scientific manuscript in microsoft word
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
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
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
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 ...
parallel computing
Threading in c/c++ programming language using pthreads
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...
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 ...
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...
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.
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
basemap
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.
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...
threading
Threading in c/c++ programming language using pthreads
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)
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)
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)
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.
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!
fourier transform
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!
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
spectrogram
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.
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.
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.
Threading in c/c++ programming language using pthreads
Using multiple threads in C for concurrent process flow
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...
clone
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
downloading-data
Getting started with obspy: downloading waveform data (codes included)
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...
loading
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...
spectral analysis
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...
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 ...
flask
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...
bandpass filter
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...
instaseis
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...
global optimization method
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...
eigen functions
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...
matlab engine api
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...
matlab compiler
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...
matlab runtime
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...
Genetic Algorithm
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
coffee cup size 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
looping
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.
while loop
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.
for loop
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.
speedup
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.
subplots in GMT
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 basics
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...
Wavelet 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!
Sea surface temperature data
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!
three-dimensional map
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
sktime
A quick look into the sktime for time-series forecasting (codes included)
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
forecasting
A quick look into the sktime for time-series forecasting (codes included)
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
NaiveForecaster
A quick look into the sktime for time-series forecasting (codes included)
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
KNeighborsRegressor
A quick look into the sktime for time-series forecasting (codes included)
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
ExponentialSmoothing
A quick look into the sktime for time-series forecasting (codes included)
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
ARIMA
A quick look into the sktime for time-series forecasting (codes included)
I used the sktime library to forecast the airline data using NaiveForecaster, KNeighborsRegressor, Statistical forecasters, and auto ARIMA model.
arrays
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...
methods and attributes
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...
big data
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
read_csv
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
PyTorch
Getting familiar with pytorch data structures (codes included)
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 ...
Tensors
Getting familiar with pytorch data structures (codes included)
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 ...
Alexnet
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
transfer learning
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
fft
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
filter
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
empirical orthogonal function analysis
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...
space-time field
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...
signals processing
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 transform
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.
Gabor transform
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.
scp in python
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.
ssh in python
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.
executable
Make your python script executable from anywhere in linux
Follow the instructions to make your python script executable from anywhere in Linux system.
http.server. webserver
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
hosting blogs
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
hosting server locally
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
inversion method
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...
get stations
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...
plotting stations map
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...
pandas
How to extract station data from the web page
Uses pandas to read the html page and extract the html data into pandas dataframe
web extract
How to extract station data from the web page
Uses pandas to read the html page and extract the html data into pandas dataframe
unsplash
How to link unsplash images directly to your jekyll blog
How can you link unsplash images directly to your blog without hosting it locally.
unsplash api
How to link unsplash images directly to your jekyll blog
How can you link unsplash images directly to your blog without hosting it locally.
jekyll blog
How to link unsplash images directly to your jekyll blog
How can you link unsplash images directly to your blog without hosting it locally.
SKS splitting
How to plot shear-wave splitting measurements using pygmt
How you can plot the shear-wave splitting measurements from splitting database using PyGMT.
plotting polarization directions
How to plot shear-wave splitting measurements using pygmt
How you can plot the shear-wave splitting measurements from splitting database using PyGMT.
Euler method
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...
shapefiles
How to overlay shapefile data on pygmt maps
We learn how to plot selected shapefile data using geopandas on top of PyGMT maps
geopandas
How to overlay shapefile data on pygmt maps
We learn how to plot selected shapefile data using geopandas on top of PyGMT maps
Taiwan
How to overlay shapefile data on pygmt maps
We learn how to plot selected shapefile data using geopandas on top of PyGMT maps
North America
How to overlay shapefile data on pygmt maps
We learn how to plot selected shapefile data using geopandas on top of PyGMT maps
librosa
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
Shooting method
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 ...
boundary value problem
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 ...
regression analysis
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...
linux
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.
arm processor
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 ...
jekyll
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 ...
M1
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 ...
install
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 ...
amazon
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
text to speech
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
ai
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
aws polly
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
audio
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
read
Read and write text files in modern fortran
We will read and write text and numeric data into a file using modern fortran.
write
Read and write text files in modern fortran
We will read and write text and numeric data into a file using modern fortran.
convolution
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...
dask
Using dask python library to read a huge global earthquake catalog file
We will see how we can read a large data file with earthquake catalog in Python using the Dask library
data
Using dask python library to read a huge global earthquake catalog file
We will see how we can read a large data file with earthquake catalog in Python using the Dask library
earthquake catalog
Using dask python library to read a huge global earthquake catalog file
We will see how we can read a large data file with earthquake catalog in Python using the Dask library
Ray
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 ...
distributed computing
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 ...
data encryption
Encrypt your data and software with python
Learn to encrypt and decrypt any files, data or software with python
decryption
Encrypt your data and software with python
Learn to encrypt and decrypt any files, data or software with python
security
Encrypt your data and software with python
Learn to encrypt and decrypt any files, data or software with python
3D
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.
publish package
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
pip
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
poetry
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
irregular data
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...
Lomb–Scargle periodogram
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...
interpolate
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...
tabular data
Back to top ↑earth model
Back to top ↑seismic earth model
Back to top ↑read text file in fortran
Back to top ↑Fortran 90
Back to top ↑Fortran 95
Back to top ↑SEMUCB
Back to top ↑gmt
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
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
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...
topographic map
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...
plot earthquakes
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...
latex
Create beamer slides using python
We will use the Python package to create beamer presentation and append existing figures to each slide
slides
Create beamer slides using python
We will use the Python package to create beamer presentation and append existing figures to each slide
beamer
Create beamer slides using python
We will use the Python package to create beamer presentation and append existing figures to each slide
pylatex
Create beamer slides using python
We will use the Python package to create beamer presentation and append existing figures to each slide
L-BFGS
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
gradient descent method
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
multithreading in C
Threading in c/c++ programming language using pthreads
Using multiple threads in C for concurrent process flow
concurrency
Threading in c/c++ programming language using pthreads
Using multiple threads in C for concurrent process flow
bash
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.
shell
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.
yml
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.
yaml
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.
loading libraries
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...
dynamic library
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...