输入：

$\mathit{Y}\in {\mathit{R}}^{L\times N},{\mathit{Y}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}}\in {\mathit{R}}^{L\times q}bySLIC,{\mathit{A}}^{\left(0\right)}\in {\mathit{R}}^{K\times N}byFCLS,$

${\mathit{A}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}}\in {\mathit{R}}^{K\times q}bySLIC,{\mathit{M}}^{\left(0\right)}\in {\mathit{R}}^{L\times K}byVCA,\alpha ,\beta$$\gamma ,\eta ,q,\mathit{d}{\mathit{M}}^{\left(0\right)}=[\mathit{d}{\mathit{M}}_1|\cdots |\mathit{d}{\mathit{M}}_q],\forall \mathit{d}{\mathit{M}}_i\in {\mathit{R}}^{L\times K}=0_{L\times K},$

${\mathit{S}}^{\left(0\right)}=[{\mathit{S}}_1|\cdots |{\mathit{S}}_q],\forall S_i\in {\mathit{R}}^{K\times K}={\mathit{I}}_K$

输出：

$\mathit{A}\in {\mathit{R}}^{K\times N},\mathit{M}\in {\mathit{R}}^{L\times K},{\mathit{S}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}}=[{\mathit{S}}_1|\cdots |{\mathit{S}}_q],\forall S_i\in {\mathit{R}}^{K\times K},$

$\mathit{d}{\mathit{M}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}}=[\mathit{d}{\mathit{M}}_1|\cdots |\mathit{d}{\mathit{M}}_q],\forall \mathit{d}{\mathit{M}}_i\in {\mathit{R}}^{L\times K}$

1.  $i\leftarrow 1$，如果没达到停止条件执行（2-6）：

2.  从$n=\mathrm{1,2},\cdots ,q$使用公式（16）~公式（20）依次求解${\mathit{a}}_n$，${\mathit{W}}_n^{(A)}$，${\mathit{T}}_n^{(A)}$，${\mathit{V}}_n^{(A)}$和${\mathit{G}}_n^{(A)}$；

3.  从$l=\mathrm{1,2},\cdots ,L$使用公式（23）~公式（25）依次求解${{\displaystyle \tilde{\mathit{m}}}}_l$，${\mathit{W}}_l^{\left(M\right)}$和${\mathit{V}}_l^{(M)}$；

4.  从$n=\mathrm{1,2},\cdots ,q$使用公式（27）公式（29）依次求解$\mathit{d}{\mathit{M}}_n$，${\mathit{W}}_n^{(dM)}$和${\mathit{V}}_n^{(dM)}$；

5.  从$n=\mathrm{1,2},\cdots ,q$使用公式（31）求解${\mathit{S}}_n$，使用公式（32），公式（33）依次求解${\mathit{W}}_n^{(S)}$和${\mathit{V}}_n^{(S)}$；

6.  $i=i+1$。

7.  利用上述步骤求得的${\mathit{A}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}},\mathit{M},\mathit{d}{\mathit{M}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}},{\mathit{S}}_{\mathrm{s}\mathrm{l}\mathrm{i}\mathrm{c}}$，在每一个超像素中，将求得的变异端元作为端元使用FCLS逐像元求解丰度$\mathit{A}$。