CC BY
29
Обчислювальні методи в радіоелектроніці
УДК 621.372.061
РЕСТАВРАЦІЯ ЗОБРАЖЕНЬ МЕТОДОМ УМОВНОЇ ДЕКОНВОЛЮЦІЇ В БАЗИСІ ПЕРЕТВОРЕННЯ RTF
Шебеста В. , д.т.н., проф.; Іванюк Н. О. , ассистент; Солохненко К. С.2, магістрант; Рибін О. І.2, д.т.н., проф.
1 Технічний університет Брно, м. Брно, Чехія,
2
Національний технічний університет України
«Київський політехнічний інститут» , м. Київ, Україна
IMAGE RESTORATION BY CONVENTIONAL DECONVOLUTION METHOD IN BASIS OF RTF TRANSFORMATION
1 2 Shebesta V. Doctor of Science (Technics), Professor; Ivaniuk N. , assistant;
Solokhnenko K.2, undergraduate student, Rybin A.2, Doctor of Science (Technics),
Professor
1Brno University of Technology, Brno, Czech Republic
2
National Technical University of Ukraine “Kyiv Polytechnic Institute ”, Kyiv, Ukraine
Вступ
Задача реставрації зображень є однією з важливих задач сучасної техніки та медицини. При розв’язанні задачі реставрації на базі обраної моделі спотворення (деградації) образу слід використовувати апріорну інформацію про причину такого спотворення. В більшості випадків використовують лінійну модель деградації, що відповідає сучасному стану розвитку систем відображення інформації [1]
ГО ГО
L (х, у)= J J g(x, y, x, y') ■ f (x', y' )dx'dy '+V(x, y), (1)
—ГО—ГО
і і
де g(x, y, x, y ) — результуюча імпульсна характеристика каналу прийому-передачі (в якому і виникають спотворення образу); f (x' , y' ) — вихідний сигнал-зображення; V(x, y) — адитивний шум; L(x, y) — отримуваний спотворений образ; x, y — відповідні просторові координати.
У випадку однорідного поля образу вираз (1) перетворюється в інтеграл Дюамеля
ГО ГО
L(x,y)= J J g(x — x,y —y) ■ f(x',y')dx'dy'+V(x,y), (2)
—ГО—ГО
Задачею реставрації є пошук оцінки f (x, y) вихідного сигналу, такої, що вона за деяким наперед заданим об’єктивним критерієм «мало» відрізняється від сигналу f (x, y) . У випадку реставрації за методом умовної де-конволюції таким критерієм (при виконанні умови «гладкості» та не пере-
Вісник Національного технічного університету України «КПІ» 25 Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
вищення енергії шуму в оцінці f (x, у) над енергією шуму в деградованому образі L (х, у)) є норма відстані між оцінкою f (х, у) та вихідним еталонним образом f (х, у).
Пошук оцінки f (х, у) за аналітичним виразом (2) при відомих енергії шуму в образі L (х, у) та імпульсній характеристиці деградації (часто параметри цієї імпульсної характеристики треба шукати «всліпу» на основі апріорної інформації про можливі джерела спотворень [1]) є складною задачею. Тому звичайно переходять від аналітичного виразу (2) до його запису в дискретному вигляді
L = G • f + V (3)
де f, L — вектори-стовпці відповідно шуканого та спотвореного образів
9 —
розміру N х 1; V — вектор-стовпець (невідомої) реалізації шуму в образі
розміру N х1; G — матричний дискретний оператор згортки [1] порядку
N 2.
Умова «гладкості (імпульсна характеристика фільтра) має вигляд
'0 1 0]
c = 1 -4 1 (4)
0 1 0
Розв’язання задачі реставрації методом умовної деконволюції в натуральних координатах при дискретних наближеннях (3) моделі (2) має вигляд
f
=T = 1 =T =
G • G +1 • c • c Я
-1
G • L,
(5)
де T — знак транспонування; Я — коефіцієнт варіації Лагранжа.
Великі порядки матриць в (5) значно обмежують використяння цього виразу при реставрації зображень. Так, наприклад, при розмірі зображення
N = 1024 порядок матриці
=Т = 1 =т =
G • G + - • c • c Я
дорівнює N2 = 1048576, а кі-
лькість операцій множення/ділення при її обертанні методом гаусса дорв-
(2 \ 3 і о
N ) , тобто 1,1529 -10 . Тому подальшим шляхом реалізації реставрації є перехід до спектральних перетворень, так, щоб радикально скоротити кількість операцій при оберненні матриць [2—9]. Серед таких перетворень найбільшого поширення набуло дискретне перетворення Фур’є [2, 3]. При цьому вираз (5) перетворюється до вигляду
-1-1
Ю
G2о +\ • c2о
Я
Go • L,
о
(6)
26
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
де Gіа та Є2а — діагональні матриці спектрів потужності матриці деградації образу та умови «гладкості» відповідно порядку N , в діагоналі яких розташовані послідовно квадрати модулів спектрів рядків двовимірного
перетворення Фур’є імпульсної характеристики деградації та «гладкості»;
—*
Ga — діагональна матриця спряженого двовимірного перетворення Фур’є імпульсної характеристики деградації порядку N2; * — знак комплексно-
го спряження. Матриця
G2a +
1
I
' с2а
є діагональною. Тому її обернення
не становить складності, і вираз (6) для кожної окремої частоти перетворюється до вигляду
F (m, n)
G 2а(m, n)
Ga(m, n) ( G 2 a( m, n )
2 1 +— I
= -----Г ' La(m n)•
C 2 a( m, П) )
2
1
(7)
де m, n — номери частот двовимірного перетворення Фур’є.
Таким чином, для реалізації методу умовної деконволюції необхідно обчислити енергію шуму в деградованому образі L(m, n) (вважаючи шум ергодичним) та використати умову «гладкості», фільтруючи образ в процесі реставрації фільтром з імпульсною характеристикою (4), спектр якої c ( m, n ).
Реставрація в базисі перетворення RTF
Фільтрацію спотвореного сигналу (7) за допомогою фільтра
= 2
Kr (m, n)
G 2a(m, n)
Gа(m, n) ( G 2 а (m, n) +1' C2 a( m, n) )
I
1
(7а)
можна представити у матричному (для простоти у одновимірному випадку) у вигляді
fa Ka " F N " sx
або для оригіналу оцінки реставрованого образу
fx = FN ' K a ' FN ' ,
(8)
де Ka — діагональна матриця коефіцієнтів передач фільтра (на різних частотах) порядку N; Fn , Fn — нормовані (діленням на VN) матричні дискретні оператори прямого та зворотного перетворень Фур’є; fx, sx — відповідно вектори-стовпці реставрованого та деградованого сигналів розміру N х 1. Як відомо, частоти для дискретного перетворення Фур’є це числа, що утворюють дискретний ряд з кроком 2л / N. В роботах [10,11] по-
Вісник Національного технічного університету України «КПІ» 27 Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
казано, що при переході від різницевих рівнянь до рівнянь в частотній області такий перехід до перетворення RTF. При цьому вираз (8) перетворюється до вигляду
_ fx = Fn ■ Кp- Fn ■ Jx, (9)
де К p — діагональна матриця коефіцієнтів передачі фільтра, в якій на ві-
дміну від К ю замість частот cok = 2л- k / N аргументами є власні частоти Рк. Ці власні частоти надаються наступним чином.
Перші два власних значення дорівнюють 0; 2 відповідно ( Р0 = 0; Р1 = 2 ). Далі Рз =j; Р4 = -у.
Усі наступні власні значення (як і чотири вище згаданих) розташовані на комплексній площині на одиночному крузі (рис 1) з центром в точці (1,0).
Так, наступні чотири власні значення знаходять діленням чотирьох дуг круга на рівні півдуги (зірочки на рис 1). Для блоку восьмого порядку вісім дуг ділять на рівні частини (кружечки на рис 1) і т.д. [11 ].
Рис. 1 Комплексна площина на Для впорядкування отриманих власних
одиничному крузі значень їх нумерацію слід проводити, по-
чинаючи від точки (0,0) за ходом годинникової стрілки. При форматі N перетворення RTF кількість власних значень, отриманих за указаним на рис.1 правилом, становить N. Причому при подальшому русі за годинниковою стрілкою (після повного оберту) номер к + n х N = к, де k,n - цілі числа ; k<N.
Значення k-го відліку (k-го власного значення) на комплексній площині можна обчислити за формулою
л л
рк = 2sin(- -
КЛ-Лк) л -j-k
k)e 2 N = 2j cos(— k)e N
(10a)
(106)
2 N
або для неперервного аргументу к
Р = 2 cos фе-*ф
л л
причому рє[-;--].
Якщо нумерацію власних значень Р вести проти часової стрілки, починаючи з точки (2,0) на рис. 1, то формула (10а) матиме вигляд
Р = \ + ejc°k (Юв)
з еквідистантним кроком cok
—k; k
N
0,1,..,N -1.
28
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
Перехід від спектрів перетворення Фур’є до спектрів перетворення RTF Більш детально описано в [12—14].
Таким чином, для організації реставрації в базисі перетворення RTF достатньо Знайти відповідні спектри RTF для Kr (m, n) в (7а) і за виразом
(7) (використовуючи для спотвореного образу спектр Фур’є L(m,n), обчислити за (9) оригінал одного з наближень (при даному коефіцієнті варіації А) та обчислити по отриманому образу енергію шуму.
Результати обчислень
Алгоритм реставрації в базисі перетворення RTF апробовано на прикладах реставрації різних геометричних образів, обраних виходячи з міркувань простоти візуальної перевірки результатів.
Так, для образу рис.2а проводилась його деградація за допомогою згортки з неточкової імпульсною характеристикою (рис. 2б та рис. 2г). На рис. 2в, рис. 2д та рис. 2е показані результати реставрації.
Рис.2а Тестовий образ
Рис. 2б Образ спотворений
неточковою імпульсною Рис. 2в Реставрований об-
характеристикою та адити- раз
вним шумом с = 1
Рис. 2г Образ спотворений
неточковою імпульсною Рис. 2д Реставрований об- Рис. 2е Реставрований об-
характеристикою та адити- раз при с = 0,9 раз при с = 0,9
вним шумом с = 0,9
Аналогічно, на рис. 3 наведемо еталон (рис. 3а), один зі спотворених образів (рис. 3б) та результати реставрації з наведеними для кожного тако-
Вісник Національного технічного університету України «КПІ» 29 Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
го результату діючими значеннями шуму в спотвореному образі (рис. 3в -Зг).
Рис.За Тестовий образ
Рис. 3б Образ спотворений неточковою імпульсною характеристикою та адитивним шумом <Т = 0,04
Рис. 3в Реставрований образ при <Т = 0,04
Рис. 3г Реставрований об- Рис. 3д Реставрований образ при <Г = 0,9 раз при <Г = 0,5
Висновки
Реставрація в базисі перетворення RTF з використанням методу умовної деконволюції дає задовільні результати, що ілюструють наведені приклади.
Певна зміна динамічного діапазону реставрованих образів може бути скомпенсована при подальшій обробці зміною шкали рівнів сірого.
Отримані позитивні результати свідчать про перспективність методу, що є підставою для подальшого його дослідження.
Перелік посилань
1. Jan Jin Qslicova filtrace, analyza a restaurace signlu / Jan Jin. — VUT v BRNE. — 1997,— 438 p.
2. Рибін О .І. Реставрація образів методом умовної деконволюції в області просторових частот / О. І. Рибін, В. Ю. Корольов // Вісник Технічного університету Поділля. — 2000. — №1. — C. 145—147.
3. Рыбин А. И. Реставрация образов в частотной области методом взвешенной фильтрации / А. И. Рыбин, В. Ю. Королев // Известия вузов. Радиоэлектроника. — 2001. — Т. 44, № 4. — С. 51—56.
30
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
4. Рибін О. І. Аналіз лінійних систем в області трансформант перетворення Уол-ша-Адамара / О. І. Рибін, А. П. Ткачук // Вісник НТУУ «КПІ». Серія Радіотехніка. Ра-діоапаратобудування. — 2006. — № 33. — С. 14—23.
5. Рыбин А. И. Анализ линейных цепей в базисе преобразований Уолша / А.И. Рыбин // Известия вузов. Радиоэлектроника — 2004. — Т. 47, №5. — С. 36—41.
6. Рибін О. І. Реставрація образів за методом умовної деконволюції в області трансформант перетворення Адамара / О. І. Рибін, Н. О. Іванюк // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2011. — № 46. — С. 51—58.
7. Рибін О. І. Алгоритм побудови матриці деградації образу з використанням дискретного перетворення Адамара / О. І. Рибін, Н. О. Іванюк // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2012. — № 49. — С. 18—32.
8. Рибін О. І. Умовна деконволюція в області трансформант Фур’є. Побудова матриці деградації образу. / О. І. Рибін, Н. О. Іванюк // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2011. — № 47. — С. 30—41.
9. Рибін О. І. Умовна деконволюція в області трансформант Фур’є. Побудова оберненої матриці деградації образу / О. І. Рибін, Н. О. Іванюк // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2012. — № 50. — С. 21—29.
10. Рибін О. І. Аналіз лінійних ситем в області трансформант кратного перетворення EIWAL / О. І. Рибін, А. П. Ткачук // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2006. — № 33. — С. 31—38.
11. Рибін О. І. Аналіз лінійних систем з використанням кратних перетворень /
О. І. Рибін, Ю. Х. Ніжебецька, І. О. Рибіна // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2010. — № 40. — С. 5—11.
12. Рибін О. І. Зв’язок між перетвореннями RTF та Фур’є / О. І. Рибін, Н. О. Іванюк // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2012. — № 48. — С. 23—37.
13. Іванюк Н. О. Зв’язок між РТФ та Фур’є перетвореннями / Н. О. Іванюк,
0. І. Рибін // Радіотехнічні поля, сигнали, апарати та системи : міжн. конф., 22 —29 лютого 2012 р. : тези доп. — м. Київ, НТУУ «КПІ» — 2012. — С. 25—26. — Режим доступу: http://conf.rtf.kpi.ua/ua/2012/2012sec1/42-ivaniuk-rybin
14. Наталенко С. С. Умовна деконволюція образу в області просторових власних частот / С. С. Наталенко, О. І. Рибін // Вісник НТУУ «КПІ». Серія Радіотехніка. Радіоапаратобудування. — 2009. — № 39. — С. 8—12.
References
1. Jan Jin (1997) Cislicova filtrace, analyza a restaurace signlu. VUT v BRNE, 438 p.
2. Rybin O.I. and Korolev V.Yu. (2000) Restavratsiia obraziv za metodom umovnoi dekonvoliutsii v oblasti prostorovykh chastot. Visnik Tekhnichnogo universytetu Podillia. No
1, pp. 145-147.
3. Rybin A.I. and Korolev V.Yu. (2001) Restavratsiia obrazov v chastotnoi oblasti metodom vzveshenoi filtratsii Uolsha [Analysis of linear circuits in the basis of Walsh transformation]. Radioelectronics and Communications Systems. Vol. 44, No 4. pp. 51-56.
4. Rybin, A. I., Tkachyk, A. P. (2006) The analisis of linear systems in basis of Walsh-Hadamard transformation. Visn. NTUU KPI, Ser. Radioteh. radioaparatobuduv. no. 33, pp. 14-23. (in Ukrainian)
5. Rybin A.I. (2004) Analiz lineinykh tsepei v bazise preobrazovania Uolsha [Analysis of linear circuits in the basis of Walsh transformation]. Radioelectronics and Communications Systems. Vol. 47, No 5, pp. 36-41.
Вісник Національного технічного університету України «КПІ» 31 Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
6. Rybin, A. I., Ivanyuk, N. A. (2011) Images restoration by theconditional deconvolution methodin basis of Hadamard transformation. Visn. NTUUKPI, Ser. Radioteh. radioapa-ratobuduv. no. 46, pp. 51-58. (in Ukrainian)
7. Rybin, A. I., Ivanyuk, N. A. (2012) Image degradation matrix formation algorithm with usage of Adamar discrete transformation. Visn. NTUU KPI, Ser. Radioteh. radioaparat-obuduv. no. 49, pp. 18-32. (in Ukrainian)
8. Rybin, A. I., Ivanyuk, N. A. (2011) Fourier transformant area deconvolution condition. Construction of the image degradation matrix. Visn. NTUU KPI, Ser. Radioteh. radi-oaparatobuduv. no. 47, pp. 30-41. (in Ukrainian)
9. Rybin, A. I., Ivanyuk, N. A. (2012) Conventional deconvolution in the Fourier transforms. The image degradation inverse matrix construction. Visn. NTUU KPI, Ser. Radioteh. radioaparatobuduv. no. 50, pp. 21-29. (in Ukrainian)
10. Rybin, A. I., Tkachyk, A. P. (2006) The analisis of linear systems in basis of EIWAL transformation. Visn. NTUU KPI, Ser. Radioteh. radioaparatobuduv. no. 33, pp. 31-38. (in Ukrainian)
11. Rybin, A. I., Nizhebetska, Y. Kh., Rybina, I. O. (2010) Analysis of the linear systems with the use of multiple transform. Visn. NTUU KPI, Ser. Radioteh. radioaparatobuduv. no. 40, pp. 5-11. (in Ukrainian)
12. Rybin, A. I., Ivanyuk, N. A. (2012) The relationship between RTF and Fourier transform. Visn. NTUU KPI, Ser. Radioteh. radioaparatobuduv. no. 48, pp. 23-37. (in Ukrainian)
13. Ivaniuk N.O. and Rybin A.I. (2012) Zviazok mizh peretvorenniamy RTF ta Furie
[Relationship between RTF and the Fourier transformation]. Radioengineering fields, signals, devices and systems : Int. Sci. and Tech. Conf. pp. 25-26. Available at:
http://conf.rtf.kpi.ua/ua/2012/2012sec1/42-ivaniuk-rybin
14. Natalenko, S. S., Rybin, A. I. (2009) Conditional deconvolution ща image in the domain of spatial own frequencies. Visn. NTUU KPI, Ser. Radioteh. radioaparatobuduv. no. 39, pp. 8-12. (in Ukrainian)
Шебеста В., Іванюк Н. О., Солохненко К. С. , Рибін О. І. Реставрація зображень методом умовної деконволюції в базисі перетворення RTF. Розглянуто метод умовної деконволюції реставрації образів, для реалізації якого потрібно по спотвореному образу обчислити енергію шуму та на базі апріорної інформації про джерела спотворення (деградації) оцінити характер та параметри спотворюючої імпульсної характеристики системи відображення. В роботі запропоновано використовувати для реставрації базис перетворення RTF, який відрізняються від перетворення Фур’є і є «калькою» різницевих методів. Приведено приклади реставрації еталонних образів з різним рівнем адитивного шуму та різними імпульсними характеристиками деградації (при цьому імпульсні характеристики не підбиралися «всліпу», а вважалися відомими).
Ключові слова: реставрація образу, умовна деконволюція, перетворення RTF.
Шебеста В., Иванюк Н. О., Солохненко К. С. , Рыбин А. И. Реставрация изображений методом условной деконволюции в базисе преобразования RTF. Рассмотрен метод условной деконволюции реставрации образов, для реализации которого необходимо по искаженному образу вычислить энергию шума и на базе априорной информации про источники искажений (деградации) оценить характер и параметры искажающей импульсной характеристки системы отображения. В работе предложено использовать для реставрации базис преобразования RTF, отличный от преобразования Фурье и являющийся «калькой» разностных методов. Приведены примеры реставрации эталонных образов с разным уровнем шума и разными характеристиками дегра-
32
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
Обчислювальні методи в радіоелектроніці
дации (при этом импульсные характеристики не подбирались «вслепую», а считались известными).
Ключевые слова: реставрация образа, условная деконволюция, преобразование RTF.
Shebesta V., Ivaniuk N., Solohnenko K., Rybin A. Image restoration by conventional deconvolution method in basis of RTF transformation.
Introduction. Conventional deconvolution method for image restoration is considered. It is need to calculate the noise energy on the base of priori information on the sources of distortion (degradation), to estimate the character and parameters of the distorting impulse response of reflection system.
The problem statement. It should be moved from natural coordinates to the discrete orthogonal transformation region using this method to avoid overly cumbersome calculations. The RTF transformation method for restoration is proposed in this paper. It differs from the Fourier transform and it is the "tracing" of difference methods.
Conclusions. Examples of reference images restoration with different levels of additive noise and different degradation impulse responses (caused not point character of the reflection system impulse response) are shown. Impulse responses are not selected "blindly" and assumed to be known.
Keywords: image restoration, conventional deconvolution, RTF transformation.
Вісник Національного технічного університету України «КПІ» 33 Серія — Радіотехніка. Радіоапаратобудування. — 2014. — №56
