Synchrosqueezing Wavelet Transform explanation?
Synchrosqueezing is a powerful reassignment method. To grasp its mechanisms, we dissect the (continuous) Wavelet Transform, and how its pitfalls can be remedied. Physical and statistical interpretations are provided.
If unfamiliar with CWT, I recommend this tutorial. SSWT is implemented in MATLAB as wsst, and in Python, ssqueezepy. (-- All answer code)

Begin with CWT of a pure tone:

A straight line in the time-frequency (rather, time-scale) plane, for our fixed-frequency sinusoid over all time - fair. ... except is it a straight line? No, it's a band of lines, seemingly centered about some maximum, likely the &quot;true scale&quot;. Zooming,

makes this more pronounced. Let's plot rows within this zoomed band, one by one:

and all superimposed, each for samples 0 to 127 (horizontal zoom):

Notice anything interesting? They all have the same frequency. It isn't particular to this sinusoid, but is how CWT works in correlating wavelets with signals.
It appears much of information &quot;repeats&quot;; there is redundancy. Can we take advantage of this? Well, if we just assume that all these adjacent bands actually stem from one and the same band, then we can merge them into one - and this, in a nutshell, is what synchrosqueezing does. Naturally it's more nuanced, but the underlying idea is that we sum components of the same instantaneous frequency to obtain a sharper, focused time-frequency representation.
Here's that same CWT, synchrosqueezed:

Now that is a straight line.

How's it work, exactly?
We have an idea, but how exactly is this mathematically formulated? Motivated by speaker identification and Empirical Mode Decomposition, SSWT builds upon the modulation model:
$$
f(t) = \sum_{k=1}^{K} A_k(t) \cos(\phi_k (t)), \tag{1}
$$
where $A_k(t)$ is the instantaneous amplitude and
$$
\omega_k(t) = \frac{d}{dt}(\phi_k(t)) \tag{2}
$$
the instantaneous frequency of component $k$, where we seek to find $K$ such &quot;components&quot; that sum to the original signal. More on this below, &quot;MM vs FT&quot;.
At this stage, we only have the CWT, $W_f(a, b)$ (a=scale, b=timeshift); how do we extract $\omega$ from it? Revisit the zoomed pure tone plots; again, the $b$-dependence preserves the original harmonic oscillations at the correct frequency, regardless of $a$. This suggests we compute, for any $(a, b)$, the instantaneous frequency via
$$
\omega(a, b) = -j[W_f(a, b)]^{-1} \frac{\partial}{\partial b}W_f(a, b), \tag{3}
$$
where we've taken the log-derivative, $f' / f$. To see why, we can show that CWT of $f(t)=A_0 \cos (\omega_0 t)$ is:
$$
W_f(a, b) = \frac{A_0}{4 \pi} \sqrt{a} \overline{\hat{\psi}(a \omega_0)} e^{j b \omega_0} \tag{4}
$$
and thus partial-diffing w.r.t. $b$, we extract $\omega_0$, and the rest in (3) gets divided out. (&quot;But what if $f$ is less nice?&quot; - see caveats).
Finally, equipped with $\omega (a, b)$, we transfer the information from the $(a, b)$-plane to a $(\omega, b)$ plane:
$$
\boxed{ S_f (\omega_l, b) = \sum_{a_k\text{ such that } |\omega(a_k, b) - w_l| \leq \Delta \omega / 2} W_f (a_k, b) a_k^{-3/2}} \tag{5}
$$
with $w_l$ spaced apart by $\Delta w$, and $a^{-3/2}$ for normalization (see &quot;Notes&quot;).
And that's about it. Essentially, take our CWT, and reassign it, intelligently.

So where are the &quot;components&quot;? -- Extracted from high-valued (ridge) curves in the SSWT plane; in the pure tone case, it's one line, and $K=1$. More examples; we select a part of the plane and invert over it as many times as needed.

Modulation Model vs Fourier Transform:
What's $(1)$ all about, and why not just use FT? Consider a pendulum oscillating with fixed period and constant damping, and its FT:
$$
s(t) = e^{-t} \cos (25t) u(t)\ \Leftrightarrow\ S(\omega) = \frac{1 + j\omega}{(1 + j\omega)^2 + 625}
$$

What does the Fourier Transform tell us? Infinitely many frequencies, but at least peaking at the pendulum's actual frequency. Is this a sensible physical description? Hardly (only in certain indirect senses); the problem is,  FT uses fixed-amplitude complex sinusoid frequencies as its building blocks (basis functions, or &quot;bases&quot;), whereas here we have a variable amplitude that cannot be easily represented by constant frequencies, so FT is forced to &quot;compensate&quot; with all these additional &quot;frequencies&quot;.
This isn't limited to amplitude modulation; the less sinusoidal or non-periodic the function, the less meaningful its FT spectrum (though not always). Simple example: 1Hz triangle wave, multiple FT frequencies. Frequency-modulation suffers likewise; more intuition here.
These are the pitfalls the Modulation Model aims to address - by decoupling amplitude and frequency over time from the global signal, rather than assuming the same (and constant!) amplitude and frequency for all time.
Meanwhile, SSWT - perfection:


Is synchrosqueezing magic?
We seem to gain a lot by ssqueezing - an apparently perfect frequency resolution, violating Heisenberg's uncertainty, and partial noise cancellation (&quot;Notes&quot;). How can this be?
A prior. We assume $f(t)$ is well-captured by the $A_k(t) \cos(\phi_k (t))$ components, e.g. based on our knowledge of the underlying physical process. In fact we assume much more than that, shown bit later, but the idea is, this works well on a subset of all possible signals:

Indeed, there are many ways synchrosqueezing can go awry, and the more the input obeys SSWT's assumptions (which aren't too restrictive, and many signals naturally comply), the better the results.

What are SSWT's assumptions? (when will it fail?)
This is a topic of its own (which I may post on later), but briefly, the formulation's as follows. Firstly note that we must somehow restrict what $A(t)$ and $\psi(t)$ can be, else, for example, $A(t)$ can simply cancel out the cosine and become any other function. More precisely, the components are to be such that:

More info in ref 2.

How would it be implemented? There's now Python code, clean &amp; commented. Regardless, worth noting:

For very small CWT coefficients, phase is unstable (just like for DFT), which we work around by zeroing all such coefficients below a given threshold.
For any frequency row/bin $w_l$ in SSWT plane, we reassign from $W_f(a, b)$ based on what's closest to $w_l$ according to $\omega (a, b)$, and for log-scaled CWT we use log-distance.


Summary:
SSWT is a time-frequency analysis tool. CWT extracts the time-frequency information, and synchrosqueezing intelligently reassigns it - providing a sparser, sharper, noise-robust, and partly denoised representation. The success of synchrosqueezing is based in and explanied by its prior; the more the input obeys assumptions, the better the results.

Notes &amp; caveats:

What if $f$ isn't nice in $\omega(a,  b)$ example? Valid question; in practice, the more the function satisfies aforementioned assumptions, the less of a problem this is, as authors demonstrate through various lemmas.
In the SSWT of damped pendulum, I cheated a little by extending signal's time to $(-2, 6)$; this is only to prevent boundary effects, which is a CWT phenomenon that can be remedied; here's directly 0 to 6.
Partial noise cancellation? Indeed; see pg 536 of ref 1.
What's the $a^{-3/2}$ in $(5)$? Synchrosqueezing effectively inverts $W_f$ onto the reassigned plane, using one-integral iCWT.
&quot;Fourier bad?&quot; My earlier comparison is prone to criticism. To be  clear, FT is the most solid and general-purpose basis that we have for a signals framework. But it's not an all-purpose-best; depending on context, other constructions are more meaningful and more useful.


Where to learn more?
The refernced papers are a good source, so are MATLAB's wsst and cwt docs and ssqueezepy's source code. I may also write further Q&amp;A's, which you can be notified of by subbing this thread.

References:

A Nonlinear Squeezing of the CWT Based on Auditory Nerve Models - I. Daubechies, S. Maes. Excellent origin paper with succinct intuitions.
Synchrosqueezed Wavelet Transforms: a tool for Empirical Mode Decomposition - I. Daubechies, J. Lu, H.T. Wu. Good followup paper with examples.
The Synchrosqueezing algorithm for time-varying spectral analysis: robustness properties and new paleoclimate applications - G. Thakur, E. Brevdo, et al. Further exploration of robustness properties and implementation details (including threshold-setting).