Topographic Differencing

The 16 April 2016 Mw 7.0 Kumamoto earthquake:

The 16 April 2016 Mw 7.0 Kumamoto, Japan, earthquake, which ruptured ~40 km of the Futagawa-Hinagu Fault Zone on Kyushu Island with an oblique strike-slip mechanism and surface offset exceeding 2 m. We use a windowed implementation of the Iterative Closest Point algorithm (ICP; Besl & McKay, 1992; Chen &Medioni, 1992) to calculate the 3-D surface deformation produced by the earthquake.

Matlab Code for performing 3D Iterative Closest Point Differencing

This code requires the following codes from the Matlab File Exchange:

• https://www.mathworks.com/matlabcentral/fileexchange/48073-lasdata
• https://www.mathworks.com/matlabcentral/fileexchange/27804-iterative-closest-point

Example pre-earthquake las topography file for topographic differencing from the 2016 M7 Kumamoto, Japan earthquake

The full dataset is available at:

http://opentopo.sdsc.edu/lidarDataset?opentopoID=OTLAS.052018.2444.2&minX=130.736&minY=32.669&maxX=131.033&maxY=32.881

Post-earthquake las topography

http://opentopo.sdsc.edu/lidarDataset?opentopoID=OTLAS.052018.2444.1&minX=130.736&minY=32.669&maxX=131.033&maxY=32.881