opera-adt / disp-s1

OPERA Displacement workflows
Other
5 stars 7 forks source link

Add "average seasonal coherence" as an input to the interface #6

Closed scottstanie closed 6 months ago

scottstanie commented 11 months ago

Charlie has given a nice easy-to-use code example to pull the Josef global coherence dataset over an area of interest

https://github.com/OPERA-Cal-Val/s1-coherence-2020/blob/main/3_Seasonal_Coherence_View.ipynb

Reference date selection

Background

For each tile, seasonal composites of the coherence at different repeat intervals and backscatter imagery were calculated. We calculated the median coherence based on all coherence estimates per tile of a given repeat interval (6, 12, 18, 24, 36, and 48) per three-month period: 1) December, January, February 2) March, April, May 3) June, July, August, and 4) September, October, November. We chose the median operation to account for outliers. In the case of the backscatter intensity products, we calculated per three-month period the average backscatter intensity in VV and VH, or HH and HV, polarization.

The decay of coherence with increasing repeat interval was modelled for each season at pixel-level with the exponential model38 𝛾𝑑(𝑑)=(1βˆ’πœŒβˆž)π‘’βˆ’π‘‘/𝜏+𝜌∞ (3)

where ρ∞ and Ο„ denote the long-term coherence and rate of coherence decay with increasing repeat interval, respectively.