Обработка рентгеновских изображений в среде MATLAB | Статья в журнале «Молодой ученый»

Авторы: ,

Рубрика: Технические науки

Опубликовано в Молодой учёный №8 (112) апрель-2 2016 г.

Дата публикации: 07.04.2016

Статья просмотрена: 346 раз

Библиографическое описание:

Куликов А. Ю., Романова В. О. Обработка рентгеновских изображений в среде MATLAB // Молодой ученый. — 2016. — №8. — С. 251-253. — URL https://moluch.ru/archive/112/28261/ (дата обращения: 23.07.2018).



Описан способ обработки рентгеновских изображений в среде MATLAB с результатами обработки реальных рентгенограмм.

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

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

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

Подобную обработку удобно проводить в среде MATLAB, которая является языком высокого уровня для выполнения технических и научных вычислений [3]. Проводятся эти преобразования согласно выражению

(1)

где L, Lres — массив значений соответственно исходного и результирующего изображений;

Lmin, Lmax — соответственно минимальное и максимальное значения исходного изображения.

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

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

,(2)

где α — коэффициент нелинейности.

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

for k=1:3;

Lmin=min(min(L(:,:,k)));Lmax=max(max(L(:,:,k)));

for i=1:640; for j=1:640;

Lvyh(i,j,k)=(((L(i,j,k)-Lmin)/(Lmax-Lmin)))^.5;

if Lvyh(i,j,k)>1;

Lvyh(i,j,k)=1;

end;

if Lvyh(i,j,k)<0;

Lvyh(i,j,k)=0;

end;

end;

end;

end;

figure, imshow(Lvyh);

На рис. 1 изображены исходные рентгеновские изображения с соответствующими гистограммами, а на рисунке 2 изображены обработанные методом растяжения динамического диапазона изображения с соответствующими гистограммами.

Рис. 1. Исходные рентгеновские изображения и их гистрограммы

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

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

Рис. 2. Обработанные методом растяжения динамического диапазона рентгеновские изображения

Результаты экспертных оценок выражены с помощью среднего балла, определяемого как:

,

где nk — число изображений, отнесенных к k-й категории,

Ck — соответствующий ей балл.

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

Литература:

  1. Смирнов А. Я. Критерии качества дискретизированных изображений // Труды ГОИ им. С. И. Вавилова. — т. 57. — вып. 191. — Л. — 1984.
  2. Гонсалес Р., Вудс Р. Цифровая обработка изображений — М.: Техносфера, 2005. — 1072 с.
  3. Гонсалес Р., Вудс Р., Эддинс С. Цифровая обработка изображений в среде МATLAB — М.: Техносфера, 2006. — 616 с.
Основные термины (генерируются автоматически): MATLAB, исходное изображение, интенсивность пикселей, рентгеновское изображение, узкий динамический диапазон.


Похожие статьи

Обработка рентгеновских изображений с применением оценки...

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

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

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

где — интенсивность пикселя с координатами (y, x) в исходном изображении.

Обработка рентгеновских изображений в среде MATLAB.

Реализация частотной фильтрации рентгеновских изображений...

Рис. 2. Рентгеновское изображение: а) — исходное, б) после фильтрации с усилением высоких частот.

Обработка рентгеновских изображений в среде MATLAB.

Особенности применения фильтров обработки изображений...

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

Обработка паспортных данных в среде Matlab | Статья в журнале...

Пусть исходное цветное изображение имеет разрешение 500х1000 пикселей.

Обработка рентгеновских изображений в среде MATLAB. Разработка программных средств синтеза и анализа весовых функций в среде MATLAB.

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

Если значение интенсивности пикселя меньше порога, то мы его относим к области объекта, в противном случае — к области фона

На рисунке 9 показано исходное изображение, промежуточные результаты и

Обработка рентгеновских изображений в среде MATLAB.

Сегментация микроскопических изображений эпителиальных клеток

[1] Алгоритм также ищет на исходном изображении «водосборные бассейны» и «линию водораздела», где белые пиксели расположены выше, а темные — ниже.

В Matlab сегментация цветных изображений методом K-средних осуществляется следующим образом

Разработка двумерных сглаживающих фильтров на основе...

Эти элементы принято называть элементами изображения или пикселями. Для цифровых изображений в оттенках серого интенсивность представляется с помощью чисел без знака, изменяющимся в диапазоне от 0 до 255, с 8-битовым

Рис. 2. Исходное изображение.

Исследование способов увеличения резкости рентгеновских...

Для получения цифровых изображений используются фоточувствительные элементы.

Производная первого порядка одномерной функции f(x) определяется как разность яркостей соседних пикселей

К вопросу о динамической устойчивости маятниковых систем.

Обсуждение

Социальные комментарии Cackle

Похожие статьи

Обработка рентгеновских изображений с применением оценки...

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

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

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

где — интенсивность пикселя с координатами (y, x) в исходном изображении.

Обработка рентгеновских изображений в среде MATLAB.

Реализация частотной фильтрации рентгеновских изображений...

Рис. 2. Рентгеновское изображение: а) — исходное, б) после фильтрации с усилением высоких частот.

Обработка рентгеновских изображений в среде MATLAB.

Особенности применения фильтров обработки изображений...

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

Обработка паспортных данных в среде Matlab | Статья в журнале...

Пусть исходное цветное изображение имеет разрешение 500х1000 пикселей.

Обработка рентгеновских изображений в среде MATLAB. Разработка программных средств синтеза и анализа весовых функций в среде MATLAB.

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

Если значение интенсивности пикселя меньше порога, то мы его относим к области объекта, в противном случае — к области фона

На рисунке 9 показано исходное изображение, промежуточные результаты и

Обработка рентгеновских изображений в среде MATLAB.

Сегментация микроскопических изображений эпителиальных клеток

[1] Алгоритм также ищет на исходном изображении «водосборные бассейны» и «линию водораздела», где белые пиксели расположены выше, а темные — ниже.

В Matlab сегментация цветных изображений методом K-средних осуществляется следующим образом

Разработка двумерных сглаживающих фильтров на основе...

Эти элементы принято называть элементами изображения или пикселями. Для цифровых изображений в оттенках серого интенсивность представляется с помощью чисел без знака, изменяющимся в диапазоне от 0 до 255, с 8-битовым

Рис. 2. Исходное изображение.

Исследование способов увеличения резкости рентгеновских...

Для получения цифровых изображений используются фоточувствительные элементы.

Производная первого порядка одномерной функции f(x) определяется как разность яркостей соседних пикселей

К вопросу о динамической устойчивости маятниковых систем.

Задать вопрос