5.4. ВЛИЯНИЕ ИЗМЕНЕНИЯ ЭНЕРГЕТИЧЕСКОГО СПЕКТРА РЕНТГЕНОВСКОГО ИЗЛУЧЕНИЯ ПРИ ЕГО ПРОХОЖДЕНИИ ЧЕРЕЗ ВЕЩЕСТВО НА КАЧЕСТВО ВОССТАНОВЛЕННЫХ ИЗОБРАЖЕНИЙ

Эффекты изменения спектрального состава рентгеновского излучения рассмотрим на примере полихроматических проекционных данных фантома головы, используя при этом выражение (4.6). Число дискретных уровней энергии равно пяти, а значения (вероятности того, что фотон, зарегистрированный при проведении калибровочных измерений, имеет энергию ) задаются табл. 5.3. В реальных условиях энергетический спектр рентгеновского излучения является сплошным; в табл. 5.3 приведены значения для дискретных уровней, аппроксимирующих сплошной спектр рентгеновского излучения, обычно используемый в реконструктивной томографии.

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

Данные о размерах и расположении элементарных объектов при использовании выражения (4.6) взяты из табл. 4.2, а величины (плотности элементарного объекта при энергии рассчитаны по данным табл. 4.1. Эти полихроматические проекционные данные вычислялись в приближении точечных источника и детекторов, в отсутствие рассеяния и без учета статистики фотонов.

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

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

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

(кликните для просмотра скана)

(кликните для просмотра скана)

Используя выражение (3.14) для рассматриваемого дискретного энергетического спектра и считая, что материалом «а» является воздух, а материалом «б» - метастаз карциномы груди, получаем следующее соотношение между полихроматической лучевой суммой и монохроматической лучевой суммой

(значения и приведены в табл. 5.3).

Полином

является «решением» уравнения (5.9) в том смысле, что если проанализировать спектр значений монохроматических проекционных данных для рассматриваемого фантома головы, то значения находятся между 0 и 4. Пусть

а

Другими словами, равно значению правой части выражения (5.9), если заменить на ту Если вычислить и для (т.е. для ), то получим

Это значит, что для того диапазона значений который нас интересует, среднеквадратичная погрешность из-за расхождения между фактическим значением и величиной, полученной из соотношений (5.10), равна всего 0,00006. С этой точностью (5.10) можно считать решением (5.9).

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

Отметим между прочим, что (5.11) подтверждает наше утверждение, приведенное в разд. 3.2, что использование в качестве корректирующей

функции полиномов низкого порядка, по всей вероятности, достаточно оправдано. Мы показали, что для объектов, состоящих из воздуха и ограниченного количества ткани, метастаза легкого, кубическая коррекция в выражении (3.15)] дает прекрасные результаты.

На рис. 5.9,б представлена реконструкция, полученная по полихроматическим проекционным данным, скорректированным при помощи (5.10), а на рис. 5.10,б - соответствующий график зависимости плотности вдоль 63-го столбца. То, что данная реконструкция хуже, чем изображение, которое получают по монохроматическим данным (рис. 5.3,б и 5.4), определяется тем, что приближение о наличии только двух видов ткани в области реконструкции является слишком грубым.

Чтобы дать меру того, насколько хороша поправка на изменение спектра рентгеновского излучения, основанная на использовании выражения (5.10), необходимо проанализировать значения лучевых сумм, полученных в рабочих измерениях для описанного выше фантома. В этом случае используется 288 положений источника и 165 положений детекторов, или полное число лучей 47 520. Пронумеруем эти лучи некоторым образом, и пусть обозначают монохроматические и полихроматические лучевые суммы соответственно для луча

Количественная мера качества коррекции определяется тем фактом, что

Следует отметить, что правая часть выражения (5.14) значительно больше правой части выражения (5.13). Теперь необходимо найти такой полином что если его использовать вместо в выражении (4.6), то правая часть этого выражения будет такой же малой, как и при использовании любой другой полиномной функции вместо Находим, что полином

удовлетворяет этому условию. Действительно,

превосходит значение, указанное в (5.14).

Интересно отметить, что данный полином в (5.15) является полиномом первой степени. Можно показать, что для указанного набора данных использование полиномов более высокой степени не приводит к существенному улучшению результата. Наилучшее, чего можно достичь при использовании полиномов пятой степени, это снизить значение правой части выражения (5.16) до 0,0381. Следует подчеркнуть, что это есть

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

На рис. 5.9,в представлена реконструкция, полученная полихроматическим данным, которые были представлены при помощи выражения (5.15), а соответствующий график зависимости плотности вдоль 63-го столбца приведен на рис. 5.10,е. В этой реконструкции числа Хаунсфилда очень точно воспроизводят реальные коэффициенты ослабления рентгеновского излучения в тканях, находящихся в центральной области. Однако при помощи линейных полиномов коррекции размазывание исключить не удается, так как при такой коррекции можно только «повысить» данные числа Хаунсфилда одновременно для всех элизов. При поверхностном сравнении нескорректированного изображения и изображения, скорректированного указанным методом при помощи полиномов, можно прийти к выводу о том, что предпочтение следует отдать первому изображению из-за его визуальных качеств. Такой парадокс связан с выбором того окна, которое используют при отображении рисунков.

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

Если получены оценки для моно- и полихроматических лучевых сумм, такие, что достаточно близко к где корректирующий полином, то можно построить следующее приближение к монохроматической лучевой сумме:

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

Пусть есть -мерный вектор компонента которого равна длине пересечения луча с элизом. Учитывая все это, определим псевдомонохроматическую лучевую сумму как

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

информацию, получить I функций для таких, что

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

Сравнивая выражения (5.18) и (2.4), а также (5.20) и (3.10), мы видим, что если дает хорошую аппроксимацию данной дискретизации изображения, которое надо восстановить, то (5.17), вероятно, будет справедливо. Правая часть выражения (5.17) дает новое значение монохроматических проекционных данных, которые затем можно использовать повторно в качестве входной информации в алгоритме реконструкции. На рис. 5.9,г представлена реконструкция, полученная по этим входным данным, а соответствующий ей график зависимости плотности вдоль 63-го столбца приведен на рис. 5.10,г. Как реконструкция, так и график показывают, насколько хорошо действует эта одноразовая итерация. Область в центре, которая уже была достаточно хорошо скорректирована, почти не затрагивается, тогда как ббльшая часть размазанной области скорректирована, а граница между черепом и мозговым веществом видна намного четче. Замечательным является тот факт, что на снимке снова видна гематома. Сравнение этой реконструкции с реконструкцией, полученной по монохроматическим проекционным данным, относительно которой мы испытываем методы коррекции, достаточно благоприятно.

Однако с улучшением качества реконструкции растет стоимость этой реконструкции, а описанный только что метод обходится достаточно дорого. Если коррекцию ведут при помощи полиномов, то требуется выполнить несколько операций умножения на лучевую сумму даже для полиномов высокого порядка. Между тем при вычислении псевдосумм приходится выполнять большое число операций на каждый луч, определять длину пересечения луча и элиза, а также возводить в степень. Время, затраченное на проведение коррекции при помощи полиномов, с одной стороны, и получения проекционных данных для одного шага итерирования — с другой, равно 50 и 1000 с соответственно. Несмотря на то что в этих программах не используют специальное для реконструкций математическое обеспечение и поэтому можно ожидать более коротких коммерческих времен, указанные значения времени дают количественное представление о разнице в стоимости этих двух методов коррекции.

Если считать, что заново полученная реконструкция (рис. 5.9,г) лучше реконструкции рис. 5.9,в, то следует ожидать, что повторное вычисление величин по этой лучшей реконструкции при помощи (5.17)

Таблица 5.4 (см. скан) Мера различия изображений, полученных по монохроматическим, полихроматическим и скорректированным полихроматическим проекционным данным

приведет к еще более высокому качеству реконструкции. Подобная реконструкция приведена на рис. 5.9,д, а соответствующая кривая зависимости плотности по 63-му столбцу — на рис. 5.10,д. По сравнению с предыдущей реконструкцией имеется лишь незначительное улучшение качества изображения на этом рисунке. С точки зрения стоимости дальнейшие итерации нецелесообразны.

В табл. 5.4 приведены меры различия изображений для указанных выше реконструкций.