\section{Introduction}
Theoretical as well as experimental studies have shown that
chiral sculptured thin films (STFs) display the so--called Bragg
phenomenon on axial excitation,
in consequence of their
periodically and unidirectionally nonhomogeneous constitution
\c{VLbook, WHL, HWKLR}. Let the
direction of nonhomogeneity of a chiral STF be parallel to the $z$ axis, while
the film completely occupies the region $0 \leq z \leq L$. When circularly
polarized, mono\-chromatic light falls normally on this film of sufficient
thickness, then it is
(i) almost perfectly reflected if the handedness of the incident
light coincides with the structural handedness of the film, and
(ii) almost perfectly transmitted if otherwise~---~provided
absorption within the film is negligible
and the free--space
wavelength $\lambda_0$ of the incident light lies within the so--called Bragg
regime.  Optical devices such as circular--polarization filters and spectral--hole filters
exploiting the
circular Bragg phenomenon have been fabricated using STF technology,
delivering performance
comparable to commercially available non--STF analogs \c{WHL,HWLM,HWTLM}.

The bandwidth of the
Bragg regime is directly proportional to the local birefringence of the
chiral STF. Although the serial bideposition technique
\c{HWKLR} yields highly birefringent chiral STFs, further enhancement
of the bandwidth can be desirable for certain applications. The objective of this communication
is to show that stepwise chirping of the periodicity of a chiral STF
can easily deliver a considerable enhancement of the bandwidth.

The stepwise chirped chiral STF is described in Section 2, and its
response to  normally incident plane waves is presented in Section 3.
Comparison is made with the response of an unchirped chiral STF.
An $\exp(-i\omega t)$ time--dependence
is implicit in this communication, with $\omega$ as the
angular frequency; vectors are underlined and dyadics are double--underlined; $\#r=x\#{u}_x+y\#{u}_y+z\#{u}_z$ is the
position vector with $\#{u}_x$, $\#{u}_y$ and $\#{u}_z$ as the cartesian
unit vectors; and  $\lambda_0$ is the free--space wavelength.

\section{Constitutive Relations}
The nonhomogeneous permittivity
dyadic of a chiral STF is expressed as follows \c{VLbook}:
\begin{equation}
\label{epsbasic}
\=\epsilon(\#r) = \epsilon_0\, \=S_z(z,h)\.\=S_y(\chi)\.
\Big[ \epsilon_a \,\#{u}_z\uz +\epsilon_b\,\#{u}_x\ux
\,  +\,\epsilon_c\,\#{u}_y\uy\Big]\.\=S_y^{-1}(\chi)\.
\=S_z^{-1}(z,h)\, , \quad 0 \leq z \leq L\, .
\end{equation}
In this equation, $\epsilon_{a,b,c}$ are three relative permittivity
scalars; $\epsilon_0 = 8.854\times 10^{-12}$~F~m$^{-1}$
is the free--space permittivity; the  tilt dyadic
\begin{equation}
\=S_y(\chi) = \#{u}_y\uy + (\#{u}_x\ux + \#{u}_z\uz) \, \cos\chi + (\#{u}_z\#{u}_x-\#{u}_x\#{u}_z)\,
\sin\chi
\end{equation}
represents the locally columnar microstructure of any chiral STF
with $\chi > 0^\circ$; while the rotation dyadic
\begin{equation}
\=S_z(z,h)=
\#{u}_z\uz + \left( \#{u}_x\ux+\#{u}_y\uy\right)\,\cos\left(\frac{\pi z}{\Omega}\right) \,
+ h\, \left( \#{u}_y\#{u}_x-\#{u}_x\#{u}_y\right)\,\sin\left( \frac{\pi z}{\Omega}\right)
\end{equation}
contains  $2\Omega$ as the structural period, with $h=1$ for
structural right--handedness and $h=-1$ for structural left--handedness.
Taking a whole number of periods to encompass the film thickness,
we choose the ratio $L/2\Omega$ to be a non--prime integer for filtering applications.

The stepwise chirped chiral STF comprises $N_p$ steps, each step having
$\ell_p$ complete periods, with $N_p \ell_p = L/2\Omega$.
The half--period of the n$^{th}$ step
is denoted by
\begin{equation}
\Omega_n = \Omega + \frac{2n-N_p-1}{2}\, \delta_\Omega \,,
\quad n \in [1,\,N_p]\,,
\end{equation}
wherein the chirping step $\delta_\Omega$ is selected to be a small fraction
of $\Omega$ to assure that  $\Omega_1 > 0$.
Thus, the expression of the
permittivity dyadic of the stepwise chirped chiral STF is the
same as \r{epsbasic}, except that the dyadic $\=S_z(z,h)$ has to be replaced
by $\=S_z^{chirp}(z,h)$, where
\begin{eqnarray}
\=S_z^{chirp}(z,h)
&=&
\#{u}_z\uz + \left( \#{u}_x\ux+\#{u}_y\uy\right)\,\cos\left[\frac{\pi (z-\zeta_n)}{\Omega_n}\right] \,
\nonumber
\\
&+& h\,\left( \#{u}_y\#{u}_x-\#{u}_x\#{u}_y\right)\,\sin\left[ \frac{\pi (z-\zeta_n)}{\Omega_n}\right]\,,
\quad \zeta_n \leq z \leq \zeta_{n+1}\,,\quad n \in [1,\,N_p]\,,
\end{eqnarray}
and
\begin{equation}
\zeta_n = 2\ell_p (n-1)\left( \Omega +\frac{n-N_p-1}{2}\,\delta_\Omega\right)\,
\end{equation}
is the total thickness of the first $n-1$ steps.
Note that the thickness $\zeta_{{N_p}+1}$ of the chirped chiral STF equals
$L$.

\section{Numerical Results and Discussion}
Let an arbitrarily polarized plane wave be normally incident
from the lower half--space $z \leq 0$  on either
of the two films described in the previous section. As a
result, a plane wave is reflected into the lower half--space and another
is transmitted into the upper half--space $z \geq L$.
he electric field phasors associated with the two plane
waves in the lower half--space are stated as
\begin{equation}
\#{E}_{inc}(\#{r},\lambda_0) =
 \left(  a_L \, \#{u}_+ + a_R \, \#{u}_- \right) \,
\exp\left( i k_0 z \right)  \, ;  \,\, z \leq 0\, ,
\end{equation}
\begin{equation}
\#{E}_{ref}(\#{r},\lambda_0) =
\left(  r_L \, \#{u}_- + r_R \, \#{u}_+ \right) \,
\exp\left( -i k_0 z \right) \, ;  \,\, z \leq 0\, ,
\end{equation}
where $\#u_\pm = (\#{u}_x \pm i \#{u}_y)/\sqrt{2}$ and
$k_0 = 2\pi/\lambda_0$ is the
free--space wavenumber. Likewise,
the electric field phasor in the upper half--space is
represented as
\begin{equation}
\#{E}_{trs}(\#{r},\lambda_0) =
 \left(  t_L \, \#{u}_+ + t_R \, \#{u}_- \right) \,
\exp\left[ i k_0 (z -L)\right]  \, ;  \,\, z \geq L \, .
\end{equation}
Here, $a_L$ and $a_R$ are the known amplitudes
of the left-- and the right--circularly polarized (LCP \& RCP)
components
of the incident plane wave;
$r_L$ and $r_R$ are the unknown amplitudes
of the reflected plane wave components; while $t_L$ and $t_R$
are the unknown amplitudes
of the transmitted plane wave components.

The reflection amplitudes $r_{L,R}$ and the transmission
amplitudes $t_{L,R}$  can be computed for specified incident amplitudes
($a_L$ and $a_R$) by following standard procedures
described elsewhere; see, e.g., \c{VLbook}.
Reflection and transmission coefficients are thus defined as the entries
in the 2$\times$2 matrixes in the following two relations:
\begin{eqnarray}
\label{eq15}
\left[ \begin{array}{cccc} r_L \\ r_R  \end{array}\right]  &=&
\left[ \begin{array}{cccc} r_{LL} & r_{LR} \\ r_{RL} & r_{RR}\end{array}\right] \,
\left[ \begin{array}{cccc} a_L \\ a_R  \end{array}\right] \, , \\
\label{eq16}
\left[ \begin{array}{cccc} t_L \\ t_R  \end{array}\right]  &=&
\left[ \begin{array}{cccc} t_{LL} & t_{LR} \\ t_{RL} & t_{RR}\end{array}\right] \,
\left[ \begin{array}{cccc} a_L \\ a_R  \end{array}\right]
\, .
\end{eqnarray}
The co--polarized transmission coefficients are denoted by $t_{LL}$ and
$t_{RR}$,
and the cross--polarized ones by $t_{LR}$ and $t_{RL}$; and similarly for the
reflection coefficients in \r{eq15}.

The reflectances $R_{LL} = \vert r_{LL}\vert^2$, etc., and the
transmittances $T_{LL} = \vert t_{LL}\vert^2$, etc., were calculated
as functions of the free--space wavelength $\lambda_0$ for a chiral STF
described by the following parameters: $h=1$, $\epsilon_a = 2.8$, $\epsilon_b = 3.4$,
$\epsilon_c = 2.9$, $\chi = 20^\circ$, $\Omega = 155$~nm, and $L = 126 \,\Omega$.
Graphs are presented in Figure 1, which clearly show the Bragg regime
delineated by $\lambda_0 \in [528,\,564]$~nm almost agrees
with the theoretical range $\lambda_0 \in \left[ 2\Omega\epsilon_c^{1/2},\,
2\Omega\tilde{\epsilon}_d^{1/2}\right]$ \c{VLbook}.
In this regime,
$R_{RR}$ and $T_{LL}$ are very high, while $R_{LL}$ and $T_{RR}$
are very low, because the chosen film is structurally right--handed.

For the stepwise chirped chiral STF, we chose $N_p=21$, $\ell_p = 3$
and $\delta_\Omega = 0.5$~nm. Thus, the smallest half--period $\Omega_1=150$~nm
and the largest half--period $\Omega_2=160$~nm, while the average
half--period is 155~nm (the same as $\Omega$ for Figure 1).
The computed reflectances
and transmittances, plotted in Figure 2 as functions of $\lambda_0$, show
that the Bragg regime now
is $\lambda_0 \in [514,\,580]$~nm.  The lower limit of the Bragg regime
is somewhat higher than $2\Omega_1\epsilon_c^{1/2}$, and the upper limit is
somewhat lower than $2\Omega_{N_p}\tilde{\epsilon}_d^{1/2}$.

Relative
to the 36~nm bandwidth in Figure 1, the Bragg regime in Figure 2
has a bandwidth of 66~nm~---~an enhancement by more than 80\%.
The enhancement comes without any noticeable reduction in the
polarization--sensitivity of the film. Furthermore, calculations show
that mild dispersion and absorption do not affect the enhancement
significantly.

Stepwise changes in the periodicity of a chiral STF are likely to be
accompanied by changes in the relative permittivity scalars $\epsilon_{a,b,c}$,
which attribute of STFs was not accommodated in the theoretical study presented.
Because the chirping step is very small compared with the mean half--period
$\left( \Sigma_{n=0}^{N_p}\,\Omega_n\right) /N$ however, those changes are expected
have insignificant effect.

To conclude, the presented theoretical study has demonstrated that
a stepwise chirping of the periodicity of a chiral sculptured thin film
considerably enhances the bandwidth of the Bragg regime,
thereby extending the frequency range of operation as a circular--polarization
filter.

\bigskip

\section{Introduction}
Theoretical as well as experimental studies have shown that
chiral sculptured thin films (STFs) display the so--called Bragg
phenomenon on axial excitation,
in consequence of their
periodically and unidirectionally nonhomogeneous constitution
\c{VLbook, WHL, HWKLR}. Let the
direction of nonhomogeneity of a chiral STF be parallel to the $z$ axis, while
the film completely occupies the region $0 \leq z \leq L$. When circularly
polarized, mono\-chromatic light falls normally on this film of sufficient
thickness, then it is
(i) almost perfectly reflected if the handedness of the incident
light coincides with the structural handedness of the film, and
(ii) almost perfectly transmitted if otherwise~---~provided
absorption within the film is negligible
and the free--space
wavelength $\lambda_0$ of the incident light lies within the so--called Bragg
regime.  Optical devices such as circular--polarization filters and spectral--hole filters
exploiting the
circular Bragg phenomenon have been fabricated using STF technology,
delivering performance
comparable to commercially available non--STF analogs \c{WHL,HWLM,HWTLM}.

The bandwidth of the
Bragg regime is directly proportional to the local birefringence of the
chiral STF. Although the serial bideposition technique
\c{HWKLR} yields highly birefringent chiral STFs, further enhancement
of the bandwidth can be desirable for certain applications. The objective of this communication
is to show that stepwise chirping of the periodicity of a chiral STF
can easily deliver a considerable enhancement of the bandwidth.

The stepwise chirped chiral STF is described in Section 2, and its
response to  normally incident plane waves is presented in Section 3.
Comparison is made with the response of an unchirped chiral STF.
An $\exp(-i\omega t)$ time--dependence
is implicit in this communication, with $\omega$ as the
angular frequency; vectors are underlined and dyadics are double--underlined; $\#r=x\#{u}_x+y\#{u}_y+z\#{u}_z$ is the
position vector with $\#{u}_x$, $\#{u}_y$ and $\#{u}_z$ as the cartesian
unit vectors; and  $\lambda_0$ is the free--space wavelength.

\section{Constitutive Relations}
The nonhomogeneous permittivity
dyadic of a chiral STF is expressed as follows \c{VLbook}:
\begin{equation}
\label{epsbasic}
\=\epsilon(\#r) = \epsilon_0\, \=S_z(z,h)\.\=S_y(\chi)\.
\Big[ \epsilon_a \,\#{u}_z\uz +\epsilon_b\,\#{u}_x\ux
\,  +\,\epsilon_c\,\#{u}_y\uy\Big]\.\=S_y^{-1}(\chi)\.
\=S_z^{-1}(z,h)\, , \quad 0 \leq z \leq L\, .
\end{equation}
In this equation, $\epsilon_{a,b,c}$ are three relative permittivity
scalars; $\epsilon_0 = 8.854\times 10^{-12}$~F~m$^{-1}$
is the free--space permittivity; the  tilt dyadic
\begin{equation}
\=S_y(\chi) = \#{u}_y\uy + (\#{u}_x\ux + \#{u}_z\uz) \, \cos\chi + (\#{u}_z\#{u}_x-\#{u}_x\#{u}_z)\,
\sin\chi
\end{equation}
represents the locally columnar microstructure of any chiral STF
with $\chi > 0^\circ$; while the rotation dyadic
\begin{equation}
\=S_z(z,h)=
\#{u}_z\uz + \left( \#{u}_x\ux+\#{u}_y\uy\right)\,\cos\left(\frac{\pi z}{\Omega}\right) \,
+ h\, \left( \#{u}_y\#{u}_x-\#{u}_x\#{u}_y\right)\,\sin\left( \frac{\pi z}{\Omega}\right)
\end{equation}
contains  $2\Omega$ as the structural period, with $h=1$ for
structural right--handedness and $h=-1$ for structural left--handedness.
Taking a whole number of periods to encompass the film thickness,
we choose the ratio $L/2\Omega$ to be a non--prime integer for filtering applications.

The stepwise chirped chiral STF comprises $N_p$ steps, each step having
$\ell_p$ complete periods, with $N_p \ell_p = L/2\Omega$.
The half--period of the n$^{th}$ step
is denoted by
\begin{equation}
\Omega_n = \Omega + \frac{2n-N_p-1}{2}\, \delta_\Omega \,,
\quad n \in [1,\,N_p]\,,
\end{equation}
wherein the chirping step $\delta_\Omega$ is selected to be a small fraction
of $\Omega$ to assure that  $\Omega_1 > 0$.
Thus, the expression of the
permittivity dyadic of the stepwise chirped chiral STF is the
same as \r{epsbasic}, except that the dyadic $\=S_z(z,h)$ has to be replaced
by $\=S_z^{chirp}(z,h)$, where
\begin{eqnarray}
\=S_z^{chirp}(z,h)
&=&
\#{u}_z\uz + \left( \#{u}_x\ux+\#{u}_y\uy\right)\,\cos\left[\frac{\pi (z-\zeta_n)}{\Omega_n}\right] \,
\nonumber
\\
&+& h\,\left( \#{u}_y\#{u}_x-\#{u}_x\#{u}_y\right)\,\sin\left[ \frac{\pi (z-\zeta_n)}{\Omega_n}\right]\,,
\quad \zeta_n \leq z \leq \zeta_{n+1}\,,\quad n \in [1,\,N_p]\,,
\end{eqnarray}
and
\begin{equation}
\zeta_n = 2\ell_p (n-1)\left( \Omega +\frac{n-N_p-1}{2}\,\delta_\Omega\right)\,
\end{equation}
is the total thickness of the first $n-1$ steps.
Note that the thickness $\zeta_{{N_p}+1}$ of the chirped chiral STF equals
$L$.

\section{Numerical Results and Discussion}
Let an arbitrarily polarized plane wave be normally incident
from the lower half--space $z \leq 0$  on either
of the two films described in the previous section. As a
result, a plane wave is reflected into the lower half--space and another
is transmitted into the upper half--space $z \geq L$.
he electric field phasors associated with the two plane
waves in the lower half--space are stated as
\begin{equation}
\#{E}_{inc}(\#{r},\lambda_0) =
 \left(  a_L \, \#{u}_+ + a_R \, \#{u}_- \right) \,
\exp\left( i k_0 z \right)  \, ;  \,\, z \leq 0\, ,
\end{equation}
\begin{equation}
\#{E}_{ref}(\#{r},\lambda_0) =
\left(  r_L \, \#{u}_- + r_R \, \#{u}_+ \right) \,
\exp\left( -i k_0 z \right) \, ;  \,\, z \leq 0\, ,
\end{equation}
where $\#u_\pm = (\#{u}_x \pm i \#{u}_y)/\sqrt{2}$ and
$k_0 = 2\pi/\lambda_0$ is the
free--space wavenumber. Likewise,
the electric field phasor in the upper half--space is
represented as
\begin{equation}
\#{E}_{trs}(\#{r},\lambda_0) =
 \left(  t_L \, \#{u}_+ + t_R \, \#{u}_- \right) \,
\exp\left[ i k_0 (z -L)\right]  \, ;  \,\, z \geq L \, .
\end{equation}
Here, $a_L$ and $a_R$ are the known amplitudes
of the left-- and the right--circularly polarized (LCP \& RCP)
components
of the incident plane wave;
$r_L$ and $r_R$ are the unknown amplitudes
of the reflected plane wave components; while $t_L$ and $t_R$
are the unknown amplitudes
of the transmitted plane wave components.

The reflection amplitudes $r_{L,R}$ and the transmission
amplitudes $t_{L,R}$  can be computed for specified incident amplitudes
($a_L$ and $a_R$) by following standard procedures
described elsewhere; see, e.g., \c{VLbook}.
Reflection and transmission coefficients are thus defined as the entries
in the 2$\times$2 matrixes in the following two relations:
\begin{eqnarray}
\label{eq15}
\left[ \begin{array}{cccc} r_L \\ r_R  \end{array}\right]  &=&
\left[ \begin{array}{cccc} r_{LL} & r_{LR} \\ r_{RL} & r_{RR}\end{array}\right] \,
\left[ \begin{array}{cccc} a_L \\ a_R  \end{array}\right] \, , \\
\label{eq16}
\left[ \begin{array}{cccc} t_L \\ t_R  \end{array}\right]  &=&
\left[ \begin{array}{cccc} t_{LL} & t_{LR} \\ t_{RL} & t_{RR}\end{array}\right] \,
\left[ \begin{array}{cccc} a_L \\ a_R  \end{array}\right]
\, .
\end{eqnarray}
The co--polarized transmission coefficients are denoted by $t_{LL}$ and
$t_{RR}$,
and the cross--polarized ones by $t_{LR}$ and $t_{RL}$; and similarly for the
reflection coefficients in \r{eq15}.

The reflectances $R_{LL} = \vert r_{LL}\vert^2$, etc., and the
transmittances $T_{LL} = \vert t_{LL}\vert^2$, etc., were calculated
as functions of the free--space wavelength $\lambda_0$ for a chiral STF
described by the following parameters: $h=1$, $\epsilon_a = 2.8$, $\epsilon_b = 3.4$,
$\epsilon_c = 2.9$, $\chi = 20^\circ$, $\Omega = 155$~nm, and $L = 126 \,\Omega$.
Graphs are presented in Figure 1, which clearly show the Bragg regime
delineated by $\lambda_0 \in [528,\,564]$~nm almost agrees
with the theoretical range $\lambda_0 \in \left[ 2\Omega\epsilon_c^{1/2},\,
2\Omega\tilde{\epsilon}_d^{1/2}\right]$ \c{VLbook}.
In this regime,
$R_{RR}$ and $T_{LL}$ are very high, while $R_{LL}$ and $T_{RR}$
are very low, because the chosen film is structurally right--handed.

For the stepwise chirped chiral STF, we chose $N_p=21$, $\ell_p = 3$
and $\delta_\Omega = 0.5$~nm. Thus, the smallest half--period $\Omega_1=150$~nm
and the largest half--period $\Omega_2=160$~nm, while the average
half--period is 155~nm (the same as $\Omega$ for Figure 1).
The computed reflectances
and transmittances, plotted in Figure 2 as functions of $\lambda_0$, show
that the Bragg regime now
is $\lambda_0 \in [514,\,580]$~nm.  The lower limit of the Bragg regime
is somewhat higher than $2\Omega_1\epsilon_c^{1/2}$, and the upper limit is
somewhat lower than $2\Omega_{N_p}\tilde{\epsilon}_d^{1/2}$.

Relative
to the 36~nm bandwidth in Figure 1, the Bragg regime in Figure 2
has a bandwidth of 66~nm~---~an enhancement by more than 80\%.
The enhancement comes without any noticeable reduction in the
polarization--sensitivity of the film. Furthermore, calculations show
that mild dispersion and absorption do not affect the enhancement
significantly.

Stepwise changes in the periodicity of a chiral STF are likely to be
accompanied by changes in the relative permittivity scalars $\epsilon_{a,b,c}$,
which attribute of STFs was not accommodated in the theoretical study presented.
Because the chirping step is very small compared with the mean half--period
$\left( \Sigma_{n=0}^{N_p}\,\Omega_n\right) /N$ however, those changes are expected
have insignificant effect.

To conclude, the presented theoretical study has demonstrated that
a stepwise chirping of the periodicity of a chiral sculptured thin film
considerably enhances the bandwidth of the Bragg regime,
thereby extending the frequency range of operation as a circular--polarization
filter.

\bigskip


\section{Introduction}

High-level computer vision tasks, including object detection and segmentation, are highly dependent on the quality of captured images. However, outdoor vision systems such as security monitoring and vision-based self-driving cars or driving assistance systems are easily influenced by severe weather or other environments. As shown in Fig.~\ref{fig:1}(a), raindrops adhered to camera lens or windshield can highly degrade the image or obstruct the visibility.

\begin{figure}
	\centering
	\includegraphics[width=0.9\linewidth]{figures/fig1_0306}
	\caption{Examples of adherent raindrops and/or mist. (a) Raindrops adhered to a camera lens (Image source: The Nature Conservancy). (b) Inhomogeneous mist adhered to a windshield. (c) Mist and raindrops adhered to a camera lens}
	\label{fig:1}
\end{figure}

Adherent raindrops (or waterdrops) exist widely in different scenes. Without enough shelter, rain streaks will result in raindrops naturally. In special environments like fishing-boats, splashing water can also form waterdrops. Raindrops can have various visual effects because of various drop sizes, different distances to the camera, and complex refraction, which makes it extremely difficult to model the raindrops manually.

Besides raindrops, another severe and commonly-observed interference, adherent mist, is rarely mentioned in computer vision. As shown in Fig.~\ref{fig:1}(b), mist adhered to the windshield can also severely hamper the vision, which is vitally dangerous for driving.

Adherent mist is often caused by the temperature difference of air between two sides of the glass window. Unlike the haze in the air, mist adhered to glass windows is usually inhomogeneous (e.g., Fig.~\ref{fig:1}(b)).

Interestingly, adherent mists can sometimes develop along with the raindrops, especially for driving, where there is a temperature difference between the warm steam inside and the cooled glass by the raindrops outside. Therefore, it is necessary to handle both adherent mist and raindrops together. Fig.~\ref{fig:1}(c) shows the overlapping of mist and raindrops adhered to a camera lens.

By modifying the formulas in \cite{chen2019gated,li2019single,qian2018attentive,ren2018gated}, the degraded image $I^{cap}$ captured by the camera can be modeled as a combination of raindrops $R$, the scattering of adherent mist $A$, and the clean background $B$ as follows:
\begin{align}
I^{cap} = ((1-M)\otimes B + R)\odot t + A \odot (1-t)
\end{align}
where M denotes to a binary mask showing the existence of raindrops, $\odot$ means element-wise multiplication, and t is the transmission map indicating how much information has passed through the adherent mist. $I^{cap}$, $M$, $B$, $R$, $t$, and $A$ are location-related maps (i.e., matrices) instead of constants.

In this study, we aim to remove the adherent mist and raindrops from a single impaired image to obtain a clean image. The proposed approach, based on convolutional neural networks (CNNs), can automatically visually remove raindrops and mist of varying degrees. It can benefit high-level vision algorithms and improve their robustness in severe weather.

As both raindrops and mist might block some objects partially or even entirely, it is challenging to restore a purely clean image. To better solve the problem, we apply smoothed dilated convolution \cite{wang2018smoothed} into our model to obtain a relatively large receptive field instead of general pooling operations that might lose spatial information. Dual attention modules consisting of channel attention and spatial attention are adopted for information selection and efficient learning. Perceptual loss is also applied in order to generate intuitively realistic images and avoid over-smoothing details \cite{johnson2016perceptual,ledig2017photo}.

Besides, to strengthen the spatial attention and handle the inhomogeneity, the classification activation map (CAM) \cite{zhou2016learning} for each captured image is then calculated based on a newly trained classification network and acts as the spatial attention input concatenated with the original RGB image. The classification network is designed to accurately judge if the collected image is clean, with raindrops only, or with both raindrops and adherent mist. The proposed CAM attention can effectively locate important regions of the input image and improve the restoration performance, without hand-crafted priors.

Thus, our contributions are four-fold as follows:
\begin{itemize}
	\item We address the visual interference problem caused by adherent mist and raindrops, which is frequent and essential but less studied. Without appropriate solutions, this problem possibly limits extensive outdoor camera-based systems.
	\item A CNN-based method is proposed to jointly remove adherent mist and raindrops. The proposed CNN is designed with smoothed dilated convolution, dual attention module, and perceptual loss to achieve state-of-the-art performance. The peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM) \cite{wang2004image} values of the test set is improved from 17.9213 and 0.6283 to 21.8073 and 0.7549, compared with the degraded input.
	\item We introduce CAM attention as an efficient spatial attention input. Except for accessible classification labels, it does not require standard assumptions to find out worthy-attention pixel locations. This technique can be easily implemented in other low-level image generation problems.
	\item The improvements of high-level algorithms using enhanced images are demonstrated, showing the importance of both the problem and the proposed method. We also provide a new public dataset, which contains image pairs of ground truth (clean) images and degraded images (with raindrops and adherent mist). Further studies could use our dataset to develop new algorithms.
\end{itemize}

\section{Related Work}

Most visibility enhancement algorithms related to rain weather focused on rain streaks \cite{fu2017removing,li2019single,li2017single,wang2019erl,zhang2018density}. However, the shapes and physical effects of raindrops exceedingly differ from those of rain streaks. Thus, these streak removal algorithms can not be directly applied to raindrop removal.

A few studies about raindrop detection have been proposed for years. Kurihata et al. \cite{kurihata2005rainy} utilize principal component analysis to learn the characteristics of raindrops, while Ito et al. \cite{ito2015adherent} employ the maximally stable extremal regions to detect raindrop candidates. However, these methods do not concentrate on raindrop removal. You et al. \cite{you2015adherent} pay attention to both raindrop detection and removal in videos. They notice the temporal change of intensity of raindrop pixels is different from that of non-raindrop pixels, which is not applicable to the single image case. Eigen et al. \cite{eigen2013restoring} may be first to address the single image raindrop removal problem. They use a straightforward idea which trains a shallow CNN with pairs of raindrop images and clean images. Unfortunately, the limited capacity of the CNN model restricts its performance, especially for large and dense raindrops.

Qian et al.'s work (DeRaindrop) in 2018 \cite{qian2018attentive} might be the first practical solution for the raindrop removal task. Specifically, they build a generative adversarial network (GAN) \cite{goodfellow2014generative} to generate raindrop-free images. Besides the GAN, DeRaindrop also includes a recurrent network to gradually locate raindrops’ location in the input image, working as the spatial attention for the GAN. The performance of DeRaindrop is related to this raindrop mask, whose ground truth is calculated by the difference between the degraded image and the clean image with a hand-crafted thresholding 30. However, it is difficult to apply this manual threshold to images with adherent mist, which contain complex variations of pixel values. Quan et al. \cite{quan2019deep} propose shape-driven attention for raindrop regions, which is based on priors about the shape of raindrops.

Few works primarily focus on the adherent mist. Some similar studies mainly concentrate on atmospheric haze (or fog) \cite{cai2016dehazenet,chen2019gated,fattal2008single,he2010single,ren2018gated}. On the one hand, they both have a scattering effect, which results in “white” color and degrades the original information. But compared to atmospheric haze, adherent mist is more likely to be inhomogeneous as it is related to the specific design of the equipment, such as which direction the warm steam comes from and how the window looks (e.g., mist adhered to the windshield near the driver’s seat in Fig.~\ref{fig:1}(b)).

Single image dehazing is initially explored using some priors-based methods. Fattal et al. \cite{fattal2008single} work on estimating the albedo of the scene for single image dehazing. He et al. \cite{he2010single} propose a dark channel prior which talks about the local minimum of the dark channel varies between haze and haze-free images.

Then some learning-based methods have been proposed for dehazing using a single image. Cai et al. \cite{cai2016dehazenet} introduce an end-to-end CNN to estimate the transmissions from images with haze. Ren et al. \cite{ren2018gated} propose multi-scale gated fusion network to improve the performance. And the smoothed dilation technique is adopted by Chen et al. \cite{chen2019gated} to remove gridding artifacts.

However, in our work, the superposition of adherent mist and raindrops results in complex local features. The influence of raindrops limits general dehazing algorithms using transmission estimation. It is also difficult to manually input an attention map similar to that in \cite{qian2018attentive} or \cite{quan2019deep}. Therefore, the proposed network is designed to efficiently utilize features and restore images for the practical scene.

\section{Methods}
\subsection{Network Architecture}

The proposed network is shown in Fig.~\ref{fig:2}. The degraded image concatenated with its CAM attention map is used as the input for a generative network to output the clean image. The CAM attention is generated through a newly trained classification network. Since the classification network can accurately distinguish the type of the degraded image, it can also inform us which regions contribute more to the classification result. The generative network applies the architecture of an autoencoder with some shortcut connections. The input is initially convolved by two convolutional layers with dual attention modules. After a convolutional layer with stride 2, the spatial size of the features is decreased. Smoothed dilated convolution modules, as well as dual attention modules, are successively applied to handle those features further. Shortcut connections are inspired by the residual learning concept \cite{he2016deep}. Finally, the feature maps are converted back to the original size by a transpose convolutional layer and tuned by two convolutional layers to generate the output. Instance normalization \cite{ulyanov2016instance} is also applied in the generative network.

\begin{figure}[t]
	\centering
	\includegraphics[width=1.0\linewidth]{figures/fig2_0306_2}
	\caption{The architecture of the proposed network. The degraded input image is firstly used to generate a CAM attention map with a classification network. Then the input image concatenated with the CAM attention map is handled by the generative network to obtain the clean output. The details of the smoothed dilated convolution module and the dual attention module are elaborated in the following figure}
	\label{fig:2}
\end{figure}

\subsubsection{CAM Attention}
It might be helpful to handle the inhomogeneous features \cite{qian2018attentive} as shown in Fig.~\ref{fig:1}, with additional spatial attention maps instead of pure CNNs. Although related works (i.e., \cite{qian2018attentive,quan2019deep}) have generated such attention maps, they use a manually set threshold or shape prior, which might be less feasible to do so for complex images with adherent mist and raindrops. Instead, we utilize CAM to adaptively generate the attention maps.

CAM is proposed by Zhou et al. \cite{zhou2016learning}, which contributes to the visual interpretation of a CNN ``black box''. For a classification network, CAM can show the importance of various spatial regions in the image to the prediction result. This mechanism is achieved by weighting different channels of the final convolutional layer's features (denoting as $f_k(x,y)$ for a spatial position $(x,y)$ of channel $k$) using the weights of the dense layer (denoting as $\omega_k$) with respect to the prediction result $t$, which can be described as \cite{zhou2016learning}:
\begin{align}
M_t(x,y) = \sum_{k}\omega_{k}^{t}f_k(x,y)
\end{align}
where $M_t$ denotes to the CAM result. The features are resized to the shape of the input before multiplication, and the activation map is generated by adding values from all channels for each spatial position and normalization.

The architecture of our simple classification network is briefly shown in Fig.~\ref{fig:2}, which is similar to the structure in \cite{he2020restoration}, with 7 convolutional layers and 1 dense layer. A global maximum pooling layer is omitted in Fig.~\ref{fig:2}, which is used to flatten 2D features to 1D feature vector before the dense layer. A softmax activation function is applied to obtain the one-hot result.

Specifically, our classification task aims to distinguish images in three situations: clean, with raindrops (and without mist), and with both adherent mist and raindrops. Once the class of the image is accurately identified, CAM can hopefully explain which regions in the input image contribute more to the prediction. Four examples of CAM attention are shown in Fig.~\ref{fig:3}.

\begin{figure}
	\centering
	\includegraphics[width=1.0\linewidth]{figures/fig3_0306}
	\caption{Examples of CAM attention. The first row shows the original RGB images, and the corresponding CAM is shown below. (a--b) Images from \cite{qian2018attentive} degraded by raindrop only. (c--d) with raindrops only. (c--d) Images collected in our dataset with both adherent mist and raindrops. All CAM results share the same color bar on the right}
	\label{fig:3}
\end{figure}

For images with raindrops only (Fig.~\ref{fig:3}(a--b)), CAMs show high activation for raindrop regions. And for images with both adherent mist and raindrops (Fig.~\ref{fig:3}(c--d)), except for some highlighted raindrops, we can visualize some highlighted edge textures in CAMs (e.g., the tall building on the left of Fig.~\ref{fig:3}(c), and the grids on the wall in Fig.~\ref{fig:3}(d)). The attention results are reasonable as the CNN might classify the image based on the existence of raindrops and the edge variation behind the mist.

Raindrops regions and edge textures (i.e., high-frequency details) are actually the parts that we mostly focus on when observing the restoration results, which means regions that are important to the classification task are probably significant for the restoration task too. Therefore, we concatenate the original RGB image with its CAM and feed them to the generative network for image restoration.

\begin{figure}
	\centering
	\includegraphics[width=0.9\linewidth]{figures/fig4_0305}
	\caption{The architecture of the smoothed dilated convolution module and the dual attention module. (a) The smoothed dilated convolution module.  The same module adopts the same dilation rate. The adding operation indicates the element-wise sum. (b) The dual attention module. The two dense layers in the channel attention part share the same weights for output from both GAP and GMP in each iteration. The multiplication operation in channel attention weights different channels, and the multiplication operation in spatial attention processes various spatial positions}
	\label{fig:4}
\end{figure}

\subsubsection{Smoothed Dilated Convolution Module}
Dilated convolution is one of the useful ways to obtain a relatively large receptive field without decreasing the feature size like pooling. The convolution kernel covers small regions larger than the kernel size with some skipped pixels. However, dilated convolution might lead to gridding artifacts as described in \cite{chen2019gated}. To handle this disadvantage, \cite{wang2018smoothed} propose smoothed dilated convolution by adding a general convolution layer before the dilated convolution operation. As a result, the result of the dilated convolution is not only relied on spaced pixels but also some continuous information.

The smoothed dilated convolution module applied in our method is shown in Fig.~\ref{fig:4}(a). A standard convolution layer is added before each dilated convolution layer. When the dilation rate of the dilated convolution kernel is $d$, the kernel size of the corresponding standard convolution kernel should be $2d-1$ in each dimension. In this way, the information from the skipped feature pixels in the input feature map could also be used in the dilated convolution to avoid gridding artifacts. Specifically, there are 6 smoothed dilated convolution modules in our network as indicated in Fig.~\ref{fig:2}. The first 3 modules apply the dilation rate of 2, and the last 3 modules adopt the dilation rate of 4.

\subsubsection{Dual Attention Module}
Considering the inhomogeneous distribution of mist and raindrops, another spatial attention mechanism with a self-attention concept besides the CAM attention is also adopted for feature selection. Besides, as introduced in \cite{hu2018squeeze}, channel attention can also be useful to improve information extraction. These two attention parts form the dual attention module. And 12 modules are inserted in the generative network as indicated in Fig.~\ref{fig:2}.

As shown in Fig.~\ref{fig:4}(b), we build the dual attention module based on \cite{woo2018cbam}. Firstly, feature maps are processed by global average pooling (GAP) and global maximum pooling (GMP), individually, and become 1D vectors with the length of the number of input channels. Then, the two vectors are fed into two shared dense layers individually before being added together. The shared dense layer means the same weights are used to the two inputs from GAP and GMP. The multiplication operation in the channel attention gives different channels in the feature tensor different weights according to the vector.

The spatial attention is followed by the channel attention. Maximum pooling along the channel axis (MPC) indicates that the maximum value at each spatial position is obtained by comparing pixels from all channels at the same position. Average pooling along the channel axis (APC) is a similar operation for average values. The multiplication operation weights different spatial positions in the input tensor according to the 1-channel attention map.

\subsection{Loss Function}
As indicated in \cite{johnson2016perceptual,ledig2017photo}, although pixel-wise mean squared error (MSE) might bring improvements on metrics like peak signal-to-noise (PSNR) and structural similarity index (SSIM) \cite{wang2004image}, the generated image is likely to be overly smoothed. Instead, a perceptual loss term can be applied to calculate the error from high-level views, which might be helpful in generating subjectively realistic images.

We calculate the perceptual loss based on the MSE of the output features from the 7th convolution layer of VGG19 \cite{ledig2017photo,simonyan2014very}. The VGG19 model is pretrained on the ImageNet dataset \cite{deng2009imagenet}, which might give a relatively high-level description of an image. Based on the restored image $I^{out}$ and the corresponding ground truth $I^{gt}$, the two output feature maps of the 7th convolutional layer are $f(I^{out})$ and $f(I^{gt})$, respectively. Then the perceptual loss is:
\begin{align}
Loss_{per} = \frac{1}{n_x \times n_y \times n_c} \sum_{x=1}^{n_x} \sum_{y=1}^{n_y} \sum_{c=1}^{n_c}(f(I^{gt})_{x,y,c} - f(I^{out})_{x,y,c})^2
\end{align}
where $x$, $y$, $c$ denote the dimensions of the feature map.

Finally, the perceptual loss is added together with the general MSE loss as:
\begin{align}
Loss_{total} = \lambda_1Loss_{mse} + \lambda_2Loss_{per}
\end{align}
where $\lambda_1$, $\lambda_2$ are the weights of items and $Loss_{mse}$ is calculated by directly comparing $I^{out}$ and $I^{gt}$. $\lambda_1$, $\lambda_2$ was set as 1.0 and 0.05, respectively.

\subsection{Dataset Preparation}
To prepare the image pairs consisting of clean and degraded images, we take photos following strategies similar to those in \cite{qian2018attentive}. Specifically, we utilize two same glass panels to simulate different cases. Firstly, we fixed the camera using a tripod. Secondly, one of the panels was randomly sprinkled with waterdrops on one side, and randomly sprayed with mist by a cosmetic sprayer on the other side. Then, we placed the blurred panel in front of the camera and took a photo as the degraded image. At last, the corresponding ground truth (i.e., clean) image was taken with the other clean glass panel.

The distributions of waterdrops and mist are random to keep good generalization. And the distance between the camera and the glass panel is also a random value between 0.5cm--4.0cm, thus the shape of waterdrops can have different appearances in the image due to the influence of camera focal length. The integrated camera of Redmi K20 Pro is used for taking pictures.

Finally, all the collected images were resized with bicubic interpolation to the shape of $640\times480$ to form the dataset. There are 572 image pairs for training, 73 image pairs for validation, and 72 image pairs for test. All the following evaluation results are based on images in the testing set.

\section{Experiments}
\subsection{Experiment Setup}

Both the classification network and the generative network are implemented on Keras framework with Tensorflow backend. Both models are trained using an Nvidia Titan RTX GPU with Adam optimizer. The learning rate of 1e-4 and the batch size of 8 were utilized for the classification network, while for the generative network, the learning rate is set as 1e-3, and the batch size is 4. During training, data augmentation operations, including flipping, rotation (180 degrees), brightness and contrast fluctuation, and gamma tuning, were applied.

In addition, only the collected dataset described in Sec. 3.3 is used to train the generative network. However, the image pairs of raindrops dataset from \cite{qian2018attentive} are mixed with our dataset to train the classification network for raindrops classification. To avoiding data leaking, the testing images of our captured data that will be used to evaluate the performance of the image restoration are not used for training the classification network.

PSNR and SSIM are adopted for evaluating our model's performance. Their values are calculated based on the luminance channel, so the RGB images are converted into YCrCb mode before evaluation.

\subsection{Classification Results}
The classification model was evaluated with 210 test images (66 from the dataset of \cite{qian2018attentive}, and 144 from our dataset). The classification network achieves the accuracy and f1-score of both about 96\%. The relatively high classification performance might contribute to generate relatively reliable CAM attention results.

\subsection{Restoration Results}
To our knowledge, there are no other comparable studies that work on the exact same problem (i.e., jointly removal of adherent mist and raindrops). Thus, we compare the proposed method with the following three strategies:
\begin{itemize}
	\item[(1)] Our method is compared with two two-step combinations. A two-step combination method consists of a pre-trained general dehazing approach (that originally focuses on atmospheric haze) and the pre-trained raindrop removal model DeRaindrop (D) by \cite{qian2018attentive} to restore the image. In other words, the degraded image is firstly processed by a dehazing model (i.e., GCANet \cite{chen2019gated} (G), or FFANet \cite{qin2019ffa} (F)) and then by D. The two combinations are named as ``G+D'' and ``F+D'', individually.
	\item[(2)] We also tried to apply D before G (or F) for comparison. These two inverse combinates are denoted as ``D+G'' and ``D+F'', respectively.
	\item[(3)] Although originally designed for other tasks, D, G, and F are advanced generative network that could be applied to our dataset. Therefore, we re-trained these three networks using our dataset. The three re-trained methods are denoted as ``re-D'', ``re-G'', and ``re-F'', respectively.
\end{itemize}
The quantitative results of the three strategies above are listed in Table~\ref{table:1}. DI stands for ``Degraded Images''.
\setlength{\tabcolsep}{2pt}
\begin{table}
	\begin{center}
		\caption{The quantitative results of different methods}
		\label{table:1}
		\begin{tabular}{p{25pt}p{30pt}p{30pt}p{30pt}p{30pt}p{30pt}p{30pt}p{30pt}p{30pt}p{30pt}}
			\hline\noalign{\smallskip}
			   & DI & G+D & F+D & D+G & D+F & re-D & re-G & re-F & Ours\\
			\noalign{\smallskip}
			\hline
			\noalign{\smallskip}
			PSNR & 17.9213 & 17.4249 & 19.4054 & 18.1926 & 19.3557 & 19.3927 & 19.4788 & 20.3022 & \it{21.8073}\\
			SSIM & 0.6283 & 0.6601 & 0.6825 & 0.6758 & 0.6831 & 0.6982 & 0.6844 & 0.7197 & \it{0.7549}\\
			\hline
		\end{tabular}
	\end{center}
\end{table}

The results of Table~\ref{table:1} show that, firstly, it is important to handle adherent mist and raindrops together. The co-existence of mist and raindrops might lead to horrible results that even worse than the degraded images (e.g., the PSNR result of G+D in Table~\ref{table:1}). Secondly, the combination of F and D outperforms the combination of G and D with both processing orders. Thirdly, the three re-trained models generally show the effectiveness of re-training, especially for G and F. But our method always outperforms other methods.

In practical scenes, the combination of adherent mist and raindrops might severely impair the image, which might block objects from the subjective view. Examples in Fig.~\ref{fig:5} and Fig.~\ref{fig:6} exhibit the awful shooting environment and the restoration results from several methods above. Due to the space limit, we select D+G and F+D methods (Fig.~\ref{fig:5}), and re-F and re-D methods (Fig.~\ref{fig:6}) as representatives for visual comparison. The intuitive restoration examples of all the above methods are provided in the supplementary material.

\begin{figure}[t]
	\centering
	\includegraphics[width=1.0\linewidth]{figures/fig5_0305}
	\caption{Examples of results from D+G, F+D, and our method. Each row corresponds to the same sample. The values below images (except ground truth images) indicate the corresponding values of PSNR / SSIM}
	\label{fig:5}
\end{figure}
\begin{figure}[t]
	\centering
	\includegraphics[width=1.0\linewidth]{figures/fig6_0305}
	\caption{Examples of results from re-F, re-D, and our method. Each row corresponds to the same sample. The values below images (except ground truth images) indicate the corresponding values of PSNR / SSIM. The black bars in three of the sub-images of the second row are used to hide the location information}
	\label{fig:6}
\end{figure}

From Fig.~\ref{fig:5} and Fig.~\ref{fig:6}, we can find that some objects are not distinguishable in the degraded image. For example, in the first row of Fig.~\ref{fig:5}, it is hard to count how many cars in the degraded image, and in the second row of Fig.~\ref{fig:6}, the bicycles in front of the car almost disappear in the degraded image, which is vital for self-driving applications. If without the correct judgment, the safety of driving algorithms is to be doubted. This problem can be effectively handled by our method, as shown in the fifth column in Fig.~\ref{fig:5} and Fig.~\ref{fig:6}. In addition, the proposed method can adaptively handle degradations with different extents. For both light degradation (e.g., the last example in Fig.~\ref{fig:6}) and severe degradation (e.g., the second example in Fig.~\ref{fig:5}), our method can always keep the reasonable brightness and fine textures.

According to the attached PSNR / SSIM values (Fig.~\ref{fig:5} and Fig.~\ref{fig:6}), which are used to measure the similarity between the given image and the ground truth, it is proved that the proposed method can obtain relatively cleaner images and fewer artifacts than the other methods. The two combination methods (D+G and F+D in Fig.~\ref{fig:5}) sometimes result in serious luminance deviations. And the two re-trained methods (re-F and re-D in Fig.~\ref{fig:6}) show more artifacts than the proposed method.

\subsection{Ablation Study}
To verify the effectiveness of various components of the proposed method, we investigate the ablation study as follows. Besides the proposed architecture, we trained a model without CAM attention (Non-CAM), a model without channel attention (Non-CA), a model without spatial attention (Non-SA), and a similar model that replaces the smoothed dilated convolution by standard dilated convolution (Non-SD). Other designs of each model remain the same.

\setlength{\tabcolsep}{7pt}
\begin{table}
	\begin{center}
		\caption{The ablation study results from different network designs}
		\label{table:2}
		\begin{tabular}{cccccc}
			\hline\noalign{\smallskip}
			& Non-CAM & Non-CA & Non-SA & Non-SD & Ours\\
			\noalign{\smallskip}
			\hline
			\noalign{\smallskip}
			PSNR & 21.2289 & 21.4883 & 21.4294 & 21.2753 & \it{21.8073}\\
			SSIM & 0.7494 & 0.7546 & 0.7514 & 0.7442 & \it{0.7549}\\
			\hline
		\end{tabular}
	\end{center}
\end{table}

The statistical results of these four models using our testing set are shown in Table~\ref{table:2}. The performance decreases if any one of the above components is removed, especially for PSNR metrics. For example, without the CAM attention mechanism, the PSNR decrease by about 0.6dB. It is demonstrated that all those techniques effectively contribute to our model's performance.

\subsection{Applications}
Besides the quantitative and qualitative evaluation in Sec. 4.2, the importance of the proposed problem and the effectiveness of our method are also validated in practical applications, such as for instance segmentation. We apply a cutting-edge segmentation algorithm, Cascade R-CNN method \cite{cai2019cascade} from Detectron2 (https://github.com/facebookresearch/detectron2) model zoo to segment the images including the ground truth images, degraded images, and restored images as shown in Fig.~\ref{fig:7}. This advanced algorithm above also detects 2D objects with bounding boxes. For a fair and representative comparison, methods with relatively high-performance such as F+D and re-F are applied in this subsection. Segmentation results based on images restored by other methods are provided in the supplementary material.

\begin{figure}[t]
	\centering
	\includegraphics[width=1.0\linewidth]{figures/fig7_0305}
	\caption{Instance segmentation results using images from different sources: (From left to right) ground truth, degraded, F+D, re-F, and our method. Each row corresponds to the same sample. Small objects are best viewed by zoom-in}
	\label{fig:7}
\end{figure}

\begin{figure}[t]
	\centering
	\includegraphics[width=1.0\linewidth]{figures/fig8_0903}
	\caption{Panoptic segmentation results using images from various sources: (From left to right) ground truth, degraded, F+D, re-F, and our method. For each scene, the segmentation results are attached below their corresponding images. The white bars in the second sample are added to hide the location information. Small objects are best viewed by zoom-in}
	\label{fig:8}
\end{figure}

As shown in the first row of Fig.~\ref{fig:7}, there is a car ``missed'' in the degraded image, the image from F+D method, and the image from re-F method. And in the second row, the segmentation algorithm can not find the bicycles in the degraded image, and F+D recovered image. Thus, it is proved that adherent mist degraded images not only impede our subjective views, but also impair the performance of computer algorithms. For the third sample, the segmentation algorithm mistakenly regards the building as a train in the degraded image, and F+D recovered image. The error possibly leads to severe consequences in practical applications.

In addition, panoptic segmentation is a newly proposed segmentation task \cite{kirillov2019panoptic} which tries to label all the pixels in the image. The results of panoptic segmentation based on various images are shown in Fig.~\ref{fig:8}. We applied R101-FPN method in Detectron2 (https://github.com/facebookresearch/detectron2) model zoo as the panoptic segmentation algorithm.

Panoptic segmentation might be more natural to show the influence of adherent mist and raindrops, as it demands all the spatial pixels of the image. In the first sample (1st and 2nd rows in Fig.~\ref{fig:8}), the segmentation algorithm can not clearly find the border between the road and the grass, as well as the border between the tree and the sky in the degraded image. In the second sample (3rd and 4th rows in Fig.~\ref{fig:8}), segmentation results of the middle three images are horrible, especially for the degraded image. Because of the impact of adherent mist and raindrops, the trees, the building, and some bicycles can not be identified by the panoptic segmentation algorithm. The segmentation result with the output from our method shows relatively distinguishable contours for the bicycles, the pavement, the trees, and the building.

\section{Conclusion}
In this work, we propose an image degradation problem caused by the co-existence of mist and raindrops adhered to a glass window or camera lens. This phenomenon shows different image characteristics compared to atmospheric haze or rain streaks. The problem is widely observed in reality but lacks investigation. Instead of using hand-crafted priors to generate spatial attention maps to handle this problem, we introduce CAM attention to generate attention maps without priors, and utilize dual attention modules as well as smoothed dilated convolution modules in the restoration CNN to generate relatively clean images. Both quantitative evaluation and intuitive examples show quality improvement with the proposed method. The restoration with our method can also benefit high-level tasks such as instance segmentation and panoptic segmentation for generating reliable results, which demonstrates a vital significance in vision applications.

\bibliographystyle{splncs04}


\section{Introduction}

Massive close binaries are progenitors of many interesting objects as X-ray binaries and binary radio pulsars \citep{1993ARA&A..31...93V},  supernovae Ib and Ic \citep{1992ApJ...391..246P}, as well as collapsars resulting in gamma ray bursts \citep{2005A&A...435..247P}
and gravitational wave sources \citep{2018MNRAS.481.1908K}.
Massive binaries trace the evolutive history of these interesting objects
and can explain their galactic populations  \citep[e.g.][]{2020arXiv200207230Z}.
In addition, massive binaries are natural laboratories to understand the
physical processes occurring during the periods of mass exchange and mass loss, that regulate
their evolution determining its final fate \citep{2014ApJ...782....7D}.

The eclipsing close binary  UU\,Cassiopeiae  (UU\,Cas, BD +60 2629, 2MASS J23503951+6054391,
$\alpha_{2000}$= 23:50:39.52, $\delta_{2000}$= +60:54:39.14, $V$ = 9.74, Spectral Type B0.5III)\footnote{http://simbad.u-strasbg.fr/simbad/}
has  a distance based on the GAIA DR2 parallax of 3256 [+392 -319] pc \citep{2018AJ....156...58B}.
It is listed in the catalogue of massive close binaries of \citet{2004A&AT...23..213P} that  compiles observable data for 176 massive close binaries with main sequence components earlier than approximately B5. In this catalogue, the object occupies the place 7th among 65
binaries with orbital period known, ordered in descending order of period length.
UU\,Cas was studied by means of photographic spectra by \citet{1934ApJ....79...84S} who found an orbital period P$_{\rm o}$=  8\fd520676.
The General Catalog of Variable Stars\footnote{http://www.sai.msu.su/gcvs/} \citep{2017ARep...61...80S} gives an orbital period of  8\fd51929 with reference to \citet{PK40}.

\citet{2002ARep...46..900P} obtained $UBVR$ differential magnitudes between 1984-1989 and
derived stellar masses of 34.5 $\pm$ 1.5 and 25.7 $\pm$ 0.6 $M_{\odot}$.
She also noticed  large magnitude deviations from the assumed over-contact binary model at some epochs
which she has explained in terms of mass flows. In particular, she noticed the large
variability of eclipse depths and overall light curve shape. The complex variability of the light curve and large deviations were ratified
nine years later with new data  by \citet{2009arXiv0907.1047K}. Using
photometric observations of UU\,Cas obtained until year 2000, the following ephemerides is given for the primary minimum by \citet{2004AcA....54..207K}\footnote{https://www.as.up.krakow.pl/o-c/index.php3}:
Min I = JD 2428751.6762 + 8\fd519281 E.

\citet{2017AstBu..72..321G}  derived stellar masses of 17.7 and 9.5  $M_{\odot}$ assuming an orbital inclination of 69$^{\circ}$ and using medium resolution ($R$ = 15\,000) spectra.
He also  presented evidence
for a disk surrounding the hotter component and a common expanding envelope.
Recent H$\alpha$ Doppler tomography revealed the importance of the gas flows in the semi-detached system to understand the
different  flux contributions.  Gas stream from the cooler star, accretion disk around the hotter star and a wind
are taken into account by   \citet{2019ApJ...883..186K} to model the orbital variability of the H$\alpha$ profiles. These authors and  \citet{2017AstBu..72..321G}  noticed that the apparent paradox of a  deeper light minimum when
the more massive star is occulted by the less massive star can be explained if the more massive star is surrounded by an accretion disk. The presence of a disk
was already suggested by  \citet{2010ASPC..435..301D} and
\citet{2010POBeo..90..159M,2011BlgAJ..15...87M}.
The larger masses obtained by \citet{2002ARep...46..900P} compared with those of \citet{2017AstBu..72..321G} are  explained
by the assumed over-contact system configuration and the absence of the disk in her model.

In this study we hope to contribute to: (1) get for the first time a physical representation of the accretion disk  surrounding the more massive star
and at the same time get reliable and consistent stellar and system parameters,
(2) shed light on the complex photometric variability of the system and
the conflicting masses obtained from previous photometric and spectroscopic studies, (3)  use
 the photometric datasets acquired over last years by all-sky surveys and satellites to advance in the understanding of this system, (4)
 improve the
 value of the orbital period and check for its variability,
and (5)  get insights on the system structure using published infrared photometric data.

The paper is organized as follows: in Section 2 we give details of the photometric observations analyzed in this paper, in Section 3
we present our results based on photometric studies, in Section 4 we present the light curve model and their physical parameters, in Section 5 we discuss our findings in the context of earlier
work on the system and in Section 6 we summarize our conclusions.

\section{Observations}

\subsection{Photometric data and methodology}

The photometric data used in this paper were collected from previously published articles and publicly available databases.
We have studied the photometric time series found in The All-Sky Automated Survey for Supernovae ASAS-SN\footnote{\url{http://www.astronomy.ohio-state.edu/asassn/index.shtml}} in $V$ and $g$,
the data from Kamogata Kiso Kyoto Wide-field Survey \citep{110009799479} in $I_{c}$ and $V$. In addition, we include
data acquired with The Optical Monitoring Camera \citep[OMC,][]{2003A&A...411L.261M} on board the high-energy INTEGRAL satellite, that provides photometry in the Johnson $V$-band within a 5 by 5 degree field of view
\citep{2010ASSP...14..493D} and it is able to detect optical sources brighter than around $V$ $\sim$ 18.
We also include  differential photometry from
\citet[][$UBV$]{2009arXiv0907.1047K} and   \citet[][$UBVR$]{2002ARep...46..900P}.

We query for infrared magnitudes through the NASA/IPAC infrared science
archive\footnote{\url{http://irsa.ipac.caltech.edu/Missions/wise.html}}.
We investigated the Wide-field Infrared Survey Explorer WISE, a NASA medium class explorer mission that conducted
an all-sky survey at mid-infrared bandpasses centered around wavelengths 3.4, 4.6, 12 and 22 $\mu$m (hereafter
$W1, W2, W3$ and $W4$; \citet{Wright10}).
The survey was conducted with a 40 cm cryogenically-cooled telescope in sun-synchronous polar orbit.
Four infrared detectors imaged the same sky field of view during 7.7 s ($W1, W2$) and 8.8 s ($W3, W4$).
We use here the data of the second-pass processing,
obtained with improved calibration and processing algorithms, superseding those obtained for the preliminary data release.

From the all aforementioned datasets we present for the first time the analysis of the INTEGRAL, ASAS-SN, KWS and WISE data.
Table~\ref{Tab:phot} gives a summary of the data sources, the bands that were used, and the time span of the observations.

We search for periods using the GLS periodogram \citep{2009A&A...496..577Z}.
This algorithm uses the principle of the \citet{1976Ap&SS..39..447L} \& \citet{1982ApJ...263..835S} periodograms with some modifications, such as the addition of a displacement in the adjustment of the fit function and the consideration of measurement errors. Compared with the classical periodogram, it gives us more accurate frequencies and a better determination of the amplitudes.  We complement this analysis with
the $O-C$ method described by \citet{2005ASPC..335....3S}, where predicted times of eclipses of a  test ephemerides are compared with observed times
in order to get an improved value of the orbital period.

Furthermore, we disentangled the light curve subtracting from the original light curve a representation of the orbital variability
constructed with a Fourier series including the orbital frequency and its harmonics. This methodology was introduced in  \citet{Mennickent2012}
and used in several past studies of close binaries \citep[e.g.][]{2019MNRAS.487.4169M}. This allowed us to search for new periodicities in the residual data.
The results of this search are given in the next section.

\begin{table}
\caption{Photometric time series analyzed in this paper. KW refers to \citet{2009arXiv0907.1047K}, KWS means  Kamogata Kiso Kyoto Wide-field Survey
and P02 refers to \citet{2002ARep...46..900P}. N means number of data points and $\sigma$ average data point error.}
\label{Tab:phot}
\centering
\begin{tabular}{l l c c r}
\hline\hline
Source & Band & $\Delta$ HJD-2400000 &N &$\sigma$ (mag)\\
\hline
KW & $U$ & 42655-46055 &955 &NA\\
KW & $B$ & 42655-46055 &955 &NA\\
KW & $V$ & 42655-46055 &955 &NA\\
INTEGRAL & $V$ & (1093-6774)* &3697 &0.013\\
ASAS-SN  & $V$ & 57009-58451&615&0.0058\\
ASAS-SN   & $g$ & 58610-58795&245&0.0025\\
KWS & $V$& 56268-58807 &748&0.027\\
KWS & $I_c$ &  56484-58807 &616&0.027\\
WISE &$W1$&(55216-55411)** &55  &0.028\\
WISE &$W2$&(55216-55411)** &55  &0.021\\
WISE &$W3$&(55216-55411)** &37  &0.023\\
WISE &$W4$&(55216-55411)** & 37 &0.163\\
P02&  $U$ & *** &47 & 0.028\\
P02& $B$ & *** &55 & 0.020\\
P02& $V$ & *** &55 & 0.018\\
P02& $R$ & *** &59 & 0.016\\
    %
\hline
\end{tabular}
\vspace{0.25cm}\\
* The INTEGRAL satellite data is in Barycentric Julian Date
(BJD-2451544.5). ** The WISE satellite data is in Mean Julian Date (MJD-2400000).  ***Authors give only phases versus differential magnitudes.
\end{table}

\begin{figure}[]
\begin{center}
\includegraphics[angle=-90,width=1\linewidth]{perorbV2.pdf}
\end{center}
\caption{Periodogram obtained with the $V$-band photometry. $\rm 2f_{0} = 2/P_{o}$, where $\rm P_{o} = 8\fd519252 \pm 0\fd000012$ d.}
\label{Operiodogram}
\end{figure}

\begin{figure*}[h!]
\begin{center}
\includegraphics[width=1.0\linewidth]{OC.pdf}
\end{center}
\caption{Observed (O) minus calculated (C) epochs for primary minima versus cycle number for 117 years of observations, according to the ephemeris given by Eq.\,1.
 Black and red dots indicate photographic and visual data, respectively.  Blue, green and magenta dots show CCD or photoelectric $V$, $I_c$ and $g$-band data, respectively.
The data points with the same X and Y coordinates are displayed with a small X offset for better visualization. See Section 3.1. }
\label{ocfit}
\end{figure*}

\begin{figure}[]
\begin{center}
\includegraphics[width=1.0\linewidth]{newper.pdf}\\
\end{center}
\caption{Generalized Lomb-Scargle periodogram obtained with $I_c$-band data after subtraction of
the orbital period. The light gray curve shows the spectral window for the data. The best 2$^{th}$ order polynomial fit is also shown.  95\% confidence and prediction bands are shown by dashed and  dashed-light areas, respectively.}
\label{perlong}
\end{figure}

\begin{figure}
\begin{center}
\includegraphics[width=1.0\linewidth]{lcycle.pdf}
\end{center}
\caption{ Upper: $I_{c}$-KWS light curve after removing the orbital period. Lower: The same but phased with the long period of 269 days. The best sinus fit is also shown. }
\label{lcycle}
\end{figure}

\section{Results}
\label{Sec:Results}

\subsection{The orbital period }

We searched for the orbital period
using the $V$-band photometry, finding an orbital period of $8\fd519252(12)$ (Fig.\ref{Operiodogram}). This value compares well with that provided by
\citet{2004AcA....54..207K}, viz.\,8\fd519281. The following ephemeris for the primary eclipse was found:

\begin{equation}
  \rm HJD_{min} = 2443132.29826 + E\times 8\fd519252(12)
 \label{eq1}
\end{equation}

In order to improve the accuracy of the orbital period, we selected the photometric data points with phases close to the primary minimum
($\rm 0.98 \leq \Phi_{o} \leq 1.02$) and performed an analysis of observed (O) minus calculated (C) eclipse times following  \citet{2005ASPC..335....3S}.
In this analysis the O-C deviations can be represented as a function of the number of cycles with a straight line whose
zero point and slope are the corrections needed for the (linear) ephemeris zero point and period, respectively. We include 48 times of primary minimum
studied by \citet{2004AcA....54..207K} and listed in the O-C Gateway\footnote{http://var2.astro.cz/ocgate/index.php?lang=en}. Additionally, we include
209 new times of minima measured from the data listed in Table\,1,   excluding satellite data because of the different time system. The dataset of minima covers about 5000 cycles i.e.
117 years of observations.
 Since some minima times are published without errors, we use a simple least square fit for our analysis.
Our result displayed in Fig.\ref{ocfit} shows that the fit can be performed with a straight line, indicating a constant period.
The new ephemeris is given by:

\begin{eqnarray}
 \rm HJD_{min} = 2443132.25696 + E\times 8\fd519296(8)
 \label{eq2}
\end{eqnarray}

\begin{figure*}[]
\begin{center}
\includegraphics[width=1.0\linewidth]{integral.pdf}
\end{center}
\caption{Integral $V$-band data phased with the ephemeris given by Eq.\,2, along with a locally weighted scatterplot smoothing curve.
Residuals are shown in the upper part with data points error bars. The orbital phase has been displaced for better visualization.}
\label{integral}
\end{figure*}

\begin{figure}[]
\begin{center}
\includegraphics[width=8.5cm]{histo.pdf}
\end{center}
\caption{Histogram of residuals shown in Fig.\ref{integral} modeled by the sum of two  Gaussians. A  non-linear fit using the
Levenberg Marquardt iteration algorithm results in two distributions (black lines) and its sum (red line). The  adjusted parameters full width at half maximum, center, area and their errors
are given for the two  Gaussians.
}
\label{histo}
\end{figure}

\begin{figure}[]
\begin{center}
\includegraphics[width=8.5cm]{Vg.pdf}
\end{center}
\caption{ASAS-SN $V$ and $g$-band magnitudes phased with the ephemeris given by Eq.\,2. The drops in brightness seen in INTEGRAL data are also present.  Phase zero corresponds to the primary (deeper) eclipse.}
\label{Vg}
\end{figure}

\subsection{Search for additional photometric periods }

We disentangled the light curve removing the orbital frequency finding a long period in the residuals  of the $I_{c}$-KWS band photometry; 268.7 $\pm$ 1.6 d with  full
amplitude 0.1096 $\pm$ 0.0066 mag. The following ephemeris  for the minimum of the long-cycle was found:

\begin{eqnarray}
 \rm HJD_{min} = 2456493.37 + E\times (268\fd7 \pm 1\fd6)
 \label{eq3}
\end{eqnarray}

No periodicity was found in the residuals of the other datasets.
The periodogram showing the peak at the long cycle frequency is shown in
Fig.\ref{perlong}. The residual light curve phased with the long cycle is shown in
Fig.\ref{lcycle}. The peak of the $I_c$-band periodogram is significant  and cannot be explained in terms of the
data sampling. \\

\subsection{Flickering in the light curve}

We observe flickering in the light curve, better represented in the INTEGRAL database (Fig.\ref{integral}). This phenomenon
consists of drops of brightness in the light curve regarding the mean level for a given orbital phase. It was previously reported as model deviations in several bands up to 0.1 mag
in earlier work with less time-resolved datasets by \citet{2002ARep...46..900P}. Larger variations of up to 0.2 mag during 9 nights between 1972 and 1973 were recorded by \citet{2009arXiv0907.1047K} but interpreted in terms of an instrumental effect. The drops in brightness observed in INTEGRAL data cannot be due to measurements errors:  individual data point errors are much smaller than the observed deviations and the expected instrumental error for each data point for the given system magnitude is less than 0.02 mag  \citep{2003A&A...411L.261M}.

 A Levenberg Marquardt nonlinear curve fit was performed with two Gaussians characterized by free parameters
position, full width at half maximum and area. The fit converged after 19 iterations with a $\chi^2$ tolerance  value of 10$^{-9}$ and R-square value of 0.977.
The best fit reveals two different distributions
separated by 0.062 $\pm$ 0.014 mag and whose heights are in the ratio 4.1 $\pm$ 0.6 (Fig.\ref{histo}). The lower
distribution is 2.4 $\pm$ 0.3 times broader than the tallest distribution. The areas of the distributions are $A_1$ = 11.5 $\pm$ 1.0 and $A_2$ = 6.8 $\pm$ 1.5 mag $\times$ counts.  The magnitudes in the half of the smaller distribution corresponding to larger residuals (fainter magnitudes) can be identified as fading events.
The corresponding ratio of areas $\eta$ = $0.5A_2/A_1$ is the relative fraction of magnitudes that can be classified as magnitude drops. We find that 29.6\% of the observations corresponds to fading events in this dataset. In addition, we find  that the drops  occurs randomly,  none periodicity was found in the residuals of INTEGRAL data.
They also appear in ASAS-SN $V$ and $g$-band data (Fig.\ref{Vg}).
The events  have not clear dependence on the orbital phase, they appear at all phases, but seems to be more abundant in quadratures. The larger drops amount to 0.3 mag in $V$-band, i.e.
a fall of roughly 30\% of the total flux.

 In order to investigate the time scale in which these fading events occur, we performed a correlation analysis measuring the temporal distance between adjacent data points
and the corresponding difference of magnitude. We find
that jumps in magnitude of the order of 0.1 mag already occur in time scale of minutes, and conclude that the time resolution of the data is not enough to resolve these events
 (Fig.\ref{residuals}).

\subsection{Infrared colors and $W1$ orbital variability}

We analyzed the $W1$ variability and compared the infrared colors with those of Double Periodic Variables (DPVs) and Be stars. DPVs
are close binaries mostly consisting of a B-type hotter star surrounded by an accretion disk fed by Roche-lobe overflow and mass transfer from a cooler giant companion
\citet{2017SerAJ.194....1M}. Be stars are rapidly rotating B-type stars surrounded by a circumstellar disk. Both types of objects usually show color excess mostly attributed to circumstellar reddening, therefore they are used here for comparison of presence of circumstellar matter. No correction for interstellar reddening was performed. However, at these wavelengths, the effect of interstellar reddening is negligible. In addition, we included as a reference the average colors for 136
main sequence stars in the range of spectral types from B1\,V to K3\,V from the Hipparcos catalogue \citep{1997A&A...323L..49P}.
For details of data selection of DPVs, Be stars and Hipparcos stars see \citet{2016MNRAS.455.1728M}. Our results are shown in Fig.\ref{wise}.

We find that the light $W$1 magnitude follows the orbital period. Our time resolution and cadence does not allow to resolve the whole orbital cycle neither the possible long-cycle.
Part of the variability occurs in longer time scales than the orbital one, this can be seen in the
data acquired in the end of the time series, that deviate to brighter magnitudes compared with the rest of the data of similar orbital phase.
In addition, the system appears in the  $W1-W2$ vs. $W2-W3$ diagram in a place populated by systems with circumstellar envelopes, far from the place
where the HIPPARCOS non-variable stars are located. This finding strongly suggests that circumstellar matter is also present in UU\,Cas.
In addition, during the orbital cycle, the system transits mostly in the $W1-W2$ axis, describing a horizontal path in the color-color diagram.

The same tendency is found when comparing the published 2MASS  $J-H$ (0.194) versus $H-K$ (0.274) colors of UU\,Cas \citep{2006AJ....131.1163S} with those of synthetic models for dwarfs and giants, along with
observed colors of Be stars and DPVs (Fig.\ref{jhk}). The UU\,Cas system is located in a similar place regarding the Double Periodic Variables, suggesting the existence of circumstellar matter and eventually an accretion disk.

\begin{figure*}[]
\begin{center}
\includegraphics[width=1.0\linewidth]{corr.pdf}
\end{center}
\caption
  The temporal distance between adjacent data points of all INTEGRAL $V$-band data is plotted against the corresponding difference of magnitude. We find that flickering occurs in time scales of minutes. Data points errors are $\approx$ 0.01 mag.}
\label{residuals}
\end{figure*}

\begin{figure*}[]
\begin{center}
\includegraphics[width=0.8\linewidth]{wise.pdf}
\end{center}
\caption{Upper: $W$1 magnitude versus orbital phase according to ephemeris given by Eq.\,2. Colors represent time intervals.
Lower: Color-color diagram for non variable HIPPARCOS stars of spectral type B1\,V - K3\,V, Be stars and Double Periodic Variables (DPVs)
according to \citet{2016MNRAS.455.1728M}.  The UU\,Cas data is added with colors  along with numbers (in some cases) showing the orbital phases.   Phase zero corresponds to the primary (deeper) eclipse.}
\label{wise}
\end{figure*}

\begin{figure}[]
\begin{center}
\includegraphics[width=0.9\linewidth]{jhk.pdf}
\end{center}
\caption{Colors for Galactic DPVs \citep{2016MNRAS.455.1728M}, Be stars \citep{2001A&A...369...99H} and synthetic stellar atmosphere models with $\log g$\ = 4.0 (squares with lines)
and $\log g$\ = 3.0 (open circles with lines) from \citet{1998A&A...333..231B}. The effective temperature of 2 selected synthetic models have been labeled along with the position of UU\,Cas.}
\label{jhk}
\end{figure}

\section{Light curve model}

In this section we model the orbital light curves with a theoretical code that solves the inverse problem considering a semi-detached
system consisting of a  less massive and cooler "donor"  star and a more massive and hotter "gainer" star surrounded by an accretion disk, that is both optically and geometrically thick
\citep{1992Ap&SS.196..267D, 1992Ap&SS.197.17D, 1996Ap&SS.240..317D}. The model includes hot and bright spots in the disk, following evidence found in previous observations of algols \citep{2004AN....325..229R}. These active regions influence the shape of the light curve during the ingress and egress of the eclipses. The temperature of the disk matches the
gainer's one in its inner edge, and decreases with a radial profile described by an exponent $a_T$. The model has been described in several papers, so we remit the interested reader to \citet{2013MNRAS.432..799M, bib:paper6, 2019MNRAS.487.4169M}.

We fix the donor temperature  and mass ratio to $\rm T_{c}$ = 22\,700\,K and $q$= 0.52, following results from the spectroscopic studies by \citet{2019ApJ...883..186K} and \citet{2017AstBu..72..321G}
and also our own work based on disentangling of spectra (work in preparation).
We also assume synchronous rotation for the  donor, as expected for a close binary that rapidly synchronize stellar spins with the orbit due to tidal forces.
On the other hand, the gainer might have been spun-up to a high rotation due to nearly tangentially infalling material \citep{1981A&A...102...17P},
hence we have assumed critical rotation for it. For $I_c$ data, we modeled the residuals after subtracting the long cycle.

\subsection{Results of the light curve model}

The parameters of the best fits are shown in Table\, \ref{TabUUCas}.
 KWV-$UBV$ and P02-$UBVR$ multi-color data was fit independently and parameters were averaged.
Examples of fits to the KWS-$I_c$ and Integral $V$-band light curves are shown in Fig.\ref{models}. We obtain a good match between observed and calculated magnitudes in all studied cases,
along with consistency  for the stellar and orbital parameters. The disk parameters show some variations, especially the bright spot position, the spots temperature and the
disk vertical thickness and radius.
Our model indicates a system seen at angle 74$^{\circ}$\ and a stellar separation of 52 ${\rm R_{\odot}}$. The stellar temperatures are 22\,700\,$K$ and 30\,200\,$K$.
The stellar masses are 9 and 17.4 ${\rm M_{\odot}}$ and surface gravities $\log g$\ =  2.94 and 3.98. The stellar radii are 16.9 and  7.0 ${\rm R_{\odot}}$. The stellar
temperatures indicate spectral types of B\,2 and B\,0  for the donor and gainer stars \citep{1988BAICz..39..329H}.

The disk radius is relatively stable around a value of 21 ${\rm R_{\odot}}$, reaching its smaller value at 18 ${\rm R_{\odot}}$ (INTEGRAL-$V$ data) and its maximum value at 22 ${\rm R_{\odot}}$ (P02 data). It has a temperature of 16\,000\,$K$ at its outer edge.
The inner disk thickness has an average value of 9 ${\rm R_{\odot}}$ and the outer vertical thickness has an average  value of 6.5 ${\rm R_{\odot}}$, being maximum at 9.4 ${\rm R_{\odot}}$
(KWS-$I_c$ data).
The hot spot is 36\% (P02 data) to 70\% (KWS-$I_c$ data) hotter than the disk and is located consistently 42$^{\circ}$\ apart from the line joining the stellar centers in direction opposite to the orbital motion,
i.e. it  is located roughly at the expected region where the gas stream impacts the disk.  The bright spot is 24\% (P02 data) to 50\% (KWS-$I_c$ data) hotter than the disk and changes its location considerably, being located between 82$^{\circ}$ (KW data) and 141$^{\circ}$ (KWS-$g$ data) from the line joining the stellar centers in direction of the orbital motion.
 At quadratures, the flux ratio between donor and disk  is $\approx$ 2 in the $I_c$-band,  while in the $V$-band this ratio is $\approx$ 3.
The gainer is  almost completely hidden by the disk at all orbital phases contributing less than 10\% to the total flux in both bands.

\begin{table*}[h!]
\small
\caption{Results of the analysis of  {$\rm UU \ Cas$} orbital light curves.
Dataset labels refers to Table\,1.
The models are obtained by solving the inverse problem for the Roche model with
an accretion disk around the more-massive (hotter) gainer in critical
non-synchronous rotation regime \citep{1992Ap&SS.196..267D, 1992Ap&SS.197.17D, 1996Ap&SS.240..317D}. Mean values are also given.}
 \label{TabUUCas}

        \[
        \begin{array}{lllllllll}

               \hline

               \noalign{\smallskip}
{\rm Quantity}&UBV& UBVR& g& V& I_c& V &{\rm Average \pm std}\\
{\rm  }&{\rm  KW}& {\rm  P02}& {\rm  ASAS-SN}& {\rm  KWS}& {\rm  KWS}&  {\rm INTEGRAL} &\\
            \noalign{\smallskip}

            \hline

            \noalign{\smallskip}

   q                               & 0.52           & 0.52          & 0.52         & 0.52        & 0.52        & 0.52      &0.52  \\
   i {\rm [^{\circ}]}              & 74.2 \pm 0.3   & 74.7 \pm 0.3  & 74.3 \pm 0.6 &74.3 \pm 0.4 &74.9 \pm 0.4 &74.6 \pm 0.5 &74.5 \pm 0.3\\
{\rm F_d}                          & 0.91 \pm 0.04  & 0.97 \pm 0.03 & 0.95 \pm 0.03& 0.88\pm 0.02&0.96 \pm 0.03&0.80 \pm 0.03&0.91 \pm 0.06\\
{\rm T_d} [{\rm K}]                & 16180\pm 390   & 14970\pm 300  &16030 \pm 300 &16320\pm 300 &16100\pm 300 &15720 \pm 300&15887 \pm 491\\
{\rm d_e} [a_{\rm orb}]            & 0.116\pm 0.006 & 0.138\pm 0.03 &0.082 \pm 0.01&0.109\pm 0.02&0.179\pm 0.01&0.121 \pm 0.02&0.12 \pm 0.03\\
{\rm d_c} [a_{\rm orb}]            & 0.171\pm 0.016 & 0.178\pm 0.01 &0.190 \pm 0.01&0.183\pm 0.01&0.107\pm0.01&0.189\pm 0.01&0.17 \pm 0.03\\
{\rm a_T}                          & 8.5  \pm 0.6   & 8.0  \pm 1.4  & 8.3  \pm 0.9 & 8.4 \pm 0.6 &8.5 \pm 0.6 &8.2  \pm 0.5 &8.3 \pm 0.2\\
{\rm f_g}                          & 11.8 \pm 0.2   &11.95 \pm 0.1  &11.91 \pm 0.2 &11.97\pm 0.2&12.0 \pm 0.3 &12.0 \pm 0.2 &11.94 \pm 0.08\\
{\rm F_h}                          & 1.00           & 1.00          & 1.00         & 1.00        & 1.00        & 1.00       &1.00 \\
{\rm T_h} [{\rm K}]                       & 30260\pm 400   & 30150\pm 400  & 30260\pm 400 &30230 \pm 400&30200 \pm 400&30220\pm 400 &30220 \pm 42 \\
{\rm T_c} [{\rm K}]                       &  22700          & 22700         & 22700        & 22700       & 22700       & 22700      &22700 \\
{\rm A_{hs}=T_{hs}/T_d}             & 1.46 \pm 0.13  & 1.36 \pm 0.12 & 1.50 \pm 0.08&1.57 \pm 0.07& 1.70\pm 0.06&1.54 \pm 0.07&1.52 \pm 0.11\\
{\rm \theta_{hs}}{\rm [^{\circ}]}    & 22.8 \pm 1.2   & 24.4 \pm 0.9  & 22.9 \pm 1.2 &22.6 \pm 1.0 &22.8 \pm 0.9 &23.2 \pm 1.1& 23.1	\pm 0.7\\
{\rm \lambda_{hs}}{\rm [^{\circ}]} & 313.0\pm 5.1   & 316.5\pm 12.8 &326.0 \pm 7.1 &313.5\pm 5.6 &320.5\pm 6.2 &319.1\pm 7.1&318.1 \pm	4.9 \\
{\rm \theta_{rad}}{\rm [^{\circ}]}   & -31.0\pm 6.1   & -30.2\pm 5.1  &-27.4 \pm 5.6 &-29.3\pm 7.0 &-31.8\pm 8.2 &-30.9\pm 9.2&-30.1	\pm 1.6 \\
{\rm A_{bs}=T_{bs}/T_d}             & 1.31 \pm 0.05  & 1.24 \pm 0.08 & 1.40 \pm 0.08&1.33 \pm 0.07&1.50 \pm 0.08&1.34 \pm 0.07&1.35 \pm 0.09\\
{\rm \theta_{bs}} {\rm [^{\circ}]}   & 52.4 \pm 3.1   & 51.8 \pm 2.2  & 52.0 \pm 2.7 &51.6 \pm 3.6 &52.5 \pm 4.0 &40.5 \pm 3.2 &50.1 \pm 4.7\\
{\rm \lambda_{bs}}{\rm [^{\circ}]} & 82.2 \pm 13.2  &124.9 \pm 12.9 & 141.2\pm 9.7 &107.9\pm 8.9 &118.4\pm 9.4 &98.2 \pm 7.5&112.1 \pm 20.8 \\
{\Omega_{\rm h}}                        & 9.467\pm 0.03  & 9.554\pm 0.02 & 9.534\pm 0.06&9.562\pm 0.05&9.592\pm 0.06&9.586\pm 0.05&9.55	\pm 0.04\\
{\Omega_{\rm c}}                        & 2.914\pm 0.02  & 2.914\pm 0.02 & 2.914\pm 0.02&2.914\pm 0.02&2.914\pm 0.02&2.914\pm 0.02&2.914 \pm  0.02\\
\cal M_{\rm_h} {[\rm M_{\odot}]}  & 17.4 \pm 0.2   & 17.4 \pm 0.2  & 17.4 \pm 0.3 &17.4 \pm 0.3 &17.4 \pm 0.3 &17.4 \pm 0.3&17.4 \pm  0.3 \\
\cal M_{\rm_c} {[\rm M_{\odot}]}  & 9.0  \pm 0.2   & 9.0  \pm 0.2  & 9.0  \pm 0.2 & 9.0 \pm 0.2 &9.0 \pm 0.2 &9.0 \pm 0.2 &9.0 \pm  0.2\\
\cal R_{\rm_h} {\rm [R_{\odot}]}   &  7.1 \pm 0.1   & 7.0  \pm 0.1 & 7.0  \pm 0.1 & 7.0 \pm 0.1& 7.0 \pm 0.1&7.0 \pm 0.1& 7.0 \pm  0.1\\
\cal R_{\rm_c} {\rm [R_{\odot}]}   & 16.9 \pm 0.1   & 16.9 \pm 0.1  & 16.9 \pm 0.1&16.9 \pm 0.1&16.9\pm 0.1&16.9\pm 0.1& 16.9 \pm  0.1\\
{\rm log} \ g_{\rm_h}                    & 3.98 \pm 0.02  & 3.98 \pm 0.02 & 3.98 \pm 0.02&3.98 \pm 0.02&3.99 \pm 0.02&3.99 \pm 0.02&3.983 \pm 0.005\\
{\rm log} \ g_{\rm_c}                    & 2.94 \pm 0.02  & 2.94 \pm 0.02 & 2.94 \pm 0.02&2.94 \pm 0.02&2.94 \pm 0.02&2.94 \pm 0.02&2.94 \pm  0.02\\
M^{\rm h}_{\rm bol}                     &-6.7  \pm 0.1   &-6.6 \pm 0.1   &-6.6  \pm 0.1 &-6.6 \pm 0.1 &-6.6\pm 0.1 &-6.6\pm 0.1&-6.62 \pm 0.04\\
M^{\rm c}_{\rm bol}                     &-7.3  \pm 0.1   &-7.3 \pm 0.1   &-7.3  \pm 0.1 &-7.3 \pm 0.1 &-7.3 \pm 0.1 &-7.3 \pm 0.1&-7.3 \pm  0.1\\
a_{\rm orb}  {\rm [R_{\odot}]}     & 52.2 \pm 0.3   & 52.2 \pm 0.2  &52.2  \pm 0.3 &52.2 \pm 0.3 &52.2 \pm 0.3 &52.2 \pm 0.3&52.2 \pm  0.3 \\
\cal{R}_{\rm d} {\rm [R_{\odot}]} & 20.7 \pm 1.1   & 22.0 \pm 0.6  & 21.6 \pm 0.3 &20.1 \pm 0.3 &21.9 \pm 0.4 &18.3 \pm 0.4 &20.8 \pm 1.4\\
\rm{d_e}  {\rm [R_{\odot}]}        & 6.0  \pm 0.4   & 7.2  \pm 1.2  & 4.3  \pm 0.7 & 5.7 \pm 0.5 &9.4  \pm 0.4 &6.3  \pm 0.4 &6.5 \pm 1.7\\
\rm{d_c}  {\rm [R_{\odot}]}        & 8.9  \pm 0.9   & 9.3  \pm 0.5  & 9.9  \pm 0.7 & 9.6 \pm 0.7 &5.6  \pm 0.6 &9.8  \pm 0.5&8.8 \pm 1.6 \\

\noalign{\smallskip}
            \hline
         \end{array}
      \]

FIXED PARAMETERS: $q={\cal M}_{\rm c}/{\cal M}_{\rm
h}=0.52$ - mass ratio of the components, ${\rm T_c=22\,700 K}$  - temperature
of the less-massive (cooler) donor, ${\rm F_c}=1.0$ -
filling factor for the critical Roche lobe of the donor,
$f{\rm _{c}}=1.00$ - non-synchronous rotation coefficients
of the donor, ${\rm F_h}=R_h/R_{zc}$ - filling factor for the
critical non-synchronous lobe of the hotter, more-massive gainer (ratio of the
stellar polar radius to the critical Roche lobe
radius along z-axis for a star in critical non-synchronous rotation regime) ,
${\rm \beta_{h,c}=0.25}$
- gravity-darkening coefficients of the components, ${\rm
A_{h,c}=1.0}$  - albedo coefficients of the components.

\smallskip \noindent Note: Origin of the photometric observations,
photometric band, $q$ - mass ratio of the components, $i$ - orbit
inclination (in arc degrees), ${\rm F_d=R_d/R_{yc}}$ - disk
dimension factor (the ratio of the disk radius to the critical Roche
lobe radius along y-axis), ${\rm T_d}$ - disk-edge temperature,
$\rm{d_e}$, $\rm{d_c}$,  - disk thicknesses (at the edge and at
the center of the disk, respectively) in the units of the distance
between the components, $a_{\rm T}$ - disk temperature
distribution coefficient, $f{\rm _g}$ - non-synchronous rotation coefficient
of the more massive gainer (in the synchronous rotation regime),
 ${\rm T_h}$  - temperature of the more-massive (hotter) gainer,
${\rm A_{hs,bs}=T_{hs,bs}/T_d}$ - hot and bright spots' temperature
coefficients, ${\rm \theta_{hs,bs}}$ and ${\rm \lambda_{hs,bs}}$ -
spots' angular dimensions and longitudes (in arc degrees, the longitude
is measured in direction of the orbital motion), ${\rm
\theta_{rad}}$ - angle between the line perpendicular to the local
disk edge surface and the direction of the hot-spot maximum
radiation, ${\Omega_{\rm h,c}}$ - dimensionless surface potentials
of the hotter gainer and cooler donor, $\cal M_{\rm_{h,c}} {[\cal
M_{\odot}]}$, $\cal R_{\rm_{h,c}} {\rm [R_{\odot}]}$ - stellar
masses and mean radii of stars in solar units, ${\rm log} \
g_{\rm_{h,c}}$ - logarithm (base 10) of the system components
effective gravity, $M^{\rm {h,c}}_{\rm bol}$ - absolute stellar
bolometric magnitudes, $a_{\rm orb}$ ${\rm [R_{\odot}]}$,
$\cal{R}_{\rm d} {\rm [R_{\odot}]}$, $\rm{d_e} {\rm [R_{\odot}]}$,
$\rm{d_c} {\rm [R_{\odot}]}$ - orbital semi-major axis, disk
radius and disk thicknesses at its edge and center, respectively,
given in solar units.

\end{table*}

\begin{figure*}[h!]
\begin{center}
\includegraphics[width=0.8\linewidth]{fUUCasI-Int.pdf}
\caption{Observed (LCO) and synthetic (LCC) light-curves of
{$\rm UU \ Cas$} obtained by analyzing {\rm KWS-$I_c$} data and INTEGRAL  $V$-band  photometric observations; final
O-C residuals between the observed and optimum synthetic light curves;
fluxes of donor, gainer and of the accretion disk, normalized to the donor flux
at phase 0.25; the views of the optimal model at orbital phases 0.20, 0.48 and 0.80
({\rm $I_c$-KWS})  and  0.20 and 0.80 (INTEGRAL), obtained with parameters estimated by the light curve analysis.}
\label{models}
\end{center}
\end{figure*}

\section{Discussion}

Our results confirm the picture summarized by \citet{2017AstBu..72..321G} and \citet{2019ApJ...883..186K}
of a semidetached binary where the less massive star feeds an accretion disk located around the more massive and partly hidden star.
Our light curve models based on 11 different datasets obtained at different epochs and with different photometric bands, indicate consistent results for the overall
system configuration, especially the orbital and stellar parameters, consolidating the overall view of a massive interacting binary of early B-type components.

Our light curve models indicate a system of B0\,IV + B2\,III stars with a relatively large total
mass of about 26 ${\rm M_{\odot}}$. Our stellar masses of 17.4 and 9 ${\rm M_{\odot}}$ are  relatively close to those reported by \citet{2017AstBu..72..321G},
namely 17.7 and 9.5 ${\rm M_{\odot}}$. The mass, radius and temperature of the more massive component fit well the relationships for B-type dwarfs reported by \citet{1988BAICz..39..329H}. However, the less massive component is much larger than expected for its temperature and indicates an evolved star, consistent with its lower surface gravity.
In addition, we find an accretion disk  around the more massive star and derive for the first time their properties.
The average disk radius of  21  ${\rm R_{\odot}}$ means ${\rm R_{d}/a}$ $\approx$ 0.40, i.e. the disk outer border is just below the tidal radius for the given mass ratio \citep{1977ApJ...216..822P, 1995CAS....28.....W}. The disk radius shows small changes during the observing epochs, while the disk is characterized by two hot and active photometrically variable regions.
The changes in these regions (position, extension, temperature) might indicate changes in mass transfer rate modulating the mass flows and disk properties.
The infrared excess obtained in our photometric study confirms the existence of circumstellar matter while the enigmatic obscuration events mostly visible in INTEGRAL $V$-data
might be related to inhomogeneities present in the flow of mass loss.  The existence of these fading events suggests a mass flow structure  more complex than provided by our simple model of disk plus two spots. In the future, new data of better accuracy could add new constraints to improve the understanding of the system.

The long period observed in the $I_c$-band and the light curve morphology suggest that UU\,Cas might be a Double Periodic Variable \citep[DPV,][]{2003A&A...399L..47M, 2017SerAJ.194....1M}, a subtype of algols showing $\beta$-Lyrae type light curves that have  an accretion disk surrounding a B-type gainer.  DPVs show long cycles of larger amplitude at red passbands \citep{2010ASPC..435..357M}. This might indicate that the long cycle arises from a light source emitting mostly
at redder wavelengths than the other sources in the system. This could explain why the long cycle is detected (with small amplitude) only in the $I_c$-band in UU\,Cas. All DPVs (except $\beta$ Lyrae) show a constant orbital period
in spite of the evidence of active mass exchange and occasional mass outflows. All these credentials are similar to those of UU\,Cas. In addition, the ratio between the long cycle length and the orbital period is 31.7 for UU\,Cas,
 relatively close to the average $\sim$ 33 for DPVs.  Since the long period is observed only in one dataset, we suggest confirmation with additional data before claiming that UU\,Cas is a DPV. If confirmed, UU\,Cas should be the  most massive DPV;  typically they have total masses just around 10 ${\rm M_{\odot}}$.

The reported evidence for an H$\alpha$ P-Cygni like profile led \citet{2017AstBu..72..321G} to conclude that a stellar wind emerges from the system. This might indicate that the system
is loosing mass. However, the constancy of the orbital period through one century suggests that this mass loss should be  mild, i.e. the angular momentum of the system remains almost unaltered,
as well as the orbital period. This is consistent with a binary in a relatively stable and mild mass transfer stage.

The error in the orbital period $\epsilon$ = 0\fd000008, given by Eq.\,2, might be interpreted as a possible drift of the orbital period between
P$_o$ - $\epsilon$ and P$_o$ + $\epsilon$ in $\Delta$t =117 years. The above implies a possible change in the orbital period roughly of less than 2$\epsilon_P$/$\Delta$t = 1.37 $\times$ 10$^{-7}$ d/yr.
For conservative mass transfer this imposes an upper limit for the mass transfer rate. The expected period change in the conservative case is \citep{1963ApJ...138..471H}:

\begin{eqnarray}
\rm \frac{\dot{P}_o}{P_o} = 3 \dot{\cal M_{\rm_c}} \left( \frac{1}{\cal M_{\rm_{c}}}-\frac{1}{\cal M_{\rm_{h}}} \right)
\end{eqnarray}

\noindent
Using the aforementioned numbers and the derived stellar masses we determine $\dot{\cal M_{\rm_c}}$ $<$ 1.0 $\times$ 10$^{-7}$  $\rm M_{\odot}$/yr.

\section{Conclusions}

Based on the analysis of multi-wavelength light curves and published spectroscopic results we have arrived to the following conclusions:

\begin{itemize}

\item We find an improved orbital period  of $\rm P_{o} =  8\fd519296(8)$ that seems to be stable during the last century. This suggests that the rate of system  mass loss or mass exchange between the stars should be small.
\item We find a long cycle of length $T$ $\sim$ 270 d in the $I_c$-band data.
\item We find a system seen at angle 74$^{\circ}$\  with a stellar separation of 52 ${\rm R_{\odot}}$.
The stellar temperatures are 22\,700\,$K$ and 30\,200\,$K$, masses 9 and 17.4 ${\rm M_{\odot}}$, radii 16.9 and 7.0 ${\rm R_{\odot}}$
and surface gravities $\log g$\ = 2.94 and 3.98.
\item The system follows an horizontal path in the $W1-W2$ vs. $W2-W3$ diagram and is located in the region of systems with circumstellar matter like DPVs and Be stars.
The location in the $J-H$ vs. $H-K$ diagram also suggests infrared excess due to circumstellar matter.
\item We find
an accretion disk around the hotter star, with a radius of 21 ${\rm R_{\odot}}$. Two active regions are found, one located roughly in the expected position where
the stream impacts the disk and the other one in the opposite side of the disk. Changes are observed in the disk and spot parameters at different datasets; this might indicate variable mass transfer rate.
\item We find that the light curves show enigmatic obscuration events consisting of rapid drops of brightness of up to $\Delta$$V$ = 0.3 mag. These events might be related to inhomogeneities present in the flow of mass loss but more studies are needed to clarify their nature.

\end{itemize}

\section*{Acknowledgements}
 We thank the anonymous referee that helped to improve the first version of this manuscript.
R. E. M.,  J. G. and M. C. acknowledge support by  BASAL Centro de Astrof{\'{i}}sica y Tecnolog{\'{i}}as Afines (CATA) PFB-06/2007 and FONDECYT 1190621. JG acknowledges ANID project 21202285.
G. D., I. V., J. P. and M. I. J.  acknowledge the financial support of the Ministry of Education, Science and Technological Development of the Republic of Serbia through the contract No 451-03-68/2020/14/20002.
M. C. and P. H. acknowledge support by  Astronomical Institute of the Czech Academy
of Sciences through the project RVO 67985815.
D. K. thanks GA 17-00871S of the Czech Science Foundation. This research has made use of:
(1) the SIMBAD database, operated at CDS, Strasbourg, France
\citep{2000A&AS..143....9W}, (2) data products from the Wide-field Infrared
Survey Explorer, which is a joint project of the University of California, Los
Angeles, and the Jet Propulsion Laboratory/California Institute of Technology,
funded by the National Aeronautics and Space Administration, (3)
the NASA/ IPAC Infrared Science Archive, which
is operated by the Jet Propulsion Laboratory, California Institute of
Technology, under contract with the National Aeronautics and Space
Administration, (4) data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation and (5) data from the OMC Archive at CAB (INTA-CSIC), pre-processed by ISDC.