PSI Precision, accuracy and validation aspects

Transcription

1 PSI Precision, accuracy and validation aspects Urs Wegmüller Gamma Remote Sensing AG, Gümligen, Switzerland,

2 Contents - Precision - Accuracy - Systematic errors - Atmospheric effects - Phase unwrapping - Reference point - Time series potential and limitations

3 Precision and accuracy Precision The precision relates to the separation of values that can be distinguished in a measurement Accuracy The accuracy of a measurements relates to how well the measurement corresponds to the truth

4 Precision in DINSAR and PSI context Precision of SAR phase measurement is a very small fraction of a phase cycle precision < 0.1mm Interferometric phase difference partial decorrelation of signal reduced coherence (0.9) higher phase noise reduced precision (e.g. 0.3 radian precision 1.3mm); this value can be improved by spatial averaging Considering the atmospheric path delay effects as an uncertainty of the measurement results in a phase standard deviations of the order of 1.5 radian precision 7mm A PSI result is based on many observations. So the estimation of an average deformation rate based on 50 observations over 5 years results in a precision < 1mm/year

5 Precision in DINSAR and PSI context (cont.) This precision of the linear rate estimation can be calculated with or without subtracting atmospheric phase screens typical precisions as small as 0.1mm/year if atmospheric phase screens are subtracted and precisions around 0.4mm/year if atmospheric phase screens are included in the statistical uncertainty estimation Precision of linear rates at outer boundary of a result and next to spatial gaps is lower. The precision of a specific observation in the time series is again different. There the conservative assumption is that it corresponds to the uncertainty of the single observation (precision of 7mm). Assuming that some of the atmosphere can be modeled results in lower values (e.g. 3mm). In any case it is significantly larger than just the coherence related signal noise unless we are sure that we can reliably model the atmospheric delay. Precision of specific observation is at beginning and end of series (and next to longer gaps in temporal coverage) is lower.

6 Accuracy in PSI context The accuracy of a measurements relates to how well the measurement corresponds to the truth. The precision of the measurement is one factor that relates to the accuracy (but often not the main one). Other relevant factors potentially affecting the accuracy of the deformation rate include: - phase unwrapping errors (see below) - overall tilts of the result (see below) - errors related to subtracted atmospheric phase screens (see below)

7 Estimation of accuracy of a PSI result Based on the data itself the precision can be estimated statistically. Based on the result the temporal and spatial consistence can be checked (e.g. to detect and discard outliers). During the TF project special processing was conducted to crosscompare results obtained by the different OSP and to compare the results against in-situ observations (see contributions / reports by Michele Crossetto). Validation with in-situ measurements can be done. In practice this is often not trivial due to non-identical locations and observation times.

8 Systematic errors (1) Large scale tilts There are causes to have large scale phase trends (e.g. phase ramps) across the image (related to SAR processing, slight errors in baseline model, ). A typical approach to scope with these is to assume that the terrain is not tilting at large scale. If that is a correct assumption then such tilts are correctly removed. Nevertheless, if the area shows a real large scale tilting (e.g. significant subsidence near the left border of the processed section) then the processing possibly removes this trend. This problem is most critical in the case of tectonic movements which are at scales larger than the size of the typical area processed. Be aware of this possibility; consider this in selection of area to be processed (ideally stable part around moving area)

9 Systematic errors (2) Underestimation of deformation rates Some people postulate a general underestimation of deformation rates in PSI results. There are clear limitations to the applicability of PSI in the case of fast deformations. Slower moving targets may be more easily identified as PS than moving targets, especially than non-uniformly moving targets. In the validation experiments we did not really identify such a general underestimation.

10 Atmospheric component (1) Scope and Approach Atmospheric path delay heterogeneity causes a phase component not related to the deformation. This atmospheric phase is spatially quite smooth. A typical approach is to estimate the uniform movement, identify then for each pair the residual phase not modeled by the uniform motion model, to spatially filter this and to subtract it as atmospheric phase. This approach permits to significantly reduce the noise on the time series considered. A limitation of this approach is that non-uniform motion at intermediate and large spatial scale is may also be removed from the result.

11 Atmospheric component (2) Height dependent atmospheric component The atmospheric path delay has a term relating to the dry atmosphere and a part relating to the water vapor. Depending on the topographic height of the terrain there is more or less atmosphere and water vapor above a site height dependent atmospheric component Simple linear model can be used to approximate height dependent atmospheric component modeling and compensation Possible errors: - height dependent atmospheric component may be interpreted as deformation -deformation that correlates with terrain height may be reduced by correction of height dependent atmospheric component

12 Phase unwrapping (1) Phase unwrapping is the most challenging step in SAR interferometry as well as in PSI Limits applicability Possible source for errors An unwrapping (ambiguity) error results in an offset of 2.8 cm (at C- band). Spatial or temporal steps in the deformation values (not rates) of 2.8 cm may be cause by an unwrapping error. For spatially and temporally smooth motion patterns and low deformation values between observations there is a high spatiotemporal consistency of the result high reliability of result, robustness of processing

13 Phase unwrapping (2) Higher probability for unwrapping errors in case of: - Faster deformation / higher deformation gradients - Coarse spatial coverage (areas with few points) - Non-uniform motion - Coarse temporal sampling (e.g. gap between 2001 and 2003) - Noisy data Characteristics of unwrapping error: - Single point, single date - Single point dates after a specific date - Entire section (=several to many points, e.g. patch or isolated region with points)

14 Reference point PSI result is relative to a spatial reference point Reference point can be changed - to transform average deformation rate result subtract reference point rate from all other values - subtract reference point values from all other point s time series

15 Time series potential and limitations PSI time series potential? What is included in a PSI time series? - terms included : linear model phase + noise + non-linear model What is not included in a PSI time series? - terms excluded : - subtracted topographic phase - atmospheric phase - orbital phase - residual linear or quadratic phase trends ( see also tilts) PSI time series provides some quality information (on reliability of average rate estimate / detection of potential unwrapping errors) PSI time series provides some information on non-uniform motion

16 Non-uniform motion potential and limitations Potential of PSI with respect to non-uniform movements - Very localized, low amplitude ( good potential with PSI) - Intermediate to large scale, low amplitude (separation from atmospheric phase is critical as overlapping spatial scale) - Intermediate to large scale, high amplitude (PSI processing very challenging, standard processing may fail gap in result) - Sensors with good spatial resolution and good temporal coverage improve potential of PSI for non-uniform motion monitoring