Our new objective model is a function of variables (u, v, I1, I2). To simultaneously solve all the unknowns, an alternating minimization framework [28] can be used. Since the optical flow terms (i.e. the first and second terms of Eq. (14)) depend on the restored images, and the image restoration terms (i.e. the second and third terms of Eq. (14)) depend on the computed optical flow, we split the objective function into two coupled modules: EO(u,v,I)=λ∫ΩΨS(∥▿u∥2+∥▿v∥2)dxdy+∫ΩΨD((I2(x+u,y+v)−I1(x,y