Вычислите полосы основных компонентов и уменьшите спектральную размерность входных данных с помощью MNF или PCA. Количество полос основного компонента, которые будут извлечены, устанавливается равным количеству членов конца, которые будут извлечены. Конечные элементы извлекаются из полос основных компонентов.

Случайным n выберите количество пиксельных спектров из уменьшенных данных в качестве начального набора конечных элементов.

Для итерации 1 обозначайте начальный набор конечных элементов как $$\left\{{\text{e}}_1^{(1)},{\text{e}}_2^{(1)},\cdots ,{\text{e}}_p^{(1)}\right\}$$.

Рассмотрим конечные элементы как вершины симплекса и вычислите объем при помощи

где $${\text{E}}^{(1)}=\left[\begin{array}{l}11\cdots 1\\ {\text{e}}_1^{(1)}{\text{e}}_2^{(1)}\cdots {\text{e}}_p^{(1)}\end{array}\right]$$.

Для итерации 2 Выберите новые пиксельные спектры r, такой что $$\text{r}\notin \left\{{\text{e}}_1^{(1)},{\text{e}}_2^{(1)},\cdots ,{\text{e}}_p^{(1)}\right\}$$.

Замените каждый конечный элемент в наборе на r и вычислите объем V симплекса (E⁽²⁾).

Замените i^th конечный элемент в наборе с r, если вычисленный объем V (E⁽²⁾) больше V (E⁽¹⁾). Это приводит к обновлению набора конечных элементов. Для примера, если i = 2, новый набор конечных элементов, выведенных в конце второй итерации, $$\left\{{\text{e}}_1^{(2)},{\text{e}}_2^{(2)}=\text{r},\cdots ,{\text{e}}_p^{(2)}\right\}$$.

Для каждой итерации выберите новые пиксельные спектры r и повторите шаги 5 и 6. Каждая итерация приводит к обновлению набора конечных элементов. Итерация заканчивается, когда общее количество итераций достигает заданного значения NumIterations.