ICP_Topo_Diff_Matlab

Matlab code and earthquake topography datasets to experiment with ICP 3D differencing

Contents: Matlab code that performs windowed ICP and plots the results. There are several Matlab files listed in the included Matlab script that must be downloaded from MathWorks.

To perform the ICP Differencing:

Download the Matlab script differencing_matlab.m and open in Matlab.
Download these functions from Mathworks:

Matlab ICP File Exchange (Jacob Wilm): %https://www.mathworks.com/matlabcentral/fileexchange/27804-iterative-closest-point

Lasdata File Exchange (Teemu Kumpumäki): %https://www.mathworks.com/matlabcentral/fileexchange/48073-lasdata

Place the downloaded files in the same directory as differencing_matlab.m.
Download already processed datasets for the Kumamoto earthquake here: https://cloud.sdsc.edu/v1/AUTH_opentopography/www/shortcourses%2FA2HRT_RCN%2FKumamoto_Lidar.zip or process your own lidar data given the links below.
Put the downloaded lidar datasets in the same directory as the Matlab scripts. The Matlab ICP differencing script (differencing_matlab.m) assumes that the pre-earthquake file is called 'pre.las' and the post-earthquake file is called 'post.las'.
Run differencing_matlab.m on Matlab. It will probably take 10-30 minutes depending your computer's processing speeds. When it completes, Matlab will display a vector field map indicating the 3D displacements for the 2016 M7 Kumamoto, Japan.
To experiment further with ICP, change the sz parameter (Line 40). This is the window size/ resolution used in the ICP calculation. Also change the grd parameter (Line 41). This is spacing between the individual ICP measurements.