909
Views
9
CrossRef citations to date
0
Altmetric
Article

Dualization and Automatic Distributed Parameter Selection of Total Generalized Variation via Bilevel Optimization

ORCID Icon, ORCID Icon, &
Pages 887-932 | Received 17 Feb 2021, Accepted 19 Apr 2022, Published online: 05 May 2022

Figures & data

Figure 1. Gaussian denoising: Typical difference between TV (piecewise constant) and TGV reconstructions (piecewise affine).

Figure 1. Gaussian denoising: Typical difference between TV (piecewise constant) and TGV reconstructions (piecewise affine).

Figure 2. Suitability of the functional F(R·) as an upper level objective. Evaluation of F(Ru) where u solves the TGV denoising problem (1.1) (T = Id), for a variety of scalar parameters (α0,α1).

Figure 2. Suitability of the functional F(R·) as an upper level objective. Evaluation of F(Ru) where u solves the TGV denoising problem (1.1) (T = Id), for a variety of scalar parameters (α0,α1).

Figure 3. Test images, resolution 256×256.

Figure 3. Test images, resolution 256×256.

Figure 4. Upper level objective values vs projected gradient iterations for the problem (PTGVh) (right) of (scalar α0, spatial α1).

Figure 4. Upper level objective values vs projected gradient iterations for the problem (PTGVh) (right) of Figure 5 (scalar α0, spatial α1).

Figure 5. First row: noisy images. Second row: bilevel TV. Third row: Best scalar TGV (SSIM). Fourth row: bilevel TGV.

Figure 5. First row: noisy images. Second row: bilevel TV. Third row: Best scalar TGV (SSIM). Fourth row: bilevel TGV.

Figure 6. Details of the reconstructions shown in .

Figure 6. Details of the reconstructions shown in Figure 5.

Figure 7. First row: the computed regularization functions α for bilevel TV. Second row: the computed regularization functions α1 for bilevel TGV.

Figure 7. First row: the computed regularization functions α for bilevel TV. Second row: the computed regularization functions α1 for bilevel TGV.

Table 1. PSNR and SSIM comparisons for the images of .

Figure 8. Experiments with optimizing over a spatially varying α0. Top row: the automatically computed scalar parameters α0, that correspond to the images of the last row of . Middle row: the automatically computed spatially varying parameters α0, where α1 has been kept fixed (last row of ). The weight α0 is adapted to piecewise constant parts having there large values and hence promoting TV like behavior, see for instance the parrot image at the last row. On the contrary α0 has low values in piecewise smooth parts promoting a TGV like behavior reducing the staircasing.

Figure 8. Experiments with optimizing over a spatially varying α0. Top row: the automatically computed scalar parameters α0, that correspond to the images of the last row of Figure 5. Middle row: the automatically computed spatially varying parameters α0, where α1 has been kept fixed (last row of Figure 7). The weight α0 is adapted to piecewise constant parts having there large values and hence promoting TV like behavior, see for instance the parrot image at the last row. On the contrary α0 has low values in piecewise smooth parts promoting a TGV like behavior reducing the staircasing.

Figure 9. Histograms of the PSNR and SSIM differences between the bilevel spatial TGV and bilevel spatial TV (top) and between the bilevel spatial TGV and bilevel scalar TGV (bottom), for reconstructions of noisy versions of the 50 images on the left. In all bilevel algorithms the ground truth-free statistics-based upper level objective was used.

Figure 9. Histograms of the PSNR and SSIM differences between the bilevel spatial TGV and bilevel spatial TV (top) and between the bilevel spatial TGV and bilevel scalar TGV (bottom), for reconstructions of noisy versions of the 50 images on the left. In all bilevel algorithms the ground truth-free statistics-based upper level objective was used.