Способ определения координат центра тяжести изображения

 

Использование: в телевидении, автоматике и измерительной технике, в системах определения координат центра тяжести изображений. Сущность изобретения: проецируется оптический сигнал объекта на фотоприемную матрицу, дискретизируют оптический сигнал, преобразуют оптический сигнал в электрический сигнал Э (m : n), где m 1 - M, а n 1 - N число элементов разложения по ортогональным осям матрицы, определяют координаты центра тяжести изображения относительно оси выбранной системы координат, при этом формируют из полученных электрических сигналов две группы дискретизированных электрических сигналов и две группы интегральных сигналов. Устройство, реализующее способ, содержит телекамеру 1, аналого-цифровой преобразователь 2, накапливающий сумматор 3, блок управления 4, который содержит инвертор 5, первый, второй, третий блоки совпадения 6,7,8, две схемы выделения заднего фронта 9,10, два счетчика 11,12 и десять выходов 13 - 22, восемь накапливающих сумматоров 23, 23, 23, 23, 24, 24, 24, 24, шесть комбинационных сумматоров 25, 25, 25, 26, 26, 26, вычислитель 27. Цель изобретения - повышение точности определения координат центра тяжести изображения. 7 ил.

Изобретение относится к автоматике и измерительной технике и может использоваться в различных контрольно-измерительных устройствах в частности в системах обработки изображений, в устройствах в частности в системах обработки изображений, в устройствах позиционирования.

Известен способ определения координат центра формы плоских неподвижных объектов в декартовой системе координат методом уравновешивания площадей (журнал Зарубежная радиоэлектроника, 1983 г. , N 11, с. 93-97), заключающийся в том, что изображение объекта проэцируют на матрицу фотоприемников, последовательно опрашивают ячейки фотоприемной матрицы и находят абсциссу Х_с = l центра формы из условия h-hh-h0, l= (1) и координату Y_c = k из условия h(i)-h(i)h(i)-h(i)0, к= (2) где ||h|| - норма матрицы, i - номер строки; j - номер столбца; n - размер матрицы; так, чтобы площади фигуры, лежащие по обе стороны от прямой Х = l и прямой Y = k были равны с заданной точностью.

Недостатком этого способа является , во-первых, то, что центр формы определяется неоднозначно, его координаты зависят от ориентации объекта относительно осей координат и не являются инвариантными относительно преобразований изображения; во-вторых, то, что для нахождения центра формы требуется большой объем вычислений и длительное время; в-третьих, то, что способ позволяет работать с изображениями, имеющими только две градации яркости, что ограничивает его функциональные возможности.

Известен способ определения координат центра тяжести изображения плоского объекта (Анисимов Б. А. , Курганов В. Д. , Злобин В. К. Распознавание и цифровая обработка изображений, М. : Высшая школа, 1983 г. , с. 81-89), заключающийся в том, что изображение объекта воспринимается фотоприемным устройством, затем вычисляют площадь фигуры через координаты точек ее контура, находят координаты центра тяжести фигуры путем деления статических моментов инерции на площадь.

Этот способ имеет следующие недостатки: для вычислений используются приближенные формулы, что снижает точность; требуется большой объем вычислений и длительное время для определения координат центра тяжести изображения, работать можно только с изображениями, имеющими две градации яркости, сложность аппаратурной реализации.

Известен способ нахождения центра изображения и центрирования изображения (В. П. Кожемяко. Оптоэлектронные-логиковременные информационно-вычислительные среды, Тбилиси, Мецниереба, 1984 г. , с. 236, рис. 6.31), заключающийся в том, что изображение проэцируют на фотоприемную матрицу, центр которой совпадает с началом декартовой системы координат ХОY, определяют разность яркостей в полуплоскостях ОY и О(-Y) и ОХ и О(-Х), на основании полученных значений определяют направление смещения, производят одновременное смещение всех точек изображения по оси ОХ и по оси ОY на [log₂(2N₁)] и [log₂ (2N₂)] , где [N] обозначает целую часть числа N, заключенного в скобки, и размеры матрицы по осям ОХ и ОY, затем определяют новое значение разностей суммарных яркостей в полуплоскостях и производят смещение изображения на вдвое меньшую величину, чем на предыдущем шаге и продолжают этот процесс до тех пор, пока суммарные яркостные составляющие относительно осей координат не станут равными с точностью до одного дискрета: при этом центр изображения совпадает с центром матрицы, а исходные координаты изображения можно определить по величинам и направлениям производящихся смещений.

Этот способ имеет следующие недостатки: необходимый центр изображения определяют неоднозначно, его положение зависит от ориентации изображения относительно осей координат и он не является инвариантным по отношению к ориентации и преобразованиям изображения (поворот изображения); вследствии инерционного характера центрирования изображения увеличивается время, необходимое для нахождения центра изображения, сложность аппаратурной реализации.

Наиболее близким по технической сущности к предлагаемому способу является способ определения координат центра тяжести изображения неподвижных объектов, заключающийся в том, что изображение объекта воспринимается фотоприемным устройством; 1) осуществляется последовательный опрос ячеек матрицы, 2) определяют массу объекта через площадь изображения как норму матрицы, в которой элементам, принадлежащим изображению, соответствует единица, а не принадлежащим - нуль, путем подсчета ячеек, принадлежащих изображению, 3) определяют статические моменты инерции относительно осей системы координат путем суммирования моментов инерции элементов изображения объекта относительно осей координат и определяют абсциссу Х_с и ординату Y_c центра тяжести путем деления соответствующих моментов инерции на норму матрицы по формулам X_c= Hh (3) Y_c= Hh (4) где ||Н_nn|| - норма матрицы, соответствующая фотоприемной H_nn матрице; i - номер строки; j - номер столбца.

Этот способ имеет следующие недостатки: требуется большой объем вычислений и соответственно длительное время определения координат центра тяжести изображения объекта, способ позволяет работать с изображениями, имеющими две градации яркости; сложность реализации.

Целью изобретения является сокращение объема вычислений/ повышение быстродействия/ расширение функциональных возможностей способа за счет обеспечения возможности работы с полутоновыми изображениями/ упрощение аппаратурной реализации способа.

Указанная цель достигается тем/ что согласно способу/ заключающемуся в проэцировании объекта на фотоприемную матрицу/ с помощью которой преобразуют совокупность оптических сигналов/ полученных в результате пространственной дискретизации в совокупность электрических сигналов Э^o(m; n), / где m = 1 ÷M/ а n = 1 ÷ N ; здесь M и N - число элементов разложения по ортогональным осям матрицы; с амплитудой/ пропорциональной интенсивности соответствующих оптических сигналов/ определении статических моментов инерции изображения относительно осей выбранной системы координат/ формируют на полученных электрических сигналов две группы новых совокупностей сигналов Э(m; n), . . . , Эⁱ_x(m; n). . . , Э) и Э(m; n), . . . , Э^j_y(m; n). . . , Э) по правилу Эⁱ_x(m; n) = (5) здесь [А] - выделение целой части А Э^j_y(m; n) = (6) одновременно формируют две группы интегральных сигналов
S. . . , Sⁱ_x, . . . , S и S. . . , S^j_y, . . . , S и вспомогательный интегральный сигнал S^о, соответствующий интегральной интенсивности изображения объекта, путем интегрирования (накапливания, суммирования) всех сигналов Э_xⁱ(m; n), Э_y^j(m; n), Э^o(m; n) соответствующих совокупностей по правилу
Sⁱ_x= Эⁱ_x(m; n) (7);
S^j_y= Э^j_y(m; n) (8);
S°= Э°(m; n) (9) координату центра тяжести изображения определяют по формулам
X_цт= (10)
Y_цт= (11) для чего известными методами суммируют с учетом весовых коэффициентов интегральные сигналы и производят деление на вспомогательный интегральный сигнал.

Анализ научно-технической и патентной литературы показывает, что применение предложенного метода формирования частных интегральных сигналов и соответствующих им частных (отфильтрованных) изображений позволил отказаться от сложной процедуры умножения отсчетов изображения на соответствующие коды координат отсчетов с целью вычисления статических моментов инерции изображения, заменить операцию умножения операцией сложения и производить операцию накапливания, сложения сигналов одновременно со считыванием, вводом информации и производить обработку по двум каналам одновременно. Такой способ получения и определения моментов изображения был ранее не известен, применение его позволяет повысить быстродействие, сократить время определения координат до времени одного кадра, упростить аппаратурную реализацию и расширить функциональные возможности.

На фиг. 1, 2, 3, приведены исходная (1) и формируемые в процессе обработки две группы совокупностей сигналов. Соответственно первая группа (см. фиг. 2, 3) и вторая группа (см. фиг. 4, 5), которые поясняют сущность способа; на фиг. 6 приведена структурная схема устройства, реализующего предлагаемый способ; на фиг. 7 - функциональная схема устройства управления.

Предлагаемый способ определения координат центра тяжести изображения заключается в следующем.

По известному способу координаты Х_ц.т и Y_ц.т центра тяжести изображения любого объекта вычисляются путем деления соответствующих моментов инерции по оси ОХ и ОY на площадь "массу" изображения. Для примера возьмем квадратную матрицу NxM = 9 элементам разложения. На фиг. 1 показано исходное изображение. (Здесь цифры внутри клеток соответствуют коду яркости светового пятна, т. е. амплитуде электрического сигнала Э^о(m; n), отсутствие цифр соответствует нулевой амплитуде). Определим центр изображения по известному способу. Обозначим статические моменты инерции по оси ОХ через М_х и по оси ОY через М_y. Интегральная яркость изображения составит
S^о = 1 + 3 + 7 + 9 + 10 + 10 + 2 + 4 + 5 + 8 + 11+ + 5 + 17 + 14 + 13 + 6 + 3 + 64 + 3 + 1 + 6 + 1 + + 2 + 1 + 2 + 8 + 7 + 9 = 232
М_х = (1 + 3 + 7 + 9 + 10 + 10)х3 + (2 + 4 + 5 + 8+ + 11 +5)х4 + (17 + 14 + 13 + 6)х5 + (3 + 64 + 3+ + 1)х6 + (6 + 1 + 2 + 1)х7 + (9 + 7 + 8 + 9)х8 = = 1214
М_y = (1 + 2)х3 + (3 + 14)х4 + (7 + 5 + 17 + 3 + 6 + +9)х5 + (9 + 8 + 14 + 64 + 1 + 7)х6 + (10 + 11 + +13 + 3 + 2 + 8 )х7 + (10 + 5 + 6 + 1 + 1 + 2)х8= = 1419
X_цт= = = 5,23;
Y_цт= = = 6,12 На фиг. 2, 3 показана, сформированная первая группа совокупностей соответственно Э_х¹(m; n); Э_х²(m; n); Э_х³(m; n); Э_х⁴(m; n). На фиг. 4, 5 показана сформированная вторая группа совокупностей электрических сигналов соответственно Э_y¹(m; n); Э_y²(m; n); Э_y³⁽m; n); Э_y⁴(m; n). Интегральные сигналы первой группы равны
S_x¹ = 1 + 3 + 7 + 9 + 10 + 10 + 17 + 14 + 13 + 6+ + 6 + 1 + 2 + 1 = 100
S_x² = 1 + 3 + 7 + 9 + 10 + 10 +3+ 64 + 3 + 1 + +6 + 1 + 2 + 1 = 121
S_x³ = 2 + 4 + 5 + 8 + 11 + 5 + 17 + 14 + 13 + 6+ + 3 +64+3+ 1 + 6 + 1 + 2 + 1 = 166
S_x⁴ = 9 + 7 + 8 + 2 = 26 Интегральные сигналы для второй группы равны
S_y¹ = 1 + 2 + 7 + 5 + 17 + 3 + 6 + 9 + 10 + +11 + 13 + 3 + 2 + 8 = 97
S_y² = 1 + 2 + 9 + 8 + 14 + 64 + 1 + 7 + 10 + +11 + 13 + 3 + 2 + 8 = 153
S_y³ = 3 + 4 + 7 + 5 + 17 + 3 + 6 + 9 + 9 + 8+ + 14 + 64 + 1 + 7 + 10 + 11 + 13 + 3 + 2 + 8 = = 204
S_y⁴ = 10 + 5 + 6 + 1 + 1 + 2 = 25 Статические моменты по обеим осям будут соответственно равны
M= ^{= 100}×¹⁺¹²¹×²⁺¹⁶⁶×⁴⁺²⁶×⁸⁼¹²¹⁴
M= ^{= 97}×¹⁺¹⁵³×²⁺²⁰⁴×⁴⁺²⁵×⁸⁼¹⁴¹⁹ И, наконец, координаты центра тяжести изображения будут равны
X_ц.т = = 5,23 Y_ц.т = = 6,12
Рассматриваемый способ предполагает выполнение необходимых операций следующим образом.

Изображение объекта проэцируется на телевизионную камеру на базе ПЗС матрицы. В камере происходит преобразование оптического сигнала в совокупности электрических сигналов и полученный видеосигнал поступает на параллельный быстродействующий аналого-цифровой преобразователь (АЦП). В первом канале посредством потактового суммирования в накапливающем сумматоре получают интегральную яркость изображения объекта S^о, соответствующую "массе" объекта, во втором - путем суммирования в накапливающих сумматорах получают интегральные сигналы первой группы совокупностей, которые складываются в конце кадра с учетом необходимого сдвига в комбнационных сумматорах, получая таким образом значение момента инерции М_х. В третьем канале путем суммирования в аналогичных накапливающих сумматорах получают интегральные сигналы второй группы совокупностей, которые складываются в конце кадра с учетом необходимого сдвига в комбинационных сумматорах, получая таким образом значение момента инерции М_y. Результат, т. е. координаты центра тяжести изображения получают на миниЭВМ путем деления по стандартной программе полученных моментов инерции на интегральную яркость изображения.

Устройство, реализующее способ, содержит телевизионную камеру на базе ПЗС матрицы (типа КТ-2-2) или на базе матрицы фотоприемников (типа МФ-14) 1, быстродействующий параллельный аналого-цифровой преобразователь (АЦП) (типа 1107ПВ1) 2, накапливающие сумматоры (см. например в кн. Применение интегральных микросхем в электронной вычислительной технике, с. 114 - справочник под ред. Б. Н. Файзулаева и Б. В. Тарабрина, М. : Радио и связь, 1986 г. ) (3,23_1-46, 24_1-4); устройство управления 4, комбинационные сумматоры (типа К155ИМ3) (25_1-3, 26_1-3) и микроЭВМ 27.

Видеовыход телекамеры 1 соединен с АЦП 2, выходы АЦП соединены с входами накапливающих сумматоров 3,23_1-4, 24_1-4, выходы накапливающего сумматора 3 соединены с входной шиной микроЭВМ 27. Выходы накапливающего сумматора 23₁, кроме выхода младшего разряда соединены с первыми входами комбинационного сумматора 25₁, причем выход с разрядом 2ⁱ соединен со входом с разрядом 2ⁱ⁺¹, на вход старшего разряда сумматора 23₁ подается логический ноль; выходы накапливающего сумматора 23₂ соединены со вторыми входами комбинационного сумматора 25₁, выходы комбинационного сумматора 25₁, кроме выхода младшего разряда соединены с первыми входами комбинационного сумматора 25₂, причем выход с разрядом 2ⁱ соединен с входом с разрядом 2ⁱ⁺¹, а выход переноса комбинационного сумматора 25₁соединен с входом старшего разряда комбинационного сумматора 25₂; выходы накапливающего сумматора 23₃ соединены со вторыми входами комбинационного сумматора 25₂, выходы комбинационного сумматора 25₂, кроме выхода младшего разряда соединены с первыми входами комбинационного сумматора 25₃, причем выход с разрядом 2ⁱ соединен с входом с разрядом 2ⁱ⁺¹, а выход переноса комбинационного сумматора 25₂соединен с входом старшего разряда комбинационного сумматора 25₃, выходы накапливающего сумматора 23₄ соединены со вторыми входами комбинационного сумматора 25₃; выход нулевого (младшего) разряда накапливающего сумматора 23₁ соединен с входом нулевого разряда входной шины канала Х микроЭВМ 27, выход младшего разряда комбинационного сумматора 25₁ соединен с входом первого разряда входной шины, выход младшего разряда комбинационного сумматора 25₂ соединен с входом второго разряда входной шины микроЭВМ 27, выходы комбинационного сумматора 25₃соединены с входами входной шины микроЭВМ 27 таким образом, что выход с разрядом 2ⁱ соединен с входом с разрядом 2ⁱ⁺³, а выход переноса соединен с старшим разрядом входной шины канала Х микроЭВМ 27. Выходы накапливающих сумматоров 24_1-4 соединены с комбинационными сумматорами 26_1-3 и с входами входной шины канала Y микроЭВМ 27 аналогично. Служебные (управляющие) выходы телевизионной камеры 1 соединены с входами устройства управления 4. Устройство управления содержит инвертор 5, девять схем совпадения 6, 7_1-4, 8_1-4, две схемы выделения заднего фронта (СВЗФ) 9, 10, и два двоичных счетчика 11 и 12.

Выход, с которого снимается строчный синхроимпульс (ССИ) телевизионной камеры 1 соединен с входом инвертора 5, суммирующим входом счетчика 12 и входом установки в ноль счетчика 11, выход инвертора 5 соединен с первым входом схемы совпадения 6, второй вход которой соединен с выходом телевизионной камеры, с которого снимается сигнал тактовой частоты F_так, этот же выход телевизионной камеры 1 соединен с тактовым входом АЦП 2, выход схемы совпадения 6 соединен с суммирующим входом счетчика 11 и СВЗФ 9, выход СВЗФ 9 соединен с первыми входами схем совпадения 7_1-4 и 8_1-4 и входом записи накапливающего сумматора 3, выход 14, вторые входы этих схем совпадения соединены с выходами счетчиков 11 и 12 соответственно выходы счетчика 11 - с входами схем совпадения 7_1-4, а выходы счетчика 12 - с входами схем совпадения 8_1-4, выход телевизионной камеры 1, с которого снимается кадровый синхроимпульс (КСИ), соединен с входом СВЗФ 10 и входом "ЗАПИСЬ" микроЭВМ 27, выход СВЗФ 10 соединен с входам установки в ноль счетчика 12 и входами "Установка в ноль" накапливающих сумматоров 23_1-4 и 24_1-4, выходы схем совпадения 7_1-4 и 8_1-4 - 14-22 соединены с входами "ЗАПИСЬ" накапливающих сумматоров 23_1-4 и 24_1-4 соответственно.

Устройство работает следующим образом: задним фронтом КСИ происходит сброс накапливающих сумматоров и счетчика управления 12, счетчик управления 11 сбрасывается передним фронтом ССИ. С началом кадра видеосигнал с телевизионной камеры 1 поступает на АЦП 2, где потактово оцифровывается и поступает на входы всех накапливающих сумматоров 3,23_1-4, 24_1-4. В накапливающем сумматоре 3 происходит потактовое суммирование всех значений оцифрованного видеосигнала и к концу кадра в нем будет храниться число, соответствующее значению вспомогательного интегрального сигнала S^о, т. е. интегральной интенсивности объекта. В накапливающих сумматорах 23_1-4 канала Х происходит суммирование значений оцифрованного видеосигнала и формирование интегральных сумм. Прохождение тактовых сигналов на входы записи накапливающих сумматоров 23_1-4регламентируется состоянием выходов счетчика управления 11, состояние которых изменяется с тактовой частотой. К концу кадра в накапливающих сумматора 23_1-4 хранятся интегральные совокупности S_xⁱ, а на входной шине канала Х микроЭВМ 27 к концу кадра в результате сложения в комбинационных сумматорах 25_1-3 получают сумму, соответствующую статическому моменту инерции по оси ОХ - М_х. Статический момент М_yвычисляется в канале Y аналогично, только прохождение тактовых сигналов на входы записи накапливающих сумматоров 24_1-4 регламентируется состоянием выходов счетчика управления 12, состояние которых изменяется с частотой строк. По переднему фронту КСИ происходит считывание значений М_хМ_y и S^о в микроЭВМ 27. В микроЭВМ выполняется операция деления, т. е. вычисление координат центра тяжести изображения, которая занимает время обратного хода кадровой развертки. С началом кадра цикл повторяется. Разрядность комбинационных сумматоров и количество накапливающих сумматоров выбирается в зависимости от предполагаемого размера объекта и количества светочувствительных ячеек фотоприемной матрицы.

Докажем, что поставленная цель достигнута. Во-первых, что касается сокращения объема вычислений, то по сравнению с прототипом, в котором требуется производить умножение сигналов отсчетов элементов изображения на координаты, а таких элементов всего N х М, в предложенном способе операция умножения отсутствует. Это приводит к тому, что не требуется задержка тактовых сигналов, синхронизирующих работу фотоприемной матрицы на время умножения, и формирование интегральных сигналов, на основе которых вычисляются статические моменты, происходит в темпе считывания видеосигнала. Кроме того, возможен вариант реализации способа, когда формирование всех частных изображений (групповых совокупностей сигналов) может происходит параллельно при оптических методах формирования. В этом случае применяются фильтры (транспаранты) и интегрирующие по пространству элементы, (линзы, фоконы и т. д. ).

Во-вторых, что касается повышения быстродействия, то время формирования всех интегральных сигналов (осуществляемое с высокой точностью, поскольку применяется предварительная оцифровка видеосигнала) идет одновременно для двух каналов и заканчивается к концу кадра, что позволяет за время обратного хода кадровой развертки (порядка 1 мсек) определить координаты центра тяжести изображения.

В-третьих, что касается упрощения аппаратурной реализации, то способ предполагает наличие только формирователей интегральных сигналов (интеграторы для аналоговых сигналов и накапливающие сумматоры для цифровых), включая простые коммутаторы сигналов для выделения по требуемым условиям нужных элементов и матрицу фоточувствительных элементов (например телекамера на ПЗС матрице). При этом не требуется сложных вычислительных устройств для умножения и управления алгоритмом обработки как в прототипе. Вычислительный блок требуется лишь на конечном этапе вычислений и может быть выполнен в виде табличного делителя ПЗУ.

В-четвертых, что касается расширения функциональных возможностей, то предлагаемый способ позволяет определить координаты центра тяжести не только бинарных, но и полутоновых изображений, с большой точностью производить обработку изображений (до долей элемента ПЗС матрицы) и учитывать фон, посредством пороговой обработки видеосигнала на АЦП. (56) Жаботинский Ю. Д. Адаптивные промышленные роботы и их применение в микроэлектронике, М. , Радио и связь, 1985 г. , с. 31.

Формула изобретения

СПОСОБ ОПРЕДЕЛЕНИЯ КООРДИНАТ ЦЕНТРА ТЯЖЕСТИ ИЗОБРАЖЕНИЯ, при котором проецируют оптический сигнал объекта на фотоприемную матрицу, дискретизируют оптический сигнал, преобразуют оптический сигнал в дискретизированный электрический сигнал Э^o(m, n), где m = 1 ÷ M, а n = 1 ÷ N - число элементов разложения по ортогональным осям матрицы, с амплитудой, пропорциональной интенсивности соответствующих оптических сигналов, определяют координаты центра тяжести изображения относительно оси выбранной системы координат, отличающийся тем, что, с целью повышения точности определения координат центра тяжести, после преобразования аналогового сигнала в дискретизированный электрический сигнал формируют из полученных электрических сигналов две группы дискретизированных электрических сигналов
Э(m; n), . . . , Эⁱ_x(m; n). . . , Э),
Э(m; n), . . . , Э^j_y(m; n). . . , Э)
из исходной совокупности по правилу
,
Э^j_y(m; n) = ,
где Э^o(m, n), Э_xⁱ(m, n), Э_yⁱ(m, n) - электрические сигналы соответственно исходной первой и второй совокупностей,
одновременно с формированием двух групп дискретизируемых электрических сигналов формируют две группы интегральных сигналов
S. . . , Sⁱ_x, . . . , S и S. . . , S^j_y, . . . , S
и вспомогательный интегральный сигнал S⁰, соответствующий интегральной интенсивности изображения объекта, путем суммирования всех сигналов соответствующих совокупностей по формуле
Sⁱ_x= Эⁱ_x(m; n) ; S^j_y= Э^j_y(m; n) ; S°= Э°(m; n),
где S^o₁, Sⁱ_x, S^j_y - интегральные сигналы соответственно из исходной первой и второй совокупностей,
а сигнал координаты центра тяжести изображения определяют по формулам
X_цт= ;
Y_цт= .

РИСУНКИ

Рисунок 1, Рисунок 2, Рисунок 3, Рисунок 4, Рисунок 5, Рисунок 6, Рисунок 7