1.2. Наблюдение и измерение цифровых изображений

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

Рис.1.1

Координаты центров пикселов в левой прямоугольной системе координат цифрового изображения о_C x_C у_C .(рис.1.1), началом которой является левый верхний угол цифрового изображения, определяются в, так называемых, пиксельных координатах (единицей измерения в этом случае является пиксел).

Пиксельные координаты центров пикселов в системе координат цифрового изображения о_C х_C у_C определяют по формулам:

. (1.1)

Для измерения координат точек цифрового изображения его визуализируют на экране дисплея. Если пиксел изображения на

экране дисплея соответствует пикселу исходного цифрового

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

Для получения подпиксельной (субпиксельной) точности можно увеличить матрицу изображения на экране монитора относительно исходного цифрового изображения. В этом случае каждый пиксел исходного изображения будет изображаться матрицей n×n пикселов, численное значение всех элементов a'_ij которой будут равны численному значению элемента матрицы исходного изображения.

Пиксельные координаты точек увеличенного изображения можно измерить с точностью до 1/n пиксела исходного изображения (рис.1.2.).

Рис. 1.2

Пиксельные координаты (в пикселах исходного изображения) элемента a'_ij увеличенного изображения определяют по формуле:

, (1.2)

в которых: i,j - номера строки и столбца элемента матрицы исходного изображения, в котором находится элемент a'_ij

увеличенного изображения:

i’,j’ - номера строки и столбца элемента a`_ij подматрицы n×n;

n – коэффициент увеличения изображения.

Например, для элемента a’₃₃ (рис.1.2) пиксельные координаты:

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

Значения физических координат определяют по формулам:

. (1.3)

Например, координаты центра пиксела, соответствующего элементу a’₃₃ (рис.1.2) при величине Δ=20 мкм будут равны х_C = 34 мкм и y_C = 34 мкм.

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

В этом случае значения пиксельных координат вычисляют по формулам:

, (1.4)

при измерениях с точностью до пиксела и по формулам:

, (1.5)

при измерениях с подпиксельной точностью.

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

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

Принцип измерения координат точек цифрового изображения по этому методу иллюстрируется на рис. 1.3 и 1.4.

Рис. 1.3 Рис. 1.4

На рисунке 1.3 представлен фрагмент исходного цифрового изображения с измерительной маркой и точкой изображения m, координаты которой необходимо измерить. Как следует из рис.1.3 центр изображения измерительной марки не совпадает с изображением точки m, причем разности значений их пиксельных координат составляют величины x _P и y _P.

Для совмещения центра изображения измерительной марки с точкой m можно создать фрагмент цифрового изображения снимка, в котором координаты начала системы координат o’_C x’_C y’_C будут иметь значения , а.

Создание такого фрагмента цифрового изображения производится следующим образом. По координатам центра каждого пиксела фрагмента изображения x’_pi,y’_piопределяют значения координат его проекцииx_pi,y_pi в системе координат о_Cх_Cу_Cисходного изображения.

Их значения определяют по формулам:

. (1.6)

Затем по значениям координат x_pi,y_piнаходят ближайшие к изображению точкиi, соответствующей центру пиксела

Рис. 1.5

создаваемого фрагмента цифрового изображения, четыре пиксела исходного цифрового изображения, например, M,K,L,N(рис.1.5)

Далее методом билинейного интерполирования определяют значения оптической плотности i-го пиксела создаваемого фрагмента изображения по формуле:

, (1.7)

в которой

.

Таким же образом формируются все элементы создаваемого фрагмента цифрового изображения.

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

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

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