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

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

Для итерации 1 укажите начальный набор конечных элементов как $$\text{{}{\text{e1}}_{\mathrm{}}^{(1)},{\text{e2}}_{\mathrm{}}^{(1)},\text{⋯},\text{ep}{\text{(}}_1^{\mathrm{})\}.}$$

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

где $${\text{E}}^{\mathrm{(}1}\text{)}\text{=}[\begin{array}{l}1\text{1}\text{⋯}\text{1e1}\text{(}\text{1}\text{)}\mathrm{}\text{e2}\text{(}\text{1}\text{)}\text{⋯}\text{ep}\text{(}\text{1}\mathrm{)].}\\ {}_{\mathrm{}}^{\mathrm{}}{}_{\mathrm{}}^{\mathrm{}}{}_{}^{\mathrm{}}\end{array}$$

Для итерации 2 выберите новый пиксельный спектр r, такой как r $$\notin \{e1{}_{\mathrm{(}}^{1\mathrm{)},}\text{e2}{}_{\mathrm{(}}^{1\mathrm{)},}\text{⋯}\text{,}ep\text{(}{\text{1}}_{}^{\mathrm{)\}.}}$$

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

Замените ^i-й конечный элемент в наборе на r, если вычисленный объем V (^{E (}2)) больше V ^(E (1)). Это приводит к обновлению набора конечных элементов. Например, если i = 2, новое множество конечных элементов, полученных в конце второй итерации, $${\text{равно}}_{\mathrm{}}^{\{\mathrm{}e1}\text{(}\text{2}{\text{)}}_{\mathrm{,}}^{\mathrm{}e2}\text{(}\text{2}\text{)}\text{=}r\text{,}\text{⋯}\text{,}{}_{}^{ep\mathrm{}(}2$$)}.

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