compare_and_find_gaps

HSTB.kluster.xarray_helpers.compare_and_find_gaps(source_dat, new_dat, max_gap_length=1.0, dimname='time')

So far, mostly used with Applanix POSPac SBETs. Converted SBET would be the new_dat and the existing navigation in Kluster would be the source_dat. You’d be interested to know if there were gaps in the sbet greater than a certain length that did not coincide with existing gaps related to stopping/starting logging or something. Here we find gaps in the new_dat of size greater than max_gap_length and trim them to the gaps found in source_dat.

Parameters

  • source_dat (Union[DataArray, Dataset]) – xarray DataArray/Dataset, object with dimname as coord that you want to use as the basis for comparison

  • new_dat (Union[DataArray, Dataset]) – xarray DataArray/Dataset that you want to find the gaps in

  • max_gap_length (float) – maximum acceptable gap

  • dimname (str) – name of the dimension you want to find the gaps in

Returns

numpy array, nx2 where n is the number of gaps found

Return type

np.array