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()