Convolutional Dictionary Learning with Spatial Mask¶
This example demonstrates the use of
cbpdndlmd.ConvBPDNMaskDictLearn
for convolutional dictionary
learning with a spatial mask, from a set of colour training images
[51]. The dictionary learning algorithm
is based on the hybrid mask decoupling / ADMM consensus dictionary
update [26].
from __future__ import print_function
from builtins import input
import pyfftw # See https://github.com/pyFFTW/pyFFTW/issues/40
import numpy as np
from sporco.admm import tvl2
from sporco.dictlrn import cbpdndl
from sporco.dictlrn import cbpdndlmd
from sporco import util
from sporco import signal
from sporco import plot
plot.config_notebook_plotting()
Load training images.
exim = util.ExampleImages(scaled=True, zoom=0.25)
S1 = exim.image('barbara.png', idxexp=np.s_[10:522, 100:612])
S2 = exim.image('kodim23.png', idxexp=np.s_[:, 60:572])
S = np.stack((S1, S2), axis=3)
Construct initial dictionary.
np.random.seed(12345)
D0 = np.random.randn(8, 8, 3, 32)
Create random mask and apply to training images.
frc = 0.5
W = signal.rndmask(S.shape[0:3] + (1,), frc, dtype=np.float32)
Sw = W * S
\(\ell_2\)-TV denoising with a spatial mask as a non-linear lowpass filter.
lmbda = 0.1
opt = tvl2.TVL2Denoise.Options({'Verbose': False, 'MaxMainIter': 200,
'DFidWeight': W, 'gEvalY': False, 'AutoRho': {'Enabled': True}})
b = tvl2.TVL2Denoise(Sw, lmbda, opt, caxis=2)
sl = b.solve()
sh = Sw - sl
CDL without a spatial mask using
dictlrn.cbpdndl.ConvBPDNDictLearn
. (Note that
prlcnscdl.ConvBPDNMaskDcplDictLearn_Consensus
solves the
same problem, but is substantially faster on a multi-core architecture.)
lmbda = 0.05
opt1 = cbpdndl.ConvBPDNDictLearn.Options({'Verbose': True,
'MaxMainIter': 200, 'AccurateDFid': True,
'CBPDN': {'rho': 50.0*lmbda + 0.5},
'CCMOD': {'rho': 3e2}}, dmethod='cns')
d1 = cbpdndl.ConvBPDNDictLearn(D0, sh, lmbda, opt1, dmethod='cns')
D1 = d1.solve()
Itn Fnc DFid ℓ1 Cnstr r_X s_X ρ_X r_D s_D ρ_D
--------------------------------------------------------------------------------------------------------
0 1.25e+03 1.27e+02 2.24e+04 8.19e-07 3.31e-01 2.40e+00 3.00e+00 4.31e-02 2.30e-01 3.00e+02
1 1.12e+03 1.53e+02 1.93e+04 7.31e-07 1.67e-01 7.28e-01 3.00e+00 5.11e-02 3.36e-01 3.00e+02
2 1.03e+03 7.68e+01 1.90e+04 6.72e-07 1.25e-01 5.48e-01 3.00e+00 3.15e-02 2.75e-01 3.00e+02
3 9.73e+02 5.52e+01 1.84e+04 5.69e-07 9.77e-02 4.53e-01 3.00e+00 2.31e-02 2.18e-01 3.00e+02
4 9.36e+02 6.36e+01 1.74e+04 6.26e-07 7.85e-02 3.87e-01 3.00e+00 1.95e-02 1.93e-01 3.00e+02
5 9.04e+02 4.99e+01 1.71e+04 8.75e-07 6.36e-02 3.38e-01 3.00e+00 1.55e-02 1.76e-01 3.00e+02
6 8.78e+02 4.73e+01 1.66e+04 8.47e-07 5.33e-02 3.09e-01 3.00e+00 1.31e-02 1.53e-01 3.00e+02
7 8.56e+02 4.75e+01 1.62e+04 8.47e-07 4.51e-02 2.81e-01 3.00e+00 1.12e-02 1.35e-01 3.00e+02
8 8.38e+02 4.53e+01 1.59e+04 8.48e-07 3.83e-02 2.56e-01 3.00e+00 9.59e-03 1.22e-01 3.00e+02
9 8.23e+02 4.36e+01 1.56e+04 8.82e-07 3.30e-02 2.36e-01 3.00e+00 8.36e-03 1.11e-01 3.00e+02
10 8.09e+02 4.30e+01 1.53e+04 6.49e-07 2.87e-02 2.19e-01 3.00e+00 7.40e-03 1.01e-01 3.00e+02
11 7.98e+02 4.24e+01 1.51e+04 8.73e-07 2.52e-02 2.04e-01 3.00e+00 6.61e-03 9.29e-02 3.00e+02
12 7.89e+02 4.18e+01 1.49e+04 6.48e-07 2.24e-02 1.90e-01 3.00e+00 5.97e-03 8.62e-02 3.00e+02
13 7.80e+02 4.13e+01 1.48e+04 8.05e-07 2.00e-02 1.78e-01 3.00e+00 5.44e-03 8.03e-02 3.00e+02
14 7.73e+02 4.07e+01 1.46e+04 8.08e-07 1.81e-02 1.67e-01 3.00e+00 5.00e-03 7.53e-02 3.00e+02
15 7.66e+02 4.03e+01 1.45e+04 8.92e-07 1.64e-02 1.57e-01 3.00e+00 4.64e-03 7.10e-02 3.00e+02
16 7.60e+02 3.99e+01 1.44e+04 7.33e-07 1.51e-02 1.49e-01 3.00e+00 4.34e-03 6.73e-02 3.00e+02
17 7.55e+02 3.95e+01 1.43e+04 7.59e-07 1.39e-02 1.41e-01 3.00e+00 4.09e-03 6.40e-02 3.00e+02
18 7.50e+02 3.91e+01 1.42e+04 7.48e-07 1.29e-02 1.34e-01 3.00e+00 3.87e-03 6.11e-02 3.00e+02
19 7.45e+02 3.88e+01 1.41e+04 7.35e-07 1.21e-02 1.28e-01 3.00e+00 3.69e-03 5.85e-02 3.00e+02
20 7.39e+02 3.81e+01 1.40e+04 7.51e-07 2.17e-02 1.20e-01 1.50e+00 5.47e-03 6.04e-02 3.00e+02
21 7.33e+02 3.65e+01 1.39e+04 7.89e-07 2.06e-02 1.10e-01 1.50e+00 4.94e-03 6.30e-02 3.00e+02
22 7.28e+02 3.69e+01 1.38e+04 7.74e-07 1.94e-02 1.03e-01 1.50e+00 5.35e-03 6.22e-02 3.00e+02
23 7.24e+02 3.76e+01 1.37e+04 6.86e-07 1.82e-02 9.59e-02 1.50e+00 4.85e-03 5.93e-02 3.00e+02
24 7.20e+02 3.73e+01 1.37e+04 7.21e-07 1.70e-02 9.00e-02 1.50e+00 4.83e-03 5.64e-02 3.00e+02
25 7.17e+02 3.70e+01 1.36e+04 6.58e-07 1.59e-02 8.50e-02 1.50e+00 4.48e-03 5.40e-02 3.00e+02
26 7.14e+02 3.69e+01 1.35e+04 7.93e-07 1.49e-02 8.05e-02 1.50e+00 4.41e-03 5.22e-02 3.00e+02
27 7.11e+02 3.67e+01 1.35e+04 6.46e-07 1.40e-02 7.63e-02 1.50e+00 4.16e-03 5.07e-02 3.00e+02
28 7.08e+02 3.65e+01 1.34e+04 8.24e-07 1.31e-02 7.24e-02 1.50e+00 4.04e-03 4.94e-02 3.00e+02
29 7.06e+02 3.63e+01 1.34e+04 8.87e-07 1.24e-02 6.90e-02 1.50e+00 3.84e-03 4.81e-02 3.00e+02
30 7.03e+02 3.62e+01 1.33e+04 6.02e-07 1.17e-02 6.58e-02 1.50e+00 3.72e-03 4.70e-02 3.00e+02
31 7.01e+02 3.59e+01 1.33e+04 7.12e-07 1.11e-02 6.29e-02 1.50e+00 3.56e-03 4.61e-02 3.00e+02
32 6.99e+02 3.57e+01 1.33e+04 7.20e-07 1.05e-02 6.03e-02 1.50e+00 3.45e-03 4.52e-02 3.00e+02
33 6.97e+02 3.56e+01 1.32e+04 7.67e-07 1.01e-02 5.79e-02 1.50e+00 3.32e-03 4.43e-02 3.00e+02
34 6.96e+02 3.56e+01 1.32e+04 7.25e-07 9.60e-03 5.56e-02 1.50e+00 3.22e-03 4.34e-02 3.00e+02
35 6.94e+02 3.54e+01 1.32e+04 8.45e-07 9.17e-03 5.36e-02 1.50e+00 3.11e-03 4.26e-02 3.00e+02
36 6.93e+02 3.53e+01 1.31e+04 7.14e-07 8.80e-03 5.17e-02 1.50e+00 3.02e-03 4.18e-02 3.00e+02
37 6.91e+02 3.52e+01 1.31e+04 6.35e-07 8.45e-03 5.00e-02 1.50e+00 2.93e-03 4.11e-02 3.00e+02
38 6.90e+02 3.50e+01 1.31e+04 7.05e-07 8.14e-03 4.85e-02 1.50e+00 2.85e-03 4.04e-02 3.00e+02
39 6.88e+02 3.49e+01 1.31e+04 8.34e-07 7.85e-03 4.71e-02 1.50e+00 2.78e-03 3.98e-02 3.00e+02
40 6.87e+02 3.48e+01 1.30e+04 7.36e-07 7.59e-03 4.57e-02 1.50e+00 2.71e-03 3.91e-02 3.00e+02
41 6.86e+02 3.47e+01 1.30e+04 8.71e-07 7.35e-03 4.45e-02 1.50e+00 2.66e-03 3.85e-02 3.00e+02
42 6.84e+02 3.47e+01 1.30e+04 8.50e-07 7.13e-03 4.33e-02 1.50e+00 2.60e-03 3.79e-02 3.00e+02
43 6.83e+02 3.46e+01 1.30e+04 6.50e-07 6.93e-03 4.23e-02 1.50e+00 2.55e-03 3.73e-02 3.00e+02
44 6.82e+02 3.45e+01 1.29e+04 7.37e-07 6.74e-03 4.13e-02 1.50e+00 2.50e-03 3.68e-02 3.00e+02
45 6.81e+02 3.45e+01 1.29e+04 8.36e-07 6.57e-03 4.04e-02 1.50e+00 2.45e-03 3.62e-02 3.00e+02
46 6.80e+02 3.44e+01 1.29e+04 7.00e-07 6.40e-03 3.96e-02 1.50e+00 2.40e-03 3.57e-02 3.00e+02
47 6.78e+02 3.43e+01 1.29e+04 7.53e-07 6.25e-03 3.87e-02 1.50e+00 2.36e-03 3.53e-02 3.00e+02
48 6.77e+02 3.42e+01 1.29e+04 8.00e-07 6.11e-03 3.80e-02 1.50e+00 2.32e-03 3.48e-02 3.00e+02
49 6.76e+02 3.42e+01 1.28e+04 7.92e-07 5.98e-03 3.73e-02 1.50e+00 2.29e-03 3.44e-02 3.00e+02
50 6.75e+02 3.41e+01 1.28e+04 8.49e-07 5.85e-03 3.67e-02 1.50e+00 2.25e-03 3.39e-02 3.00e+02
51 6.74e+02 3.40e+01 1.28e+04 6.60e-07 5.74e-03 3.61e-02 1.50e+00 2.22e-03 3.35e-02 3.00e+02
52 6.73e+02 3.39e+01 1.28e+04 6.14e-07 5.64e-03 3.55e-02 1.50e+00 2.19e-03 3.31e-02 3.00e+02
53 6.72e+02 3.39e+01 1.28e+04 8.21e-07 5.54e-03 3.49e-02 1.50e+00 2.16e-03 3.27e-02 3.00e+02
54 6.71e+02 3.38e+01 1.27e+04 7.55e-07 5.45e-03 3.44e-02 1.50e+00 2.14e-03 3.24e-02 3.00e+02
55 6.70e+02 3.38e+01 1.27e+04 9.59e-07 5.36e-03 3.39e-02 1.50e+00 2.12e-03 3.20e-02 3.00e+02
56 6.69e+02 3.37e+01 1.27e+04 7.69e-07 5.27e-03 3.35e-02 1.50e+00 2.09e-03 3.17e-02 3.00e+02
57 6.68e+02 3.37e+01 1.27e+04 6.65e-07 5.19e-03 3.30e-02 1.50e+00 2.07e-03 3.14e-02 3.00e+02
58 6.67e+02 3.36e+01 1.27e+04 9.44e-07 5.12e-03 3.26e-02 1.50e+00 2.04e-03 3.10e-02 3.00e+02
59 6.67e+02 3.36e+01 1.27e+04 7.42e-07 5.05e-03 3.22e-02 1.50e+00 2.02e-03 3.07e-02 3.00e+02
60 6.66e+02 3.35e+01 1.26e+04 7.61e-07 4.98e-03 3.18e-02 1.50e+00 2.00e-03 3.04e-02 3.00e+02
61 6.65e+02 3.35e+01 1.26e+04 9.00e-07 4.91e-03 3.14e-02 1.50e+00 1.98e-03 3.01e-02 3.00e+02
62 6.64e+02 3.34e+01 1.26e+04 6.54e-07 4.85e-03 3.11e-02 1.50e+00 1.96e-03 2.98e-02 3.00e+02
63 6.63e+02 3.34e+01 1.26e+04 6.31e-07 4.79e-03 3.07e-02 1.50e+00 1.94e-03 2.95e-02 3.00e+02
64 6.62e+02 3.34e+01 1.26e+04 8.93e-07 4.73e-03 3.04e-02 1.50e+00 1.92e-03 2.92e-02 3.00e+02
65 6.61e+02 3.33e+01 1.26e+04 7.64e-07 4.68e-03 3.01e-02 1.50e+00 1.90e-03 2.89e-02 3.00e+02
66 6.60e+02 3.33e+01 1.25e+04 5.99e-07 4.62e-03 2.98e-02 1.50e+00 1.89e-03 2.87e-02 3.00e+02
67 6.60e+02 3.32e+01 1.25e+04 9.00e-07 4.58e-03 2.95e-02 1.50e+00 1.87e-03 2.84e-02 3.00e+02
68 6.59e+02 3.32e+01 1.25e+04 7.91e-07 4.53e-03 2.92e-02 1.50e+00 1.86e-03 2.82e-02 3.00e+02
69 6.58e+02 3.31e+01 1.25e+04 8.29e-07 4.48e-03 2.90e-02 1.50e+00 1.84e-03 2.80e-02 3.00e+02
70 6.57e+02 3.31e+01 1.25e+04 8.68e-07 4.44e-03 2.87e-02 1.50e+00 1.83e-03 2.78e-02 3.00e+02
71 6.56e+02 3.30e+01 1.25e+04 8.44e-07 4.40e-03 2.85e-02 1.50e+00 1.81e-03 2.75e-02 3.00e+02
72 6.55e+02 3.30e+01 1.24e+04 8.40e-07 4.36e-03 2.82e-02 1.50e+00 1.80e-03 2.73e-02 3.00e+02
73 6.54e+02 3.29e+01 1.24e+04 6.78e-07 4.33e-03 2.80e-02 1.50e+00 1.79e-03 2.71e-02 3.00e+02
74 6.54e+02 3.29e+01 1.24e+04 8.79e-07 4.29e-03 2.78e-02 1.50e+00 1.78e-03 2.69e-02 3.00e+02
75 6.53e+02 3.29e+01 1.24e+04 6.40e-07 4.26e-03 2.76e-02 1.50e+00 1.77e-03 2.67e-02 3.00e+02
76 6.52e+02 3.28e+01 1.24e+04 8.52e-07 4.23e-03 2.74e-02 1.50e+00 1.75e-03 2.65e-02 3.00e+02
77 6.51e+02 3.28e+01 1.24e+04 8.01e-07 4.19e-03 2.72e-02 1.50e+00 1.74e-03 2.63e-02 3.00e+02
78 6.50e+02 3.28e+01 1.24e+04 8.34e-07 4.16e-03 2.70e-02 1.50e+00 1.73e-03 2.61e-02 3.00e+02
79 6.49e+02 3.27e+01 1.23e+04 7.41e-07 4.13e-03 2.68e-02 1.50e+00 1.72e-03 2.59e-02 3.00e+02
80 6.49e+02 3.27e+01 1.23e+04 6.76e-07 4.10e-03 2.66e-02 1.50e+00 1.71e-03 2.57e-02 3.00e+02
81 6.48e+02 3.27e+01 1.23e+04 7.73e-07 4.08e-03 2.65e-02 1.50e+00 1.70e-03 2.55e-02 3.00e+02
82 6.47e+02 3.26e+01 1.23e+04 7.76e-07 4.05e-03 2.63e-02 1.50e+00 1.69e-03 2.54e-02 3.00e+02
83 6.46e+02 3.26e+01 1.23e+04 7.91e-07 4.02e-03 2.61e-02 1.50e+00 1.68e-03 2.52e-02 3.00e+02
84 6.45e+02 3.25e+01 1.23e+04 7.07e-07 4.00e-03 2.60e-02 1.50e+00 1.67e-03 2.50e-02 3.00e+02
85 6.45e+02 3.25e+01 1.22e+04 8.10e-07 3.98e-03 2.59e-02 1.50e+00 1.66e-03 2.49e-02 3.00e+02
86 6.44e+02 3.25e+01 1.22e+04 8.95e-07 3.95e-03 2.57e-02 1.50e+00 1.65e-03 2.47e-02 3.00e+02
87 6.43e+02 3.24e+01 1.22e+04 7.76e-07 3.93e-03 2.56e-02 1.50e+00 1.64e-03 2.46e-02 3.00e+02
88 6.42e+02 3.24e+01 1.22e+04 6.87e-07 3.91e-03 2.55e-02 1.50e+00 1.64e-03 2.44e-02 3.00e+02
89 6.41e+02 3.24e+01 1.22e+04 7.81e-07 3.90e-03 2.54e-02 1.50e+00 1.63e-03 2.43e-02 3.00e+02
90 6.41e+02 3.23e+01 1.22e+04 7.10e-07 3.88e-03 2.52e-02 1.50e+00 1.62e-03 2.41e-02 3.00e+02
91 6.40e+02 3.23e+01 1.22e+04 7.54e-07 3.86e-03 2.51e-02 1.50e+00 1.61e-03 2.40e-02 3.00e+02
92 6.39e+02 3.23e+01 1.21e+04 7.26e-07 3.84e-03 2.49e-02 1.50e+00 1.61e-03 2.39e-02 3.00e+02
93 6.38e+02 3.22e+01 1.21e+04 7.28e-07 3.82e-03 2.48e-02 1.50e+00 1.60e-03 2.37e-02 3.00e+02
94 6.37e+02 3.22e+01 1.21e+04 8.94e-07 3.81e-03 2.47e-02 1.50e+00 1.59e-03 2.36e-02 3.00e+02
95 6.37e+02 3.22e+01 1.21e+04 7.58e-07 3.79e-03 2.46e-02 1.50e+00 1.58e-03 2.35e-02 3.00e+02
96 6.36e+02 3.21e+01 1.21e+04 7.55e-07 3.77e-03 2.45e-02 1.50e+00 1.58e-03 2.33e-02 3.00e+02
97 6.35e+02 3.21e+01 1.21e+04 7.17e-07 3.76e-03 2.44e-02 1.50e+00 1.57e-03 2.32e-02 3.00e+02
98 6.34e+02 3.21e+01 1.20e+04 8.54e-07 3.74e-03 2.43e-02 1.50e+00 1.56e-03 2.31e-02 3.00e+02
99 6.34e+02 3.20e+01 1.20e+04 5.65e-07 3.73e-03 2.42e-02 1.50e+00 1.56e-03 2.30e-02 3.00e+02
100 6.33e+02 3.20e+01 1.20e+04 8.57e-07 3.71e-03 2.41e-02 1.50e+00 1.55e-03 2.28e-02 3.00e+02
101 6.32e+02 3.20e+01 1.20e+04 7.44e-07 3.69e-03 2.40e-02 1.50e+00 1.55e-03 2.27e-02 3.00e+02
102 6.31e+02 3.20e+01 1.20e+04 7.59e-07 3.68e-03 2.39e-02 1.50e+00 1.54e-03 2.26e-02 3.00e+02
103 6.30e+02 3.19e+01 1.20e+04 6.93e-07 3.67e-03 2.38e-02 1.50e+00 1.53e-03 2.25e-02 3.00e+02
104 6.30e+02 3.19e+01 1.20e+04 7.23e-07 3.65e-03 2.37e-02 1.50e+00 1.53e-03 2.23e-02 3.00e+02
105 6.29e+02 3.19e+01 1.19e+04 8.91e-07 3.64e-03 2.36e-02 1.50e+00 1.52e-03 2.22e-02 3.00e+02
106 6.28e+02 3.19e+01 1.19e+04 7.23e-07 3.62e-03 2.35e-02 1.50e+00 1.51e-03 2.21e-02 3.00e+02
107 6.28e+02 3.18e+01 1.19e+04 6.34e-07 3.61e-03 2.34e-02 1.50e+00 1.51e-03 2.20e-02 3.00e+02
108 6.27e+02 3.18e+01 1.19e+04 8.35e-07 3.60e-03 2.33e-02 1.50e+00 1.50e-03 2.19e-02 3.00e+02
109 6.26e+02 3.18e+01 1.19e+04 6.77e-07 3.59e-03 2.32e-02 1.50e+00 1.49e-03 2.18e-02 3.00e+02
110 6.25e+02 3.18e+01 1.19e+04 6.98e-07 3.57e-03 2.32e-02 1.50e+00 1.49e-03 2.17e-02 3.00e+02
111 6.25e+02 3.17e+01 1.19e+04 8.49e-07 3.56e-03 2.31e-02 1.50e+00 1.48e-03 2.16e-02 3.00e+02
112 6.24e+02 3.17e+01 1.18e+04 7.45e-07 3.55e-03 2.30e-02 1.50e+00 1.48e-03 2.15e-02 3.00e+02
113 6.23e+02 3.17e+01 1.18e+04 1.06e-06 3.55e-03 2.30e-02 1.50e+00 1.47e-03 2.14e-02 3.00e+02
114 6.22e+02 3.17e+01 1.18e+04 8.15e-07 3.54e-03 2.29e-02 1.50e+00 1.47e-03 2.13e-02 3.00e+02
115 6.22e+02 3.16e+01 1.18e+04 8.43e-07 3.53e-03 2.28e-02 1.50e+00 1.46e-03 2.12e-02 3.00e+02
116 6.21e+02 3.16e+01 1.18e+04 5.96e-07 3.52e-03 2.28e-02 1.50e+00 1.46e-03 2.11e-02 3.00e+02
117 6.20e+02 3.16e+01 1.18e+04 7.34e-07 3.51e-03 2.27e-02 1.50e+00 1.45e-03 2.10e-02 3.00e+02
118 6.20e+02 3.16e+01 1.18e+04 6.75e-07 3.50e-03 2.27e-02 1.50e+00 1.45e-03 2.09e-02 3.00e+02
119 6.19e+02 3.15e+01 1.17e+04 7.70e-07 3.49e-03 2.26e-02 1.50e+00 1.44e-03 2.08e-02 3.00e+02
120 6.18e+02 3.15e+01 1.17e+04 7.89e-07 3.48e-03 2.25e-02 1.50e+00 1.44e-03 2.07e-02 3.00e+02
121 6.17e+02 3.15e+01 1.17e+04 6.94e-07 3.48e-03 2.25e-02 1.50e+00 1.44e-03 2.06e-02 3.00e+02
122 6.17e+02 3.14e+01 1.17e+04 8.16e-07 3.47e-03 2.24e-02 1.50e+00 1.43e-03 2.05e-02 3.00e+02
123 6.16e+02 3.14e+01 1.17e+04 7.76e-07 3.46e-03 2.24e-02 1.50e+00 1.43e-03 2.04e-02 3.00e+02
124 6.15e+02 3.14e+01 1.17e+04 7.72e-07 3.46e-03 2.23e-02 1.50e+00 1.43e-03 2.03e-02 3.00e+02
125 6.15e+02 3.14e+01 1.17e+04 7.69e-07 3.45e-03 2.22e-02 1.50e+00 1.42e-03 2.02e-02 3.00e+02
126 6.14e+02 3.14e+01 1.16e+04 8.37e-07 3.45e-03 2.22e-02 1.50e+00 1.42e-03 2.01e-02 3.00e+02
127 6.13e+02 3.14e+01 1.16e+04 9.47e-07 3.44e-03 2.21e-02 1.50e+00 1.41e-03 2.01e-02 3.00e+02
128 6.12e+02 3.13e+01 1.16e+04 8.30e-07 3.43e-03 2.20e-02 1.50e+00 1.41e-03 2.00e-02 3.00e+02
129 6.12e+02 3.13e+01 1.16e+04 8.60e-07 3.42e-03 2.19e-02 1.50e+00 1.40e-03 1.99e-02 3.00e+02
130 6.11e+02 3.13e+01 1.16e+04 9.22e-07 3.41e-03 2.19e-02 1.50e+00 1.40e-03 1.98e-02 3.00e+02
131 6.10e+02 3.13e+01 1.16e+04 8.86e-07 3.40e-03 2.18e-02 1.50e+00 1.40e-03 1.97e-02 3.00e+02
132 6.10e+02 3.13e+01 1.16e+04 8.01e-07 3.39e-03 2.17e-02 1.50e+00 1.39e-03 1.96e-02 3.00e+02
133 6.09e+02 3.12e+01 1.16e+04 7.71e-07 3.39e-03 2.17e-02 1.50e+00 1.39e-03 1.95e-02 3.00e+02
134 6.08e+02 3.12e+01 1.15e+04 9.10e-07 3.38e-03 2.16e-02 1.50e+00 1.38e-03 1.94e-02 3.00e+02
135 6.08e+02 3.12e+01 1.15e+04 8.38e-07 3.37e-03 2.16e-02 1.50e+00 1.37e-03 1.93e-02 3.00e+02
136 6.07e+02 3.12e+01 1.15e+04 7.42e-07 3.37e-03 2.15e-02 1.50e+00 1.37e-03 1.92e-02 3.00e+02
137 6.06e+02 3.12e+01 1.15e+04 7.45e-07 3.36e-03 2.15e-02 1.50e+00 1.37e-03 1.91e-02 3.00e+02
138 6.06e+02 3.12e+01 1.15e+04 6.48e-07 3.35e-03 2.14e-02 1.50e+00 1.36e-03 1.90e-02 3.00e+02
139 6.05e+02 3.11e+01 1.15e+04 7.51e-07 3.34e-03 2.14e-02 1.50e+00 1.36e-03 1.89e-02 3.00e+02
140 6.04e+02 3.11e+01 1.15e+04 8.45e-07 3.34e-03 2.13e-02 1.50e+00 1.35e-03 1.89e-02 3.00e+02
141 6.04e+02 3.11e+01 1.15e+04 8.34e-07 3.33e-03 2.13e-02 1.50e+00 1.35e-03 1.88e-02 3.00e+02
142 6.03e+02 3.11e+01 1.14e+04 6.15e-07 3.32e-03 2.12e-02 1.50e+00 1.34e-03 1.87e-02 3.00e+02
143 6.02e+02 3.11e+01 1.14e+04 9.42e-07 3.32e-03 2.12e-02 1.50e+00 1.34e-03 1.87e-02 3.00e+02
144 6.02e+02 3.11e+01 1.14e+04 7.36e-07 3.31e-03 2.11e-02 1.50e+00 1.34e-03 1.86e-02 3.00e+02
145 6.01e+02 3.10e+01 1.14e+04 7.76e-07 3.31e-03 2.11e-02 1.50e+00 1.34e-03 1.85e-02 3.00e+02
146 6.00e+02 3.10e+01 1.14e+04 7.35e-07 3.30e-03 2.10e-02 1.50e+00 1.33e-03 1.85e-02 3.00e+02
147 6.00e+02 3.10e+01 1.14e+04 8.02e-07 3.30e-03 2.10e-02 1.50e+00 1.33e-03 1.84e-02 3.00e+02
148 5.99e+02 3.10e+01 1.14e+04 9.67e-07 3.30e-03 2.10e-02 1.50e+00 1.33e-03 1.83e-02 3.00e+02
149 5.98e+02 3.10e+01 1.13e+04 8.39e-07 3.29e-03 2.10e-02 1.50e+00 1.32e-03 1.82e-02 3.00e+02
150 5.98e+02 3.09e+01 1.13e+04 1.03e-06 3.29e-03 2.09e-02 1.50e+00 1.32e-03 1.82e-02 3.00e+02
151 5.97e+02 3.09e+01 1.13e+04 7.95e-07 3.29e-03 2.09e-02 1.50e+00 1.32e-03 1.81e-02 3.00e+02
152 5.96e+02 3.09e+01 1.13e+04 7.22e-07 3.29e-03 2.09e-02 1.50e+00 1.31e-03 1.80e-02 3.00e+02
153 5.96e+02 3.09e+01 1.13e+04 9.42e-07 3.29e-03 2.09e-02 1.50e+00 1.31e-03 1.79e-02 3.00e+02
154 5.95e+02 3.09e+01 1.13e+04 9.80e-07 3.28e-03 2.09e-02 1.50e+00 1.31e-03 1.79e-02 3.00e+02
155 5.95e+02 3.09e+01 1.13e+04 9.19e-07 3.28e-03 2.08e-02 1.50e+00 1.30e-03 1.78e-02 3.00e+02
156 5.94e+02 3.09e+01 1.13e+04 6.74e-07 3.28e-03 2.08e-02 1.50e+00 1.30e-03 1.78e-02 3.00e+02
157 5.93e+02 3.09e+01 1.12e+04 8.17e-07 3.28e-03 2.08e-02 1.50e+00 1.30e-03 1.77e-02 3.00e+02
158 5.93e+02 3.09e+01 1.12e+04 7.41e-07 3.27e-03 2.08e-02 1.50e+00 1.30e-03 1.77e-02 3.00e+02
159 5.92e+02 3.08e+01 1.12e+04 7.07e-07 3.27e-03 2.07e-02 1.50e+00 1.29e-03 1.76e-02 3.00e+02
160 5.91e+02 3.08e+01 1.12e+04 6.39e-07 3.26e-03 2.07e-02 1.50e+00 1.29e-03 1.76e-02 3.00e+02
161 5.91e+02 3.08e+01 1.12e+04 9.40e-07 3.26e-03 2.06e-02 1.50e+00 1.29e-03 1.75e-02 3.00e+02
162 5.90e+02 3.08e+01 1.12e+04 7.25e-07 3.26e-03 2.06e-02 1.50e+00 1.29e-03 1.75e-02 3.00e+02
163 5.90e+02 3.08e+01 1.12e+04 6.49e-07 3.25e-03 2.06e-02 1.50e+00 1.29e-03 1.74e-02 3.00e+02
164 5.89e+02 3.08e+01 1.12e+04 7.05e-07 3.25e-03 2.05e-02 1.50e+00 1.28e-03 1.74e-02 3.00e+02
165 5.88e+02 3.08e+01 1.12e+04 8.28e-07 3.25e-03 2.05e-02 1.50e+00 1.28e-03 1.73e-02 3.00e+02
166 5.88e+02 3.08e+01 1.11e+04 7.10e-07 3.25e-03 2.05e-02 1.50e+00 1.28e-03 1.73e-02 3.00e+02
167 5.87e+02 3.08e+01 1.11e+04 8.04e-07 3.24e-03 2.05e-02 1.50e+00 1.28e-03 1.72e-02 3.00e+02
168 5.86e+02 3.08e+01 1.11e+04 7.37e-07 3.24e-03 2.04e-02 1.50e+00 1.27e-03 1.72e-02 3.00e+02
169 5.86e+02 3.07e+01 1.11e+04 6.23e-07 3.24e-03 2.04e-02 1.50e+00 1.27e-03 1.71e-02 3.00e+02
170 5.85e+02 3.07e+01 1.11e+04 5.73e-07 3.23e-03 2.03e-02 1.50e+00 1.27e-03 1.71e-02 3.00e+02
171 5.85e+02 3.07e+01 1.11e+04 8.24e-07 3.23e-03 2.03e-02 1.50e+00 1.27e-03 1.70e-02 3.00e+02
172 5.84e+02 3.07e+01 1.11e+04 7.28e-07 3.23e-03 2.03e-02 1.50e+00 1.27e-03 1.70e-02 3.00e+02
173 5.83e+02 3.07e+01 1.11e+04 8.42e-07 3.23e-03 2.03e-02 1.50e+00 1.26e-03 1.69e-02 3.00e+02
174 5.83e+02 3.07e+01 1.10e+04 7.32e-07 3.22e-03 2.02e-02 1.50e+00 1.26e-03 1.69e-02 3.00e+02
175 5.82e+02 3.07e+01 1.10e+04 8.24e-07 3.22e-03 2.02e-02 1.50e+00 1.26e-03 1.68e-02 3.00e+02
176 5.82e+02 3.07e+01 1.10e+04 8.34e-07 3.21e-03 2.02e-02 1.50e+00 1.26e-03 1.68e-02 3.00e+02
177 5.81e+02 3.07e+01 1.10e+04 6.00e-07 3.21e-03 2.01e-02 1.50e+00 1.26e-03 1.67e-02 3.00e+02
178 5.80e+02 3.06e+01 1.10e+04 8.71e-07 3.21e-03 2.01e-02 1.50e+00 1.25e-03 1.67e-02 3.00e+02
179 5.80e+02 3.06e+01 1.10e+04 8.56e-07 3.21e-03 2.01e-02 1.50e+00 1.25e-03 1.66e-02 3.00e+02
180 5.79e+02 3.06e+01 1.10e+04 8.70e-07 3.21e-03 2.01e-02 1.50e+00 1.25e-03 1.66e-02 3.00e+02
181 5.79e+02 3.06e+01 1.10e+04 6.90e-07 3.21e-03 2.01e-02 1.50e+00 1.25e-03 1.66e-02 3.00e+02
182 5.78e+02 3.06e+01 1.09e+04 8.65e-07 3.21e-03 2.01e-02 1.50e+00 1.25e-03 1.65e-02 3.00e+02
183 5.78e+02 3.06e+01 1.09e+04 7.93e-07 3.22e-03 2.01e-02 1.50e+00 1.25e-03 1.65e-02 3.00e+02
184 5.77e+02 3.06e+01 1.09e+04 5.19e-07 3.22e-03 2.01e-02 1.50e+00 1.25e-03 1.65e-02 3.00e+02
185 5.76e+02 3.06e+01 1.09e+04 8.18e-07 3.22e-03 2.01e-02 1.50e+00 1.25e-03 1.64e-02 3.00e+02
186 5.76e+02 3.06e+01 1.09e+04 7.48e-07 3.23e-03 2.01e-02 1.50e+00 1.24e-03 1.64e-02 3.00e+02
187 5.75e+02 3.06e+01 1.09e+04 8.47e-07 3.23e-03 2.02e-02 1.50e+00 1.24e-03 1.64e-02 3.00e+02
188 5.75e+02 3.06e+01 1.09e+04 8.91e-07 3.24e-03 2.02e-02 1.50e+00 1.24e-03 1.63e-02 3.00e+02
189 5.74e+02 3.06e+01 1.09e+04 7.38e-07 3.24e-03 2.02e-02 1.50e+00 1.24e-03 1.63e-02 3.00e+02
190 5.73e+02 3.06e+01 1.09e+04 6.98e-07 3.24e-03 2.02e-02 1.50e+00 1.24e-03 1.63e-02 3.00e+02
191 5.73e+02 3.06e+01 1.08e+04 8.36e-07 3.25e-03 2.02e-02 1.50e+00 1.24e-03 1.62e-02 3.00e+02
192 5.72e+02 3.05e+01 1.08e+04 5.53e-07 3.25e-03 2.02e-02 1.50e+00 1.24e-03 1.62e-02 3.00e+02
193 5.72e+02 3.05e+01 1.08e+04 7.05e-07 3.25e-03 2.02e-02 1.50e+00 1.24e-03 1.62e-02 3.00e+02
194 5.71e+02 3.05e+01 1.08e+04 7.46e-07 3.26e-03 2.02e-02 1.50e+00 1.24e-03 1.62e-02 3.00e+02
195 5.71e+02 3.05e+01 1.08e+04 9.17e-07 3.26e-03 2.02e-02 1.50e+00 1.24e-03 1.61e-02 3.00e+02
196 5.70e+02 3.05e+01 1.08e+04 7.59e-07 3.26e-03 2.02e-02 1.50e+00 1.24e-03 1.61e-02 3.00e+02
197 5.69e+02 3.05e+01 1.08e+04 9.11e-07 3.26e-03 2.02e-02 1.50e+00 1.24e-03 1.61e-02 3.00e+02
198 5.69e+02 3.05e+01 1.08e+04 7.52e-07 3.26e-03 2.01e-02 1.50e+00 1.24e-03 1.60e-02 3.00e+02
199 5.68e+02 3.05e+01 1.08e+04 7.81e-07 3.26e-03 2.01e-02 1.50e+00 1.24e-03 1.60e-02 3.00e+02
--------------------------------------------------------------------------------------------------------
Reconstruct from the CDL solution without a spatial mask.
sr1 = d1.reconstruct().squeeze() + sl
CDL with a spatial mask using
cbpdndlmd.ConvBPDNMaskDictLearn
.
opt2 = cbpdndlmd.ConvBPDNMaskDictLearn.Options({'Verbose': True,
'MaxMainIter': 200, 'AccurateDFid': True,
'CBPDN': {'rho': 50.0*lmbda + 0.5},
'CCMOD': {'rho': 1.0}}, dmethod='cns')
d2 = cbpdndlmd.ConvBPDNMaskDictLearn(D0, sh, lmbda, W, opt2,
dmethod='cns')
D2 = d2.solve()
Itn Fnc DFid ℓ1 Cnstr r_X s_X ρ_X r_D s_D ρ_D
--------------------------------------------------------------------------------------------------------
0 1.70e+03 2.19e+02 2.96e+04 6.48e-07 1.65e-01 1.01e+00 3.00e+00 2.00e-01 1.85e+03 1.00e+00
1 4.99e+02 3.98e+01 9.18e+03 8.30e-07 1.37e-01 9.73e-01 3.00e+00 6.31e-02 2.47e+02 1.00e+00
2 9.29e+02 7.14e+01 1.72e+04 5.93e-07 1.02e-01 9.97e-01 3.00e+00 1.03e-01 8.04e+02 1.00e+00
3 5.67e+02 4.68e+01 1.04e+04 8.08e-07 7.86e-02 9.77e-01 3.00e+00 9.78e-02 4.87e+02 1.00e+00
4 6.02e+02 1.87e+01 1.17e+04 6.45e-07 6.77e-02 9.71e-01 3.00e+00 1.33e-01 5.03e+02 1.00e+00
5 5.48e+02 1.44e+01 1.07e+04 7.61e-07 5.45e-02 9.52e-01 3.00e+00 1.01e-01 3.55e+02 1.00e+00
6 4.87e+02 1.34e+01 9.47e+03 9.34e-07 4.49e-02 9.30e-01 3.00e+00 8.78e-02 3.37e+02 1.00e+00
7 4.52e+02 1.13e+01 8.82e+03 7.67e-07 3.95e-02 9.27e-01 3.00e+00 9.91e-02 3.58e+02 1.00e+00
8 4.08e+02 1.18e+01 7.92e+03 8.98e-07 3.39e-02 9.20e-01 3.00e+00 1.78e-01 3.97e+02 1.00e+00
9 3.76e+02 1.29e+01 7.25e+03 9.44e-07 2.97e-02 9.03e-01 3.00e+00 6.68e-02 8.83e+01 1.00e+00
10 3.48e+02 1.34e+01 6.70e+03 8.16e-07 2.69e-02 8.76e-01 3.00e+00 9.28e-02 9.62e+01 1.00e+00
11 3.16e+02 1.37e+01 6.04e+03 7.47e-07 2.47e-02 8.47e-01 3.00e+00 5.54e-02 7.89e+01 1.00e+00
12 2.85e+02 1.44e+01 5.41e+03 6.43e-07 2.28e-02 8.29e-01 3.00e+00 6.39e-02 7.68e+01 1.00e+00
13 2.60e+02 1.49e+01 4.91e+03 7.24e-07 2.12e-02 8.14e-01 3.00e+00 5.65e-02 5.25e+01 1.00e+00
14 2.39e+02 1.54e+01 4.47e+03 8.97e-07 1.98e-02 7.90e-01 3.00e+00 5.59e-02 4.89e+01 1.00e+00
15 2.18e+02 1.58e+01 4.05e+03 8.40e-07 1.85e-02 7.61e-01 3.00e+00 6.98e-02 2.42e+01 1.00e+00
16 2.00e+02 1.62e+01 3.67e+03 7.92e-07 1.76e-02 7.36e-01 3.00e+00 8.37e-02 1.99e+01 1.00e+00
17 1.84e+02 1.67e+01 3.34e+03 9.33e-07 1.67e-02 7.11e-01 3.00e+00 7.98e-02 2.41e+01 1.00e+00
18 1.69e+02 1.72e+01 3.03e+03 8.12e-07 1.59e-02 6.82e-01 3.00e+00 8.80e-02 2.04e+01 1.00e+00
19 1.55e+02 1.76e+01 2.74e+03 8.09e-07 1.52e-02 6.53e-01 3.00e+00 8.99e-02 1.99e+01 1.00e+00
20 1.42e+02 1.79e+01 2.49e+03 8.08e-07 1.45e-02 6.28e-01 3.00e+00 9.73e-02 1.69e+01 1.00e+00
21 1.32e+02 1.82e+01 2.27e+03 8.05e-07 1.38e-02 6.06e-01 3.00e+00 9.94e-02 1.60e+01 1.00e+00
22 1.22e+02 1.85e+01 2.07e+03 8.61e-07 1.33e-02 5.83e-01 3.00e+00 1.04e-01 1.46e+01 1.00e+00
23 1.14e+02 1.88e+01 1.89e+03 8.19e-07 1.28e-02 5.61e-01 3.00e+00 1.07e-01 1.37e+01 1.00e+00
24 1.06e+02 1.90e+01 1.74e+03 8.90e-07 1.23e-02 5.41e-01 3.00e+00 9.75e-02 1.27e+01 1.00e+00
25 9.92e+01 1.91e+01 1.60e+03 6.97e-07 1.19e-02 5.23e-01 3.00e+00 8.84e-02 1.17e+01 1.00e+00
26 9.34e+01 1.93e+01 1.48e+03 8.29e-07 1.14e-02 5.09e-01 3.00e+00 8.19e-02 1.06e+01 1.00e+00
27 8.81e+01 1.94e+01 1.38e+03 8.52e-07 1.11e-02 4.96e-01 3.00e+00 7.67e-02 9.64e+00 1.00e+00
28 8.34e+01 1.93e+01 1.28e+03 8.11e-07 1.09e-02 4.85e-01 3.00e+00 7.33e-02 8.58e+00 1.00e+00
29 7.90e+01 1.91e+01 1.20e+03 5.46e-07 1.07e-02 4.80e-01 3.00e+00 7.12e-02 7.30e+00 1.00e+00
30 7.51e+01 1.90e+01 1.12e+03 7.39e-07 1.08e-02 4.83e-01 3.00e+00 7.01e-02 5.72e+00 1.00e+00
31 7.22e+01 1.92e+01 1.06e+03 6.51e-07 1.13e-02 5.00e-01 3.00e+00 7.04e-02 3.39e+00 1.00e+00
32 7.15e+01 2.11e+01 1.01e+03 7.80e-07 1.28e-02 5.43e-01 3.00e+00 7.22e-02 1.74e+00 1.00e+00
33 7.75e+01 2.82e+01 9.85e+02 7.01e-07 1.66e-02 6.29e-01 3.00e+00 7.78e-02 3.27e+00 1.00e+00
34 9.15e+01 3.84e+01 1.06e+03 7.55e-07 2.70e-02 7.65e-01 3.00e+00 1.06e-01 4.20e+00 1.00e+00
35 1.01e+02 4.59e+01 1.11e+03 8.35e-07 4.51e-02 8.65e-01 3.00e+00 2.13e-01 4.91e+00 1.00e+00
36 8.47e+01 4.52e+01 7.89e+02 7.21e-07 5.24e-02 8.33e-01 3.00e+00 3.73e-01 3.70e+00 1.00e+00
37 6.66e+01 3.66e+01 6.01e+02 7.10e-07 4.93e-02 7.76e-01 3.00e+00 5.27e-01 1.81e+00 1.00e+00
38 5.78e+01 2.79e+01 5.97e+02 8.00e-07 4.22e-02 6.92e-01 3.00e+00 5.44e-01 1.20e+00 1.00e+00
39 5.26e+01 2.10e+01 6.31e+02 7.79e-07 3.48e-02 6.01e-01 3.00e+00 4.45e-01 1.02e+00 1.00e+00
40 4.91e+01 1.81e+01 6.20e+02 5.99e-07 2.89e-02 5.34e-01 3.00e+00 3.64e-01 9.01e-01 1.00e+00
41 4.69e+01 1.81e+01 5.75e+02 6.71e-07 2.52e-02 4.98e-01 3.00e+00 3.04e-01 9.00e-01 1.00e+00
42 4.55e+01 1.84e+01 5.42e+02 7.37e-07 2.26e-02 4.78e-01 3.00e+00 2.58e-01 9.19e-01 1.00e+00
43 4.46e+01 1.84e+01 5.24e+02 7.26e-07 2.04e-02 4.55e-01 3.00e+00 2.14e-01 9.21e-01 1.00e+00
44 4.33e+01 1.76e+01 5.15e+02 8.69e-07 1.81e-02 4.25e-01 3.00e+00 1.74e-01 8.88e-01 1.00e+00
45 4.22e+01 1.67e+01 5.11e+02 7.09e-07 1.64e-02 3.94e-01 3.00e+00 1.55e-01 9.09e-01 1.00e+00
46 4.16e+01 1.59e+01 5.15e+02 9.39e-07 1.50e-02 3.70e-01 3.00e+00 1.41e-01 9.11e-01 1.00e+00
47 4.10e+01 1.50e+01 5.20e+02 7.45e-07 1.34e-02 3.51e-01 3.00e+00 1.18e-01 8.63e-01 1.00e+00
48 4.03e+01 1.44e+01 5.19e+02 8.35e-07 1.22e-02 3.34e-01 3.00e+00 1.10e-01 8.26e-01 1.00e+00
49 3.98e+01 1.40e+01 5.16e+02 8.65e-07 1.12e-02 3.22e-01 3.00e+00 1.04e-01 8.23e-01 1.00e+00
50 3.95e+01 1.40e+01 5.10e+02 7.87e-07 1.06e-02 3.13e-01 3.00e+00 9.82e-02 8.26e-01 1.00e+00
51 3.94e+01 1.42e+01 5.04e+02 7.72e-07 1.00e-02 3.04e-01 3.00e+00 8.97e-02 8.37e-01 1.00e+00
52 3.93e+01 1.43e+01 4.99e+02 1.03e-06 9.57e-03 2.97e-01 3.00e+00 8.44e-02 8.43e-01 1.00e+00
53 3.90e+01 1.42e+01 4.98e+02 9.20e-07 9.06e-03 2.90e-01 3.00e+00 7.97e-02 8.27e-01 1.00e+00
54 3.86e+01 1.36e+01 5.01e+02 6.66e-07 8.44e-03 2.83e-01 3.00e+00 7.79e-02 7.85e-01 1.00e+00
55 3.82e+01 1.28e+01 5.06e+02 6.93e-07 8.00e-03 2.76e-01 3.00e+00 7.59e-02 7.37e-01 1.00e+00
56 3.79e+01 1.24e+01 5.12e+02 8.02e-07 7.77e-03 2.70e-01 3.00e+00 7.38e-02 7.27e-01 1.00e+00
57 3.80e+01 1.23e+01 5.13e+02 8.58e-07 7.79e-03 2.69e-01 3.00e+00 7.10e-02 7.58e-01 1.00e+00
58 3.80e+01 1.26e+01 5.09e+02 6.90e-07 7.79e-03 2.70e-01 3.00e+00 6.73e-02 7.88e-01 1.00e+00
59 3.79e+01 1.28e+01 5.03e+02 6.89e-07 7.52e-03 2.68e-01 3.00e+00 6.48e-02 7.98e-01 1.00e+00
60 3.78e+01 1.29e+01 4.98e+02 5.58e-07 7.24e-03 2.67e-01 3.00e+00 6.36e-02 7.83e-01 1.00e+00
61 3.75e+01 1.27e+01 4.95e+02 7.50e-07 6.99e-03 2.68e-01 3.00e+00 6.21e-02 7.45e-01 1.00e+00
62 3.72e+01 1.25e+01 4.93e+02 7.85e-07 6.80e-03 2.68e-01 3.00e+00 6.01e-02 7.01e-01 1.00e+00
63 3.70e+01 1.24e+01 4.91e+02 5.69e-07 6.67e-03 2.68e-01 3.00e+00 5.78e-02 6.82e-01 1.00e+00
64 3.68e+01 1.24e+01 4.89e+02 6.56e-07 6.60e-03 2.65e-01 3.00e+00 5.50e-02 6.98e-01 1.00e+00
65 3.68e+01 1.25e+01 4.86e+02 7.75e-07 6.49e-03 2.59e-01 3.00e+00 5.24e-02 7.22e-01 1.00e+00
66 3.68e+01 1.25e+01 4.85e+02 8.27e-07 6.38e-03 2.52e-01 3.00e+00 5.06e-02 7.32e-01 1.00e+00
67 3.67e+01 1.24e+01 4.86e+02 8.32e-07 6.23e-03 2.48e-01 3.00e+00 4.98e-02 7.28e-01 1.00e+00
68 3.66e+01 1.22e+01 4.88e+02 5.96e-07 5.98e-03 2.47e-01 3.00e+00 4.99e-02 7.12e-01 1.00e+00
69 3.64e+01 1.20e+01 4.89e+02 7.20e-07 5.77e-03 2.47e-01 3.00e+00 5.04e-02 6.90e-01 1.00e+00
70 3.63e+01 1.18e+01 4.89e+02 7.26e-07 5.57e-03 2.48e-01 3.00e+00 5.11e-02 6.69e-01 1.00e+00
71 3.61e+01 1.17e+01 4.88e+02 5.66e-07 5.48e-03 2.47e-01 3.00e+00 5.09e-02 6.56e-01 1.00e+00
72 3.60e+01 1.17e+01 4.87e+02 7.61e-07 5.45e-03 2.43e-01 3.00e+00 4.96e-02 6.55e-01 1.00e+00
73 3.60e+01 1.17e+01 4.86e+02 7.10e-07 5.45e-03 2.38e-01 3.00e+00 4.83e-02 6.62e-01 1.00e+00
74 3.60e+01 1.18e+01 4.84e+02 7.35e-07 5.39e-03 2.34e-01 3.00e+00 4.71e-02 6.67e-01 1.00e+00
75 3.60e+01 1.18e+01 4.83e+02 7.18e-07 5.28e-03 2.31e-01 3.00e+00 4.61e-02 6.66e-01 1.00e+00
76 3.59e+01 1.18e+01 4.82e+02 7.95e-07 5.16e-03 2.31e-01 3.00e+00 4.53e-02 6.61e-01 1.00e+00
77 3.58e+01 1.17e+01 4.82e+02 7.86e-07 5.04e-03 2.33e-01 3.00e+00 4.47e-02 6.52e-01 1.00e+00
78 3.57e+01 1.17e+01 4.81e+02 7.11e-07 4.92e-03 2.34e-01 3.00e+00 4.35e-02 6.41e-01 1.00e+00
79 3.56e+01 1.15e+01 4.81e+02 8.21e-07 4.82e-03 2.34e-01 3.00e+00 4.17e-02 6.28e-01 1.00e+00
80 3.55e+01 1.14e+01 4.82e+02 7.76e-07 4.77e-03 2.30e-01 3.00e+00 4.01e-02 6.15e-01 1.00e+00
81 3.55e+01 1.14e+01 4.81e+02 7.41e-07 4.74e-03 2.23e-01 3.00e+00 3.89e-02 6.03e-01 1.00e+00
82 3.54e+01 1.14e+01 4.80e+02 5.68e-07 4.72e-03 2.18e-01 3.00e+00 3.79e-02 5.96e-01 1.00e+00
83 3.54e+01 1.14e+01 4.79e+02 6.44e-07 4.68e-03 2.16e-01 3.00e+00 3.73e-02 5.95e-01 1.00e+00
84 3.53e+01 1.14e+01 4.78e+02 7.65e-07 4.62e-03 2.18e-01 3.00e+00 3.69e-02 5.98e-01 1.00e+00
85 3.53e+01 1.14e+01 4.77e+02 1.12e-06 4.54e-03 2.22e-01 3.00e+00 3.67e-02 6.05e-01 1.00e+00
86 3.52e+01 1.14e+01 4.76e+02 7.22e-07 4.44e-03 2.23e-01 3.00e+00 3.67e-02 6.10e-01 1.00e+00
87 3.52e+01 1.14e+01 4.75e+02 7.18e-07 4.36e-03 2.21e-01 3.00e+00 3.68e-02 6.10e-01 1.00e+00
88 3.51e+01 1.14e+01 4.74e+02 8.42e-07 4.27e-03 2.17e-01 3.00e+00 3.71e-02 6.02e-01 1.00e+00
89 3.51e+01 1.14e+01 4.74e+02 6.53e-07 4.20e-03 2.12e-01 3.00e+00 3.67e-02 5.87e-01 1.00e+00
90 3.50e+01 1.13e+01 4.73e+02 7.39e-07 4.16e-03 2.08e-01 3.00e+00 3.64e-02 5.66e-01 1.00e+00
91 3.49e+01 1.13e+01 4.73e+02 8.80e-07 4.12e-03 2.04e-01 3.00e+00 3.62e-02 5.47e-01 1.00e+00
92 3.49e+01 1.13e+01 4.72e+02 7.64e-07 4.08e-03 2.02e-01 3.00e+00 3.58e-02 5.36e-01 1.00e+00
93 3.48e+01 1.13e+01 4.72e+02 6.88e-07 4.04e-03 2.01e-01 3.00e+00 3.52e-02 5.36e-01 1.00e+00
94 3.48e+01 1.13e+01 4.71e+02 8.15e-07 3.99e-03 2.01e-01 3.00e+00 3.46e-02 5.45e-01 1.00e+00
95 3.48e+01 1.12e+01 4.71e+02 7.85e-07 3.94e-03 2.01e-01 3.00e+00 3.40e-02 5.55e-01 1.00e+00
96 3.48e+01 1.12e+01 4.72e+02 8.67e-07 3.88e-03 2.01e-01 3.00e+00 3.33e-02 5.63e-01 1.00e+00
97 3.48e+01 1.11e+01 4.72e+02 6.86e-07 3.82e-03 2.00e-01 3.00e+00 3.24e-02 5.63e-01 1.00e+00
98 3.47e+01 1.11e+01 4.73e+02 8.18e-07 3.75e-03 1.97e-01 3.00e+00 3.14e-02 5.56e-01 1.00e+00
99 3.47e+01 1.10e+01 4.73e+02 7.21e-07 3.69e-03 1.94e-01 3.00e+00 3.07e-02 5.43e-01 1.00e+00
100 3.46e+01 1.10e+01 4.73e+02 9.35e-07 3.63e-03 1.90e-01 3.00e+00 3.01e-02 5.26e-01 1.00e+00
101 3.46e+01 1.09e+01 4.73e+02 6.68e-07 3.59e-03 1.87e-01 3.00e+00 2.97e-02 5.11e-01 1.00e+00
102 3.45e+01 1.09e+01 4.71e+02 9.50e-07 3.55e-03 1.84e-01 3.00e+00 2.94e-02 5.01e-01 1.00e+00
103 3.45e+01 1.10e+01 4.70e+02 6.84e-07 3.53e-03 1.83e-01 3.00e+00 2.91e-02 4.98e-01 1.00e+00
104 3.44e+01 1.10e+01 4.68e+02 6.39e-07 3.49e-03 1.83e-01 3.00e+00 2.90e-02 5.00e-01 1.00e+00
105 3.44e+01 1.11e+01 4.67e+02 7.51e-07 3.46e-03 1.84e-01 3.00e+00 2.88e-02 5.05e-01 1.00e+00
106 3.44e+01 1.11e+01 4.66e+02 8.45e-07 3.41e-03 1.84e-01 3.00e+00 2.88e-02 5.09e-01 1.00e+00
107 3.44e+01 1.11e+01 4.66e+02 8.90e-07 3.36e-03 1.83e-01 3.00e+00 2.88e-02 5.10e-01 1.00e+00
108 3.44e+01 1.11e+01 4.66e+02 7.98e-07 3.30e-03 1.80e-01 3.00e+00 2.88e-02 5.08e-01 1.00e+00
109 3.43e+01 1.10e+01 4.66e+02 8.66e-07 3.24e-03 1.76e-01 3.00e+00 2.86e-02 5.02e-01 1.00e+00
110 3.43e+01 1.10e+01 4.66e+02 6.33e-07 3.18e-03 1.72e-01 3.00e+00 2.83e-02 4.96e-01 1.00e+00
111 3.43e+01 1.10e+01 4.66e+02 9.04e-07 3.12e-03 1.69e-01 3.00e+00 2.81e-02 4.88e-01 1.00e+00
112 3.42e+01 1.09e+01 4.66e+02 8.30e-07 3.08e-03 1.66e-01 3.00e+00 2.78e-02 4.81e-01 1.00e+00
113 3.42e+01 1.09e+01 4.66e+02 6.63e-07 3.05e-03 1.65e-01 3.00e+00 2.75e-02 4.75e-01 1.00e+00
114 3.42e+01 1.09e+01 4.66e+02 7.47e-07 3.03e-03 1.64e-01 3.00e+00 2.73e-02 4.68e-01 1.00e+00
115 3.41e+01 1.08e+01 4.66e+02 7.00e-07 3.00e-03 1.63e-01 3.00e+00 2.70e-02 4.62e-01 1.00e+00
116 3.41e+01 1.08e+01 4.66e+02 7.65e-07 2.97e-03 1.63e-01 3.00e+00 2.68e-02 4.57e-01 1.00e+00
117 3.41e+01 1.08e+01 4.66e+02 6.93e-07 2.94e-03 1.62e-01 3.00e+00 2.65e-02 4.55e-01 1.00e+00
118 3.41e+01 1.08e+01 4.65e+02 7.52e-07 2.91e-03 1.60e-01 3.00e+00 2.61e-02 4.55e-01 1.00e+00
119 3.41e+01 1.08e+01 4.65e+02 6.75e-07 2.87e-03 1.59e-01 3.00e+00 2.56e-02 4.57e-01 1.00e+00
120 3.41e+01 1.08e+01 4.65e+02 8.27e-07 2.83e-03 1.56e-01 3.00e+00 2.52e-02 4.61e-01 1.00e+00
121 3.40e+01 1.08e+01 4.64e+02 8.18e-07 2.79e-03 1.54e-01 3.00e+00 2.48e-02 4.63e-01 1.00e+00
122 3.40e+01 1.08e+01 4.64e+02 7.49e-07 2.76e-03 1.52e-01 3.00e+00 2.45e-02 4.62e-01 1.00e+00
123 3.40e+01 1.08e+01 4.64e+02 7.50e-07 2.73e-03 1.51e-01 3.00e+00 2.42e-02 4.57e-01 1.00e+00
124 3.40e+01 1.08e+01 4.63e+02 8.03e-07 2.69e-03 1.50e-01 3.00e+00 2.40e-02 4.48e-01 1.00e+00
125 3.39e+01 1.08e+01 4.63e+02 9.57e-07 2.66e-03 1.50e-01 3.00e+00 2.39e-02 4.36e-01 1.00e+00
126 3.39e+01 1.08e+01 4.63e+02 7.06e-07 2.62e-03 1.48e-01 3.00e+00 2.38e-02 4.25e-01 1.00e+00
127 3.39e+01 1.08e+01 4.63e+02 7.90e-07 2.60e-03 1.47e-01 3.00e+00 2.37e-02 4.17e-01 1.00e+00
128 3.39e+01 1.08e+01 4.62e+02 6.06e-07 2.57e-03 1.44e-01 3.00e+00 2.35e-02 4.13e-01 1.00e+00
129 3.39e+01 1.07e+01 4.62e+02 7.10e-07 2.56e-03 1.42e-01 3.00e+00 2.33e-02 4.15e-01 1.00e+00
130 3.38e+01 1.07e+01 4.62e+02 6.87e-07 2.54e-03 1.40e-01 3.00e+00 2.30e-02 4.20e-01 1.00e+00
131 3.38e+01 1.07e+01 4.62e+02 7.03e-07 2.51e-03 1.39e-01 3.00e+00 2.28e-02 4.25e-01 1.00e+00
132 3.38e+01 1.07e+01 4.62e+02 5.88e-07 2.49e-03 1.39e-01 3.00e+00 2.25e-02 4.28e-01 1.00e+00
133 3.38e+01 1.07e+01 4.62e+02 8.58e-07 2.46e-03 1.39e-01 3.00e+00 2.23e-02 4.28e-01 1.00e+00
134 3.38e+01 1.07e+01 4.62e+02 7.04e-07 2.43e-03 1.38e-01 3.00e+00 2.22e-02 4.24e-01 1.00e+00
135 3.38e+01 1.07e+01 4.62e+02 8.68e-07 2.40e-03 1.37e-01 3.00e+00 2.21e-02 4.17e-01 1.00e+00
136 3.38e+01 1.07e+01 4.61e+02 6.46e-07 2.38e-03 1.36e-01 3.00e+00 2.21e-02 4.09e-01 1.00e+00
137 3.37e+01 1.07e+01 4.61e+02 6.63e-07 2.35e-03 1.34e-01 3.00e+00 2.20e-02 4.00e-01 1.00e+00
138 3.37e+01 1.07e+01 4.61e+02 7.96e-07 2.32e-03 1.32e-01 3.00e+00 2.19e-02 3.94e-01 1.00e+00
139 3.37e+01 1.07e+01 4.61e+02 7.96e-07 2.29e-03 1.31e-01 3.00e+00 2.18e-02 3.90e-01 1.00e+00
140 3.37e+01 1.07e+01 4.61e+02 6.15e-07 2.27e-03 1.29e-01 3.00e+00 2.16e-02 3.88e-01 1.00e+00
141 3.37e+01 1.07e+01 4.60e+02 7.20e-07 2.25e-03 1.28e-01 3.00e+00 2.14e-02 3.88e-01 1.00e+00
142 3.37e+01 1.06e+01 4.60e+02 8.81e-07 2.23e-03 1.27e-01 3.00e+00 2.12e-02 3.89e-01 1.00e+00
143 3.37e+01 1.06e+01 4.60e+02 7.09e-07 2.21e-03 1.26e-01 3.00e+00 2.10e-02 3.89e-01 1.00e+00
144 3.36e+01 1.06e+01 4.60e+02 7.77e-07 2.19e-03 1.25e-01 3.00e+00 2.08e-02 3.89e-01 1.00e+00
145 3.36e+01 1.06e+01 4.60e+02 6.66e-07 2.17e-03 1.24e-01 3.00e+00 2.07e-02 3.87e-01 1.00e+00
146 3.36e+01 1.06e+01 4.60e+02 7.74e-07 2.15e-03 1.24e-01 3.00e+00 2.05e-02 3.85e-01 1.00e+00
147 3.36e+01 1.06e+01 4.60e+02 9.70e-07 2.13e-03 1.23e-01 3.00e+00 2.03e-02 3.82e-01 1.00e+00
148 3.36e+01 1.06e+01 4.59e+02 8.02e-07 2.11e-03 1.22e-01 3.00e+00 2.00e-02 3.79e-01 1.00e+00
149 3.36e+01 1.06e+01 4.59e+02 7.37e-07 2.09e-03 1.20e-01 3.00e+00 1.97e-02 3.76e-01 1.00e+00
150 3.36e+01 1.06e+01 4.59e+02 8.97e-07 2.07e-03 1.18e-01 3.00e+00 1.95e-02 3.73e-01 1.00e+00
151 3.36e+01 1.06e+01 4.59e+02 8.56e-07 2.06e-03 1.17e-01 3.00e+00 1.92e-02 3.69e-01 1.00e+00
152 3.36e+01 1.06e+01 4.59e+02 6.31e-07 2.04e-03 1.15e-01 3.00e+00 1.90e-02 3.65e-01 1.00e+00
153 3.35e+01 1.06e+01 4.59e+02 6.84e-07 2.02e-03 1.15e-01 3.00e+00 1.89e-02 3.61e-01 1.00e+00
154 3.35e+01 1.06e+01 4.59e+02 7.95e-07 2.01e-03 1.14e-01 3.00e+00 1.87e-02 3.57e-01 1.00e+00
155 3.35e+01 1.06e+01 4.59e+02 7.86e-07 1.99e-03 1.15e-01 3.00e+00 1.86e-02 3.55e-01 1.00e+00
156 3.35e+01 1.06e+01 4.59e+02 7.42e-07 1.98e-03 1.15e-01 3.00e+00 1.85e-02 3.54e-01 1.00e+00
157 3.35e+01 1.06e+01 4.58e+02 6.67e-07 1.96e-03 1.14e-01 3.00e+00 1.84e-02 3.54e-01 1.00e+00
158 3.35e+01 1.06e+01 4.58e+02 6.47e-07 1.94e-03 1.12e-01 3.00e+00 1.83e-02 3.55e-01 1.00e+00
159 3.35e+01 1.06e+01 4.58e+02 7.59e-07 1.92e-03 1.10e-01 3.00e+00 1.83e-02 3.55e-01 1.00e+00
160 3.35e+01 1.06e+01 4.58e+02 6.58e-07 1.91e-03 1.09e-01 3.00e+00 1.82e-02 3.55e-01 1.00e+00
161 3.35e+01 1.06e+01 4.58e+02 6.79e-07 1.89e-03 1.08e-01 3.00e+00 1.82e-02 3.53e-01 1.00e+00
162 3.35e+01 1.06e+01 4.58e+02 8.56e-07 1.87e-03 1.07e-01 3.00e+00 1.81e-02 3.50e-01 1.00e+00
163 3.34e+01 1.06e+01 4.58e+02 9.02e-07 1.85e-03 1.07e-01 3.00e+00 1.80e-02 3.46e-01 1.00e+00
164 3.34e+01 1.05e+01 4.58e+02 9.41e-07 1.83e-03 1.07e-01 3.00e+00 1.79e-02 3.41e-01 1.00e+00
165 3.34e+01 1.05e+01 4.58e+02 6.28e-07 1.81e-03 1.06e-01 3.00e+00 1.78e-02 3.36e-01 1.00e+00
166 3.34e+01 1.05e+01 4.58e+02 7.22e-07 1.79e-03 1.05e-01 3.00e+00 1.77e-02 3.32e-01 1.00e+00
167 3.34e+01 1.05e+01 4.58e+02 9.56e-07 1.78e-03 1.03e-01 3.00e+00 1.76e-02 3.29e-01 1.00e+00
168 3.34e+01 1.05e+01 4.58e+02 7.84e-07 1.77e-03 1.01e-01 3.00e+00 1.75e-02 3.28e-01 1.00e+00
169 3.34e+01 1.05e+01 4.58e+02 6.52e-07 1.76e-03 1.00e-01 3.00e+00 1.74e-02 3.28e-01 1.00e+00
170 3.34e+01 1.05e+01 4.57e+02 7.33e-07 1.74e-03 9.97e-02 3.00e+00 1.72e-02 3.29e-01 1.00e+00
171 3.34e+01 1.05e+01 4.57e+02 8.15e-07 1.73e-03 9.99e-02 3.00e+00 1.70e-02 3.29e-01 1.00e+00
172 3.34e+01 1.05e+01 4.57e+02 7.00e-07 1.71e-03 1.00e-01 3.00e+00 1.68e-02 3.29e-01 1.00e+00
173 3.34e+01 1.05e+01 4.57e+02 8.28e-07 1.69e-03 9.99e-02 3.00e+00 1.66e-02 3.28e-01 1.00e+00
174 3.34e+01 1.05e+01 4.57e+02 5.61e-07 1.68e-03 9.90e-02 3.00e+00 1.64e-02 3.25e-01 1.00e+00
175 3.33e+01 1.05e+01 4.57e+02 9.26e-07 1.66e-03 9.73e-02 3.00e+00 1.63e-02 3.23e-01 1.00e+00
176 3.33e+01 1.05e+01 4.57e+02 9.06e-07 1.65e-03 9.53e-02 3.00e+00 1.61e-02 3.19e-01 1.00e+00
177 3.33e+01 1.05e+01 4.57e+02 1.07e-06 1.64e-03 9.38e-02 3.00e+00 1.60e-02 3.16e-01 1.00e+00
178 3.33e+01 1.05e+01 4.57e+02 7.71e-07 1.63e-03 9.32e-02 3.00e+00 1.59e-02 3.13e-01 1.00e+00
179 3.33e+01 1.05e+01 4.57e+02 9.56e-07 1.62e-03 9.33e-02 3.00e+00 1.58e-02 3.11e-01 1.00e+00
180 3.33e+01 1.05e+01 4.56e+02 7.34e-07 1.61e-03 9.39e-02 3.00e+00 1.57e-02 3.08e-01 1.00e+00
181 3.33e+01 1.05e+01 4.56e+02 8.88e-07 1.59e-03 9.41e-02 3.00e+00 1.57e-02 3.06e-01 1.00e+00
182 3.33e+01 1.05e+01 4.56e+02 6.57e-07 1.56e-03 9.36e-02 3.00e+00 1.56e-02 3.05e-01 1.00e+00
183 3.33e+01 1.05e+01 4.56e+02 8.26e-07 1.54e-03 9.21e-02 3.00e+00 1.56e-02 3.04e-01 1.00e+00
184 3.33e+01 1.05e+01 4.56e+02 8.69e-07 1.53e-03 9.01e-02 3.00e+00 1.55e-02 3.03e-01 1.00e+00
185 3.33e+01 1.05e+01 4.56e+02 9.41e-07 1.52e-03 8.80e-02 3.00e+00 1.55e-02 3.02e-01 1.00e+00
186 3.33e+01 1.05e+01 4.56e+02 8.06e-07 1.51e-03 8.66e-02 3.00e+00 1.54e-02 3.02e-01 1.00e+00
187 3.33e+01 1.05e+01 4.56e+02 7.64e-07 1.51e-03 8.61e-02 3.00e+00 1.53e-02 3.02e-01 1.00e+00
188 3.33e+01 1.04e+01 4.56e+02 6.35e-07 1.50e-03 8.65e-02 3.00e+00 1.52e-02 3.01e-01 1.00e+00
189 3.33e+01 1.04e+01 4.56e+02 8.06e-07 1.49e-03 8.73e-02 3.00e+00 1.51e-02 3.00e-01 1.00e+00
190 3.32e+01 1.04e+01 4.56e+02 6.73e-07 1.48e-03 8.79e-02 3.00e+00 1.50e-02 2.97e-01 1.00e+00
191 3.32e+01 1.04e+01 4.56e+02 7.96e-07 1.46e-03 8.77e-02 3.00e+00 1.49e-02 2.94e-01 1.00e+00
192 3.32e+01 1.05e+01 4.56e+02 6.80e-07 1.44e-03 8.67e-02 3.00e+00 1.49e-02 2.90e-01 1.00e+00
193 3.32e+01 1.04e+01 4.56e+02 6.28e-07 1.42e-03 8.50e-02 3.00e+00 1.48e-02 2.86e-01 1.00e+00
194 3.32e+01 1.04e+01 4.56e+02 6.54e-07 1.41e-03 8.32e-02 3.00e+00 1.47e-02 2.83e-01 1.00e+00
195 3.32e+01 1.04e+01 4.56e+02 7.27e-07 1.40e-03 8.16e-02 3.00e+00 1.46e-02 2.81e-01 1.00e+00
196 3.32e+01 1.04e+01 4.56e+02 7.41e-07 1.39e-03 8.06e-02 3.00e+00 1.45e-02 2.81e-01 1.00e+00
197 3.32e+01 1.04e+01 4.56e+02 9.70e-07 1.38e-03 8.02e-02 3.00e+00 1.44e-02 2.81e-01 1.00e+00
198 3.32e+01 1.04e+01 4.56e+02 6.79e-07 1.38e-03 8.02e-02 3.00e+00 1.43e-02 2.82e-01 1.00e+00
199 3.32e+01 1.04e+01 4.55e+02 8.67e-07 1.37e-03 8.03e-02 3.00e+00 1.41e-02 2.83e-01 1.00e+00
--------------------------------------------------------------------------------------------------------
Reconstruct from the CDL solution with a spatial mask.
sr2 = d2.reconstruct().squeeze() + sl
Compare dictionaries.
fig = plot.figure(figsize=(14, 7))
plot.subplot(1, 2, 1)
plot.imview(util.tiledict(D1.squeeze()), fig=fig,
title='Without Mask Decoupling')
plot.subplot(1, 2, 2)
plot.imview(util.tiledict(D2.squeeze()), fig=fig,
title='With Mask Decoupling')
fig.show()
Display reference and training images.
fig = plot.figure(figsize=(14, 14))
plot.subplot(2, 2, 1)
plot.imview(S[...,0], title='Reference', fig=fig)
plot.subplot(2, 2, 2)
plot.imview(Sw[...,0], title='Test', fig=fig)
plot.subplot(2, 2, 3)
plot.imview(S[...,1], title='Reference', fig=fig)
plot.subplot(2, 2, 4)
plot.imview(Sw[...,1], title='Test', fig=fig)
fig.show()
Compare reconstructed images.
fig = plot.figure(figsize=(14, 14))
plot.subplot(2, 2, 1)
plot.imview(sr1[...,0], title='Without Mask Decoupling', fig=fig)
plot.subplot(2, 2, 2)
plot.imview(sr2[...,0], title='With Mask Decoupling', fig=fig)
plot.subplot(2, 2, 3)
plot.imview(sr1[...,1], title='Without Mask Decoupling', fig=fig)
plot.subplot(2, 2, 4)
plot.imview(sr2[...,1], title='With Mask Decoupling', fig=fig)
fig.show()