CORRECTION OF SPATIAL LOSSES IN TELEVISION REGISTRATION SYSTEM WORKING IN CERTAIN CONDITIONS

Full Text

Введение. Различные задачи производственного и исследовательского характера на предприятиях, использующих радиоактивные материалы, решаются в радиационно-защитных камерах. Учитывая высокую ответственность, а зачастую и уникальность выполняемых работ, все технологические операции документируются с целью дальнейшего анализа. Для этой цели используются телевизионные системы регистрации (ТВ-системы), которые располагают перед смотровой системой в зоне, защищенной от радиации, как показано на рис. 1. Особые условия, которые обусловлены наличием смотровой системы толщиной несколько сотен миллиметров между объектом и камерным объективом, вносят специфику в формирование изображений такой ТВ-системой. Автором выполнен анализ ТВ-системы, при котором тракт формирования и передачи информации представлялся состоящим из ряда преобразовательных и усилительных звеньев, как показано на рис. 2. Рис. 1. Функциональная схема системы регистрации Рис. 2. Представление ТВ-системы для анализа Ограниченные пространственные разрешающие способности смотровой системы, камерного объектива, ПЗС камеры и воспроизводящего монитора не позволяют регистрировать мелкие детали изображения с достаточным числом градаций, а начиная с некоторых пространственных частот эти детали вообще перестанут регистрироваться системой. Эти ограничения описываются относительной двумерной частотноградационной характеристикой системы (ЧГХ) [1], которая определяется отношением числа градаций /иэ(П) , различимых на пространственной частоте Q, к числу градаций тоэ, различимых на пространственной частоте, близкой к нулю: 8ТВС(0) = тэ(0.)/тоэ. ЧГХ является убывающей функцией пространственной частоты (см. рис. 3): 5тс (Q) = 1, при Q = 0; 1 1 > 8ТВС (П) > 0, при П > 0.J Если в целях упрощения предположить, что число градаций различимых на экране монитора ЭВМ в крупных деталях тоэ и относительная ЧГХ 8^(0.) не изменяются по полю изображения и воспользоваться их усредненными значе «Инфокоммуникационные технологии» Том 11, № 1, 2013 Никитин К.А. 89 ниями тоэ и STBC(Q), то количество информации на выходе ТВ-системы, снимаемое с экрана монитора, будет определяться выражением [1] Np 1твс j*^§2[^03^rac(^) = о Np — 3Nр log2 tnоэ + 3 J log2 \_STBC (Г2) +1 і ] d£l9 где Np - число элементов разложения в растре; Q - двумерная пространственная частота. Второе слагаемое в этом выражении определяет пространственные потери и при тоэ > 1 является отрицательным. Рис. 3. Относительная ЧГХ системы регистрации Выполненный анализ ТВ-системы показал, что с достаточной для практики точностью зависимость <5^(0) можно заменить двумерной частотноконтрастной характеристикой T{Nx,Ny), которую можно построить инструментальными средствами. Несмотря на то что ЧКХ не учитывает зависимость числа воспроизводимых градаций от шумов, она позволяет проследить за качественными изменениями пространственных потерь. Возникающие в ТВ-системе пространственные потери можно разделить на две следующие группы. 1. Потери, вызванные хроматическими аберрациями в смотровой системе (хроматизм положения и хроматизм увеличения). Такие потери приводят к ухудшению различимости мелких деталей изображения и цветовым искажениям. Причины их возникновения и предлагаемый метод коррекции рассмотрены автором в [2]. 2. Потери, вызванные монохроматическими аберрациями (сферическая аберрация, кома, астигматизм), а также неидеальными частотными характеристиками преобразовательных и усилительных звеньев. Они приводят к снижению пространственной разрешающей способности системы и выражаются в ухудшении различимости мелких деталей изображения или их полной потере. Однако если мелкие детали потеряны не полностью, то их передачу можно улучшить, используя методы коррекции и улучшения изображений. Коррекция частотно-контрастной характеристики Целью коррекции ЧКХ является такая коррекция характеристик системы, которая позволяет получить на выходе изображение, соответствующее идеальной системе. Если характеристики объектива, ТВ-камеры и монитора выбраны высокоразрешающими, то в качестве «идеальной» можно рассматривать систему регистрации, в которой отсутствует звено радиационно-защитной смотровой системы, вносящее максимальный вклад в потери пространственной информации. Будем считать такую систему эталонной. Структура эталонной системы показана на рис. 4а). Структура реальной системы регистрации отличается от эталонной наличием между объектом и камерным объективом смотровой системы и показана на рис. 4б). Рис. 4. Структура ТВ-систем регистрации а) эталонной; б) реальной Несмотря на наличие аберраций комы, и астигматизма, которые делают систему регистрации неизопланарной в целом, поле изображения можно разбить на зоны, в пределах которых условия изопланарности будут соблюдаться с определенной точностью [3]. Поэтому в первом приближении такие системы регистрации можно считать линейными, инвариантными к сдвигу и рассматривать их как линейные двумерные фильтры пространственных частот. Процесс образования изображения в этом случае можно описать математической зависимостью между характеристиками объекта, системы регистрации и получаемого изображения. Для области пространственных частот можно записать [4-5] Ё{Их^у) = Й^х,Ну)-Ь{Ых,Иу), (1) «Инфокоммуникационные технологии» Том 11, № 1, 2013 90 Никитин К.А. где E(Nx,Ny) = J JE(x,y)e-i2^x+N^dx'dy' _ —00 —00 двумерный спектр распределения яркостей в плоскости изображений; L{NX, Ny ) = J J L(x, y)e~l2*iN*x+Nyy)dxdy - дву--00 -00 мерный спектр распределения яркостей в плоскости объекта; ЩК,МУ) - двумерная передаточная функция системы регистрации; Nx, Ny - пространственные частоты по направлениям х и у в плоскостях изображений и объекта. Известно, что передаточная функция Й{Их,Иу) является преобразованием Фурье реакции системы на входное воздействие Н(х\у') в виде двумерной S -функции. Следовательно, ее можно представить в виде H(Nx,Ny)= ] ]H(x,y)e-n^I+N^dx'dy'. равен нулю, а передаточная функция H(Nx,Ny) будет являться вещественной и характеризовать пространственный фильтр с нулевым сдвигом фазы. Модуль этой двумерной функции |Я(^,Лд| является двумерной частотно-контрастной характеристикой системы и позволяет установить соответствие между двумерными спектрами распределения яркостей в объекте и его изображении T(Nx,Ny)= H(Nx,N) . (3) Исходя из изложенного, для эталонной и реальной систем регистрации, показанных на рис. 4, запишем КЖ,КУ) = \H3m(Nx,Ny)\ ■ L(Nx,Ny) ; (4) E(Nx,Ny) = \&(Nx,Ny)\-L{Nx,Ny). (5) Выполнив преобразование этого выражения с помощью формулы Эйлера е~’ = cos сох -/sin ах, где (o = 2nN, представим передаточную функцию в виде вещественной и мнимой частей: H{Nx,Ny) = 00 оо 1 !н<*' ,y)cos2n{Nxx +Nyy)dxdy -оо —оо зо оо ,y)s\n2n(Nxx +Nyy)dxdy - -оо -оо = ReH(Nx,Ny) + \mH(Nx,Ny). (2) Используя (2), запишем передаточную функцию в виде модуля T(Nx,Ny) и аргумента <p(Nx,Ny), которые зависят от пространственных частот Nx,Ny и представляют собой двумерные амплитудно-частотную и фазо-частотную характеристики вида i<p(Nx+Ny) _ 'у/ Р J4>(NX+Ny) H{Nx,Ny) = \H{Nx,Ny)\e = T(Nx,Ny)eц При этом ПК ,Ny) = ^Re H{NX, Ny )]2 + \MH(NX, Ny )]2 1тЩМх,М) <p(Nx,N ) = arctg ReH(Nx,N ) В большинстве случаев импульсный отклик системы Н(х ,у ) симметричен по направлениям x и у, следовательно, является четной функцией. Поэтому второй интеграл в выражении (2) будет Если целью коррекции считать получение на выходе реальной системы изображения, совпадающего с изображением на выходе эталонной: E3m(Nx,Ny) = np(Nx,Ny), то ЧКХ реальной системы следует умножить на некоторую двумерную корректирующую функцию |ЯК(Л^)| таким образом, чтобы это равенство выполнялось: \Hp(Nx,Ny)\\AK(Nx,Ny)\L(Nx,Ny) = = \^(Nx,Ny)\-L(Nx,Ny). После несложных преобразований получим \ll r.Y .V 1І W,) (6) ^ |h,(W„W,)| Tr(N„NJ ■ ( ) Поскольку в анализируемой системе регистрации шумовая составляющая очень мала, задача коррекции сводится к отысканию линейного корректирующего фильтра с передаточной функцией, определяемой выражением (6). Такой фильтр включается в систему регистрации как отдельное звено и называется обратным фильтром. Обобщенный алгоритм такого способа коррекции предполагает получение Фурье-образа исходного изображения, его преобразования с помощью фильтра с необходимой передаточной функцией и последующего обращения результата для получения обработанного изображения [4-5]. Коррекцию следует выполнить до синтеза изображения на экране монитора. В то же время логично вклю- «Инфокоммуникационные технологии» Том 11, № 1, 2013 Никитин К.А. 91 чить ее в ту часть системы регистрации, где полученные изображения объекта будут оцифрованы и храниться в памяти в виде массивов удобных для выполнения операций, предусмотренных алгоритмом. Поскольку в анализируемой системе используется цветная ТВ-камера, изображение объекта наблюдения, после ввода его в ЭВМ будет храниться в памяти в виде полноцветного RGB-изображения. Такое оцифрованное изображение представляет собой трехмерный массив, фактически состоящий из трех прямоугольных матриц размером MxN, каждая из которых содержит индексированные переменные со значениями, соответствующими относительной яркости красного R, зеленого G и синего B цветов. Элемент массива 1{г,с,п) определяется значениями индексов, где r и с - номер строки и столбца, которые задают положение пиксела, n - номер матрицы, задает доступ к матрице нужного цвета. Введем обозначения: [/(r,c,l)] = [*(r,c)] = R - матрица красной составляющей; R = Д(0,0) Д(0,1) . . . 5(0, W-1) 5(1,0) 5(1,1) • • • R&N-1) _R(M-1,0) Д(М-1,1) . . ,R(M-l,N-l)_ \l{r,c,2j\ = \G(r,cj\ = G - матрица зеленой составляющей; G(0,0) G(0,1) . . . G(0,N-l) G(1,0) G(l,l) . . . G(1,N-1) G(M-1,0) G(M-1,1) . . ,G(M-1,N-1) [/(r,c,3)] = [5(r,c)] = B - матрица синей со ставляющеи; 5(0,0) ( 5(1,0) в ' 5(0,1) . . . B(0,N-1) 5(1,1) . . . B&N-1) В(М-1,0) В(М-1,1) . . .B(M-\,N-\) r = 0;l;2... М-1 - номер строки матрицы; c = 0;l;2... N-1 - номер столбца матрицы; М-число отсчетов цифрового изображения по вертикали (поперек строк); N - число отсчетов цифрового изображения по горизонтали (вдоль строк). Для перехода в частотную область к каждой компоненте изображения следует применить прямое двумерное дискретное преобразование Фурье (ДПФ). Например, для красной составляющей изображения можно записать 1 М-1N-1 S^y,Nx) = — fdfiR(r,c)W1Nyrw^c. (7) М./V г=0 с=о Поскольку комплексное ядро этого преобразования разделимо W(Ny,Nx;r,c) = W^rW?*c ; Щ=ехр(-І^\ N є{0,1, ...М-1} ; ^ М) у ґ 2п\ W2 = exp -i- N и удовлетворяет условию ортогональности [5], выражение (7) можно записать в следующей матричной форме Sr = Ас R Аг, (8) где SH - матрица спектральных коэффициентов ДПФ для R-компоненты; R - матрица значений R-компоненты изображения; Ас - матрица преобразования размером МхМ: J. \[м А. = 4n 1 •К л w° . . • wl и£ w\ wf . . . 4M_1) wj2 w* . . . W?M-V преобразования размером Д w2 w2 w°2 . . . w°2 w° w\ wl . . . Wj w2 w2 . . . wl(N~x) W2 w<2~ V"-4. Таким образом, матрица спектральных коэффициентов двумерного ДПФ для красной компоненты RGB изображения SR получается как результат последовательного выполнения двух одномерных ДПФ: сначала выполняется одномерное преобразование по всем столбцам матрицы отсчетов R (умножение на матрицу Ас слева), затем выполняется второе одномерное преобразование по всем строкам полученной матрицы (умножение на матрицу Аг справа). Размер матрицы получившихся коэффициентов совпадает с «Инфокоммуникационные технологии» Том 11, № 1, 2013 92 Никитин К.А. размером матрицы исходного изображения R и равен MxN. Аналогичные матрицы отсчетов двумерного ДПФ получаются для зеленой G: и синей В: SG=AcGAr SB=AcBAr компонент полноцветного изображения. Поскольку матрица ядра ДПФ симметрична по отношению к ее главной диагонали, матрицы преобразования по столбцам Ас и строкам Аг являются унитарными и по определению Аё* = А*т = Ас; а;1 = а;т = Аг. Действительные унитарные матрицы являются ортогональными и для них справедливы соотношения Aj1 = А^ = Ас; А,1 = А, = Аг. Отсюда следует, что обратный переход от частотного представления к пространственному можно представить в виде матричных выражений R = Ас SR Аг; G “ Ас SG Аг; В = АС SB Аг. (9) Практическая реализация и результаты эксперимента Как следует из выражения (6), для вычисления передаточной функции корректирующего фильтра следует знать двумерные ЧКХ эталонной и реальной систем. Различие методов измерения ЧКХ связано в первую очередь с различием используемых для этой цели тест-объектов. В [6] автором изложен метод вычисления ЧКХ с использованием тестового изображения шумового поля. Следуя изложенным в [6] принципам, вычислим частотно-контрастные характеристики эталонной и реальной систем регистрации. Для этого оцифруем и введем в память ЭВМ эталонное тестовое изображение шумового поля, зарегистрированное без смотровой системы (см. рис. 5). Рис. 5. Эталонное изображение тестового шумового поля размером 512x512 Получим трехмерный массив полноцветного эталонного изображения \_1шэт{г,с,п) |, состоящий из трех прямоугольных матриц R^ , GM3m, Вшэт • Как показали экспериментальные исследования, частотно-контрастные характеристики для красного, зеленого и синего каналов практически совпадают. Это позволяет минимизировать вычисления и использовать общую ЧКХ, рассчитанную по любой из компонент изображения шумового поля, например по красной RMi3m . Вычислим ЧКХ эталонной системы регистрации в вертикальном направлении (поперек строк), выполняя следующие шаги: - выполним одномерное ДПФ по всем столбцам матрицы R,^ в соответствии с (8): с _ A R ШЭТПС С ШЭТП 9 полученная матрица спектральных коэффициентов будет иметь размер MxN и в общем случае состоит из комплексных значений. Для перехода к вещественным значениям вычислим модули спектральных коэффициентов и получим матрицу |®ШЭ(ЯС | > - полученную матрицу модулей спектральных коэффициентов |SM3m(,| можно рассматривать как совокупность из N векторов-столбцов размером Мх 1: 8шэис(с), где c = 0;l..../V-l - номер столбца матрицы. Каждый вектор-столбец представляет собой спектральное распределение реализации соответствующего столбца шумового поля. Вычислим усредненный спектр шумового поля в вертикальном направлении. Для этого выполним усреднение по всем столбцам «Инфокоммуникационные технологии» Том 11, № 1, 2013 Никитин К.А. 93 I N-1 ^шэтс -XJ- ^®шэтс(^) 5 ■А» с=о полученный вектор-столбец будет соответствовать двум полупериодам ЧКХ эталонной системы регистрации в вертикальном направлении TJNy). Аналогичные действия выполним для вычисления ЧКХ по горизонтали: - выполним одномерное ДПФ по всем строкам матрицы Rшэт в соответствии с (8): S — D А шэт г шэт г ’ перейдем от комплексных значений к вещественным и получим матрицу |8ШЭШ). |; - рассматривая матрицу модулей спектральных коэффициентов как совокупность из М векторов-строк размером \xN: 8шэтг (г), где г = 0;1...М-1 - номер строки матрицы, выполним усреднение спектра по всем строкам шумового поля: 1 М-1 с =_N ' с (уЛ шэт г а ж- / j шэт г V /’ МТҐо полученный вектор-строка $шзтг соответствует двум полупериодам ЧКХ эталонной системы регистрации по горизонтали Тэт(Их). Отметим, что вычисление ЧКХ эталонной системы можно выполнить заранее на этапе калибровки. Рис. 6. Реальное изображение тестового шумового поля размером 512x512 Для вычисления ЧКХ реальной системы регистрации оцифруем и введем в память ЭВМ тестовое изображение шумового поля, зарегистрированное в реальных условиях (см. рис. 6). 01-1-1-1-1-1-1-1-1-1-и О 50 100 150 200 250 300 350 400 450 500 Рис. 7. ЧКХ эталонной и реальной систем регистрации в вертикальном направлении Рассматривая трехмерный массив тестового изображения \1ш{г,с,п)\, состоящий из трех матриц , Gm , Вш, и выполняя описанную выше последовательность действий над любой из них, например над Rw , получим вектор-столбец 8ШС и вектор-строку 8ШГ, которые соответствуют двум полупериодам ЧКХ реальной системы регистрации в вертикальном Tp(Ny) и в горизонтальном Tp(Nx) направлениях. Вид полученных нормированных и сглаженных характеристик показан на рис. 7-8. Рис. 8. ЧКХ эталонной и реальной систем регистрации в горизонтальном направлении От одномерных ЧКХ легко перейти к двумерным характеристикам, если найти матричное произведение Т ҐИ N ) = s ¥ =Т • эт\ х, у' шэт с шэт г эт 5 Вид этих характеристик показан на рис. 9-10. Для вычисления передаточной функции корректирующего фильтра в соответствии с полученным ранее выражением (6) следует выполнить деление ЧКХ эталонной системы на ЧКХ «Инфокоммуникационные технологии» Том 11, № 1, 2013 94 Никитин К.А. реальной системы для каждой пространственной частоты. Известно, что операция деления матриц в общем случае не определена. Однако в нашем случае получившиеся матрицы можно рассматривать как упорядоченный набор данных, идентифицируемый с помощью индексов, то есть как массивы данных. Рис.9. Двумерная ЧКХ эталонной системы Рис. 10. Двумерная ЧКХ реальной системы Для массивов одинаковых размеров существует операция поэлементного деления •/ и поэлементного умножения • *, которые и следует использовать в подобных случаях. Таким образом, получим ht =s •/¥ ;ht = s ■/¥ ;Ht=T -/T , kc шэтс uic s кг тэт г ' uiruiy > к эт р 9 где и hkr - корректирующие функции по вертикали и по горизонтали, а Н* - двумерная передаточная функция корректирующего фильтра, показанная на рис. 11. После того как передаточная функция корректирующего фильтра сформирована, остается выполнить собственно коррекцию исходного изображения, например такого, как показано на рис. 12. Видно, что изображение слегка размыто из-за влияния монохроматических аберраций. Коррекгирующая функция Рис. 11. Двумерная передаточная функция корректирующего фильтра Рис. 12. Изображение тестовой таблицы до фильтрации Проведем коррекцию в следующей последовательности: - выполним ДПФ каждой цветовой компоненты исходного изображения в соответствии с (8) и получим три матрицы спектральных коэффициентов S„ , SG и SB ; - выполним поэлементное умножение полученных матриц на матрицу двумерной передаточной функции корректирующего фильтра и получим три матрицы скорректированных спектральных коэффициентов Sr = SR ■ *Н*; SG = SG • *Н*; SB = SB ■ *Ht; - для перехода в пространственную область над полученными матрицами выполним обратное ДПФ в соответствии с (9) и объединим полученные R , G , В' компоненты в цветное скорректированное изображение, которое показано на рис. 13. «Инфокоммуникационные технологии» Том 11, № 1, 2013 Никитин К.А. 95 Заключение Рассмотренный метод коррекции монохроматических аберраций оказался достаточно эффективным для анализируемой системы регистрации. Однако если в данной системе шумы практически не оказывают влияние на изображение, то в системах с высоким уровнем шумов результаты будут отличаться от приведенных. Это является известным недостатком обратной фильтрации. Из-за ограниченного объема в статье не рассмотрены вопросы дополнения матриц спектральных коэффициентов нулями, что необходимо выполнять с целью избавления от краевых искажений [4]. Рассмотренную коррекцию можно выполнить и в пространственной области, если по вычисленной передаточной характеристике корректирующего фильтра синтезировать соответствующий пространственный фильтр[4].

About the authors

K. A Nikitin

Email: nikitin@ist.psati.ru

References

Supplementary files