Применение локального метода обнаружения краев изображений для восстановления профиля земляного полотна the application of a local method of a images borders detection for recovering a roadbed profile

У

правление техническими системами и технологическими процессами

УДК 519,254

ПРИМЕНЕНИЕ ЛОКАЛЬНОГО МЕТОДА ОБНАРУЖЕНИЯ КРАЕН ИЗОБРАЖЕНИЙ ДЛЯ ВОССТАНОВЛЕНИЯ ПРОФИЛЯ ЗЕМЛЯНОГО ПОЛОТНА1

А.Г. Броневич, А.Н. Каркищенко, В.И. Уманский, Д.А. Якушев

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

Ключевые слова: земляное полотно, поперечный профиль, сканированные данные, случайные последовательности, отношение правдоподобия.

ВВЕДЕНИЕ

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

К числу наиболее простых и широко применяемых методов мониторинга относятся так называ-

1 Работа выполнена при финансовой поддержке РФФИ, проект № 11-07-13125-офи-м-2011-РЖД.

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

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

Эффективность диагностики земляного полотна во многом определяется такими свойствами, как объективность, достоверность, простота проведения, оперативность, дешевизна. Требуется

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

— ненормативная ширина плеча балластной призмы;

— ненормативная ширина обочины земляного полотна;

— места с завышенной крутизной откосов;

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

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

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

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

В настоящей работе использовались реальные данные сканирования, полученные с помощью технических комплексов наземного лазерного сканирования, устанавливаемых на специально оборудованном поезде [5]. Технологически это реализуется путем установки двух лазеров на задней

Рис. 1. Схема измерений земляного полотна

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

1. ОПИСАНИЕ И ОБЩАЯ ПОСТАНОВКА ЗАДАЧИ

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

— совокупность измерений, называемых в дальнейшем замерами, представляет собой конечное множество точек, число которых определяется частотой излучения лазера и может варьироваться в зависимости от его технических возможностей;

— за один проход лазеров получается совокупность замеров по общему углу охвата в 360° градусов (по 180° на каждый лазер);

— распределение точек замеров существенно зависит от высоты расположения лазера над поверхностью земли; поскольку лазеры устанавливаются на задней стенке последнего вагона, то считается, что высота подвеса лазеров равна 4—5 м над рельсами;

— характер и совокупность данных зависят от того, перемещается ли поезд в открытом пространстве (насыпь), полуоткрытом (выемка) или закрытом (например, тоннель);

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

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

— данные реальных измерений оказываются сильно зашумленными; ошибки измерений воз-

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

Далее будем считать, что выполняются предположения:

• ошибка каждого измерения в одном цикле появляется в результате действия большого числа независимых или слабо зависимых факторов (случайных величин), поэтому она может быть достаточно точно описана случайной величиной, подчиненной нормальному закону распределения;

• случайные величины имеют короткий «промежуток зависимости». Иными словами, случайные возмущения слабо связаны с их предысторией, и поэтому статистическая зависимость двух измерений быстро уменьшается с ростом расстояния между точками замеров. Перечисленные выше факторы, влияющие на возникновение ошибок в измерениях, позволяют считать данное предположение вполне реалистичным. Оно гарантирует, что получаемые на основе таких измерений стандартные статистические оценки являются асимптотически гауссовскими в смысле центральной предельной теоремы [6].

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

2. МОДЕЛЬ КРАЯ

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

1

описывается гауссовой функцией

2/ол2

-х /2А

л/2П а

Пусть идеальный край описывается кусочно-линейной функцией

Г а 1 х, х < 0,

£(х а2) =

[(а2 + а1)х, х> 0.

Рис. 2. Сглаженная модель края

Тогда сглаженную модель края можно получить

с помощью интеграла свертки

+ >2

+ х ----

1 2 А

/(х|а, а2) = [ £(х - 0 —=- е Л = *' л/2 пА

-х у

х -—

= Г (02 + аХх - О —1— е 2А Ж + Дп А

-х

+ ?

+ х --2

+ Г а (х — /) —1— е 2 А л = > 1 л/2п А

А -

= ах + х Ф(х/А) + —— е

дп

-х2/2А2

где Ф — функция Лапласа. Отметим, что если в последней формуле а2 = 0, то /(х|а1, а2) = ^(х|а1, а2) = = а1х. Кроме того, если А ^ 0, т. е. сглаживание края отсутствует, то также /(х|а1, а2) = g(х|а1, а2). Мы получим общую модель края, если с помощью функции /(х) будем описывать различные положения краевой точки х0 вдоль координаты х, и высоты точки с. Тогда

/

А х2/2а2

где ф(х) = хФ(х/А) + -е-х . Сглаженная мо-

дп

дель края для с = 1, х0 = 1, а1 = 1, а2 = —1 и А = 0,1 изображена на рис. 2.

3. ВЕРОЯТНОСТНАЯ МОДЕЛЬ КРАЯ. ДИСКРИМИНАНТНАЯ ФУНКЦИЯ

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

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

(у-я*»2

Ку|х) = -к- е 2°2 ,

л/2 па

где /(х) — модель края. Далее дискретные значения высот в поперечном сечении профиля будем описывать с помощью случайной последовательности вида (У1(х1), 72(х2), ..., У^хм)), в которой последовательно записаны значения высот, определенные в точках х1, х2, ..., хм. С учетом этого, с помощью

2

(У - f( x,-) ) -

функции h(y) = —1— e 2ст будем описывать лДлСТ

плотность вероятности случайной величины Y.. При анализе конечных участков профиля возможна ситуация, когда этот участок описывается линейной функцией, т. е. на нем нет краевых точек. В этом случае значения высот можно описать случайной функцией f (x) = b + ax + s, и на этом участке случайные величины можно охарактеризовать

(y-b-ax,)2

плотностью h.(y) = —1— e 2 ст .

J2n а

Рассмотрим идею локального обнаружения краевых точек в последовательности (Y1(x1), Y2(x2), ..., YN(xN)). Внутри отрезка [x. - m, ..., x. + m] строится функция правдоподобия гипотезы «точка x; является точкой края», против альтернативной гипотезы «внутри указанного отрезка нет краевых точек». Отметим, что это один из подходов выбора диск-риминантной функции, которая по определению должна принимать большие значения вблизи краевых точек и относительно небольшие значения на участках, где нет краевых точек.

Найдем аналитическое выражение для дискри-минантной функции, основанной на модели края (1). Будем считать, что поперечный профиль задан последовательностью точек {(x1, y1), ..., (xN, yN)}. В этом функция правдоподобия первой гипотезы определяется равенством

R1(xJx,- - m, ..., xi + m, yi - m, .", yi + m) =

= max

a^, «2, "

П

k = i- m

[yk - c - a 1 (xk - x¡) - a2 y ( xk - x¡) ]

2 ст

-У2П

Также нетрудно получить функцию правдоподобия для второй гипотезы

^2(хгК' - т, ..., хг + т, Уг - т, ..., Уг + т) _

[ у к - Ь - а хк ] 2 г + т ^ 2

П1 2 а

-е .

Ьа к = г-т*/2Ла

Далее дискриминантная функция строится как логарифм отношения функций правдоподобия гипотез:

р

¿(х,) = 1п -1 = Р2

a i, a2,

k = i - m

2 ст'

2

. г + m [y b-a

a, b ^ 2_2

k = i-m 2CT

Можно показать [8], что рассмотренная функция, основанная на отношении правдоподобия гипотез, обладает всеми требуемыми свойствами, необходимыми для дискриминантных функций. Ми-нимизационные задачи в правой части последнего равенства имеют явное решение, которое можно получить, решая системы линейных алгебраических уравнений:

относительно ax, a2 и с:

i + m i + m

c(2m + 1) + ax V (xk — x) + a2 V ф(хк — x) =

k = i-m k = i-m

i + m

= V y*

k = i-m

c Z (xk - x) + ai Z (xk - x) +

k = i-m i + m

k = i-m

+ a2 Z (xk - x>(xk - x) = Z (xk - x)yk'

k = i-m

k = i-m

Z ^(xk - x/) + ai Z (xk - x>(xk - x) +

k = i-m

k = i-m

2

+ a2 z ф (xk - x/) = Z ^(xk - x)yk;

k = i-m k = i-m

относительно a и b:

b(2m + 1) + a Z (xk - x) = Z yk'

k = i-m

k = i-m

г + m г + m

2

ПСТ

b Z (xk- xi)+a Z (xk- x) = Z (xk- x)yk-

k = i-m

k = i-m

k = i-m

С

e

4. ОЦЕНИВАНИЕ ПОЛОЖЕНИИ КРАЕВЫХ ТОЧЕК. ТРАНЗИТИВНОЕ ЗАМЫКАНИЕ ДИСКРИМИНАНТНОИ ФУНКЦИИ

Положение краевых точек можно оценить, анализируя локальные максимумы дискриминантной функции. Более достоверные оценки получаются, если воспользоваться функцией, которая называется транзитивным замыканием дискриминантной функции. Эту функцию можно вычислить по формуле

рЦ.) = тах Е[х^ ^

х( е [хк хт]

где Е[хк, хт] — среднее значение дискриминантной функции на отрезке [хк, хт], а максимум в формуле (3) берется по всем отрезкам [хк, хт], содержащим точку хг. Таким образом, по определению

1 Хт

Е хт] = х""1^ 1 Х(х)^х.

т хк

Отметим, что поскольку значение Ь известно, только в дискретных точках хт + 1, ..., х^ _ т интеграл вычисляется по приближенным формулам. В частности, если точки х; расположены равномерно с шагом АА, т. е. хг = хг _ 1 + АА, то приближенно значение Е [хк, хт] можно рассчитать по формуле

Е [хи хт] =

1

т - к + 1

I Ь(х).

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

5. РЕАЛИЗАЦИЯ МЕТОДА ДЛЯ ПРОФИЛЯ, ЗАДАННОГО АНАЛИТИЧЕСКИ

Ранее уже отмечалось, что профиль в результате предварительной обработки может быть задан в виде кусочно-линейной аппроксимации, которая получается, если мы воспользуемся триангуляцией. В этом случае анализируемый профиль / можно задать последовательностью значений /(хк) в точках хк, к = 1, ..., N причем следует выбрать равномерное распределение точек хк на оси х, т. е. х; = х;_ 1 + АА, / = 2, ..., N. Другая альтернатива заключается в обобщении дискриминантной функ-

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

Ь

а^ «2, с ) 1

х - я

х + я

— а2ф(и — х))2Аи + тт Г (/(и) — Ь -

2 а, Ь Л

а(и — х)) Аи.

Отметим, что здесь в качестве окрестности точки х, по которой вычисляется дискриминантная функция, выбирается отрезок [х — л, х + л]. Данная дискриминантная функция является обобщением дискриминантной функции, рассмотренной в § 3, так как фактически при расчете дискриминант-ной функции мы аппроксимируем анализируемый профиль на участке [х; _ т, х; + т] по методу наименьших квадратов, подбирая параметры функций /(ха2, с, х;) = с + а,(х — х;) + а2ф(х — х;) и /(х |а, Ь) = Ь + ах, т. е. формула (4) является интегральным аналогом формулы (2). При расчете дискриминантной функции (4) решение оптимизационных задач также связано с решением систем линейных уравнений. В данном случае значения а,, а2 и с находятся из системы:

х + я х + я

с | Аи+а, | (и—х)Аи +а2 | ф(и—х)Аи = | /(и)Аи,

х - я х - я

с | (и — х)Аи + а, | (и — х) Аи +

х - я х + я

+ а2 | (и — х)ф(и — х)Ам = | (и — х)/(м)Ам,

х + я

с | ф(и — х)Ам + а, | (и — х)ф(и — х)Ам +

х - я х - я

х + я х + я

+ а2 | ф2(и — х)Ам = | ф(и — х)/(м)Ам.

х - я х - я

Данная система после несложных преобразований приводится к виду:

я я

с2л + а2 |ф(м)Ам = |/(и + х)Ам,

-я -я

я я

а,(2/3)л + а2 | Иф(м)Ам = | м/(м + х)Аи

т

я

£

с |ф(и)Си + а1 |иф(и)Си + а2 |ф2(и)Си

-5 -5

5

= | ф(ы)/(ы + х)Си.

Отметим, что поскольку функция

Ф(х) - 0,5х = х(Ф(х/А) - 0,5) + -А. е-*2/2д2

л/2л

является четной, то

, то | и(ф(и) — 0,5и)Си = 0. Поэтому

|иф(и)Си = 0,5 11?с1и = (1/3)£3.

Аналогичным образом находим а и Ь из системы уравнений

X + 5 X + 5 X + 5

Ь | Си + а | (и — х)Си = | /(и)Си,

X + 5 X + 5 X + 5

2

Ь | (и — х)Си + а | (и — х) сСи = | (и — х)/(и)Си.

1 5

Отсюда находим, что Ь = — Г/(и — х)Си и

2 £ ^

3

а = —- | и/(и + х)Си.

6. ВОССТАНОВЛЕНИЕ ПРОФИЛЯ ЗЕМЛЯНОГО ПОЛОТНА

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

g(х|с, а1, ..., ак, х1, ..., хк) = к

= с + а^х + £ а,(х — х)х(х — х),

г = 1

где х(х) =

0, х< 0,

1, х> 0.

Сглаженный профиль земляного полотна можно получить с помощью оператора свертки с ядром

1 ^2/2А2

е :

Т2ПА

g(х|с, а1, ..., ак, х1, ..., хк) =

1 2 А 2

[ g(х — /|с, а1, ..., ак, х1, ..., хк)-—- е С/1 =

л/2лА

к

= с + а0х + £ а,.ф(х — хг).

г = 1

Будем считать положения точек переключения известными, а для оценки параметров с, а1, ... , ак воспользуемся методом наименьших квадратов в интегральной форме. Пусть / — наблюдаемый профиль на отрезке [х0, хк + 1]. Тогда нахождение оценок параметров с, а1, ..., ак сводится к решению следующей оптимизационной задачи:

О = тт Г (/(х) — g (х|с, а1, ..., а

с, а!,..., ак *

к

2

х1, ..., хк)) Сх.

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

с(хк + 1 — ххо) + ао0,5( хк + 1 — х2) +

к ^ +1 а

xk + 1

+ £ а, | ф(х — х,)Сх = | /(х)Сх.

г = 1

с0,5(хк + 1 х0 ) + а0 3(хк + 1 х0 ) +

+ £ аг | хф(х — х)Сх = | х/(х)Сх,

г = 1 x0 x0

xk + 1 xk + 1 с | ф(х — х)Сх + а0 | хф(х — х.)Сх +

к xk + 1

+ £ аг | ф(х — х,.)ф(х — х)Сх = | ф(х — хД/(х)Сх,

г = 1

] = 1, ..., к.

7. ЭКСПЕРИМЕНТАЛЬНЫЕ РЕЗУЛЬТАТЫ

Проверка работоспособности метода восстановления профиля земляного полотна проводилась на модельных и экспериментальных данных. Для генерации модельных данных были выбраны профили земляного полотна, отвечающие всем нормативным характеристикам в масштабе 1:1, высота подвеса лазерных сканирующих устройств варьировалась от 4 до 6 м. Статистические данные представляли собой последовательность замеров {(фг, рг)} в полярной системе координат с полюсом в точке подвеса сканера, причем предполагалось равномерное изменение угла, т. е. фг = фг _ 1 + Аф, а шумовая составляющая, распределенная по нормальному закону с среднеквадратическим отклонением а, добавлялась к значению полярного радиуса. Было установлено, что метод дает достаточно устойчивые результаты при а = 0,1 и Аф = 0,001 рад. Значение

5

5

5

5

X

k + 1

5

X

0

X

X

0

0

X

k + 1

-5

X

X

0

0

X

k + 1

X

X

0

0

+

у

-да

Рис. 3. Идеальный профиль

Рис. 4. Дискриминантная функция и ее транзитивное замыкание

Рис. 5. Статистические данные и восстановленный профиль £

Аф = 0,001 рад соответствует тому, что на анализируемых профилях получается порядка 1000 точек. Подчеркнем, что метод начинает плохо работать, если часть профиля попадает в тень, т. е. высота подвеса оказывается недостаточной. В этом случае можно восстановить работоспособность метода, выявляя затененные участки. Результаты моделирования представлены на рис. 3—5.

ЗАКЛЮЧЕНИЕ

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

рассчитывать «усредненный» профиль внутри полосы, параллельной поперечному сечению земляного полотна. Учитывая характеристики современных сканеров, ширина полосы должна быть в пределах 2—3 м. Другая проблема заключается в наличии теней, когда требуется выделение областей, где восстановление профиля невозможно. Дальнейшее развитие данного метода связано с сегментацией поверхности земляного полотна, т. е. с выделением участков, где поверхность может быть достаточно хорошо аппроксимирована плоскостью. Это позволит выделять участки земляного полотна, имеющие на всей своей протяженности приблизительно одинаковый поперечный профиль, а также участки, где такая аппроксимация невозможна, т. е. участки, где наблюдается разрушение земляного полотна.

ЛИТЕРАТУРА

1. Сазонов В., Кармалин Б., Лебедев А. и др. Современные направления диагностики и мониторинга земляного полотна // Путь и путевое хозяйство. — 2009. — № 6. — С. 34—37.

2. Ашпиз Е.С. Мониторинг эксплуатируемого земляного полотна: Теоретические основы и практические решения: автореф. дис. ... д-ра техн. наук. — М., 2002.

3. Дулина Н.Г., Уманский В.И. Структуризация проблемы улучшения пространственной согласованности баз геоданных // Сообщения по прикладной математике ВЦ РАН. — 2009. — С. 3—22.

4. Дулин С.К., Розенберг И.Н., Уманский В.И. Методы кластеризации в исследовании массивов геоданных // Системы и средства информатики. — 2009. — Доп. вып. — С. 86—114.

5. Уманский В.И. Технология построения трехмерных моделей железнодорожного полотна в высокоточном координатном пространстве // Сб. докл. 6-й междунар. науч.-практ. конф. «Геопространственные технологии и сферы их применения». — М., 2010. — С. 66—67.

6. Бриллинджер Д. Временные ряды. Обработка данных и теория. — М.: Мир, 1980.

7. Bronevich A.G, Itenberg I.I, Karkishchenko A.N. The application of a local method for edge detection // Proc. of the 6th Intern. Conf. on Control, Robotics and Vision, Singapore, 2000.

8. Броневич А.Г., Каркищенко А.Н. Вероятностные и возмож-ностные модели классификации случайных последовательностей. — Таганрог: ТРТУ, 1996.

Статья представлена к публикации членом редколлегии А.С. Манделем.

Броневич Андрей Георгиевич — д-р физ.-мат. наук, гл. науч. сотрудник, ОАО «Научно-исследовательский и проектно-конструкторский институт информатизации, автоматизации и связи на железнодорожном транспорте», г. Москва, Ш (495) 967-77-06, доб. 3-29, И brone@mail.ru, Каркищенко Александр Николаевич — д-р физ.-мат. наук, руководитель центра, ОАО «Научно-исследовательский и проектно-конструкторский институт информатизации, автоматизации и связи на железнодорожном транспорте», г. Москва, Ш (495) 967-77-06, доб. 3-94, Уманский Владимир Ильич — канд. техн. наук, ген. директор, ЗАО «ИнтехГеоТрана>, г. Москва, Ш (495) 967-77-06, доб. 109, И umanvi8@yandex.ru, Якушев Дмитрий Алексеевич — нач. отдела, ЗАО «ИнтехГеоТрано, г. Москва, Ш (495) 967-77-06, доб. 6-13, И D.Yakushev@gismps.ru.