МЕТОД АВТОМАТИЧЕСКОЙ ЛОКАЛИЗАЦИИ ПРОТЯЖЕННЫХ ГЕОПРОСТРАНСТВЕННЫХ ОБЪЕКТОВ НА КОСМИЧЕСКИХ СНИМКАХ
Метод автоматической локализации протяженных геопространственных объектов 17
УДК 004.93"1, 004.932
С. Ю. МИРОШНИЧЕНКО, В. С. ТИТОВ, А. А. ЯЩЕНКО
МЕТОД АВТОМАТИЧЕСКОЙ ЛОКАЛИЗАЦИИ ПРОТЯЖЕННЫХ ГЕОПРОСТРАНСТВЕННЫХ ОБЪЕКТОВ
НА КОСМИЧЕСКИХ СНИМКАХ
Разработан метод локализации протяженных геопространственных объектов, позволяющий повысить точность автоматизированного формирования тематических слоев электронных карт геоинформационных систем на основе космических снимков. Повышение точности достигается за счет построения фильтра с конечной импульсной характеристикой, учитывающего характерные признаки объектов дорожной сети.
Ключевые слова: геоинформационные системы, космические снимки, объекты дорожной сети, КИХ-фильтр.
Для актуализации дорожных карт и управления городским транспортом в условиях быстрого увеличения транспортных потоков и активного развития сети автомобильных дорог в настоящее время широко используются картографические данные, полученные путем анализа и дешифрировки космических снимков.
Выделение на космических снимках участков дорожной сети, относящихся к классу протяженных геопространственных объектов, является нетривиальной задачей, поскольку их геометрические и яркостные характеристики могут различаться в пределах одного изображения. Распознавание объектов дорожной сети затруднено их перекрытием другими объектами, такими как растительность, здания, тени, атмосферные явления (туман, облачность), транспортные средства. Приведенные факторы значительно усложняют локализацию объектов дорожной сети и построение их векторных представлений как в ручном, так и в автоматическом режиме. Более того, методы выделения контуров и построения кластеров [1], применяемые при автоматической сегментации космических снимков на однородные области, характеризуются невысокой точностью и эффективностью [2].
Повышение точности автоматического выделения заданной области (целевого слоя) дорожной сети на космических снимках достигается путем их сопоставления с векторными топографическими картами местности [2], что возможно лишь при наличии географической привязки, помимо того, возникают ошибки дешифровки, обусловленные погрешностью привязки и инструментальными погрешностями векторизации.
Таким образом, актуальной научно-технической задачей является повышение точности автоматического выделения объектов дорожной сети на космических снимках. Для ее решения разработан метод выделения протяженных геопространственных объектов на космических снимках, основанный на поиске областей изображения, соответствующих пространственным характеристикам объектов целевого слоя [3] с учетом их основных признаков и свойств.
К основным признакам объектов дорожной сети относятся: 1) значительная протяженность на прямолинейном участке; 2) одинаковая ширина участка; 3) равномерное распределение яркости в пределах объекта; 4) четкие контуры дорожного полотна. Для локализации объектов дорожной сети с учетом описанных свойств разработан фильтр с конечной импульсной характеристикой (КИХ-фильтр) [4]. Выбор КИХ-фильтра обусловлен высокой производительностью и простотой реализации. КИХ-фильтр предназначен
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
18 С. Ю. Мирошниченко, В. С. Титов, А. А. Ященко
для аппроксимации „идеального“ импульсного отклика, определяющего принадлежность точек космического снимка к участкам дорожной сети, с применением указанных выше признаков.
Характерные признаки дорожной сети математически описываются следующим образом. Исходное изображение I представлено как дискретное двумерное поле яркости f(·) [1]:
I = f (x, y), x = 1, K , y = 1, N ,
(1)
где x, y — координаты отсчета яркости вдоль осей абсцисс и ординат соответственно; K — ширина, N — высота изображения I.
В пределах поля изображения I прямолинейный участок объекта дорожной сети представляет собой некоторую прямоугольную область d(a,b,α), где b — ширина, a — длина, α — угол поворота дорожного полотна относительно его центра.
Для области d(a,b,α) должны выполняться следующие условия: 1) значительная протяженность, т.е. ширина области d(a,b,α) много меньше ее длины
b ∇f (x, y) .
(5)
Таким образом, область d(a,b,α) считается принадлежащей дорожной сети, если выпол-
нены все необходимые условия (2)—(5). При поиске областей определяется набор значений
функции отклика КИХ-фильтра, полученный путем вращения окна КИХ-фильтра относи-
тельно исследуемой точки изображения, взятой в качестве начала координат.
Значение угла поворота окна фильтра вычисляется по формуле
ϕ
=
180° Nϕ
,
(6)
где Nϕ — параметр, определяющий число секторов, на которое разбивается диапазон значений угла поворота окна фильтра [0°; 180°).
Для нахождения значения яркости точек, принадлежащих окну фильтра с заданным уг-
лом поворота ϕ , необходимо определить их координаты в узлах дискретной сетки исходного
космического изображения. Координаты рассчитываются путем объединения уравнения пря-
мой, которую образует окно фильтра в выбранной системе координат, и уравнения окружно-
сти, описываемой окном фильтра при вращении относительно исследуемой точки
yh = ± xh = ±
kh (1+ k
2
)
,⎫⎪ ⎪
h
⎬ ,⎪
(1+ k 2 )
⎪ ⎭
(7)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
Метод автоматической локализации протяженных геопространственных объектов 19
где (xh, yh) — координаты h-й точки, принадлежащей окну фильтра; k — угловой коэффициент прямой, которую образует окно фильтра в выбранной системе координат; индекс h оп-
ределяет положение точки в окне фильтра относительно исследуемой точки.
В случае, если в результате с помощью уравнений (7) получены вещественные значения
координат, яркость точек окна фильтра вычисляется следующим образом. Для точки с коор-
динатами (xh, yh) в узлах дискретной сетки изображения определяются четыре ближайшие точки hj с координатами (xhj, yhj), где j — индекс смежного узла сетки (рис. 1). Далее вычис-
ляются расстояния вдоль осей абсцисс и ординат ∆xj и ∆yj от смежных узлов дискретной сетки до точки h:
∆x j = xh − xhj , ⎪⎫ ∆y j = yh − yhj ,⎪⎭⎬
(8)
а также значения весовых коэффициентов ωhj(xhj,yhj) узловых точек ωhj (xhj , yhj ) = (1− ∆x j )(1− ∆y j ) .
(9)
Значение яркости h-й точки окна фильтра определяется как
xh yh
∑ ∑fh (xh , yh ) =
ωhj (xhj , yhj ) f j (x j , y j ) ,
jj
где fj(xj,yj) — значение яркости hj-го узла дискретной сетки.
(10)
Пространство точек изображения
h1(xh1, yh1) h2 (xh2, yh2) h(xh, yh)
h3(xh3, yh3) h4 (xh4, yh4)
Направление поворота окна фильтра
Рис. 1
Так как окно фильтра симметрично относительно исследуемой точки, отдельно рассматриваются первая половина окна (рис. 2), к которой относятся точки с положительным смещением координат относительно оси ординат, и вторая половина, к которой относятся точки с отрицательным смещением. Для первой и второй половины окна фильтра рассчитываются значения СКО яркости точек σ1 и σ2 соответственно [5].
Среди группы рассчитанных для различных углов поворота окна КИХ-фильтра σ1 и σ2 определяются минимальные значения σ1min и σ2min и соответствующие им углы ϕ1 и ϕ2.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
20 С. Ю. Мирошниченко, В. С. Титов, А. А. Ященко
Согласно (4), для проверки принадлежности исследуемой точки участку дорожной сети необходимо выполнение следующей совокупности условий:
σ1min < T1, σ2min < T1, σ1min → 0, σ2min → 0,
⎫ ⎪ ⎪ ⎪⎪ ⎬ ⎪
(11)
σ1min − σ2 min
→ 0,⎪ ⎪
ϕ1 ≈ ϕ2 ,
⎪⎭
где T1 — эмпирически выбираемое пороговое значение, определяющее максимально допус-
тимые значения σ1min и σ2min, при которых точка считается принадлежащей участку дорожной
сети.
Учитывая (5), сформулируем дополнительное условие
σ1 − σ1min > T2 , ⎫⎪ σ2 − σ2min > T2 ,⎬⎭⎪
(12)
где σ1 — средняя величина СКО для группы вычисленных значений σ1 различных углов поворота окна фильтра относительно исследуемой точки, σ2 — средняя величина СКО для группы вычисленных значений σ2, Т2 — эмпирически устанавливаемое пороговое значение.
Первая половина
окна фильтра
ϕ3 ϕ2
ϕ1
X
Вторая половина
окна фильтра Y
Рис. 2
Метод выделения протяженных геопространственных объектов на аэрокосмических изображениях заключается в поэлементном обходе изображения в порядке возрастания координат. Текущая точка изображения принимается за начало координат и проверяется на принадлежность объекту целевого слоя, с этой целью в точке определяется значение отклика КИХ-фильтра для каждого из углов его поворота (6). С учетом условий (2) и (3) в качестве окна КИХ-фильтра выбран прямолинейный отрезок, ширина которого равна одной точке растра, длина L устанавливается согласно масштабу аэрокосмического изображения [6, 7]. Количество положений окна фильтра относительно исследуемой точки равно значению параметра Nϕ. Схема поворота окна КИХ-фильтра относительно исследуемой точки представлена на рис. 2.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
Метод автоматической локализации протяженных геопространственных объектов 21
Таким образом, для текущего положения окна фильтра определяется величина ϕ, после чего для каждой точки окна фильтра по формуле (10) вычисляется значение яркости. Если
ϕ = 90°, формулу (10) применить невозможно, так как тангенс угла девяносто градусов не определен. В этом случае окно фильтра параллельно оси ординат и значения яркости точек в пределах окна совпадают с узлами дискретной сетки. Для каждого угла поворота фильтра от-
носительно исследуемой точки вычисляются значения σ1 и σ2 для первой и второй половины его окна соответственно. Также для обеих групп значений σ1 и σ2 определяются σ1min и σ2min, соответствующие им величины ϕ1 и ϕ2, а также σ1 и σ2 . Если для исследуемой точки одно-
временно выполняются условия (11) и (12), значение функции отклика КИХ-фильтра считается положительным.
Условие (11) позволяет определить объекты с равномерным распределением функции яркости. Условие (12) необходимо для исключения площадных объектов с одинаковым уровнем яркости (водоемы, поля, большие здания), не принадлежащих целевому слою, но удовлетворяющих условию (11). Параметры Т1 и Т2 выбираются с учетом яркости и контрастности исходного изображения.
В случае положительного отклика функции КИХ-фильтра текущая точка считается при-
надлежащей объекту целевого слоя. Для нее сохраняется минимальное значение СКО σ′min, равное меньшему из значений σ1min и σ2min, и соответствующее ему значение угла поворота ϕ′. Значения σ′min для всех точек изображения с положительным откликом КИХ-фильтра образуют карту минимальных СКО, величины углов поворота ϕ′ — карту углов поворота окна фильтра.
Построение объектов дорожной сети осуществляется с помощью полученной карты углов поворота окна фильтра с использованием метода дополнения объектов по направлению угла поворота. Для каждой точки карты углов поворота окна фильтра в направлении сохраненного для нее значения угла происходит поиск точек на исходном изображении с яркостью, близкой к яркости текущей точки. Если такие точки найдены, они добавляются на карту углов поворота окна фильтра, уточняя участки дорожной сети.
На рис. 3, а представлен исследуемый фрагмент изображения космического снимка раз-
мером 400×400 точек; б — карта минимальных значений СКО; в — результат применения метода дополнения объектов по направлению угла поворота. Для вычислений использованы
следующие значения параметров: L=10, Nϕ =20, T1=4, T2=4.
а) б)
в)
Рис. 3
На рис. 4, а представлен исследуемый фрагмент изображения аэрокосмического снимка размером 1580×870 точек; б — карта значений углов поворота окна фильтра; в — результат применения метода дополнения объектов по направлению угла поворота. Для вычислений использованы следующие значения параметров: L=40, Nϕ =20, T1=5, T2=2. Данное изображение
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
22 С. Ю. Мирошниченко, В. С. Титов, А. А. Ященко
представляет собой фрагмент космического снимка города и характеризуется наличием большого количества зданий, а также присутствием на дорогах транспортных средств.
а) б)
в)
Рис. 4
Таким образом, предложенный метод позволяет выделять протяженные геопространственные объекты на космических снимках, он обеспечивает автоматическую локализацию прямолинейных участков объектов дорожной сети. В дальнейшем метод будет использован для построения векторных представлений объектов целевого слоя дорожной сети.
СПИСОК ЛИТЕРАТУРЫ
1. Гонсалес Р., Вудс Р. Цифровая обработка изображений. М.: Техносфера, 2005. 1072 с.
2. Лабутина И. А. Дешифрирование аэрокосмических снимков: Учеб. пособие. М.: Аспект Пресс, 2004. 184 с.
3. Кашкин В. Б., Сухинин А. И. Дистанционное зондирование Земли из космоса. Цифровая обработка изображений: Учеб. пособие. М.: Логос, 2001. 264 с.
4. Емельянов С. Г., Кудряшов Е. А., Мирошниченко С. Ю., Титов В. С. Методы и системы цифровой обработки аэрокосмических изображений. Новосибирск: Наука, 2012. 175 с.
5. Кочетков П. А. Краткий курс теории вероятностей и математической статистики: Учеб. пособие. М.: МГИУ, 1999. 51 c.
6. Ященко А. А. Алгоритм локализации объектов дорожной сети // „Распознавание — 2012“: Сб. матер. Х Междунар. науч.-техн. конф. Курск: ЮЗГУ, 2012. С. 278.
7. Мирошниченко С. Ю., Титов В. С. Метод локализации теней геопространственных объектов на цифровых аэрокосмических изображениях // Изв. вузов. Приборостроение. 2010. Т. 53, № 9. С. 7—13.
Сергей Юрьевич Мирошниченко Виталий Семенович Титов Анна Александровна Ященко
Сведения об авторах — канд. техн. наук; Юго-Западный государственный университет, ка-
федра вычислительной техники, Курск; старший преподаватель; E-mail: oldguy7@rambler.ru — д-р техн. наук, профессор; Юго-Западный государственный университет, кафедра вычислительной техники, Курск; заведующий кафедрой; E-mail: titov-kstu @rambler.ru — аспирант; Юго-Западный государственный университет, кафедра вычислительной техники, Курск; E-mail: annajaschenko@rambler.ru
Рекомендовано Юго-Западным государственным университетом
Поступила в редакцию 18.02.2013
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
УДК 004.93"1, 004.932
С. Ю. МИРОШНИЧЕНКО, В. С. ТИТОВ, А. А. ЯЩЕНКО
МЕТОД АВТОМАТИЧЕСКОЙ ЛОКАЛИЗАЦИИ ПРОТЯЖЕННЫХ ГЕОПРОСТРАНСТВЕННЫХ ОБЪЕКТОВ
НА КОСМИЧЕСКИХ СНИМКАХ
Разработан метод локализации протяженных геопространственных объектов, позволяющий повысить точность автоматизированного формирования тематических слоев электронных карт геоинформационных систем на основе космических снимков. Повышение точности достигается за счет построения фильтра с конечной импульсной характеристикой, учитывающего характерные признаки объектов дорожной сети.
Ключевые слова: геоинформационные системы, космические снимки, объекты дорожной сети, КИХ-фильтр.
Для актуализации дорожных карт и управления городским транспортом в условиях быстрого увеличения транспортных потоков и активного развития сети автомобильных дорог в настоящее время широко используются картографические данные, полученные путем анализа и дешифрировки космических снимков.
Выделение на космических снимках участков дорожной сети, относящихся к классу протяженных геопространственных объектов, является нетривиальной задачей, поскольку их геометрические и яркостные характеристики могут различаться в пределах одного изображения. Распознавание объектов дорожной сети затруднено их перекрытием другими объектами, такими как растительность, здания, тени, атмосферные явления (туман, облачность), транспортные средства. Приведенные факторы значительно усложняют локализацию объектов дорожной сети и построение их векторных представлений как в ручном, так и в автоматическом режиме. Более того, методы выделения контуров и построения кластеров [1], применяемые при автоматической сегментации космических снимков на однородные области, характеризуются невысокой точностью и эффективностью [2].
Повышение точности автоматического выделения заданной области (целевого слоя) дорожной сети на космических снимках достигается путем их сопоставления с векторными топографическими картами местности [2], что возможно лишь при наличии географической привязки, помимо того, возникают ошибки дешифровки, обусловленные погрешностью привязки и инструментальными погрешностями векторизации.
Таким образом, актуальной научно-технической задачей является повышение точности автоматического выделения объектов дорожной сети на космических снимках. Для ее решения разработан метод выделения протяженных геопространственных объектов на космических снимках, основанный на поиске областей изображения, соответствующих пространственным характеристикам объектов целевого слоя [3] с учетом их основных признаков и свойств.
К основным признакам объектов дорожной сети относятся: 1) значительная протяженность на прямолинейном участке; 2) одинаковая ширина участка; 3) равномерное распределение яркости в пределах объекта; 4) четкие контуры дорожного полотна. Для локализации объектов дорожной сети с учетом описанных свойств разработан фильтр с конечной импульсной характеристикой (КИХ-фильтр) [4]. Выбор КИХ-фильтра обусловлен высокой производительностью и простотой реализации. КИХ-фильтр предназначен
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
18 С. Ю. Мирошниченко, В. С. Титов, А. А. Ященко
для аппроксимации „идеального“ импульсного отклика, определяющего принадлежность точек космического снимка к участкам дорожной сети, с применением указанных выше признаков.
Характерные признаки дорожной сети математически описываются следующим образом. Исходное изображение I представлено как дискретное двумерное поле яркости f(·) [1]:
I = f (x, y), x = 1, K , y = 1, N ,
(1)
где x, y — координаты отсчета яркости вдоль осей абсцисс и ординат соответственно; K — ширина, N — высота изображения I.
В пределах поля изображения I прямолинейный участок объекта дорожной сети представляет собой некоторую прямоугольную область d(a,b,α), где b — ширина, a — длина, α — угол поворота дорожного полотна относительно его центра.
Для области d(a,b,α) должны выполняться следующие условия: 1) значительная протяженность, т.е. ширина области d(a,b,α) много меньше ее длины
b ∇f (x, y) .
(5)
Таким образом, область d(a,b,α) считается принадлежащей дорожной сети, если выпол-
нены все необходимые условия (2)—(5). При поиске областей определяется набор значений
функции отклика КИХ-фильтра, полученный путем вращения окна КИХ-фильтра относи-
тельно исследуемой точки изображения, взятой в качестве начала координат.
Значение угла поворота окна фильтра вычисляется по формуле
ϕ
=
180° Nϕ
,
(6)
где Nϕ — параметр, определяющий число секторов, на которое разбивается диапазон значений угла поворота окна фильтра [0°; 180°).
Для нахождения значения яркости точек, принадлежащих окну фильтра с заданным уг-
лом поворота ϕ , необходимо определить их координаты в узлах дискретной сетки исходного
космического изображения. Координаты рассчитываются путем объединения уравнения пря-
мой, которую образует окно фильтра в выбранной системе координат, и уравнения окружно-
сти, описываемой окном фильтра при вращении относительно исследуемой точки
yh = ± xh = ±
kh (1+ k
2
)
,⎫⎪ ⎪
h
⎬ ,⎪
(1+ k 2 )
⎪ ⎭
(7)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
Метод автоматической локализации протяженных геопространственных объектов 19
где (xh, yh) — координаты h-й точки, принадлежащей окну фильтра; k — угловой коэффициент прямой, которую образует окно фильтра в выбранной системе координат; индекс h оп-
ределяет положение точки в окне фильтра относительно исследуемой точки.
В случае, если в результате с помощью уравнений (7) получены вещественные значения
координат, яркость точек окна фильтра вычисляется следующим образом. Для точки с коор-
динатами (xh, yh) в узлах дискретной сетки изображения определяются четыре ближайшие точки hj с координатами (xhj, yhj), где j — индекс смежного узла сетки (рис. 1). Далее вычис-
ляются расстояния вдоль осей абсцисс и ординат ∆xj и ∆yj от смежных узлов дискретной сетки до точки h:
∆x j = xh − xhj , ⎪⎫ ∆y j = yh − yhj ,⎪⎭⎬
(8)
а также значения весовых коэффициентов ωhj(xhj,yhj) узловых точек ωhj (xhj , yhj ) = (1− ∆x j )(1− ∆y j ) .
(9)
Значение яркости h-й точки окна фильтра определяется как
xh yh
∑ ∑fh (xh , yh ) =
ωhj (xhj , yhj ) f j (x j , y j ) ,
jj
где fj(xj,yj) — значение яркости hj-го узла дискретной сетки.
(10)
Пространство точек изображения
h1(xh1, yh1) h2 (xh2, yh2) h(xh, yh)
h3(xh3, yh3) h4 (xh4, yh4)
Направление поворота окна фильтра
Рис. 1
Так как окно фильтра симметрично относительно исследуемой точки, отдельно рассматриваются первая половина окна (рис. 2), к которой относятся точки с положительным смещением координат относительно оси ординат, и вторая половина, к которой относятся точки с отрицательным смещением. Для первой и второй половины окна фильтра рассчитываются значения СКО яркости точек σ1 и σ2 соответственно [5].
Среди группы рассчитанных для различных углов поворота окна КИХ-фильтра σ1 и σ2 определяются минимальные значения σ1min и σ2min и соответствующие им углы ϕ1 и ϕ2.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
20 С. Ю. Мирошниченко, В. С. Титов, А. А. Ященко
Согласно (4), для проверки принадлежности исследуемой точки участку дорожной сети необходимо выполнение следующей совокупности условий:
σ1min < T1, σ2min < T1, σ1min → 0, σ2min → 0,
⎫ ⎪ ⎪ ⎪⎪ ⎬ ⎪
(11)
σ1min − σ2 min
→ 0,⎪ ⎪
ϕ1 ≈ ϕ2 ,
⎪⎭
где T1 — эмпирически выбираемое пороговое значение, определяющее максимально допус-
тимые значения σ1min и σ2min, при которых точка считается принадлежащей участку дорожной
сети.
Учитывая (5), сформулируем дополнительное условие
σ1 − σ1min > T2 , ⎫⎪ σ2 − σ2min > T2 ,⎬⎭⎪
(12)
где σ1 — средняя величина СКО для группы вычисленных значений σ1 различных углов поворота окна фильтра относительно исследуемой точки, σ2 — средняя величина СКО для группы вычисленных значений σ2, Т2 — эмпирически устанавливаемое пороговое значение.
Первая половина
окна фильтра
ϕ3 ϕ2
ϕ1
X
Вторая половина
окна фильтра Y
Рис. 2
Метод выделения протяженных геопространственных объектов на аэрокосмических изображениях заключается в поэлементном обходе изображения в порядке возрастания координат. Текущая точка изображения принимается за начало координат и проверяется на принадлежность объекту целевого слоя, с этой целью в точке определяется значение отклика КИХ-фильтра для каждого из углов его поворота (6). С учетом условий (2) и (3) в качестве окна КИХ-фильтра выбран прямолинейный отрезок, ширина которого равна одной точке растра, длина L устанавливается согласно масштабу аэрокосмического изображения [6, 7]. Количество положений окна фильтра относительно исследуемой точки равно значению параметра Nϕ. Схема поворота окна КИХ-фильтра относительно исследуемой точки представлена на рис. 2.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
Метод автоматической локализации протяженных геопространственных объектов 21
Таким образом, для текущего положения окна фильтра определяется величина ϕ, после чего для каждой точки окна фильтра по формуле (10) вычисляется значение яркости. Если
ϕ = 90°, формулу (10) применить невозможно, так как тангенс угла девяносто градусов не определен. В этом случае окно фильтра параллельно оси ординат и значения яркости точек в пределах окна совпадают с узлами дискретной сетки. Для каждого угла поворота фильтра от-
носительно исследуемой точки вычисляются значения σ1 и σ2 для первой и второй половины его окна соответственно. Также для обеих групп значений σ1 и σ2 определяются σ1min и σ2min, соответствующие им величины ϕ1 и ϕ2, а также σ1 и σ2 . Если для исследуемой точки одно-
временно выполняются условия (11) и (12), значение функции отклика КИХ-фильтра считается положительным.
Условие (11) позволяет определить объекты с равномерным распределением функции яркости. Условие (12) необходимо для исключения площадных объектов с одинаковым уровнем яркости (водоемы, поля, большие здания), не принадлежащих целевому слою, но удовлетворяющих условию (11). Параметры Т1 и Т2 выбираются с учетом яркости и контрастности исходного изображения.
В случае положительного отклика функции КИХ-фильтра текущая точка считается при-
надлежащей объекту целевого слоя. Для нее сохраняется минимальное значение СКО σ′min, равное меньшему из значений σ1min и σ2min, и соответствующее ему значение угла поворота ϕ′. Значения σ′min для всех точек изображения с положительным откликом КИХ-фильтра образуют карту минимальных СКО, величины углов поворота ϕ′ — карту углов поворота окна фильтра.
Построение объектов дорожной сети осуществляется с помощью полученной карты углов поворота окна фильтра с использованием метода дополнения объектов по направлению угла поворота. Для каждой точки карты углов поворота окна фильтра в направлении сохраненного для нее значения угла происходит поиск точек на исходном изображении с яркостью, близкой к яркости текущей точки. Если такие точки найдены, они добавляются на карту углов поворота окна фильтра, уточняя участки дорожной сети.
На рис. 3, а представлен исследуемый фрагмент изображения космического снимка раз-
мером 400×400 точек; б — карта минимальных значений СКО; в — результат применения метода дополнения объектов по направлению угла поворота. Для вычислений использованы
следующие значения параметров: L=10, Nϕ =20, T1=4, T2=4.
а) б)
в)
Рис. 3
На рис. 4, а представлен исследуемый фрагмент изображения аэрокосмического снимка размером 1580×870 точек; б — карта значений углов поворота окна фильтра; в — результат применения метода дополнения объектов по направлению угла поворота. Для вычислений использованы следующие значения параметров: L=40, Nϕ =20, T1=5, T2=2. Данное изображение
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6
22 С. Ю. Мирошниченко, В. С. Титов, А. А. Ященко
представляет собой фрагмент космического снимка города и характеризуется наличием большого количества зданий, а также присутствием на дорогах транспортных средств.
а) б)
в)
Рис. 4
Таким образом, предложенный метод позволяет выделять протяженные геопространственные объекты на космических снимках, он обеспечивает автоматическую локализацию прямолинейных участков объектов дорожной сети. В дальнейшем метод будет использован для построения векторных представлений объектов целевого слоя дорожной сети.
СПИСОК ЛИТЕРАТУРЫ
1. Гонсалес Р., Вудс Р. Цифровая обработка изображений. М.: Техносфера, 2005. 1072 с.
2. Лабутина И. А. Дешифрирование аэрокосмических снимков: Учеб. пособие. М.: Аспект Пресс, 2004. 184 с.
3. Кашкин В. Б., Сухинин А. И. Дистанционное зондирование Земли из космоса. Цифровая обработка изображений: Учеб. пособие. М.: Логос, 2001. 264 с.
4. Емельянов С. Г., Кудряшов Е. А., Мирошниченко С. Ю., Титов В. С. Методы и системы цифровой обработки аэрокосмических изображений. Новосибирск: Наука, 2012. 175 с.
5. Кочетков П. А. Краткий курс теории вероятностей и математической статистики: Учеб. пособие. М.: МГИУ, 1999. 51 c.
6. Ященко А. А. Алгоритм локализации объектов дорожной сети // „Распознавание — 2012“: Сб. матер. Х Междунар. науч.-техн. конф. Курск: ЮЗГУ, 2012. С. 278.
7. Мирошниченко С. Ю., Титов В. С. Метод локализации теней геопространственных объектов на цифровых аэрокосмических изображениях // Изв. вузов. Приборостроение. 2010. Т. 53, № 9. С. 7—13.
Сергей Юрьевич Мирошниченко Виталий Семенович Титов Анна Александровна Ященко
Сведения об авторах — канд. техн. наук; Юго-Западный государственный университет, ка-
федра вычислительной техники, Курск; старший преподаватель; E-mail: oldguy7@rambler.ru — д-р техн. наук, профессор; Юго-Западный государственный университет, кафедра вычислительной техники, Курск; заведующий кафедрой; E-mail: titov-kstu @rambler.ru — аспирант; Юго-Западный государственный университет, кафедра вычислительной техники, Курск; E-mail: annajaschenko@rambler.ru
Рекомендовано Юго-Западным государственным университетом
Поступила в редакцию 18.02.2013
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2013. Т. 56, № 6