“…Therefore, the finite difference method [ 19 , 20 ]] is used as the discrete differential operator to solve Eq (10) , and the discrete expression of image enhancement is obtained as shown in Eq (11) : Where, ∇ N ( I t ) = I x,y−1 − I x,y , ∇ S ( I t ) = I x,y+1 − I x,y , ∇ E ( I t ) = I x−1,y+1 − I x,y , ∇ W ( I t ) = I x+1,y − I x,y is four divergence formulas, whose values are the partial derivatives of the current pixel in the directions of southeast, northwest, and northwest, respectively. cN x,y = exp(−‖∇ N ( I t )‖ 2 / k 2 ), cS x,y = exp(−‖∇ S ( I t )‖ 2 / k 2 ), cE x,y = exp(−‖∇ E ( I t )‖ 2 / k 2 ), cW x,y = exp(−‖∇ W ( I t )‖ 2 / k 2 ) is the thermal conductivity in the four directions of southeast, northwest, I is the current image.…”