Précédent : Echantillonnage Suivant : Algorithme d'estimation |
|
On utilise un développement
semi-quadratique de :
variable auxiliaire b, images auxiliaires Bx et By (gradients sur x et y)
échantillonnage de triplets (X, Bx, By) de manière alternée sur X et Bx, By
les pixels de Bx et By sont indépendants
et la proba. de X sachant Bx, By est gaussienne
le tirage de X s'effectue en une seule fois (contrairement à Gibbs et Metropolis),
le traitement est indépendant de la taille du voisinage,
ce qui en fait un algorithme rapide pour échantillonner la distribution a posteriori
|
![]() On calcule d'abord la transformée de Fourier F[h4 ] et DCT[Y], ainsi que W : |
![]() Tous les pixels b de Bx et By sont indépendants : |
![]() la distribution de X sachant B est gaussienne car |
|
![]() |
![]() |
![]() |
![]() |
|
|
![]() ![]() |
|
![]() |
![]() |
|
|
![]() |
|
![]() |
![]() |
|
|
![]() ![]() |
|
![]() |
![]() |
|
|
![]() |
|
|
![]() |
|
Distribution a priori
|
![]() On calcule d'abord la transformée de Fourier F[h4 ], ainsi que W0 : |
![]() Tous les pixels b de Bx et By sont indépendants : |
![]() la distribution de X sachant B est gaussienne car |
|
![]() |
|
![]() |
![]() |
|
|
![]() ![]() |
|
![]() |
![]() |
|
|
![]() |
|
![]() |
![]() |
|
|
![]() ![]() |
|
![]() |
![]() |
|
|
![]() |
|
|
![]() |
|