Обработка видеосигналов для измерения скорости движения железнодорожных вагонов на сортировочной горке

ОБРАБОТКА ВИДЕОСИГНАЛОВ ДЛЯ ИЗМЕРЕНИЯ СКОРОСТИ ДВИЖЕНИЯ ЖЕЛЕЗНОДОРОЖНЫХ ВАГОНОВ НА СОРТИРОВОЧНОЙ ГОРКЕ

Васин Н.Н., Куринский В.Ю. Поволжская государственная академия телекоммуникаций и информатики

Аннотация

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

Постановка задачи Измерение скорости движения железнодорожных вагонов на сортировочной горке предполагает размещение средства регистрации изображения подвижного состава (цифровой видеокамеры) (рис.1). Ширина видимой части пространства

Рис. 1. Геометрическая модель съемки

Образ перемещаемого объекта (вагона или состава), находясь в положении Л1, регистрируется на цифровом изображении так, что точка объекта А1 отображается в точку изображения В1. Через определенное время А/ производится повторная регистрация объекта. В этот момент объект находится в положении Л2, и та же точка объекта, отмеченная на рисунке как А2, отображается на цифровом изображении в точку В2. Цифровые изображения подвижного состава, получаемые устройством регистрации Я, последовательно поступают на вычислительный блок С, осуществляющий автоматический расчет параметров движения состава (рис.1). Результаты вычислений отправляются на удаленное вычислительное устройство (ПЭВМ), цель которого - определение необходимости торможения и степени воздействия на подвижный состав. Видеокамера или отдельная ПЗС-матрица обеспечивают непосредственную регистрацию цифрового изображения, которое передается в дальнейшем для обработки. В качестве изображения

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

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

- промышленную видеокамеру,

- ПЭВМ со специализированной платой ввода видеосигнала.

Основная содержательная часть обработки данных, полученных с видеокамеры, включает в себя решение следующих задач:

- определение наличия состава (вагона) на текущем изображении,

- определение величины смещения состава на смежной паре изображений,

- определение величины смещения состава в реальных координатах,

- определение скорости движения состава.

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

1. Физико-геометрическое описание работы системы

В качестве объекта рассмотрен небольшой прямоугольник площадью Д5"ь нормальный к оси объектива телевизионной камеры (далее для простоты -линза), помещенный на расстояние -а от нее (рис.2). Будем считать, что среда слева и справа от линзы имеет одинаковый показатель преломления п=1. Изображением Д£1 будет прямоугольник площадью Д£2, расположенный на расстоянии Ь от линзы. В параксиальном приближении:

Рис. 2. Объемная геометрическая модель

Если яркость объекта равна Ьч, то через линзу пройдет световой поток ДФу, значение которого определяется выражением:

где Д^2 - телесный угол, под которым линза видна из места расположения объекта. Если обозначить диаметр линзы через ё, то:

пё2

ДФ„ = Lv AS2 AQ2,

AQ„ =-

4b

За счет поглощения света веществом линзы и частичного отражения на ее поверхности поток света за линзой ДФможет только уменьшиться, т.е.

ДФ = к ДФ „,

где к - коэффициент потерь (меньше единицы).

Поток за линзой будет распространяться в пределах телесного угла Д^, под которым линза видна с площадки, являющейся изображением. Из рисунка видно, что

пё2

AQ, =-

4a1

Обозначим через L'v яркость светового потока, формирующего изображение, выражение для потока ДФ'v можно записать следующим образом:

nd2

ДФ' = L' AS,.

v v 4a2 1

Или

L^AS2 = kL' ^AS,.

v 4b2 2 v 4a2 1

L'v = kLv.

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

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

А1

Рис. 3. Горизонтальное сечение

Как видно треугольники д А1А20 и д В1В20 подобны между собой, поэтому очевидно следующее равенство:

Ь = Ь

I а

Из выражения 9 вытекает следующее соотношение между L и l:

l = L-

b

где I - ширина изображения в пикселях;

Ь - длина видимого участка железнодорожного пути, зависит от угла обзора камеры;

а - фокусное расстояние до ПЗС-матрицы;

Ь - расстояние до лицевой (передней) стороны вагона.

Обозначим через д1 единичное значение длины изображения, т. е. сколько приходится на один пиксель единиц длины. Единица длины может быть установлена в зависимости от дальности съемки. Тогда получим:

д1 = Ь —,

т  Ь

где т - это ширина изображения в пикселях.

2. Алгоритмическое описание работы системы

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

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

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

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

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

Преобразование изображения в бинарный вид

— г

Фильтрация полученного изображения

— г

Поиск объекта на изображении по заданным критериям

г

Вычисление смещения вагонов

~ т

Вычисление скорости движения вагонов

Г

Подсистема управления торможением

—Цветное • изображение?

кет

Формирование битовой матрицы изображения

Формирование битовой матрицы изображения

Рис. 4. Алгоритм работы системы

3. Определение величины смещения состава

по смежной паре изображений Движение состава происходит непрерывно во времени и в пространстве. Для сохранения непрерывности происходит квазинепрерывная съемка телевизионной камерой. Непрерывность получается за счет инертности человеческого глаза. Оптимальной скоростью следования считается 25 кадров/с, а время воспроизведения каждого кадра - 40мс. При съемке движения вагона (или другого транспортного средства) изображения двух смежных кадров п и п+1 будут отличаться (рис.5, 6). Это отличие позволило использовать метод межкадровой разности [1, 2] для определения величины смещения.

Суть данного метода заключается в нахождении разности между единичными элементами изображений п и п+1, т.е.:

/х, У)А

где 4(х, у) - значение цвета пикселя на кадре п с координатами (х, у);

^(х, у) - значение цвета пикселя на кадре п+1 с координатами (х, У).

Рис. 6. Кадр п+1

Съемка камерой производится вблизи ж/д пути, что откладывает свой отпечаток на принимаемые изображения (рис. 5, 6). Поэтому одним из существенных факторов является качество снятых кадров. После передачи данных с телевизионной камеры на обрабатывающее устройство производится анализ на наличие объекта исследования. После обнаружения последовательности с необходимыми данными информационные кадры подвергаются обработке с помощью математического алгоритма, который производит поэлементное преобразование. Таким образом, происходит выделение объекта исследования из общей массы предметов, присутствующих на изображении.

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

Область изображения, представляющая наибольший интерес, обладает спектром более светлых от-

Преобразование в черно-белое изображение

Преобразование числовой матрицы в изображение

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

Рис. 7. Разность кадров п и п+1

Рис.8. Поиск объекта с наибольшей площадью на межразностном кадре

4. Анализ межкадровой разности При анализе, не трудно заметить, что на полученном изображении присутствует большое количество криволинейных объектов не нулевого уровня яркости. Среди них находятся объекты, несущие информацию о скорости движения вагона. Чтобы выделить из общей массы объектов наиболее информативный нужно произвести выборку, например, по максимальной площади объекта. Произвести выборку однозначным образом за один шаг практически невозможно. Поэтому следует разделить поиск максимальной площади на два этапа:

- поиск некоторого количества участков максимальной длины хтах 1 и единичной высоты. Здесь под максимальной длиной участка следует понимать непрерывную последовательность единичных битов, ограниченных слева и справа нулевыми;

- поиск соответствующего количества участков максимальной высоты утах 1 и единичной длины. Причем следует понимать, что каждому хтах 1 соответствует единственный утах , имеющий общую точку пересечения.

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

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

Поиск параметров необходим для оценки площади криволинейных объектов. Именно это значащее свойство (площадь) объекта и позволяет сделать вывод о скорости движения вагонов. Конечно, этот вывод будет косвенным, так как значение площади криволинейной трапеции имеет квадратичный характер, который нельзя сопоставлять значению скорости движения, имеющий линейный. Поэтому, чтобы перейти к линейной зависимости наиболее рационально получить интегральное значение ширины выбранной криволинейной трапеции 1ср. Как видно 1ср является изображением смещения вагонов 5. Это расстояние, пройденное вагонами за интервал времени воспроизведения двух соседних кадров А/. Поэтому скорость движения вагонов равна:

V = -

Д/

Заключение

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

Литература

1. Методы компьютерной обработки изображений. Под ред. В.А. Сойфера. М.: Физмалит, 2004.

2. Еремин С.Н., Малыгин Л.Л., Михайлов А.Е., Царев В.А. Опыт использования технологии обработки изображений при проектировании интеллектуальных транспортных систем // Искусственный интеллект в XXI веке: Труды международного конгресса 1СА1 ' 2001. М.: Физмалит, 2001. Т. 2. с. 676 - 690.

5