This method is an analytic scheme of fractional steps for convection operator and for diffusion operator respectively.

Experiments also show that selecting suitable symmetrization normalization techniques while forming the diffusion operator is important to hyperspectral imagery representation.

Based on the skew direction difference operator, the implicit skew difference schemes were constructed, and the boundary condition was unified, to solve the diffusion equation.