## LEAST-SQUARES METHOD IN GEOSCIENCES

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

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

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

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

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

Numerical tests to understand the degrees of freedom in time series

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

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

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

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...

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 ...

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

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

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

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

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 fit data to non-linear model

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

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

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...

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

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

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

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

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

Generating GMT style shaded relief map in Python

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

Plotting 1 arc-minute topographic map in python

Plotting publishable quality maps using generic mapping tools (GMT)

Plotting publishable quality maps using generic mapping tools (GMT)

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

With the increasing amount of data, parallel computing is quickly becoming a necessity. Modern computers comes with more than one process and we most often o...

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

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

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

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

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

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

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

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

Shortcut code for quickly logging temperature in Apple health app

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

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

Using Randomization to test the disprove the null hypothesis

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 ...

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

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

Introduction to the time and frequency analysis in MATLAB

How to plot the track or trajectory of a hurricane on a map?