Interferogram formation (InSAR operator)

Computation of interferogram and computation/removal of the flat-earth phase

This operator computes (complex) interferogram, with or without subtraction of the flat-earth (reference) phase. The reference phase is subtracted using a 2d-polynomial that is also estimated in this operator.

If the orbits for interferometric pair are known, the flat-earth phase is estimated using the orbital and metadata information and subtracted from the complex interferogram. The flat-earth phase is the phase present in the interferometric signal due to the curvature of the reference surface. The geometric reference system of the reference surface is defined by the reference system of satellite orbits (for now only WGS84 supported, which the reference system used by all space-borne SAR systems).

The flat-earth phase is computed in a number of points distributed over the total image, after which a 2d-polynomial is estimated (using least squares) fitting these 'observations', (e.g. plane can be fitted by setting the degree to 1.)

A polynomial of degree 5 normally is sufficient to model the reference phase for a full SAR scene (approx 100x100km). While, a lower degree might be selected for smaller images, and higher degree for 'long-swath' scenes. Note that the higher order terms of the flat-earth polynomial are usually small, because the polynomial describes a smooth, long wave body (ellipsoid). To recommended polynomial degree, that should ensure the smooth surface for most image sizes and areas of the world is 5th degree.

In order to reduce the noise, as the post-processing step, you can perform multilooking (with Multilook Operator). Multilooking has to be performed separately on 'virtual' bands phase or intensity. In future releases complex Multilook operator will be released. Note that in case of ESA's ERS and Envisat sensors, the factor 5:1 (azimuth:range) or similar ratio between the factors is chosen to obtain approximately square pixels (20x20 m^2 for factors 5 and 1). Of course the resolution decreases if multilooking is applied.

Note: The reference phase polynomial is not estimated nor subtracted from the computed interferogram, if the flag for NOT performing the reference-phase estimation/subtraction is being marked. Note that if the subtraction of the flat-earth is being skipped the formed interferogram will still have the fringes caused by the earth curvature, and could hamper further interferometric processing and analysis. The intention of computing interferogram without the reference-phase being subtracted is only for demonstration and educational purposes.

Operator parameters:

   The following parameters are used by this operator:
  1. Number of points to compute reference phase for least square estimation. Default value is 501, and sufficient for 100x100km SAR scenes. For smaller/or larger scenes this number can be adapted.
  2. The degree of 2D flat-earth polynomial. Recommended degree, appropriate for most usage cases is 5th degree.
  3. Orbit interpolation method. Defaults to a polynomial of degree (number of state vectors)-1, but smaller than degree 5. Optionally, the degree for the orbit interpolation can be declared. The specified degree has to be smaller or equal to (number of state vectors)-1. The positions of state vectors (x,y,z) are independently interpolated, the velocities are estimated from the position. The annotated velocities, if available, will be used in the interpolation and not computed from the positions. NOTE: It is not recommended to use a degree smaller than the maximum possible, except if it gets too large to avoid oscillations. However, depending on the temporal posting of the state vectors and their accuracy different interpolation strategies can give better results.
  4. Flag for skipping estimation and subtraction of the reference-phase.
  5. Flag for subtracting topographic phase
  6. DEM to be used for subtracting topographic phase
  7. Extension of tile (%) for DEM simulation (optimization parameter)
  8. Flag for outputting elevation (if topographic phase subtraction is chosen)
  9. Flag for outputting orthorectified latitudes and longitudes (if topographic phase subtraction is chosen)
  10. Flag for including coherence estimation
  11. Choice of using square pixel or not (for coherence estimation)
  12. Coherence range and azimuth window sizes

Source bands:

The source bands are the set of, usually coregistered, bands of the complex product