Выбор технических средств для АРМ ВР (автоматизированное рабочее место врача-рентгенолога) во многом зависит от типа решаемых задач. Обычно в качестве технической базы для АРМ обработки изображений используют графические станции или персональные компьютеры. Графические станции, созданные прежде всего для решений задач машинной графики, оборудованы специальными графическими процессорами, ускоряющими процедуры построения графических примитивов (особенно трехмерных). Для задач обработки и анализа изображений более существенна скорость обработки видеоданных. Поэтому в качестве технической базы АРМ ВР использована широко распространенная и дешевая ПЭВМ типа IBM PC/AT.
[pic]
Рис.4 Блок-схема технических средств АРМ ВР.
1-негатоскоп; 2-телевизионная камера; 3-ПЭВМ; 4-фрейм- граббер; 5-телемонитор.
Практическая работа показала, что производительность персонального компьютера во многих случаях достаточна, чтобы решать задачи обработки видеоданных в реальном времени врача. Кроме того ПЭВМ имеют мощные технические и программные средства для организации «оконного» человеко- машинного диалога.
При использовании изображений, записанных в аналоговом виде, например рентгенограмм, необходимо устройство для ввода и визуализации их в ЭВМ. В качестве такого устройства удобно использовать фрейм-граббер конструктивно оформленный в виде платы, расположенной в корпусе ПЭВМ. Также необходимо иметь телекамеру с объективом, световой стол для подсветки рентгенограмм (негатоскоп) и телемонитор для визуализации изображений (рис.4). Устройство цифрового ввода и визуализации изображений должно обеспечивать высокое качество представления медицинских изображений, чтобы при их использовании не терялась важная диагностическая информация. [№ 7]
3.2. Области применения и преимущества цифровых систем.
К преимуществам цифровых рентгенографических систем относятся следующие четыре фактора: цифровое отображение изображения; пониженная доза облучения; цифровая обработка изображений; цифровое хранение и улучшение качества изображений.
Рассмотрим первое преимущество, связанное с отображением цифровой информации. Разложение изображения по уровням яркости на экране становится в полной мере доступным для пользователя. Весь диапазон оптических яркостей может быть использован для отображения лишь одного участка изображения, что приводит к повышению контраста в интересующей области. В распоряжении оператора имеются алгоритмы для аналоговой обработки изображения с целью оптимального использования возможностей систем отображения. [№ 5, стр. 103]
Это свойство цифровой рентгенографии также дает возможность снизить лучевую нагрузку на пациента путем уменьшения количества рентгенограмм для получения диагностической информации (той же полезности).
Цифровое отображение при его компьютерной обработке позволяет извлечь количественную и качественную информацию и таким образом перейти от интуитивно-эмпирического способа изображения к объективно измеренному.
Существенным преимуществам цифровой рентгенографии перед экранно- пленочным процессом являются простота и скорость получения изображения. Изображение становится доступным анализу врачом-рентгенологом в момент окончания экспозиции. [№ 6]
Второе преимущество цифровой рентгенологии - возможность снижения дозы облучения. Если в обычной рентгенологии доза облучения зависит от чувствительности приемника изображения и динамического диапазона пленки, то в цифровой рентгенологии оба этих показателя могут оказаться несущественными. Снижения дозы можно достичь установкой экспозиции, при которой поддерживается требуемый уровень шума в изображении. Дальнейшее уменьшение дозы возможно путем подбора такой длины волны рентгеновского излучения, которая обеспечивала бы минимальную дозу при данном отношении сигнал/шум, а также путем ликвидации любых потерь контраста с помощью описанных выше методов отображения цифровых изображений.
Третье преимущество цифровой рентгенологии - это возможность цифровой обработки изображений. Рентгенолог должен выявить аномальные образования на осложненной фоном нормальной структуре биоткани. Он может не заметить мелких деталей в изображении, которые система разрешает, или пропустить слабоконтрастную структуру, видимую на фоне шумов изображения, из-за сложного строения окружающих (или сверхлежащих) тканей. Субстракционный метод в рентгенографии позволяет устранить большую часть паразитной фоновой структуры и тем самым увеличить вероятность выявления важных деталей на рентгенограмме. Компьютерную томографию можно рассматривать как частный случай метода субстракционной рентгенографии, в котором из обычных проекционных изображений устраняется информация о вышележащих структурах. [№ 5, стр.103-104]
Особенная ценность применения цифровой рентгенографии заключается в возможности полного отказа от рентгеновской пленки и связанного с ней фотохимического процесса. Это делает рентгенологическое исследование экологически чище, а хранение информации в цифровом виде позволяет создать легкодоступные рентгеновские архивы. Новые количественные формы обработки информации открывают широкие возможности стандартизации получения изображений, приведения их к стандарту качества в момент получения и при отсроченных повторных исследованиях. Немаловажна открывающаяся возможность передачи изображения на любые расстояния при помощи средств компьютерных коммуникаций.
Приведенные соображения с достаточной наглядностью демонстрируют прогрессивность внедрения в практику цифровой рентгенографии, которая сможет перевести диагностическую рентгенологию на новый более высокий технологический уровень. Отказ от дорогостоящих расходных материалов обнаруживает и ее высокую экономическую эффективность, что в сочетании с возможностью уменьшения лучевых нагрузок на пациентов делает ее применение в практике особенно привлекательным. [№ 6]
3.3. Цифровая рентгенография с экрана электронно-оптического преобразователя (ЭОП).
Система рентгенографии с экрана ЭОП (рис. 5) состоит, как и обычная система электронно-оптического преобразования для просвечивания, из ЭОП, телевизионного тракта с высоким разрешением, рентгеновского высоковольтного генератора и рентгеновского излучателя. Сюда же входит штатив для исследования, цифровой преобразователь изображения и другие компоненты.
При обычной методике рентгенографии с экрана ЭОП с помощью 100 мм фотокамеры или кинокамеры переснимается оптическое изображение на выходном экране преобразователя.
В цифровой же системе сигнал, поступающий с видеокамеры, аналого- цифровым преобразователем трансформируется в набор цифровых данных и передается в накопительное устройство. Затем эти данные, в соответствии с выбранными исследователем параметрами, компьютерное устройство переводит в видимое изображение. [pic] Рис.5 Цифровая рентгенография с экрана ЭОП
1-генератор; 2-рентгеновская трубка; 3-пациент; 4-ЭОП; 5- видеокамера;
6-аналого-цифровой преобразователь; 7-накопитель изображений;
8-видеопроцессор; 9-сеть; 10-цифро-аналоговый преобразователь;
11-монитор; 12-снимок; 13-рентгенолог.
3.4. Цифровая люминесцентная рентгенография (ЦЛР).
Применяемые в ЦЛР (рис.6) пластины-приемники изображения после их экспонирования рентгеновским излучением последовательно, точка за точкой, сканируются специальным лазерным устройством, а возникающий в процессе лазерного сканирования световой пучок трансформируется в цифровой сигнал. После цифрового усиления контуров и контрастности элементов изображения оно лазерным принтером печатается на пленке или воспроизводится на телевизионном мониторе рабочей консоли. Люминесцентные пластины-накопители выпускаются в стандартных формах рентгеновской пленки, помещаются вместо обычных комплектов «пленка-усиливающий экран» в кассету и применяются в обычных рентгеновских аппаратах.
Такая пластина обладает значительно большей экспозиционной широтой, чем общепринятые комбинации пленка-экран, благодаря чему значительно расширяется интервал между недо- и переэкспонированием. Этим способом можно получать достаточно контрастные изображения даже при резко сниженной экспозиционной дозе, нижним пределом которой является лишь уровень квантового шума. Поэтому даже при рентгенографии в палате у постели больного методика ЦЛР гарантирует получения качественного снимка.
При ЦЛР используются цифровые преобразователи, пространственное разрешение которых выше, чем у большинства используемых в настоящее время для обычной рентгенографии комбинаций экран-пленка. Все же особым преимуществом ЦЛР является передача малоконтрастных деталей, тогда как передача очень мелких деталей, таких, например, как микрокальценаты в молочной железе, остается прерогативой рентгенографии на рентгеновской пленке. [pic] Рис. 6 Цифровая люминесцентная рентгенография.
1-генератор; 2-рентгеновская трубка; 3-пациент; 4- запоминающая пластина; 5-транспортирующее устройство; 6-аналого- цифровой преобразователь; 7-накопитель изображений;8- видеопроцессор; 9-сеть;
10-цифро-аналоговый преобразователь; 11-монитор; 12-снимок;
13-рентгенолог.
3.5. Селеновая рентгенография.
Селеновые детекторы представляют собой новейшую систему цифровой рентгенографии (рис. 7). Основной частью такого устройства служит детектор в виде барабана, покрытого слоем аморфного селена. Селеновая рентгенография в настоящее время используется только в системах рентгенографии грудной клетки. Характерная для снимков грудной клетки высокая контрастность между легочными полями и областью средостения при цифровой обработке сглаживается, не уменьшая при этом контрастности деталей изображения. Другим преимуществом селенового детектора является высокий коэффициент отношения сигнал/шум. [pic]
Рис.5 Цифровая селеновая рентгенография.
1-генератор; 2-рентгеновская трубка; 3-пациент; 4-селеновый барабан;
5-сканирующие электроды+усилитель; 6-аналого-цифровой преобразо- ватель; 7-накопитель изображений; 8-видеопроцессор; 9- сеть;
4. Математические основы компьютерной томографии
Исследования внутренней структуры объектов с помощью рентгеновского излучения широко распространены и хорошо известны. Ослабление рентгеновского излучения вдоль луча, соединяющего источник и приемник, является интегральной характеристикой плотности исследуемого объекта. С математической точки зрения речь идет о задаче восстановления функции по ее интегральным значениям вдоль некоторого семейства лучей. Различные лучи соответствуют различным (относительно объекта) положениям источника и приемника излучения. Такая модель является простейшей, но во многих случаях хорошо отражает реальную ситуацию и подтверждается исследованиям реальных тестовых объектов. Плотность реальных объектов является функцией трех пространственных координат. Однако в классической компьютерной томографии трехмерный объект представляют в виде набора тонких срезов. Внутри каждого среза плотность считают функцией только двух переменных. При исследовании фиксированного среза систему источник-приемник устраивают таким образом, что регистрируются данные только по лучам, лежащим в тонком слое относительно центральной плоскости среза. Таким образом приходят к задаче восстановления функции двух переменных по ее интегральным значениям вдоль некоторого семейства лучей Для регистрации в веерной схеме, чаще встречающейся в реальных томографах, используется линейка детекторов, различные положения источника относительно объекта обеспечиваются вращением системы регистрации или объекта.
4.1. Математическая постановка задачи рентгеновской компьютерной томографии, преобразование Радона и формулы обращения. В компьютерной рентгеновской томографии трехмерный объект представляется обычно в виде набора тонких срезов. Для восстановления плотности среза решается задача обращения двумерного преобразования Радона. Преобразованием Радона функции f(x, y) называется функция, [pic]определяемая равенством [pic]. Обычно для восстановления функции двух переменных по ее интегралам вдоль прямых используется метод свертки и обратного проецирования. В этом методе формула обращения преобразования Радона записывается без явного использования обобщенных функций. Однако наиболее общий и естественный вид формулы обращения преобразования Радона приобретают при использовании аппарата обобщенных функций. Далее будет рассмотрено соотношение между методом обобщенных функций и методом свертки и обратного проецирования. Перед изложением собственно численного алгоритма будет дан вывод формулы обращения, позволяющий естественным образом перейти к построению алгоритма. В силу равенства [pic] функция [pic]при любом фиксированном p определяется своими значениями при [pic]. Это позволяет нам перейти к функции [pic]. Здесь L(r, ?) - прямая, ортогональная лучу, имеющему угол ? ? положительным направлением оси X, и отстоящая от начала координат на расстояние r (r[pic] 0), при r < 0 L(r, ?) - прямая, симметричная относительно начала координат прямой L(|r|, ?). Выразим f(x, y) через I(r, ?). Поскольку [pic], где [pic]- преобразование Фурье функции f, то, переходя к полярным координатам после элементарных преобразований интеграла по ? на интервале [?, 2?], ?олучаем [pic]. Введем функцию S(z, ?), полагая [pic]. При фиксированном ? функция S(z, ?) ?сть обратное одномерное преобразование Фурье от произведения [pic]и |r|. Для [pic]справедливо равенство [pic]. Обратное преобразование Фурье от |r| есть обобщенная функция v1/?z2. Переходя от преобразования Фурье произведения к свертке, получаем S(z,?) = I(z,?)[pic](v1/?z2). Используя регуляризацию функции 1/z2 [19] приходим к выражению [pic]. (1.5.1) Таким образом, для f(x, y) справедлива формула [pic], (1.5.2) позволяющая выразить искомую функцию через наблюдаемые данные. Прежде чем перейти к дискретному варианту сделаем ряд замечаний, связанных с обоснованием корректности рассматриваемых алгоритмов в реальных ситуациях. Обобщенные функции являются функционалами над пространством бесконечно дифференцируемых быстро убывающих функций. Однако при построении аппроксимаций исходных реальных данных по отсчетам, заданным в дискретных точках, желательно иметь менее жесткие требования к гладкости аппроксимирующих функций. Свертка с обобщенными функциями, в частности, с функцией 1/z2, может быть определена для значительно менее гладких функций, это очень важно при доказательстве корректности применения численных алгоритмов, получаемых с помощью аппарата обобщенных функций, к реальным данным. Перейдем к дискретному варианту. Будем предполагать, что f(x, y) = 0 вне круга радиуса R с центром в нуле. Исходными данными являются величины I(ri, ?i), здесь ri v отсчеты в интервале [-R, R], 1 ? i ? M - отсчеты в интервал [0, ?], 1 ? j ? N. Если теперь при заданных значениях функции I(r, ?) ? отсчетах (ri, ?i) построить аппроксимацию I(r, ?) так, что для S(z,?) ?ыполняется равенство (1.5.1), то используя (1.5.1) и (1.5.2) можно получить приближение к f(x, y). В дальнейшем будем предполагать, что отсчеты на осях r и ? являются равноотстоящими. При каждом фиксированном ?j определим [pic]следующим образом.
Страницы: 1, 2, 3, 4, 5, 6
