Dear SNAP developers (and users),
For some reason, I need to compute the Doppler Centroid on my image. Original image contains metadata to compute it. The thing is, it is a bit boring because many steps are involved.
Summary
The Doppler Frequency in TOPSAR SLC images is given by
The Doppler centroid rate k_t is a function of the Doppler rate k_a and the Doppler centroid rate in raw data k_{rot}. The Doppler centroid rate is given by :
And k_{rot} is a function of the satellite velocity v_s, the wavelength \lambda and the steering rate of the sensor \omega_r :
v_s is computed from orbits. Here we took the central velocity estimate are assumed it as valid for the whole scene. If needed, we can draw a polynomial through all velocity estimates and interpolate the velocity at the the very precise mid-burst time.
- v_ s : â€¦/product/generalAnnotation/orbitList/velocity[x y z]
Which gives us the absolute velocity of 7566.93 m/s.
- f_c : â€¦/product/generalAnnotation/productInformation/radarFrequency
\longrightarrow 5 405 000 459 Hz.
- \omega_r : â€¦/product/generalAnnotation/productInformation/azimuthSteeringRate
\longrightarrow 0.0349 rad/sec.
- c : 299 792 458 m/s
Putting all together in equation \ref{eq:carotte2}, we find k_{rot} equal to 9522.48 Hz./sec
This allows to compute the Doppler rate k_a. It is range variant in TOPSAR acquisition. It is computed from a polynom given in the metadata.
Where i denotes the range dimension (in pixel) and \tau (i) is its related slant-range acquisition time.
- c_i : â€¦/product/generalAnnotation/azimuthFmRateList/azimuthFmRate/ \ azimutFmRatePolynomial[c_0 c_1 c_2]
Note that the value of these coefficients doesnâ€™t change a lot between bursts (relative difference of about \approx e-5).
- \tau_0 : â€¦/product/imageAnnotation/imageInformation/slantRangeTime
\longrightarrow 0.00505182323 sec.
- \tau(i) : given in the Tie-Point Grids in SNAP
Putting all together, k_a is a range varying variable going from -2455.28 Hz/sec (near range) to -2318.39 Hz/sec (far range). This is coherent with Nuno Mirandaâ€™s results (Sentinel-1 Data Quality Manager).
Then, we can compute the Doppler centroid rate in the focused SLC TOPSAR data is given by the formula
The Doppler centroid rate in the focused SLC TOPSAR data is given by the formula
Since k_a is range varying and k_{rot} is almost constant, the Doppler centroid rate is also range varying ang here goes from 1953 Hz/sec (near range) to 1861 Hz/sec (far range). This is coherent with Nuno Mirandaâ€™s results. The doppler Frequency is simply computed from first equation.
However, as soon as a product is create by SNAP, the .data/.dim â€ścomboâ€ť, Abstracted metada are created, from which we can find some parameters not present in the original metadata. More specificaly, the Doppler Centroid Coefficients are directly present for each burst as a second degree polynomial.
How should I use them to directly get the f_{DC} ?
Thanks in advance, it would greatly help me gain some time.
Quentin