Изображение Deblur с помощью метода Люси-Ричардсона
J = deconvlucy(I,psf)
J = deconvlucy(I,psf,iter)
J = deconvlucy(I,psf,iter,dampar)
J = deconvlucy(I,psf,iter,dampar,weight)
J = deconvlucy(I,psf,iter,dampar,weight,readout)
J = deconvlucy(I,psf,iter,dampar,weight,readout,subsample)
восстановления отображают J
= deconvlucy(I
,psf
)I
, который был ухудшен сверткой с функцией рассеяния точки (PSF), psf
, и возможно аддитивным шумом. Алгоритм основан на максимизации вероятности, что получившееся изображение J
является экземпляром оригинального изображения I
под статистикой Пуассона.
Чтобы улучшить восстановление, deconvlucy
поддерживает несколько дополнительных параметров, описанных ниже. Используйте []
в качестве заполнителя, если вы не задаете промежуточный параметр.
задает, какими пикселями во входном изображении J
= deconvlucy(I
,psf
,iter
,dampar
,weight
)I
рассматривается в восстановлении. Значение элемента в массиве weight
определяет, каким количеством рассматривается пиксель в соответствующем положении во входном изображении. Например, чтобы исключить пиксель из фактора, присвойте его значение 0
в массиве weight
. Можно настроить значение веса, присвоенное каждому пикселю согласно объему исправления плоского поля.
Можно использовать deconvlucy
, чтобы выполнить развертку, которая запускается где предыдущая остановленная развертка. Чтобы использовать эту функцию, передайте входное изображение I
и PSF, psf
, как массивы ячеек: {I}
и {psf}
. Когда вы делаете, функция deconvlucy
возвращает выходное изображение J
как массив ячеек, который может затем быть передан как входной массив в следующий вызов deconvblind
. Выходной массив ячеек J
содержит четыре элемента:
J{1}
содержит I
, оригинальное изображение.
J{2}
содержит результат последней итерации.
J{3}
содержит результат предпоследней итерации.
J{4}
является массивом, сгенерированным итеративным алгоритмом.
Выходное изображение J
могло показать вызов, введенный дискретным преобразованием Фурье, используемым в алгоритме. Чтобы уменьшать вызов, используйте I = edgetaper(I,psf)
прежде, чем вызвать deconvlucy
.
deconvlucy
преобразовывает PSF в double
без нормализации.
[1] Округ Южная Каролина. Четырехрядные ячмени и М. Эндрюс, Ускорение итеративных алгоритмов восстановления изображений, Прикладной Оптики, Издания 36, № 8, 1997.
[2] Р.Дж. Хэниш, Р.Л. Вайт, и Р.Л. Джиллилэнд, Развертки Изображений Космического телескопа Хабблa и Спектров, Развертка Изображений и Спектров, Эда. П.А. Янссон, 2-й редактор, Academic Press, CA, 1997.