Вход: растр, два спектра, две области интереса.
Выход: растр.
Растр преобразовывается таким образом, чтобы спектры указанных областей интереса совпали с указанными спектрами.
В растре должна быть спектральная калибровка. Число каналов и значения длин волн в растре и в спектрах должны совпадать (рекомендуется использовать интерполяцию спектров, если длины волн в спектрах не соответствуют в точности длинам волн в растре). Единицы измерения двух спектров также должны быть одинаковыми.
Процедура вычисляет выходной растр путём применения к каждому каналу этого растра линейного преобразования:
new_value(z) = a(z) * old_value(z) + b(z)
где z - номер канала, old_value и new_value - старое и новое значения некоторого пикселя в этом канале, a(z) и b(z) - коэффициенты. Пара коэффициентов a и b вычисляется для каждого канала с помощью эталонных спектров и спектров областей интереса.
Коэффициенты вычисляются таким образом, чтобы линейное преобразование с этими коэффициентами переводило спектр первой области интереса в первый эталонный спектр и спектр второй области интереса во второй эталонный спектр:
etalon_spectrum_value_1(z) = a(z) * selection_spectrum_value_1(z) + b(z) etalon_spectrum_value_2(z) = a(z) * selection_spectrum_value_2(z) + b(z)
где etalon_spectrum_value_1(z), etalon_spectrum_value_2(z), selection_spectrum_value_1(z), selection_spectrum_value_2(z) - известные величины, а именно, значения эталонных спектров и спектров областей интереса в канале z.
Если значения спектров областей интереса совпадают в каком-то канале, то для этого канала будут выставлены нулевые коэффициенты и будет выдано предупреждение об этом.