Sometimes, simulation of polarimetric synthetic aperture radar (PolSAR) data is very useful for theory or algorithm verification. In this blog, I provide one python implementation of the simulation algorithm, which is first proposed by JongSon Lee [1]. You can refer to [1] for detail.

The procedure is simply reviewed as fellows:
1. Compute \(C^{1/2}\) for a given covariance \(C\), where \(C = C^{1/2}(C^{1/2})^{*T}\)
2. Simulate a complex random vector \(v\) that is complex normal distributed.
3. Compute the complex single-look vector by \(u = C^{1/2}v\).
4. Compute the n-look covariance matrix by \(C_{n} = \frac{1}{n}\sum_{1}^{n}uu^{*T}\).
Sometimes, simulation of polarimetric synthetic aperture radar (PolSAR) data is very useful for theory or algorithm verification. In this blog, I provide one python implimentation of the simulation algorithm, which is first proposed by JongSon Lee [1]. You can refer to [1] for detail.