Эталонная коррекция по двум спектрам

Вход: растр, два спектра, две области интереса.

Выход: растр.

Растр преобразовывается таким образом, чтобы спектры указанных областей интереса совпали с указанными спектрами.

Параметры

В растре должна быть спектральная калибровка. Число каналов и значения длин волн в растре и в спектрах должны совпадать (рекомендуется использовать интерполяцию спектров, если длины волн в спектрах не соответствуют в точности длинам волн в растре). Единицы измерения двух спектров также должны быть одинаковыми.

Работа процедуры

Процедура вычисляет выходной растр путём применения к каждому каналу этого растра линейного преобразования:

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.

Если значения спектров областей интереса совпадают в каком-то канале, то для этого канала будут выставлены нулевые коэффициенты и будет выдано предупреждение об этом.

Оглавление

Предыдущий раздел

Умножение каналов растра на коэффициенты

Следующий раздел

Анализ данных