Skip to content
3D differencing in OpenTopography - Author: Chelsea Scott
Branch: master
Clone or download

Latest commit

Fetching latest commit…
Cannot retrieve the latest commit at this time.


Type Name Latest commit message Commit time
Failed to load latest commit information.
3D_differencing_python Update Mar 13, 2020

Measuring change at the Earth’s surface: On-Demand vertical and 3D topographic differencing implemented in OpenTopography

Chelsea Scott: author)

Minh Phan, Viswanath Nandigam, Christopher Crosby, Ramon Arrowsmith

A windowed implementation of the iterative closest point (ICP) algorithm is used to calculate displacement and rotation fields from topographic point clouds that span a geologic event of interest. This technique resolves surface deformation along and adjacent to active faults where other geodetic datasets commonly lack spatial resolution.

We provide two sets of codes to perform 3D ICP differencing. This option uses Matlab scripts. The Matlab script is relatively easy to set-up, although it will take a while to complete. This option uses Python and c++ functions to perform ICP differencing. The set-up is more involved than for the Matlab option, but the scripts will run much faster.

Pre- and post- earthquake topographic datasets for the M7 Kumamoto, Japan, earthquake are available here.

Pre earthquake topography:

Post earthquake topography:

More information about ICP differencing applied to earthquakes is available in these publications: Scott, C. P., Arrowsmith, J. R., Nissen, E., Lajoie, L., Maruyama, T., & Chiba, T. (2018). The M7 2016 Kumamoto, Japan, Earthquake: 3-D Deformation Along the Fault and Within the Damage Zone Constrained From Differential Lidar Topography. Journal of Geophysical Research: Solid Earth.

Scott, C., Champenois, J., Klinger, Y., Nissen, E., Maruyama, T., Chiba, T., & Arrowsmith, R. (2019). 2016 M7 Kumamoto, Japan, Earthquake Slip Field Derived From a Joint Inversion of Differential Lidar Topography, Optical Correlation, and InSAR Surface Displacements. Geophysical Research Letters.

You can’t perform that action at this time.