# Video Superresolution via Parameter-Optimized Particle Swarm Optimization

Authors: Yunyi Yan; Yujie He; Yingying Hu; Baolong Guo
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Video superresolution (VSR) aims to reconstruct a high-resolution video sequence from a low-resolution sequence. We propose a novel particle swarm optimization algorithm named as parameter-optimized multiple swarms PSO (POMS-PSO). We assessed the optimization performance of POMS-PSO by four standard benchmark functions. To reconstruct high-resolution video, we build an imaging degradation model. In view of optimization, VSR is converted to an optimization computation problem. And we take POMS-PSO as an optimization method to solve the VSR problem, which overcomes the poor effect, low accuracy, and large calculation cost in other VSR algorithms. The proposed VSR method does not require exact movement estimation and does not need the computation of movement vectors. In terms of peak signal-to-noise ratio (PSNR), sharpness, and entropy, the proposed VSR method based POMS-PSO showed better objective performance. Besides objective standard, experimental results also proved the proposed method could reconstruct high-resolution video sequence with better subjective quality.


---

## Body

## 1. Introduction

Extended and expanded from the definition of still image superresolution reconstruction proposed by Harris [1] and Goodman [2] in 1960s, video superresolution (VSR) was investigated quite intensively in recent years [3–7]. The two common types resolution algorithms to VSR are based on multiframe complementary information [8] and based on motion estimation [9, 10]. The former utilizes the redundancy information of different frames to reconstruct a high-resolution video sequence. So it is noneffective when object moves fast or no redundancy information is available. On the other hand, the latter describes movement rule of moving objects by motion vector, whose performance depends on the accuracy in motion estimation. Unfortunately, the accuracy is not very high and even very low in some cases, and yet the process is highly calculation-costed [11].Particle swarm optimization (PSO) is an established method for optimization with the advantage of simplicity of implementation [12, 13] and it has been applied in many fields including scientific research and engineering problems [12, 14–18]. Compared to other heuristic methods, like genetic algorithm (GA) and ant colony optimization (ACO), PSO is proven to have higher computational accuracy [19] and lower time cost [20] in many complicated applications.However, the parameter selection will affect the performance deeply [21], but there are no strategies proven to work well for all problems [21–24]. In many cases, we have to adjust the parameters in PSO several times to obtain satisfactory precision. The accession of subswarms with different evolutionary strategies ensured the quality both intensively and extensively [25–27]. In this paper, we introduced one swarm to determine those configuration parameters of PSO and proposed the parameter-optimized multiple swarms PSO (POMS-PSO) which could avoid the limitation of fixed parameters and achieve dynamical adjustment of parameters throughout the whole optimized process.The goal of video superresolution reconstruction (VSR) is to obtain a video with higher resolution that meets certain required specifications, for example, sharpness, succession, details, and so forth. Hence, in view of optimization, VSR can be thought as an optimal procedure, which looks for the optimal result satisfying objective and subjective quality requirements. If imaging degeneration model built, optimization technique can be used to search for the global minimum points according to proper fitness function.In this paper, we propose a novel video superresolution reconstruction(VSR) scheme which is based on image degeneration model solution via POMS-PSO. The paper is organized as follows. In Section2, the 3D video degeneration model and its solution are introduced. In Section 3, the parameter-optimized multiple swarms PSO and its performance are presented. In Section 4, we implement the video reconstruction using POMS-PSO algorithm and evaluate the result by several criteria, and finally the conclusion of this paper is given in Section 5.


## 2. Video Superresolution Model

### 2.1. 3D Video Degeneration Model

Space-time dynamic scene can be presented by coordinate system like(x,y,z,t). However in most real imaging cases, space-time dynamic scene can be presented by three-coordinate system (x,y,t) under some suitable conditions, such that (1) the scene is at one identical plane and the movement takes place within this plane; (2) the distance between various cameras is much smaller than the distance between the camera and the scene, which is to say, there is almost no parallax between two cameras. In this situation, arrange all the elements into a vector L→ according to the lexicographical order (x-y-t). In order to simplify calculations, only gray value is used to represent the characterization of pixels in images.The value onpl in a low-resolution sequence Sl can be deemed as a mapping of corresponding value on ph in a high-resolution sequence S as follows:
(1)Sl(pl)=∭p=(x,y,t)∈sup⁡⁡(Bh)‍S(p)Bh(p-ph)dp+n,
where Sl is a low-resolution sequence, pl=(xl,yl,tl) is a point in the low-resolution sequence, ph=(xh,yh,th) is the mapping of corresponding value on pl in the high-resolution sequence, and Bh is the space-time position-dependent blur kernel function in high-resolution coordinates. It consists of spatial domain point spread effect, which can be simulated by kind of Gaussian function model, and time domain fuzzy caused by integration effect by the reason of exposure time during imaging; n is white noise which is subjected to Gaussian distribution.Obviously, the range of support set of the space-time blur kernel function determines the space-time resolution relationship between the high-resolution sequence and low-resolution sequence. It determines whether the reconstruction occurs in time domain and/or spatial domain(s) and amplification factor used during reconstruction.Equation (1) is defined in continuous space which presents the map relationship between high-resolution scenes and the low-resolution ones in real world. But in digital processing system, like computers or digital signal processors (DSPs), all data need to be stored and processed in digital format. This means we have to transfer the continuous space to discrete space.The unknown continuous sceneS can be thought as the unknown high-resolution sequence Sh, and the relationship between low-resolution and high-resolution sequence can be presented by matrix computing. Hence, we obtain the discrete form of (1) as follows, one giant system composed of linear equations carrying the relation between high- and low-resolution sequence,
(2)L→=MH→+N→,
where L→ is a vector composed of all the elements in the low-resolution sequence arranging into lexicographical order, H→ is a vector composed of all the elements in the high-resolution sequence arranging into lexicographical order, N→ is a noise vector, and M is a sparse Toeplitz matrix corresponding to the space-time blur kernel function.


### 2.2. Model Solution

Like other ill-posed problems, in superresolution reconstruction, the known elements in low-resolution sequence are much less than the unknown elements in high-resolution sequence. To solve this kind of problems, extra proper criteria or constraints should be set reasonably. In this paper, a visual criteria were taken as the object optimization function, and the solution of the VSR problem can be shown as follows:(3)H*→=argmin⁡f(H→),
where H*→ is a vector composed of all the elements in rebuilt high-resolution sequence arranging into lexicographical order and f(H→) is the object optimization function known as the fitness function in r/KPSO with visual criteria taken into consideration, and it is shown as below:
(4)f(H→)=∥MH→-L→∥2+Sm(H→)x,y,t+Dt(H→)x,y.∥
M
H
→
-
L
→
∥
2 reflects fidelity, and Sm(H→)x,y,t reflects smoothness in directions of x-y-t, and more specifically, it is defined by the following equation:
(5)Sm(H→)x,y,t=λx∥WxGx(H→)∥2+λy∥WyGy(H→)∥2+λt∥WtGt(H→)∥2.λ
i
(
i
=
x
,
y
,
t
) is the image smooth coefficients in i direction, Wi(i=x,y,t) is the weight matrix in i direction, and Gi(i=x,y,t) is the second-order differential operator in i direction, and Dt(H→)x,y reflects details of the retention, which can be computed as follows:
(6)Dt(H→)x,y=(∥Dx(H→)∥2+∥Dy(H→)∥2+ɛ)-1.D
i
(
H
→
)
(
i
=
x
,
y
) is the gradient operator in i direction; ɛ is a small positive constant.The schematic diagram of the reconstruction procedure in this paper is given in Figure1. The input low definition video is transformed into a sequence, and then all elements are arranged in lexicographical order to compose a matrix. Based on the video degeneration model, there is a map between the low-resolution sequence and the high-resolution sequence.Figure 1
Schematic diagram of a reconstruction procedure.In Figure1, POMS-PSO is used to search the best reconstructed results which is represented as H*→ in (3). In fact, we code the reconstructed high definition video as a particle of subswarm in the inner loop of Figure 2. And the f(H→) in (4) is used as the fitness function during the optimization process. The best solution to (4) will be thought as reconstructed video. The utilization of POMS-PSO algorithm changes the reconstruction procedure into a problem that search a solution satisfying the video degeneration process in best.Figure 2
Schematic diagram of POMS-PSO.


## 2.1. 3D Video Degeneration Model

Space-time dynamic scene can be presented by coordinate system like(x,y,z,t). However in most real imaging cases, space-time dynamic scene can be presented by three-coordinate system (x,y,t) under some suitable conditions, such that (1) the scene is at one identical plane and the movement takes place within this plane; (2) the distance between various cameras is much smaller than the distance between the camera and the scene, which is to say, there is almost no parallax between two cameras. In this situation, arrange all the elements into a vector L→ according to the lexicographical order (x-y-t). In order to simplify calculations, only gray value is used to represent the characterization of pixels in images.The value onpl in a low-resolution sequence Sl can be deemed as a mapping of corresponding value on ph in a high-resolution sequence S as follows:
(1)Sl(pl)=∭p=(x,y,t)∈sup⁡⁡(Bh)‍S(p)Bh(p-ph)dp+n,
where Sl is a low-resolution sequence, pl=(xl,yl,tl) is a point in the low-resolution sequence, ph=(xh,yh,th) is the mapping of corresponding value on pl in the high-resolution sequence, and Bh is the space-time position-dependent blur kernel function in high-resolution coordinates. It consists of spatial domain point spread effect, which can be simulated by kind of Gaussian function model, and time domain fuzzy caused by integration effect by the reason of exposure time during imaging; n is white noise which is subjected to Gaussian distribution.Obviously, the range of support set of the space-time blur kernel function determines the space-time resolution relationship between the high-resolution sequence and low-resolution sequence. It determines whether the reconstruction occurs in time domain and/or spatial domain(s) and amplification factor used during reconstruction.Equation (1) is defined in continuous space which presents the map relationship between high-resolution scenes and the low-resolution ones in real world. But in digital processing system, like computers or digital signal processors (DSPs), all data need to be stored and processed in digital format. This means we have to transfer the continuous space to discrete space.The unknown continuous sceneS can be thought as the unknown high-resolution sequence Sh, and the relationship between low-resolution and high-resolution sequence can be presented by matrix computing. Hence, we obtain the discrete form of (1) as follows, one giant system composed of linear equations carrying the relation between high- and low-resolution sequence,
(2)L→=MH→+N→,
where L→ is a vector composed of all the elements in the low-resolution sequence arranging into lexicographical order, H→ is a vector composed of all the elements in the high-resolution sequence arranging into lexicographical order, N→ is a noise vector, and M is a sparse Toeplitz matrix corresponding to the space-time blur kernel function.


## 2.2. Model Solution

Like other ill-posed problems, in superresolution reconstruction, the known elements in low-resolution sequence are much less than the unknown elements in high-resolution sequence. To solve this kind of problems, extra proper criteria or constraints should be set reasonably. In this paper, a visual criteria were taken as the object optimization function, and the solution of the VSR problem can be shown as follows:(3)H*→=argmin⁡f(H→),
where H*→ is a vector composed of all the elements in rebuilt high-resolution sequence arranging into lexicographical order and f(H→) is the object optimization function known as the fitness function in r/KPSO with visual criteria taken into consideration, and it is shown as below:
(4)f(H→)=∥MH→-L→∥2+Sm(H→)x,y,t+Dt(H→)x,y.∥
M
H
→
-
L
→
∥
2 reflects fidelity, and Sm(H→)x,y,t reflects smoothness in directions of x-y-t, and more specifically, it is defined by the following equation:
(5)Sm(H→)x,y,t=λx∥WxGx(H→)∥2+λy∥WyGy(H→)∥2+λt∥WtGt(H→)∥2.λ
i
(
i
=
x
,
y
,
t
) is the image smooth coefficients in i direction, Wi(i=x,y,t) is the weight matrix in i direction, and Gi(i=x,y,t) is the second-order differential operator in i direction, and Dt(H→)x,y reflects details of the retention, which can be computed as follows:
(6)Dt(H→)x,y=(∥Dx(H→)∥2+∥Dy(H→)∥2+ɛ)-1.D
i
(
H
→
)
(
i
=
x
,
y
) is the gradient operator in i direction; ɛ is a small positive constant.The schematic diagram of the reconstruction procedure in this paper is given in Figure1. The input low definition video is transformed into a sequence, and then all elements are arranged in lexicographical order to compose a matrix. Based on the video degeneration model, there is a map between the low-resolution sequence and the high-resolution sequence.Figure 1
Schematic diagram of a reconstruction procedure.In Figure1, POMS-PSO is used to search the best reconstructed results which is represented as H*→ in (3). In fact, we code the reconstructed high definition video as a particle of subswarm in the inner loop of Figure 2. And the f(H→) in (4) is used as the fitness function during the optimization process. The best solution to (4) will be thought as reconstructed video. The utilization of POMS-PSO algorithm changes the reconstruction procedure into a problem that search a solution satisfying the video degeneration process in best.Figure 2
Schematic diagram of POMS-PSO.


## 3. Parameter-Optimized Multiple Swarms PSO (POMS-PSO)

Particle swarm optimization (PSO) has surprising ability of handle optimization problems with multiple local optima and its simplicity of implementation. We combined the concepts of multiswarm and parameter autooptimizing, and then we proposed the improvement of multiswarm PSO named as parameter-optimized PSO.In fact, many improved PSO algorithms [15, 28–30] have been given with fixed parameters [22, 31]. The parameter selection can affect the performance of PSO [21], but there are no strategies proven to work well for all problems [21–24]. In many cases, when facing real engineering problems, we have to adjust the parameters in PSO several times to obtain satisfactory precision. The concept of parameter-optimizing is to determine those free parameters of the subswarms via PSO program. In this paper, the concept of multi-subswarm is based on an ecological approach [25, 26], in which all the particles are divided into two different subswarms depending on their fitness. By adjusting parameters dynamically throughout the whole optimized process, the method can avoid the limitation of fixed parameters. At the same time, those parameters can automatically change when problem to be solved is changed. By introducing different evolutionary strategies using on different subswarms [27], the new framework ensures quality both intensively and extensively. All the advantages above over current methods result in better optimization performance than parameter-fixed PSO methods [32, 33].


### 3.1. Algorithm Introduction

Initialization of the particles which are used to search for the global optimum is spanned as follows:(7)(x11⋯x1nv11⋯v1nx21⋯x2nv11⋯v2n⋮⋱⋮⋮⋱⋮xN1⋯xNnvN1⋯vNn),
where N and n are the size and the dimension of particle swarm, respectively. And {xij,i=1,2,…,N;j=1,2,…,n} refers to the ith particle with position j, and its corresponding speed is vij. General speaking, one particle’s iteration and update of PSO algorithm can be presented by the following two equations:
(8)vidn+1=ωvidn+c1r1(pidn-xidn)+c2r2(pgdn-xidn),xidn+1=xidn+vidn+1,
where d is dimensions index of particle’s position vector; i is particle index; n notes for iteration number; c1 and c2 are positive constants calledcognitive andsocial factor, respectively; ω is calledinertia weight; r1 and r2 are random numbers; pi is the best position visited by this particle; and pg is the global best position over the whole swarm. More details on PSO model and its implementation can be found in [28].We proposed a novel optimization method with a feature that the parameters in PSO can be determined and optimized by its self, which was named as POMS-PSO. To accomplish this goal, we introduced two layer loops into the iteration process. In outer loop, we placed one swarm to optimize the parameter setting for the two subswarms in inner loop. And the inner loop is used to find the optimal resolution to given problem. The proposed method is named as POMS-PSO, which comes from parameter-optimized multiple swarms PSO. The schematic diagram of PMOS-PSO is shown in Figure2.In the view of particle swarm optimization, a particle means a potential solution to problem. The PSO algorithm is inspired by the bird or fish swarms’ seeking food procedure in nature. In PSO algorithm, a bird is modelled as a particle, and a bird position is presented by a vector. In Figure2, POMS-PSO employs double layer loops and three swarms. The swarm in outer loop is named as C-swarm and the particles in C-swam is called C-particle which is based on the fact that the C-particle’s valued is used to configure the two subswarms’ parameter selection (in inner loop) and the inner loop’s computing is performed with a certain C-particle’s configuration to find the best solution to given problem. In the proposed POMS-PSO, a C-particle means a parameter setting for the inner loop and a particle in inner loop means a solution to given problem.As shown in Figure2, the first step of optimization initializes the particles with random numbers normally distributed in the search space. The C-swarm’s search space is different from the r- and K-subswarm. For example, the search space for the four benchmark functions in Section 3.2 is [-100,100]D, D=30, but the C-swarm’s is recommended to set in [0,10]D, D=11 due to empirical information. In the step of fitness evaluation, all particles in inner loop are evaluated by given fitness function and are sorted in descend. The top (1-Pr)×100% particles are allocated to K-subswarm and the other to r-subswarm. For example, we totally have 20 particles in inner loop, and Pr is set to 0.1, then 18 out of 20 particles belong to K-subswarm, and 2 out of 20 belong to r-subswarm. The r- and K-subswarms will perform r- and K-selection strategies, respectively, and update their positions.In most of all optimization based on iteration, the stop criterion can be defined as the iteration could satisfy one of the following two conditions. (1) The maximum generation is reached and (2) the error or the fitness value is good enough to satisfy our requirement. To our experimental experience, when POMS-PSO is used to solve complicated problems (e.g., VSR), the error threshold is difficult to estimate and we often just set it to zero. In this case, the inner loop is often stopped by the maximum iteration generation reached unless our luck is good enough to get a solution with zero error.For every C-particle, all of its various positions mean various configuration of inner loop. The inner loop deploys two subswarms, the r-subswarm and the K-subswarm. In ecology, r-selection is termed for those species that breed many offspring and live in unstable environments and K-selection for those species that produce few offspring and live in stable environments. The r-selection can be characterized as quantitative, little parent care, large growth rate, and rapid development and K-selection can be characterized as qualitative, much parent care, small growth rate, and slow development. K-selection is performed for those particles with high fitness, and so K-subswarm only can produce relatively fewer progenies but the progenies are nurtured delicately with much parent care. On the other hand, r-selection is performed for those particles with relatively lower fitness. Due to little parent care, r-subswarm can produce a large number of progenies; the progenies have to compete with other members in r-subswarm for survival and only the best ones can survive. The main task of r-subswarm is to explore the search space as possible as they can, and K-subswarm should try to keep the current optimum solutions and exploit the space as they can.In fact, (8) are the core procedure in most of all PSO algorithms. The r-selection and K-selection are defined in ecology to describe various species’ idiosyncrasy; however, in our POMS-PSO algorithm, the two strategies are defined by different parameter setting.Judging from (8), we can find there are three parameters that could affect optimization performance, theinertia weightω, thecognitive factor c1, and thesocial factor c2. Due to their importance to PSO, these three parameters of r-subswarm together with those of K-subswarm are selected to be part of the C-particle. And POMS-PSO will optimize these six parameters.For the inner loop, we divide the whole swarm into r-subswarm and K-subswarm as shown in Figure2. That is why we call them not swarm but subswarm. Now we have to determine how many particles divided into r-subswarm. In other words, the proportion of r-subswarm Pr should be in the range of [0,1], but it needs to be determined. And on the other hand, the proportion of K-subswarm will be (1-Pr).As the r-selection is likely to produce more progenies than K-selection, the fertility rate of r-subswarmρr is larger than that of K-subswarm ρK. So we have
(9)ρr>ρK≥1.Three steps in Figure2, the r- and K-selection and update particles, are presented in details in Figure 3. Suppose the total particle number of r- and K-subswarms is N as shown in Figure 3; after one evolutionary generation, the total number of progenies will be
(10)Np=ρK×N×Pr+ρr×N×(1-Pr).Figure 3
Schematic diagram of K-selection and r-selection.Considering (9), we can draw a conclusion that NP≥N. To keep stability of r- and K-subswarms size, we only select the best N out of Np to survive to next iteration in inner loop and the others are abandoned. And we have to notice that like the classic PSO, C-particles only have one progeny.Another important difference between r- and K-subswarm is that the velocity in r-subswarm may be much higher than the K-subswarm, because r-subswarm is designed to explore the search space to the greatest extent. The maximum velocity of r-subswarm,Tr, should be larger than that of K-subswarm TK. These two parameters, Tr and TK, should be optimized by C-swarm in POMS-PSO. Totally we have eleven parameters for r-subswarm and K-subswarm to be determined. Table 1 lists these parameters and their meanings. All of these parameters will be optimized in outer loop of POMS-PSO. The optimization procedure will be running with these eleven parameters’ configuration. Although the parameter setting of r- and K-subswarm needs to be determined by the C-swarm in POMS-PSO, we have to notice that C-swarm perform the iteration procedure using the classic setting as given in [28].Table 1
Eleven parameters optimized by outer loop in POMS-PSO.
Number
Symbols
Meanings
1
ω
K
I
n
e
r
t
i
a weight of K-subswarm
2
c
K
1
C
o
g
n
i
t
i
v
e factor of K-subswarm
3
c
K
2
S
o
c
i
a
l factor of K-subswarm
4
ω
r
I
n
e
r
t
i
a
w
e
i
g
h
t of r-subswarm
5
c
r
1
C
o
g
n
i
t
i
v
e factor of r-subswarm
6
c
r
2
S
o
c
i
a
l factor of r-subswarm
7
ρ
K
Fertility rate of K-subswarm
8
ρ
r
Fertility rate of r-subswarm
9
P
r
Proportion of r-subswarm
10
T
r
Maximum velocity of r-subswarm
11
T
K
Maximum velocity of K-subswarmAs the inner loop is run under the configuration of a certain C-particle, the best fitness of inner loop got can be used to evaluate the corresponding C-particle. As shown in Figure2, we take the best particle inner loop as the validation of C-particle. By evolution and update in the outer loop, POMS-PSO can find the best C-particle. And for inner loop, r- and K-subswarms can evolute in best parameter settings.


### 3.2. Performance Evaluation of POMS-PSO

We assessed the performance of parameter-optimized PSO method by employing a suite of benchmark functions and compared the result to that of the standard PSO (SPSO) and the constriction type PSO (CPSO) methods. The benchmark suite consists of two unimodal functions (Tablet and Quadric) and two multimodal functions (Rastrigin and Schaffer) [28]. All functions except for the 2D Schaffer function are optimized in 30 dimensional spaces.The four functions are presented as follows:Schaffer:(11)f(x)=0.5+sin2x12+x22-0.51+10-3(x12+x22)2,x∈[-100,100]2.Tablet:(12)f(x)=106x12+∑i=2D‍xi2,x∈[-100,100]D,D=30.Rastrigin:(13)f(x)=10D+∑i=1D‍(xi2-10cos⁡(2πxi)),x∈[-100,100]D,D=30.Quadric:(14)f(x)=∑i=1D‍(∑j=1i‍xj)2,x∈[-100,100]D,D=30.
The four functions are often used in optimization performance evaluation [28]. To test parameter-optimized PSO’s optimal ability, we used the search space [-100,100]D which was larger than those in [28]. The Tablet and Quadric functions are evolved from spherical function, which have only one minimum point in the search space. And the Schaffer and the Rastrigin are multimodal, which means they have multiple local minimum points in the space. All of the benchmark functions have the global minimum f(x)=0 at xi=0. Unimodal and multimodal should be both evaluated for complete comparison performance assessment. General speaking, multimodal functions are much more difficult to find the optimal point.As mentioned in Section3.1, POMS-PSO employs 3 swarms in two loops, that is, the r-subswam, the K-subswarm (in inner loop), and the C-swarm (in outer loop). One task of parameter-optimized PSO is to find the optimal swarm parameters setting for r- and K-subswarms, which is performed by the C-swarm. In our experiments, one population of 20 particles was chosen for r- and K-subswarms and 10 for C-swarm. The maximum number of iterations of outer loop is set to 40 and inner loop is set to 20. Table 2 lists the optimal results of C-particles of the four benchmark functions.Table 2
Optimized swarm parameters for four benchmark functions.
Sch
Tab
Ras
Qua
ω
K
1
2.97
2.64
1.09
1.13
c
K
1
3.35
2.75
3.65
2.78
c
K
2
0.14
2.75
2.74
1.65
ω
r
1
1.05
0.97
1.81
2.74
c
r
1
3.86
2.54
2.87
2.70
c
r
2
1.94
1.92
3.38
2.26
ρ
K
2.83
3.31
1.24
1.36
ρ
r
7.76
3.79
8.69
9.21
P
r
0.74
0.01
0.01
0.01
T
r
0.56
0.77
1.00
1.00
T
K
0.07
0.36
0.01
0.01The parameter-optimized PSO can perform optimization for every benchmark function’s parameter selection to guarantee various functions can get the optimal or better solution under various configurations instead of the common used settings. The parameter-optimized PSO is a problem-oriented PSO solution. In fact, we noticed that even we run the outer loop only once, we could obtain the results comparable to SPSO and CPSO in our experiments. In Table2, all of the eleven parameters are shown for Sch (Schaffer), Tab (Tablet), Ras (Rastrigin), and Qua (Quadric), respectively. It is obvious that various functions have various configuration.The performance comparison between parameter-optimized PSO, SPSO, and CPSO methods is given in Table3 for four benchmark functions. The optimization performance comparison is based on 200 trials. We compared the minimum, maximum, mean, and standard deviation of those three methods. The statistical results are listed when the SPSO, CPSO, and parameter-optimized PSO run 200 times. The size of all the populations is 20. The maximum evolution generation of SPSO and CPSO is 800, while that of POMS-PSO is also 40×20=800 (outer loop: 40, inner loop: 20).Table 3
Performance comparison of SPSO, CPSO, and POMS-PSO implementation. Maximum evolution generation of SPSO and CPSO: 800. Maximum evolution generation of POMS-PSO outer loop: 40, inner loop: 20.
Functions
Schaffer
Tablet
Rastrigin
Quadric
Standard PSO
Mean
8.61
E
-
04
1.57
E
+
03
1.21
E
+
03
3.29
E
+
02
Standard deviation
8.56
E
-
04
9.10
E
+
02
3.84
E
+
02
3.08
E
+
02
Minimum
4.88
E
-
15
1.84
E
+
02
4.54
E
+
02
3.51
E
+
01
Median
9.85
E
-
04
1.39
E
+
03
1.18
E
+
03
1.17
E
+
02
Maximum
8.76
E
-
03
6.11
E
+
03
2.92
E
+
03
2.13
E
+
03
CPSO
Mean
8.43
E
-
04
2.52
E
+
03
1.93
E
+
03
9.69
E
+
02
Standard deviation
7.60
E
-
04
1.19
E
+
03
6.66
E
+
02
6.50
E
+
02
Minimum
0.00
E
+
00
4.01
E
+
02
7.99
E
+
02
1.06
E
+
02
Median
9.85
E
-
04
2.35
E
+
03
1.84
E
+
03
7.95
E
+
02
Maximum
8.76
E
-
03
7.38
E
+
03
5.24
E
+
03
4.71
E
+
03
POMS-PSO
Mean
2.14
E
-
05
8.61
E
+
02
4.40
E
+
02
1.25
E
+
02
Standard deviation
2.80
E
-
05
3.30
E
+
03
2.39
E
+
02
2.20
E
+
02
Minimum
9.01
E
-
12
5.59
E
-
04
4.86
E
+
01
3.20
E
-
03
Median
1.08
E
-
05
6.72
E
+
01
4.09
E
+
02
1.59
E
+
01
Maximum
1.96
E
-
04
3.80
E
+
04
9.93
E
+
02
9.45
E
+
02From the values listed in Table3, we can see that at most time, the POMS-PSO can outperform the other two, no matter the minimum, maximum, mean, or standard deviation values, although POMS-PSO has much smaller iteration times. The performance of the POMS-PSO method decreases the errors by one order of magnitude compared with the other two methods.As shown in Figure4, multiruns experiments show POMS-PSO results vary in much smaller range intervals which means the proposed methods can give more stable optimum. Due to POMS-PSO’s precision improvement under multi-subswarms strategy framework, it is a reasonable choice to adapt it into video superresolution reconstruction by solving the video degeneration model showed in (2), (3), and (4).The convergence procedure of four functions between POMS-PSO, SPSO, and CPSO methods (The corresponding evolution generation is outer loop generation for parameter-optimized PSO andGeneration× 20 for SPSO and CPSO method.).
(a)
Schaffer
(b)
Tablet
(c)
Rastrigin
(d)
Quadric


## 3.1. Algorithm Introduction

Initialization of the particles which are used to search for the global optimum is spanned as follows:(7)(x11⋯x1nv11⋯v1nx21⋯x2nv11⋯v2n⋮⋱⋮⋮⋱⋮xN1⋯xNnvN1⋯vNn),
where N and n are the size and the dimension of particle swarm, respectively. And {xij,i=1,2,…,N;j=1,2,…,n} refers to the ith particle with position j, and its corresponding speed is vij. General speaking, one particle’s iteration and update of PSO algorithm can be presented by the following two equations:
(8)vidn+1=ωvidn+c1r1(pidn-xidn)+c2r2(pgdn-xidn),xidn+1=xidn+vidn+1,
where d is dimensions index of particle’s position vector; i is particle index; n notes for iteration number; c1 and c2 are positive constants calledcognitive andsocial factor, respectively; ω is calledinertia weight; r1 and r2 are random numbers; pi is the best position visited by this particle; and pg is the global best position over the whole swarm. More details on PSO model and its implementation can be found in [28].We proposed a novel optimization method with a feature that the parameters in PSO can be determined and optimized by its self, which was named as POMS-PSO. To accomplish this goal, we introduced two layer loops into the iteration process. In outer loop, we placed one swarm to optimize the parameter setting for the two subswarms in inner loop. And the inner loop is used to find the optimal resolution to given problem. The proposed method is named as POMS-PSO, which comes from parameter-optimized multiple swarms PSO. The schematic diagram of PMOS-PSO is shown in Figure2.In the view of particle swarm optimization, a particle means a potential solution to problem. The PSO algorithm is inspired by the bird or fish swarms’ seeking food procedure in nature. In PSO algorithm, a bird is modelled as a particle, and a bird position is presented by a vector. In Figure2, POMS-PSO employs double layer loops and three swarms. The swarm in outer loop is named as C-swarm and the particles in C-swam is called C-particle which is based on the fact that the C-particle’s valued is used to configure the two subswarms’ parameter selection (in inner loop) and the inner loop’s computing is performed with a certain C-particle’s configuration to find the best solution to given problem. In the proposed POMS-PSO, a C-particle means a parameter setting for the inner loop and a particle in inner loop means a solution to given problem.As shown in Figure2, the first step of optimization initializes the particles with random numbers normally distributed in the search space. The C-swarm’s search space is different from the r- and K-subswarm. For example, the search space for the four benchmark functions in Section 3.2 is [-100,100]D, D=30, but the C-swarm’s is recommended to set in [0,10]D, D=11 due to empirical information. In the step of fitness evaluation, all particles in inner loop are evaluated by given fitness function and are sorted in descend. The top (1-Pr)×100% particles are allocated to K-subswarm and the other to r-subswarm. For example, we totally have 20 particles in inner loop, and Pr is set to 0.1, then 18 out of 20 particles belong to K-subswarm, and 2 out of 20 belong to r-subswarm. The r- and K-subswarms will perform r- and K-selection strategies, respectively, and update their positions.In most of all optimization based on iteration, the stop criterion can be defined as the iteration could satisfy one of the following two conditions. (1) The maximum generation is reached and (2) the error or the fitness value is good enough to satisfy our requirement. To our experimental experience, when POMS-PSO is used to solve complicated problems (e.g., VSR), the error threshold is difficult to estimate and we often just set it to zero. In this case, the inner loop is often stopped by the maximum iteration generation reached unless our luck is good enough to get a solution with zero error.For every C-particle, all of its various positions mean various configuration of inner loop. The inner loop deploys two subswarms, the r-subswarm and the K-subswarm. In ecology, r-selection is termed for those species that breed many offspring and live in unstable environments and K-selection for those species that produce few offspring and live in stable environments. The r-selection can be characterized as quantitative, little parent care, large growth rate, and rapid development and K-selection can be characterized as qualitative, much parent care, small growth rate, and slow development. K-selection is performed for those particles with high fitness, and so K-subswarm only can produce relatively fewer progenies but the progenies are nurtured delicately with much parent care. On the other hand, r-selection is performed for those particles with relatively lower fitness. Due to little parent care, r-subswarm can produce a large number of progenies; the progenies have to compete with other members in r-subswarm for survival and only the best ones can survive. The main task of r-subswarm is to explore the search space as possible as they can, and K-subswarm should try to keep the current optimum solutions and exploit the space as they can.In fact, (8) are the core procedure in most of all PSO algorithms. The r-selection and K-selection are defined in ecology to describe various species’ idiosyncrasy; however, in our POMS-PSO algorithm, the two strategies are defined by different parameter setting.Judging from (8), we can find there are three parameters that could affect optimization performance, theinertia weightω, thecognitive factor c1, and thesocial factor c2. Due to their importance to PSO, these three parameters of r-subswarm together with those of K-subswarm are selected to be part of the C-particle. And POMS-PSO will optimize these six parameters.For the inner loop, we divide the whole swarm into r-subswarm and K-subswarm as shown in Figure2. That is why we call them not swarm but subswarm. Now we have to determine how many particles divided into r-subswarm. In other words, the proportion of r-subswarm Pr should be in the range of [0,1], but it needs to be determined. And on the other hand, the proportion of K-subswarm will be (1-Pr).As the r-selection is likely to produce more progenies than K-selection, the fertility rate of r-subswarmρr is larger than that of K-subswarm ρK. So we have
(9)ρr>ρK≥1.Three steps in Figure2, the r- and K-selection and update particles, are presented in details in Figure 3. Suppose the total particle number of r- and K-subswarms is N as shown in Figure 3; after one evolutionary generation, the total number of progenies will be
(10)Np=ρK×N×Pr+ρr×N×(1-Pr).Figure 3
Schematic diagram of K-selection and r-selection.Considering (9), we can draw a conclusion that NP≥N. To keep stability of r- and K-subswarms size, we only select the best N out of Np to survive to next iteration in inner loop and the others are abandoned. And we have to notice that like the classic PSO, C-particles only have one progeny.Another important difference between r- and K-subswarm is that the velocity in r-subswarm may be much higher than the K-subswarm, because r-subswarm is designed to explore the search space to the greatest extent. The maximum velocity of r-subswarm,Tr, should be larger than that of K-subswarm TK. These two parameters, Tr and TK, should be optimized by C-swarm in POMS-PSO. Totally we have eleven parameters for r-subswarm and K-subswarm to be determined. Table 1 lists these parameters and their meanings. All of these parameters will be optimized in outer loop of POMS-PSO. The optimization procedure will be running with these eleven parameters’ configuration. Although the parameter setting of r- and K-subswarm needs to be determined by the C-swarm in POMS-PSO, we have to notice that C-swarm perform the iteration procedure using the classic setting as given in [28].Table 1
Eleven parameters optimized by outer loop in POMS-PSO.
Number
Symbols
Meanings
1
ω
K
I
n
e
r
t
i
a weight of K-subswarm
2
c
K
1
C
o
g
n
i
t
i
v
e factor of K-subswarm
3
c
K
2
S
o
c
i
a
l factor of K-subswarm
4
ω
r
I
n
e
r
t
i
a
w
e
i
g
h
t of r-subswarm
5
c
r
1
C
o
g
n
i
t
i
v
e factor of r-subswarm
6
c
r
2
S
o
c
i
a
l factor of r-subswarm
7
ρ
K
Fertility rate of K-subswarm
8
ρ
r
Fertility rate of r-subswarm
9
P
r
Proportion of r-subswarm
10
T
r
Maximum velocity of r-subswarm
11
T
K
Maximum velocity of K-subswarmAs the inner loop is run under the configuration of a certain C-particle, the best fitness of inner loop got can be used to evaluate the corresponding C-particle. As shown in Figure2, we take the best particle inner loop as the validation of C-particle. By evolution and update in the outer loop, POMS-PSO can find the best C-particle. And for inner loop, r- and K-subswarms can evolute in best parameter settings.


## 3.2. Performance Evaluation of POMS-PSO

We assessed the performance of parameter-optimized PSO method by employing a suite of benchmark functions and compared the result to that of the standard PSO (SPSO) and the constriction type PSO (CPSO) methods. The benchmark suite consists of two unimodal functions (Tablet and Quadric) and two multimodal functions (Rastrigin and Schaffer) [28]. All functions except for the 2D Schaffer function are optimized in 30 dimensional spaces.The four functions are presented as follows:Schaffer:(11)f(x)=0.5+sin2x12+x22-0.51+10-3(x12+x22)2,x∈[-100,100]2.Tablet:(12)f(x)=106x12+∑i=2D‍xi2,x∈[-100,100]D,D=30.Rastrigin:(13)f(x)=10D+∑i=1D‍(xi2-10cos⁡(2πxi)),x∈[-100,100]D,D=30.Quadric:(14)f(x)=∑i=1D‍(∑j=1i‍xj)2,x∈[-100,100]D,D=30.
The four functions are often used in optimization performance evaluation [28]. To test parameter-optimized PSO’s optimal ability, we used the search space [-100,100]D which was larger than those in [28]. The Tablet and Quadric functions are evolved from spherical function, which have only one minimum point in the search space. And the Schaffer and the Rastrigin are multimodal, which means they have multiple local minimum points in the space. All of the benchmark functions have the global minimum f(x)=0 at xi=0. Unimodal and multimodal should be both evaluated for complete comparison performance assessment. General speaking, multimodal functions are much more difficult to find the optimal point.As mentioned in Section3.1, POMS-PSO employs 3 swarms in two loops, that is, the r-subswam, the K-subswarm (in inner loop), and the C-swarm (in outer loop). One task of parameter-optimized PSO is to find the optimal swarm parameters setting for r- and K-subswarms, which is performed by the C-swarm. In our experiments, one population of 20 particles was chosen for r- and K-subswarms and 10 for C-swarm. The maximum number of iterations of outer loop is set to 40 and inner loop is set to 20. Table 2 lists the optimal results of C-particles of the four benchmark functions.Table 2
Optimized swarm parameters for four benchmark functions.
Sch
Tab
Ras
Qua
ω
K
1
2.97
2.64
1.09
1.13
c
K
1
3.35
2.75
3.65
2.78
c
K
2
0.14
2.75
2.74
1.65
ω
r
1
1.05
0.97
1.81
2.74
c
r
1
3.86
2.54
2.87
2.70
c
r
2
1.94
1.92
3.38
2.26
ρ
K
2.83
3.31
1.24
1.36
ρ
r
7.76
3.79
8.69
9.21
P
r
0.74
0.01
0.01
0.01
T
r
0.56
0.77
1.00
1.00
T
K
0.07
0.36
0.01
0.01The parameter-optimized PSO can perform optimization for every benchmark function’s parameter selection to guarantee various functions can get the optimal or better solution under various configurations instead of the common used settings. The parameter-optimized PSO is a problem-oriented PSO solution. In fact, we noticed that even we run the outer loop only once, we could obtain the results comparable to SPSO and CPSO in our experiments. In Table2, all of the eleven parameters are shown for Sch (Schaffer), Tab (Tablet), Ras (Rastrigin), and Qua (Quadric), respectively. It is obvious that various functions have various configuration.The performance comparison between parameter-optimized PSO, SPSO, and CPSO methods is given in Table3 for four benchmark functions. The optimization performance comparison is based on 200 trials. We compared the minimum, maximum, mean, and standard deviation of those three methods. The statistical results are listed when the SPSO, CPSO, and parameter-optimized PSO run 200 times. The size of all the populations is 20. The maximum evolution generation of SPSO and CPSO is 800, while that of POMS-PSO is also 40×20=800 (outer loop: 40, inner loop: 20).Table 3
Performance comparison of SPSO, CPSO, and POMS-PSO implementation. Maximum evolution generation of SPSO and CPSO: 800. Maximum evolution generation of POMS-PSO outer loop: 40, inner loop: 20.
Functions
Schaffer
Tablet
Rastrigin
Quadric
Standard PSO
Mean
8.61
E
-
04
1.57
E
+
03
1.21
E
+
03
3.29
E
+
02
Standard deviation
8.56
E
-
04
9.10
E
+
02
3.84
E
+
02
3.08
E
+
02
Minimum
4.88
E
-
15
1.84
E
+
02
4.54
E
+
02
3.51
E
+
01
Median
9.85
E
-
04
1.39
E
+
03
1.18
E
+
03
1.17
E
+
02
Maximum
8.76
E
-
03
6.11
E
+
03
2.92
E
+
03
2.13
E
+
03
CPSO
Mean
8.43
E
-
04
2.52
E
+
03
1.93
E
+
03
9.69
E
+
02
Standard deviation
7.60
E
-
04
1.19
E
+
03
6.66
E
+
02
6.50
E
+
02
Minimum
0.00
E
+
00
4.01
E
+
02
7.99
E
+
02
1.06
E
+
02
Median
9.85
E
-
04
2.35
E
+
03
1.84
E
+
03
7.95
E
+
02
Maximum
8.76
E
-
03
7.38
E
+
03
5.24
E
+
03
4.71
E
+
03
POMS-PSO
Mean
2.14
E
-
05
8.61
E
+
02
4.40
E
+
02
1.25
E
+
02
Standard deviation
2.80
E
-
05
3.30
E
+
03
2.39
E
+
02
2.20
E
+
02
Minimum
9.01
E
-
12
5.59
E
-
04
4.86
E
+
01
3.20
E
-
03
Median
1.08
E
-
05
6.72
E
+
01
4.09
E
+
02
1.59
E
+
01
Maximum
1.96
E
-
04
3.80
E
+
04
9.93
E
+
02
9.45
E
+
02From the values listed in Table3, we can see that at most time, the POMS-PSO can outperform the other two, no matter the minimum, maximum, mean, or standard deviation values, although POMS-PSO has much smaller iteration times. The performance of the POMS-PSO method decreases the errors by one order of magnitude compared with the other two methods.As shown in Figure4, multiruns experiments show POMS-PSO results vary in much smaller range intervals which means the proposed methods can give more stable optimum. Due to POMS-PSO’s precision improvement under multi-subswarms strategy framework, it is a reasonable choice to adapt it into video superresolution reconstruction by solving the video degeneration model showed in (2), (3), and (4).The convergence procedure of four functions between POMS-PSO, SPSO, and CPSO methods (The corresponding evolution generation is outer loop generation for parameter-optimized PSO andGeneration× 20 for SPSO and CPSO method.).
(a)
Schaffer
(b)
Tablet
(c)
Rastrigin
(d)
Quadric


## 4. Video Superresolution by POMS-PSO

### 4.1. Reconstruction Quality Evaluation Criteria

#### 4.1.1. Similarity-Based Evaluation Criteria

The PSNR is an engineering term for the ratio between the maximum possible power of a signal and the power of corrupting noise that affects the fidelity of its representation. The PSNR is most commonly used as a measure of quality of reconstruction in image compression. It is defined via the mean squared error (MSE) as follows:(15)MSE=1mn∑i=0m-1‍∑j=0n-1‍∥I(i,j)-K(i,j)∥2.I and K are two m×n monochrome images where one of the images is considered a noisy approximation of the other.The PSNR is defined as follows:(16)PSNR=10log⁡(2552MSE)(dB).


#### 4.1.2. Detail-Based Evaluation Criteria

Sharpness is a photograph term that reflects image’s detail as an integral part of its appeal. Sharpness is defined by the boundaries between zones of different tones or colors. It is illustrated by the bar pattern of increasing spatial frequency.As shown in Figure5, the top portion is sharp and its boundaries are crisp steps, not gradual. The bottom portion illustrates how the degraded pattern is blurred.Figure 5
Bar pattern: original (top); with degradation (bottom).The sharpness is most commonly defined as the following equation for onem×n monochrome image I:
(17)sharpness=22(m-1)(n-1)1111111111×∑i=1m-1‍∑j=1n-1‍(Ii,j-Ii+1,j)2+(Ii,j-Ii,j+1)2.
The sharpness lines indicate the average gradient values of all frames in the video, which express whether the variation of scene is large or not in adjacent pixels.


#### 4.1.3. Information-Based Evaluation Criteria

Image entropy reflects the average amount of information for a given image and the entropy of grayscale image indicates the amount of information contained in aggregation feature of gray distribution for the given image and it can be computed as follows:(18)Entropy=-∑i=0M-1‍Pilog⁡2(Pi),
where M is the number of gray levels and Pi is the probability associated with gray level i.


### 4.2. Reconstruction Simulation and Experimental Results

We assessed the performance of the video superresolution process using POMS-PSO algorithm on several low-resolution video sequences and compared the result to that of the bilinear interpolation method. The resolution of the reconstructed video is twice as large as incoming video, in both height and width. To compare the objective criteria, the input video used to perform superresolution procedure was obtained by 2-downsampling the ground truth video. Due to minimize the errors and noise effect on the reconstructed sequence, we introduced an Gaussian filter with size3×3.All of the reconstructed video should be compared with the ground truth video sequence. Taking 10 frames as experimental cases, we achieved the PSNR, sharpness, and entropy results as shown in Figure6. The maximum iteration of POMS-PSO is set to 15. The size of whole populations is set to 10 particles.Criteria comparison on VSR performance for 10 frames video.
(a)
PSNR
(b)
Sharpness
(c)
EntropyFrom PSNR lines in Figure6, we can see that the result using POMS-PSO algorithm achieved PSNR value approximating about 47 dB and bilinear interpolation only got 18 dB. That is to say, the proposed method has a much better objective image quality than that of bilinear interpolation method. The higher sharpness value means the more details the frame shows. We can see that the sharpness values of POMS-PSO algorithm is much higher than those of bilinear interpolation method by 1.5 in average and even better than those of ground truth by 0.5 in average. The entropy of proposed method is lower than the bilinear interpolation, but it is very close to the ground truth. The close entropy means the reconstructed sequence has similar pixel value distribution to ground truth, while the higher entropy got by interpolation could be seen as artificial effect.The results of the reconstruction using bilinear interpolation and POMS-PSO are given as Figure7. The images given are the 9th and the 10th frames of the corresponding reconstructed video. The images in Figure 7 show that video reconstruction using POMS-PSO could obtain high-quality results with a small number of iterations. The video is continuous smoothly.Results of the 9th (left) and 10th (right) frames of the video.
(a)
The input
(b)
VSR results of bilinear interpolation
(c)
VSR results of POMS-PSO
(d)
Ground truthFor more clear details, the helicopter and the railings from the 9th and 10th frames are zoomed and shown in Figure8. Obviously, the helicopter is of great interest in the whole image, and the railings are the most difficult to recognize fine details part. For the bilinear interpolation, the helicopter and the railings are blurred seriously. Judging from Figure 8, we can draw a conclusion that the objective vision quality of POMS-PSO result is much better than that of bilinear interpolation, especially in the regions of interest. Considering other parts of the images, say the sky and the trees, VSR by POMS-PSO can also achieve better performance.Zoomed parts of the 9th and 10th frames.
(a)
Zoomed helicopter of the 9th frame, bilinear (left), proposed method (middle), and ground truth (right)
(b)
Zoomed railings of the 9th frame, bilinear (left), proposed method (middle), and ground truth (right)
(c)
Zoomed helicopter of the 10th frame, bilinear (left), proposed method (middle), and ground truth (right)
(d)
Zoomed railings of the 10th frame, bilinear (left), proposed method (middle), and ground truth (right)


## 4.1. Reconstruction Quality Evaluation Criteria

### 4.1.1. Similarity-Based Evaluation Criteria

The PSNR is an engineering term for the ratio between the maximum possible power of a signal and the power of corrupting noise that affects the fidelity of its representation. The PSNR is most commonly used as a measure of quality of reconstruction in image compression. It is defined via the mean squared error (MSE) as follows:(15)MSE=1mn∑i=0m-1‍∑j=0n-1‍∥I(i,j)-K(i,j)∥2.I and K are two m×n monochrome images where one of the images is considered a noisy approximation of the other.The PSNR is defined as follows:(16)PSNR=10log⁡(2552MSE)(dB).


### 4.1.2. Detail-Based Evaluation Criteria

Sharpness is a photograph term that reflects image’s detail as an integral part of its appeal. Sharpness is defined by the boundaries between zones of different tones or colors. It is illustrated by the bar pattern of increasing spatial frequency.As shown in Figure5, the top portion is sharp and its boundaries are crisp steps, not gradual. The bottom portion illustrates how the degraded pattern is blurred.Figure 5
Bar pattern: original (top); with degradation (bottom).The sharpness is most commonly defined as the following equation for onem×n monochrome image I:
(17)sharpness=22(m-1)(n-1)1111111111×∑i=1m-1‍∑j=1n-1‍(Ii,j-Ii+1,j)2+(Ii,j-Ii,j+1)2.
The sharpness lines indicate the average gradient values of all frames in the video, which express whether the variation of scene is large or not in adjacent pixels.


### 4.1.3. Information-Based Evaluation Criteria

Image entropy reflects the average amount of information for a given image and the entropy of grayscale image indicates the amount of information contained in aggregation feature of gray distribution for the given image and it can be computed as follows:(18)Entropy=-∑i=0M-1‍Pilog⁡2(Pi),
where M is the number of gray levels and Pi is the probability associated with gray level i.


## 4.1.1. Similarity-Based Evaluation Criteria

The PSNR is an engineering term for the ratio between the maximum possible power of a signal and the power of corrupting noise that affects the fidelity of its representation. The PSNR is most commonly used as a measure of quality of reconstruction in image compression. It is defined via the mean squared error (MSE) as follows:(15)MSE=1mn∑i=0m-1‍∑j=0n-1‍∥I(i,j)-K(i,j)∥2.I and K are two m×n monochrome images where one of the images is considered a noisy approximation of the other.The PSNR is defined as follows:(16)PSNR=10log⁡(2552MSE)(dB).


## 4.1.2. Detail-Based Evaluation Criteria

Sharpness is a photograph term that reflects image’s detail as an integral part of its appeal. Sharpness is defined by the boundaries between zones of different tones or colors. It is illustrated by the bar pattern of increasing spatial frequency.As shown in Figure5, the top portion is sharp and its boundaries are crisp steps, not gradual. The bottom portion illustrates how the degraded pattern is blurred.Figure 5
Bar pattern: original (top); with degradation (bottom).The sharpness is most commonly defined as the following equation for onem×n monochrome image I:
(17)sharpness=22(m-1)(n-1)1111111111×∑i=1m-1‍∑j=1n-1‍(Ii,j-Ii+1,j)2+(Ii,j-Ii,j+1)2.
The sharpness lines indicate the average gradient values of all frames in the video, which express whether the variation of scene is large or not in adjacent pixels.


## 4.1.3. Information-Based Evaluation Criteria

Image entropy reflects the average amount of information for a given image and the entropy of grayscale image indicates the amount of information contained in aggregation feature of gray distribution for the given image and it can be computed as follows:(18)Entropy=-∑i=0M-1‍Pilog⁡2(Pi),
where M is the number of gray levels and Pi is the probability associated with gray level i.


## 4.2. Reconstruction Simulation and Experimental Results

We assessed the performance of the video superresolution process using POMS-PSO algorithm on several low-resolution video sequences and compared the result to that of the bilinear interpolation method. The resolution of the reconstructed video is twice as large as incoming video, in both height and width. To compare the objective criteria, the input video used to perform superresolution procedure was obtained by 2-downsampling the ground truth video. Due to minimize the errors and noise effect on the reconstructed sequence, we introduced an Gaussian filter with size3×3.All of the reconstructed video should be compared with the ground truth video sequence. Taking 10 frames as experimental cases, we achieved the PSNR, sharpness, and entropy results as shown in Figure6. The maximum iteration of POMS-PSO is set to 15. The size of whole populations is set to 10 particles.Criteria comparison on VSR performance for 10 frames video.
(a)
PSNR
(b)
Sharpness
(c)
EntropyFrom PSNR lines in Figure6, we can see that the result using POMS-PSO algorithm achieved PSNR value approximating about 47 dB and bilinear interpolation only got 18 dB. That is to say, the proposed method has a much better objective image quality than that of bilinear interpolation method. The higher sharpness value means the more details the frame shows. We can see that the sharpness values of POMS-PSO algorithm is much higher than those of bilinear interpolation method by 1.5 in average and even better than those of ground truth by 0.5 in average. The entropy of proposed method is lower than the bilinear interpolation, but it is very close to the ground truth. The close entropy means the reconstructed sequence has similar pixel value distribution to ground truth, while the higher entropy got by interpolation could be seen as artificial effect.The results of the reconstruction using bilinear interpolation and POMS-PSO are given as Figure7. The images given are the 9th and the 10th frames of the corresponding reconstructed video. The images in Figure 7 show that video reconstruction using POMS-PSO could obtain high-quality results with a small number of iterations. The video is continuous smoothly.Results of the 9th (left) and 10th (right) frames of the video.
(a)
The input
(b)
VSR results of bilinear interpolation
(c)
VSR results of POMS-PSO
(d)
Ground truthFor more clear details, the helicopter and the railings from the 9th and 10th frames are zoomed and shown in Figure8. Obviously, the helicopter is of great interest in the whole image, and the railings are the most difficult to recognize fine details part. For the bilinear interpolation, the helicopter and the railings are blurred seriously. Judging from Figure 8, we can draw a conclusion that the objective vision quality of POMS-PSO result is much better than that of bilinear interpolation, especially in the regions of interest. Considering other parts of the images, say the sky and the trees, VSR by POMS-PSO can also achieve better performance.Zoomed parts of the 9th and 10th frames.
(a)
Zoomed helicopter of the 9th frame, bilinear (left), proposed method (middle), and ground truth (right)
(b)
Zoomed railings of the 9th frame, bilinear (left), proposed method (middle), and ground truth (right)
(c)
Zoomed helicopter of the 10th frame, bilinear (left), proposed method (middle), and ground truth (right)
(d)
Zoomed railings of the 10th frame, bilinear (left), proposed method (middle), and ground truth (right)


## 5. Conclusions and Future Works

In this paper, we proposed POMS-PSO which could help to find out better swarm configuration for given problem. The POMS-PSO employed two layer loops and three subswarms. Performance comparison on four standard benchmark functions showed POMS-PSO could achieve higher accuracy in unimodal and multimodal functions.The model of imaging degeneration is very important during reconstruction. To some extent, proper model determination is the basis of finding the solution. Based on a imaging degeneration model and POMS-PSO, we proposed a novel VSR method in view of optimization computation. In view of computation intelligence, the video sequence with higher resolution is thought as the optimal resolution for swarm optimization. Experimental results showed that the proposed novel method could obtain high objective and subjective quality results.In the future, more efforts should be focused on the perfect degeneration model building and the fast implementation of multiple swarms based adaptive optimization algorithms.



---

*Source: 373425-2014-08-28.xml*