Convert to displacement

Most DAS interrogators provides strain-rate. For sufficiently rectilinear deployments it is more natural to convert strain-rate to displacement. Read this article for more details.

First open some data and convert it to xarray format.

strain_rate = xd.open("sample.nc")
strain_rate.plot(yincrease=False, vmin=-0.5, vmax=0.5);

Then convert strain rate to deformation and then to displacement.

import xdas.signal as xs

strain = xs.integrate(strain_rate, dim="time")
deformation = xs.integrate(strain, dim="distance")
displacement = xs.sliding_mean_removal(deformation, wlen=2000.0, dim="distance")
displacement.plot(yincrease=False, vmin=-0.5, vmax=0.5);