close

Вход

Забыли?

вход по аккаунту

?

Патент BY16426

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(46) 2012.10.30
(12)
(51) МПК
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
BY (11) 16426
(13) C1
(19)
G 01C 5/00
G 01C 3/08
(2006.01)
(2006.01)
СПОСОБ ОПРЕДЕЛЕНИЯ ДАЛЬНОСТИ ДО ОБЪЕКТА,
ЕГО ВЫСОТЫ И ШИРИНЫ
(21) Номер заявки: a 20110038
(22) 2011.01.10
(43) 2011.08.30
(71) Заявитель: Белорусский государственный университет (BY)
(72) Авторы: Козлов Владимир Леонидович; Кузьмичев Илья Романович
(BY)
(73) Патентообладатель: Белорусский государственный университет (BY)
(56) BY 12740 C1, 2009.
RU 94039372 A1, 1996.
JP 1173811 A, 1989.
ЛОБАНОВ А.Н. Большая советская
энциклопедия. - М.: Советская энциклопедия, 1977. - Т. 27. - С. 571-572,
597-600.
BY 16426 C1 2012.10.30
(57)
Способ определения дальности до объекта, его высоты и ширины, заключающийся в
том, что получают два цифровых фотографических изображения измеряемого объекта из
двух точек, разнесенных в пространстве по горизонтали на известное расстояние, на полученных изображениях формируют по окну сканирования с одинаковыми координатами,
размер которых выбирают таким образом, чтобы разность расстояний до отдельных фрагментов измеряемого объекта была меньше инструментального разрешения по дальности,
осуществляют сканирование одного окна относительно другого с шагом в один пиксель
фотоприемной матрицы используемой фотокамеры по горизонтали и вертикали, определяют точку расположения максимума корреляционной функции, в области которой строят
BY 16426 C1 2012.10.30
сетку ячеек, соответствующих узлам интерполяции, и осуществляют повторное сканирование с шагом h = 1 ⁄ k пикселя, где k - целое число от 5 до 20, причем область выбирают
шириной и высотой не менее двух пикселей, осуществляют уточнение расположения максимума корреляционной функции в субпиксельном диапазоне в соответствии с выражением:
I(x + ih, y + jh ) = (1 − ih )(1 − jh ) ⋅ I(x, y ) + (1 − jh )ih ⋅ I(x + 1, y ) +
+ (1 − ih ) jh ⋅ I(x, y + 1) + ijh 2 ⋅ I(x + 1, y + 1),
где I(x + ih, y + jh) - интенсивность узла сетки;
x, y - значения координат точки расположения максимума корреляционной функции
по горизонтали и вертикали;
i, j - горизонтальный и вертикальный индексы узла;
I(x, y) - значения интенсивности в точке расположения максимума корреляционной
функции;
I(x + 1, y), I(x, y + 1), I(x + 1, y + 1) - значения интенсивности корреляционной функции в ближайших точках к точке расположения максимума корреляционной функции;
затем осуществляют локализацию максимума между узлом сетки с наибольшим значением корреляционной функции и его соседними узлами в соответствии с выражением:
1 x 2m +1 − x 2m (I m − I m −1 ) + x 2m − x 2m −1 (I m − I m +1 )
∆x =
,
2 (x m +1 − x m )(I m − I m −1 ) + (x m − x m −1 )(I m − I m +1 )
где ∆x - сдвиг между изображениями объекта по горизонтали;
xm + 1, xm-1, Im + 1, Im-1 - координаты и интенсивности узлов сетки, соседствующих с узлом сетки с максимальным значением интенсивности;
xm, Im, - координата и интенсивность узла сетки с максимальным значением интенсивности;
а дальность R до объекта определяют из выражения:
f ⋅L
,
R=
∆x
где f - фокусное расстояние используемой фотокамеры;
L - расстояние между точками фотографирования;
а высоту H и ширину D объекта определяют из выражений:
R ⋅ y′
H=
,
f
R ⋅ x′
D=
,
f
где y′, x′ - размеры объекта на фотоприемной матрице по вертикали и горизонтали.
(
)
(
)
Изобретение относится к области информационно-измерительных систем и предназначается для решения задач измерения дальности и линейных размеров объектов по их
цифровым изображениям.
Известен способ измерения дальности и линейных размеров объектов, основанный на
формировании стереопары цифровых изображений измеряемого объекта [1], при этом
дальность до объекта определяется по сдвигу между изображениями. Недостатком такой
системы является невысокая точность измерений.
Наиболее близким к предлагаемому изобретению является способ измерения расстояний до объекта [2] с использованием двух цифровых фотокамер, разнесенных по горизонтали на известное расстояние. Дальность до объекта определяется по сдвигу между
изображениями, при этом необходимо знать расстояние между фотокамерами и фокусы
объективов камер. Недостатком этого способа является низкая точность измерения даль2
BY 16426 C1 2012.10.30
ности, обусловленная тем, что сдвиг между изображениями на фотоприемной матрице
определяется только с точностью до одного пикселя.
Задача изобретения - повышение точности измерений расстояний. Решение поставленной задачи позволит использовать изобретение для измерения расстояний и размеров
объектов как для решения задач криминалистики, так и для решения задач геодезии, строительства и т.п.
Поставленная задача решается в способе измерения расстояний на цифровой фотокамере, заключающемся в том, что получают два цифровых фотографических изображения
измеряемого объекта из двух точек, разнесенных в пространстве по горизонтали на известное расстояние, на полученных изображениях формируют по окну сканирования с
одинаковыми координатами, размер которых выбирают таким образом, чтобы разность
расстояний до отдельных фрагментов измеряемого объекта была меньше инструментального разрешения по дальности, осуществляют сканирование одного окна относительно
другого с шагом в один пиксель фотоприемной матрицы используемой фотокамеры по
горизонтали и вертикали, определяют точку расположения максимума корреляционной
функции, в области которой строят сетку ячеек, соответствующих узлам интерполяции, и
осуществляют повторное сканирование с шагом, равным h = 1 / k пикселя, где k - целое
число от 5 до 20, причем область выбирают шириной и высотой не менее двух пикселей,
осуществляют уточнение положения максимума корреляционной функции в субпиксельном диапазоне в соответствии с выражением:
I(x + ih, y + jh ) = (1 − ih )(1 − jh ) ⋅ I(x, y ) + (1 − jh )ih ⋅ I(x + 1, y ) +
+ (1 − ih ) jh ⋅ I(x, y + 1) + ijh 2 ⋅ I(x + 1, y + 1),
где I(x + ih, y + jh) - интенсивность узла сетки; x, y - координаты максимума корреляционной функции по горизонтали и вертикали; i, j - горизонтальный и вертикальный индексы
узла; I(x, y) - значение интенсивности в точке расположения максимума корреляционной
функции; I(x + 1, y), I(x, y + 1), I(x + 1, y + 1) - значения интенсивности корреляционной
функции в ближайших точках к точке расположения максимума корреляционной функции;
затем осуществляют локализацию максимума между узлом сетки с наибольшим значением корреляционной функции и его соседними узлами в соответствии с выражением:
1 x 2m +1 − x 2m (I m − I m −1 ) + x 2m − x 2m −1 (I m − Im +1 )
∆x =
,
2 (x m +1 − x m )(I m − I m −1 ) + (x m − x m −1 )(I m − I m +1 )
где ∆x - сдвиг между изображениями объекта по горизонтали; xm + 1, xm-1, Im + 1, Im-1 - координаты и интенсивности узлов сетки, соседствующих с узлом сетки с максимальным значением интенсивности; xm, Im - координата и интенсивность узла сетки с максимальным
значением интенсивности;
а дальность R до объекта определяют из выражения:
f ⋅L
R=
,
∆x
где f - фокусное расстояние используемой фотокамеры; L - расстояние между точками фотографирования;
а высоту H и ширину D объекта определяют из выражений:
R ⋅ y′
R ⋅ x′
,
,
H=
D=
f
f
где x', y' - размеры объекта на фотоприемной матрице по горизонтали и вертикали.
Свойство, появляющееся у заявляемого объекта, - это повышение точности измерения
расстояний, обусловленное тем, что в прототипе сдвиг между изображениями определяется с точностью в один пиксель, а в предлагаемом изобретении разрешение при определении сдвига составляет сотые доли пикселя, что повышает точность измерения расстояний.
(
)
(
3
)
BY 16426 C1 2012.10.30
Сущность способа измерений поясняется с помощью фигуры, на которой представлена функциональная схема измерителя расстояний на основе цифрового фотоаппарата. Система содержит измеряемый объект 1, два цифровых фотоаппарата, состоящих из линз 2 и
фотоприемных матриц 3, 4 соответственно, вычислительный блок 5.
Измеритель расстояний работает следующим образом. С помощью цифровых фотоаппаратов на фотоприемных матрицах 3 и 4 реализуются цифровые изображения измеряемых объектов. Полученные два цифровых изображения одних и тех же объектов
поступают в вычислительный блок, где производится измерение расстояний по следующему алгоритму. На первом изображении формируется окно сканирования, размер которого выбирают таким образом, чтобы разность расстояний до отдельных фрагментов
выделенного объекта была меньше инструментального разрешения по дальности. Если
поверхность объекта является плоскостью, перпендикулярной горизонтали и оси наблюдения, то все точки плоскости будут находиться на одинаковом расстоянии от объектива.
Если предмет имеет объемную форму, то необходимо на нем выделить плоскую поверхность. Автоматически окно сканирования с такими же координатами формируется и на
втором снимке. Затем осуществляется сканирование одного окна относительно другого по
горизонтали и вертикали. Каждую итерацию окно сдвигается на пиксель вдоль осей в системе координат изображения и вычисляется соответствующее значение корреляционной
функции. По положению максимального значения нормированной корреляционной функции определяют сдвиг между изображениями ∆x = x2-x1. Так как все точки объекта в выделенном окне находятся на одинаковом расстоянии и осуществляется нормировка по
величине среднего значения сигнала, то при полном совпадении изображений корреляционная функция будет близка к единице. Дальность R до выделенной области объекта
определяется из выражения:
f ⋅L
R=
,
(1)
∆x
где L - расстояние между фотокамерами, f - фокусное расстояние фотокамер.
Аналогичным образом определяются расстояния Ri до всех объектов, попавших в поле
зрения фотокамеры.
Система также позволяет измерять высоту и линейные размеры объектов. Это осуществляется следующим образом. Определив значение расстояния до объекта R и размеры этого объекта на фотоприемной матрице x' по горизонтали и y' по вертикали, ширину
объекта Di и высоту Hi определяют из выражений:
R ⋅ y′
R ⋅ x′
H=
, D=
,
(2)
f
f
Вычислительный блок с помощью геометрии полученных изображений объектов может также осуществить измерение расстояний между объектами на цифровой фотографии.
Так как сканирование одного окна относительно другого осуществляется со сдвигом
на один пиксель вдоль осей координат, то положение максимума корреляционной функции и, следовательно, сдвиг между изображениями определяются также с точностью в
один пиксель. Этот фактор значительно ограничивает точность измерения расстояний.
Для повышения точности необходимо определять положение максимума корреляционной
функции (сдвига ∆x) в пределах одного пикселя.
Для уточнения величины сдвига ∆x в субпиксельном диапазоне используется билинейная интерполяция сигнала с заданным шагом, алгоритм которой заключается в следующем. Обычно область уточнения выбирается шириной и высотой в два пикселя. В точке
максимального значения корреляционной функции I(x, y), полученного на этапе грубой
оценки с точностью в один пиксель, строится сетка ячеек, соответствующих узлам интерполяции, и ведется повторное сканирование с шагом, равным h = 1/k пикселя, где k - целое
число (k = 5…20). Таким образом, вычисляется значение корреляционной функции для k2
узлов сетки. Интенсивность I(x + ih, y + jh) узла сетки вычисляется из выражения:
4
BY 16426 C1 2012.10.30
I(x + ih, y + jh ) = (1 − ih )(1 − jh ) ⋅ I(x, y ) + (1 − jh )ih ⋅ I(x + 1, y ) +
(3)
+ (1 − ih ) jh ⋅ I(x, y + 1) + ijh 2 ⋅ I(x + 1, y + 1),
где h = 1 / k - шаг сетки; i, j - горизонтальный и вертикальный индексы узла соответственно; I(x, y) - максимальное значение корреляционной функции, I(x + 1, y), I(x, y + 1),
I(x + 1, y + 1) - значения интенсивности корреляционной функции в точках, ближайших к
максимальному.
Затем производится коррекция ошибки при помощи параболической интерполяции
между узлом сетки с наибольшим значением корреляционной функции и его соседними
узлами, локализирующая максимум параболы в заданном интервале. Допустим, найден
узел сетки с максимальным значением интенсивности Im и координатой xm. А его соседние
узлы с координатами xm + 1 и xm-1 имеют значения интенсивности Im + 1 и Im-1 причем
Im + 1 < Im и Im-1 < Im. Для уточнения величины смещения на заданном интервале проводится локализация максимума путем параболической интерполяции в соответствии с выражением:
1 x 2m +1 − x 2m (I m − I m −1 ) + x 2m − x 2m −1 (I m − I m +1 )
∆x =
.
(4)
2 (x m +1 − x m )(I m − I m −1 ) + (x m − x m −1 )(I m − I m +1 )
Очевидно, что билинейная интерполяция дает разрешение в пределах одного пикселя,
равное 1/k, а параболическая интерполяция имеет относительную погрешность 17 %. Следовательно, если взять k = 10…15, то разрешение при определении сдвига составит величину в сотую долю пикселя.
Таким образом, за счет использования предложенной методики субпиксельной интерполяции достигается повышение точности измерения расстояний в результате анализа
двух цифровых фотографических изображений, полученных с двух точек, разнесенных в
пространстве по горизонтали.
(
)
(
)
Источники информации:
1. Грузман И.С. и др. Цифровая обработка изображений в информационных системах. Новосибирск: НГТУ, 2000. - 168 с.
2. US 5432594, МПК G 01С 3/00, 1995.
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
5
Документ
Категория
Без категории
Просмотров
1
Размер файла
151 Кб
Теги
by16426, патент
1/--страниц
Пожаловаться на содержимое документа