Gts.find_offsets_ivel

class pyacs.gts.Gts.Gts(code=None, lat=None, lon=None, h=None, X0=None, Y0=None, Z0=None, t0=None, data=None, data_xyz=None, data_corr_neu=None, data_corr_xyz=None, offsets_dates=[], offsets_values=None, outliers=[], annual=None, semi_annual=None, velocity=None, ifile=None, log=None, metadata=None)[source]

find_offsets_ivel(ivelts=None, lcomponent='EN', threshold_offset=100, offsets_file=None)

Find offsets using l1trend filtering and instantaneous velocity

parameters:

ivelts: Sgts: ivel time series
lcomponent: string: component to use for the analysis
threshold_offset: float: threshold for the offset in mm/yr
offsets_file: string: name of the file to append the offsets. If None, no file is written

notes:

The function uses the l1trend of the time series. It then detects large ivel values and check whether they last for more than one day. If so, an offset is estimated by integrating ivel over two successive dates and offset is applied to the time series. Format of offsets_file is: code date (datetime.isoformat()) lon lat offset_e offset_n offset_u std_offset_e std_offset_n std_offset_u comment