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.
Utpal Kumar
Geophysicist | Geodesist | Seismologist | Open-source Developer
I am a geophysicist with a background in computational geophysics, currently working as a postdoctoral researcher at UC Berkeley. My research focuses on seismic data analysis, structural health monitoring, and understanding deep Earth structures. I have had the opportunity to work on diverse projects, from investigating building characteristics using smartphone data to developing 3D models of the Earth's mantle beneath the Yellowstone hotspot.
In addition to my research, I have experience in cloud computing, high-performance computing, and single-board computers, which I have applied in various projects. This includes working with platforms like AWS, Docker, and Kubernetes, as well as supercomputing environments such as STAMPEDE2, ANVIL, Savio and PERLMUTTER (and CORI). My work involves developing innovative solutions for structural health monitoring and advancing real-time seismic response analysis. I am committed to applying these skills to further research in computational seismology and structural health monitoring.
Empirical Orthogonal Function analysis to inspect the spatial coherency in the geospatial data
Empirical Orthogonal Functions analysis decomposes the continuous space-time field into a set of orthogonal spatial patterns along with a set of associated uncorrelated time series or principal components. Introductory concepts of EOF analysis
Signal denoising using Fourier Analysis in Python
Fourier analysis is based on the idea that any time series can be decomposed into a sum of integral of harmonic waves of different frequencies. Hence, theoretically, we can employ a number of harmonic waves to generate any signal. The…
Wavelet Analysis Applied to Real Geophysical Datasets: A Quick and Easy Approach
Explore the application of wavelet analysis on real-world geophysical datasets, like ENSO sea surface temperature and Indian monsoon rainfall. This post provides a comparison to Fourier analysis and includes Python code for easy replication.
Genetic Algorithm: a highly robust inversion scheme for geophysical applications
An introduction to the basics of genetic algorithm along with a simple numerical example and solution of an earthquake location problem
Speeding Up Your Code with Parallel Computing in Python
Parallel computing is essential for handling large datasets efficiently. In this post, we explore Python's threading, multiprocessing, and joblib libraries to speed up code execution. Learn the differences between threading and multiprocessing, and understand how to use joblib for optimized parallel processing, especially with NumPy arrays.
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
Automatically Plotting a Record Section for an Earthquake in a Given Time Range Using Python
This tutorial demonstrates how to automatically generate a seismic record section for the largest earthquake within a specified time range and geographic area using Python. The process includes downloading seismic data, applying filters, and visualizing the waveforms for interpretation.
Topographic map clipped by coastlines in Python
Learn how to create visually compelling relief maps in Python by clipping topographic data to specific geographic boundaries, such as coastlines. This tutorial walks you through the process of using matplotlib and Basemap to focus your visualizations on land features within the Taiwan coastline, or alternatively, highlight bathymetric features by masking the land.
Plotting the geospatial data clipped by coastlines in Python
In this post, we explore how to perform geospatial interpolation using the Ordinary Kriging method and visualize the results within the coastline borders of Taiwan. Learn how to efficiently handle irregularly distributed geospatial data and automate the process of clipping data outside geographic boundaries for more accurate and meaningful visualizations.