ВИЗУАЛИЗАЦИЯ С ПОМОЩЬЮ ДИНАМИЧЕСКОЙ КОМПЬЮТЕРНОЙ ТОМОГРАФИИ
VISUALIZATION USING DYNAMIC COMPUTED TOMOGRAPHY

Настоящее изобретение относится к визуализации с помощью компьютерной томографии (КТ). Она находит применение в области кардиологии и в других областях, где существует необходимость в улучшении визуализации движущихся структур.
Недавние разработки в области компьютерной томографии (КТ), такие как увеличение скорости вращения гентри, и широкое распространение сканнеров с технологией множественного среза открыли новые возможности по применению визуализации с помощью компьютерной томографии. Одна из таких областей применения относится к визуализации сердечной деятельности, в частности к визуализации коронарных артерий.
Одним из наиболее сложных аспектов визуализации сердечной деятельности является обеспечение волюметрических данных с временным разрешением, необходимым для эффективной визуализации бьющегося сердца. См., например, Blobel, et al., Optimization of Temporal and Spatial Resolution for Cardiac CT Diagnostics (2004); Shoepf, et al., CT of Coronary Artery Disease, Radiology 2004: 232:18:32 (2004).
Одной из технологий получения синхронизированной (по ЭКГ) реконструкции сердечной деятельности является способ расширенной реконструкции сердечной деятельности (ERC-extended cardiac reconstruction). См. Grass, et al., Helical Cardiac Cone Beam Reconstruction Using Retrospective ECG Gating, Phys. Med. Biol. 48 (2003) 3069-3084. Способ расширенной реконструкции сердечной деятельности представляет собой способ приближенной реконструкции со спиральным сканированием конусным пучком, основанный на отфильтрованной 3D обратной проекции. Здесь, ретроспективная схема кардиосинхронизации ограничивает временную информацию конкретным, представляющим интерес, состоянием сердечной деятельности. Избыточность данных достигается путем использования режима реконструкции с более мелким шагом спирали.
Более подробно, технология расширенной реконструкции сердечной деятельности (ECR) основана на клиновом методе реконструкции со спиральным сканированием конусным пучком, который использует пересчет распределения дозы в каждом бине даваемой веерообразным пучком к таковой для параллельного пучка перед взвешиванием с использованием косинуса в качестве весовой функции, фильтрованием и 3D обратной проекции данных, полученных с использованием конусного пучка. Дополнительно к геометрическим взвешивающим коэффициентам, свойственным геометрии реконструкции, ECR способ также включает в себя взвешивающую функцию, основанную на окне освещения, которая является результатом системы наблюдения. Способ также включает в себя взвешивающую функцию сердечной деятельности, используемую для определения части проекционных данных во временном домене, который используется для реконструкции изображения объемного изображения для требуемой фазы сердечной деятельности. Перед обратной проекцией взвешивающую функцию освещения и взвешивающую функцию сердечной деятельности объединяют для каждого вокселя, используя метод нормализации.
В технологии эффективной реконструкции сердечной деятельности (Effective Cardiac Reconstruction) для визуализации каждой фазы сердечной деятельности ECR способу требуется новая обратная проекция. Поскольку этап обратной проекции требует больших денежных затрат, ECR способ не является подходящим для использования в тех областях, где требуется относительно быстрая, интерактивная визуализация волюметрических данных множества фаз сердечной деятельности.
В другом способе, поток получаемых данных разделяют на сегменты данных, имеющие угловые размеры, которые меньше размера, требуемого для формирования полного томографического набора данных. В свою очередь, сегменты данных реконструируют для формирования потока сегментных изображений. Затем производится временная повторная сортировка данных по глубине, которые затем добавляют для формирования волюметрических данных. См. Bruder, et al., Dynamic Cardiac CT Imaging Using Detectors with Large Cone Angle, Fully 3D Image Reconstruction in Radiology and Nuclear Medicine (2005).
Несмотря на то, что в этом подходе не нужно производить перерасчет объема изображения для каждой выбранной пользователем фазы сердечной деятельности, фаза сердечной деятельности может быть скорректирована только в пределах дискретной временной сетки, которая представляет собой функцию углового размера сегментов изображения. Более того, добавляя сегменты изображения для получения полного волюметрического изображения сердца для заново выбранной фазы, требуется запись большого количества данных на жесткий диск компьютера. Следовательно, данная технология не обеспечивает возможность для быстрой, интерактивной визуализации волюметрических данных для случайно выбранной фазы сердечной деятельности.
Аспекты настоящего изобретения решают упомянутые выше и другие задачи.
Согласно первому аспекту настоящего изобретения, способ включает в себя использование проекционных данных, полученных во время томографического обследования объекта с периодическим движением для реконструкции множества уровней изображения путем объединения уровней изображения для генерации первых данных изображения, характерных для объекта на случайно выбранной первой фазе движения, и отображения воспринимаемого человеком изображения, характерного для первых данных изображения. Уровни изображения генерируются из сегментов проекционных данных, имеющих угловой размер, который меньше размера, требуемого для обеспечения полного набора томографических данных.
Согласно другому аспекту изобретения, устройство включает в себя средство использования проекционных данных, собранных во время томографического обследования объекта с периодическим движением для реконструкции множества уровней изображения; средство объединения уровней изображения для генерации первых данных изображения, характерных для объекта в случайно выбранной первой фазе движения; и средство отображения воспринимаемого человеком изображения, характерного для первых данных изображения. Уровни изображения генерируются из сегментов проекционных данных, имеющих угловой размер, который меньше размера, требуемого для обеспечения полного набора томографических данных.
Согласно другому аспекту изобретения, способ включает в себя использование проекционных данных, собранных во время томографического обследования объекта с периодическим движением для реконструкции множества уровней изображения; объединение уровней изображения согласно первой взвешивающей функции для генерации первых данных изображения, характерных для объекта в случайно выбранной первой фазе движения; отображение воспринимаемого человеком изображения, характерного для первых данных изображения; объединение уровней изображения согласно второй взвешивающей функции для генерации вторых данных изображения, характерных для объекта на второй фазе движения; и генерацию воспринимаемого человеком изображения, характерного для объекта на второй фазе движения. Уровни изображения генерируются из множества смещенных под углом сегментов проекционных данных. Уровни изображения также смещены в результате первого углового смещения, при этом первая и вторая фазы движения разделены угловым расстоянием, которое меньше, чем первое угловое смещение.
Согласно другому аспекту, читаемые компьютером носители данных содержат инструкции, которые при их исполнении компьютером приводят к тому, что компьютер реализует способ, включающий в себя этапы, на которых используют проекционные данные, собранные во время томографического обследования объекта с периодическим движением для генерации множества уровней изображения, генерируют первую функцию взвешивания относительно углового смещения, выполняют взвешивание уровней изображения таким образом, чтобы примененное к уровням изображения взвешивание аппроксимировало первую взвешивающую функцию, объединяют уровни изображения для генерации первых волюметрических данных, характерных для объекта, и генерируют воспринимаемое человеком изображение, характерное для первых волюметрических данных.
Другие преимущества и выгоды будут очевидны специалистами в данной области техники после ознакомления с прилагаемым описанием.
На Фиг. 1 изображен КТ сканер.
На Фиг. 2 изображены этапы интерактивного способа визуализации.
На Фиг. 3А, 3В и 3С изображены аспекты способа генерации уровней изображения.
На Фиг.4 изображены этапы генерации уровней изображения.
На Фиг. 5 изображены аспекты способа генерации взвешивающей функции для уровней изображения.
На Фиг. 6 изображены этапы генерации взвешивающей функции для уровней изображения.
Ссылаясь на Фиг.1, КТ сканер (10) включает в себя вращающуюся гентри 18, которая вращается вокруг оси Z. Гентри 18 поддерживает источник 12 рентгеновского излучения, такой как рентгеновская трубка, которая генерирует обычно конический пучок излучения. Гентри 18 также поддерживает чувствительный к рентгеновскому излучению детектор 20, который находится напротив угловой дуги на противоположной стороне области 14 осмотра. Детектор 20 предпочтительно представляет собой детектор технологии множественного среза, который включает в себя множество рядов или слоев элементов детектора, расположенных вдоль оси Z, и множество колонок элементов детектора, расположенных в поперечном направлении. Детектор 20 генерирует выходные сигналы, свидетельствующие об излучении, получаемом от множества лучей. Также можно использовать плоский или другой детектор 20, а также геометрию системы сбора данных четвертого поколения или другую геометрию системы сбора данных.
Блок 28 электрокардиограммы (ЭКГ) генерирует данные, показывающие фазу сердечной деятельности пациента, проходящего обследование.
Опора 16 для пациента, такая как кресло, удерживает пациента в области 14 обследования. Опора 16 для пациента предпочтительно изготовлена с возможностью перемещения по оси Z. Контроллер 28 координирует различные параметры сканирования, которые являются необходимыми для выполнения требуемого протокола сканирования, включая параметры источника 12 рентгеновского излучения, такие как напряжение и ток в трубке. Перемещение опоры 16 для пациента предпочтительно координируется вращением гентри так, чтобы генерировать в основном путь сканирования по спирали.
Средство 30 реконструкции реконструирует поток проекционных данных для генерации волюметрических данных, показывающих внутреннюю анатомию пациента. Как более подробно описано ниже, средство 30 реконструкции включает в себя средство 34 реконструкции уровня изображения, которое реконструирует множество сегментов или уровней изображения.
Процессор 42 интерактивного взвешивания обрабатывает волюметрические данные изображения, сгенерированные средством 30 реконструкции для отображения в воспринимаемом человеком виде.
Компьютер общего назначения служит в качестве консоли 44 оператора. Консоль 44 включает в себя устройство вывода воспринимаемой человеком информации, такое как монитор или дисплей, и устройство ввода информации, такое как клавиатура и мышь. Консоль 44 также включает в себя программную память и память для хранения данных, включая дисковую память, относительно большую по объему, но имеющую более низкую скорость, и сравнительно меньшую по объему, но более быстродействующую оперативную память (RAM). Резидентное программное обеспечение на консоли позволяет оператору управлять работой сканера 10 путем установки требуемого протокола сканирования, запуская и завершая сканирование, наблюдая и манипулируя иным образом волюметрическими данными изображения и взаимодействуя, иным образом, со сканером 10. Следует отметить, что процессор 42 интерактивного взвешивания и консоль оператора могут быть реализованы на том же компьютере общего назначения.
На Фиг. 2 изображены этапы способа визуализации, который можно выполнять при помощи сканера 10 и который хорошо подходит для интерактивной визуализации объектов с периодическим движением, таких как коронарные артерии сердца.
На этапе 202 выполняют КТ сканирование пациента. Данные электрокардиограммы (ЭКГ) пациента записывают вместе с проекционными данными.
На этапе 204 устройство 30 реконструкции реконструирует проекционные данные для генерации волюметрических данных, показывающих представляющий интерес объем или область пациента, например область, которая включает в себя сердце. Данные ЭКГ преимущественно используются для выполнения ретроспективной синхронизованной по ЭКГ реконструкции требуемой фазы сердечной деятельности пациента. Поскольку волюметрические данные будут использованы на следующем этапе для определения представляющей интерес подобласти, параметры реконструкции могут быть установлены для относительно быстрой генерации изображений, хотя и относительно более низкого качества. В качестве альтернативы, параметры реконструкции могут быть установлены для генерации диагностических данных с более высоким качеством. В одном из вариантов осуществления реконструкцию выполняют при помощи ECR способа, хотя можно использовать другие подходящие способы реконструкции.
На этапе 206 на консоли 44 оператора отображается еще одно воспринимаемое человеком изображение, характерное для волюметрических данных. В одном из вариантов осуществления область, которая включает в себя сердце, представлена в виде одного или нескольких 3D изображений.
На этапе 208 пользователь выбирает представляющую интерес подобласть. В примере с сердцем, пользователь обычно выбирает представляющую интерес подобласть, которая содержит одну или несколько коронарных артерий. Как более подробно описано ниже, максимальный размер выбираемой пользователем подобласти предпочтительно является ограниченным.
На этапе 210 устройство 34 реконструкции уровней изображения выполняет несинхронизованную реконструкцию проекционных данных во множество сегментов или уровней изображения. Более конкретно, проекционные данные, получаемые во время каждого поворота вдоль пути по спирали, разделены на множество проекционных сегментов, имеющих угловое смещение Δθ:
Уравнение 1
,
seg
seg
где N представляет собой число сегментов на единицу вращения. Следует отметить, что N также может подбираться пользователем. Как будет описано ниже, каждый сегмент подвергается свертке и обратной проекции в предоставляющей интерес области или подобласти для формирования сегмента или уровня изображения.
Необходимое число уровней изображения, NS, может быть оценено следующим образом.
Уравнение 2
,
где ST представляет собой время сканирования, а TS представляет собой время сдвига между проекционными сегментами, которые используются для генерации последовательных уровней изображения.
Размер подобласти и число уровней изображения предпочтительно устанавливают согласно следующему отношению.
Уравнение 3
,
voxels
voxels
mem
где N представляет собой количество вокселей в подобласти, S представляет собой количество байт, используемых для представления объема каждого вокселя в подобласти, а S представляет собой размер оперативной памяти консоли. Такая конфигурация позволяет загрузить в оперативную память консоли 44 представляющую интерес подобласть для каждого уровня изображения. По сравнению со случаями, когда на дисковой памяти консоли 44 хранятся некоторые или все уровни изображения, такая конфигурация обычно позволяет значительно более быструю обработку уровней изображения, если по существу не мгновенную с точки зрения пользователя-человека.
voxel
6
8
В качестве примера, устанавливая TS меньше примерно 20 миллисекунд (мс), можно ожидать получение изображений с приемлемым качеством для применения в кардиологии. За обычное время сканирования, равное примерно 8 секундам (с), можно получить примерно 400 уровней изображения. Если значения вокселей хранятся в виде чисел с плавающей запятой, S обычно равно 4. Следовательно, ограничение размера подобласти до 10 вокселей (т.е. в диапазоне вокселей, равном 64*64*64), что является обоснованным для получения изображения коронарных артерий, потребует оперативной памяти, равной примерно 1,6 Гбайт. Очевидно, что эти требования к памяти существенно ниже, чем те, которые потребовались бы для представляющей интерес области, включающей уровень полного изображения, который обычно составляет порядок, равный 10 вокселям (т.е. 512*512*512 вокселей), и находится в пределах оперативной памяти, доступной для обычного компьютера.
На этапе 212, также описанном ниже, вычисляют веса для уровней изображения для требуемой фазы сердечной деятельности.
На этапе 214 уровни изображения объединяют или суммируют согласно вычисленному взвешиванию для генерации волюметрических данных изображения, характерных для представляющей интерес подобласти, в требуемой фазе сердечной деятельности. Если представляющая интерес подобласть для каждого из требуемых уровней изображения не была загружена в оперативную память перед вычислением весов 212 для изображения, информация загружается в оперативную память перед обработкой уровней изображения.
На этапе 216 волюметрические данные изображения, характерные для представляющей интерес подобласти, в требуемой фазе, отображаются в воспринимаемой человеком форме на консоли 44 для просмотра и/или дальнейшей обработки, производимой пользователем.
При необходимости, на этапе 218 пользователь может выбрать другую фазу для просмотра, и обработка вернется на этап 212. Особое преимущество приведенной ниже конфигурации заключается в сохранении в оперативной памяти только представляющих интерес подобластей различных уровней изображения и в отсутствии необходимости в проведении последующей обратной проекции объема при каждом изменении фазы сердечной деятельности, обработка уровней изображения обычно может выполняться практически мгновенно, с точки зрения пользователя-человека. Следовательно, выбор фазы и ополчаемая в результате визуализация могут выполняться интерактивно.
1
n
Один из способов генерации и взвешивания уровней изображения согласно клиновому способу описан со ссылкой на Фиг.3 и 4. Ссылаясь на Фиг.3А, проекционные данные включают в себя проекционные данные 302...302, получаемые в каждом из множества позиций траектории 304 сканирования, проводимого обычно по спирали.
1
m
На этапе 402 и со ссылкой на Фиг.3В к данным проекции 302 применяется повторная сортировка веерообразного пучка в параллельный пучок по глубине таким образом, чтобы сгенерировать множество наборов сегментов проекционных данных 306...306, содержащих равноудаленные параллельные проекции. Каждый сегмент 306 проекционных данных характеризуется начальным угловым положением  и угловой длиной Δθ, выбираемой так, как описано выше. Следует отметить, что смежные сегменты данных 306 частично перекрываются, чтобы облегчить переходное взвешивание между последовательными уровнями изображения.
1
m
1
m
На этапе 404 со ссылкой на Фиг.3С выполняют операцию взвешивания для соответствующих сегментов проекционных данных 306...306, свертку с ramp-фильтром и обратную проекцию для генерации уровней изображения 308...308. Каждый слой характеризуется начальным угловым положением  и угловым смещением Δθ, которые соответствуют начальному угловому положению и угловому смещению соответствующего сегмента 306 проекционных данных.
Как показано, для каждого уровня 308 изображения выполняют взвешивание согласно функции трапецеидального взвешивания. Угловое перекрытие и переходное взвешивание преимущественно выбирают таким образом, чтобы в перекрывающихся областях 312 сумма взвешиваний, примененных к смежным уровням изображения, оставалась постоянной.
На этапе 406 взвешенные уровни изображения затем суммируются таким образом, чтобы можно было сгенерировать волюметрические данные изображения, характерные для представляющей интерес подобласти. Следует отметить, что для получения максимального временного разрешения в сгенерированном изображении суммирование предпочтительно ограничивают теми сегментами 308 изображения, которые обеспечивают полный набор КТ данных (например, при объединении угловая длина равна примерно π). Суммирование может выполняться для сегментов 308 изображения, получаемых примерно в одной и той же точке в повторяющемся цикле сердечной деятельности. Особым преимуществом этого способа является то, что изменение требуемой фазы требует смещения по сегментам 308 изображения, которые объединены для формирования волюметрических данных. По сравнению с временной повторной сортировкой проекционных данных по глубине, в этом способе отсутствует необходимость в выполнении на компьютере дорогостоящей операции обратной проекции для каждой фазы сердечной деятельности, выбираемой пользователем. Более того, ограничение уровней изображения в подобластях до размера, который может быть загружен в относительно более быстродействующую оперативную память консоли 44, облегчает интерактивную визуализацию различных фаз сердечной деятельности.
Одним из недостатков вышеприведенного способа технологии является то, что фаза сердечной деятельности может сдвигаться только с угловыми приращениями Δθ (или устанавливаться другим способом согласно приращению по времени). Способ, который обеспечивает дополнительную гибкость при выборе фазы, описан ниже со ссылкой на Фиг.5 и 6. Задачей является присвоение веса уровням изображения таким образом, чтобы весовые профили, зависящие от проекции, для каждого вокселя в представляющей интерес подобласти, полученного с учетом веса каждого уровня изображения и весов проекций, используемых для конструирования каждого уровня изображения, аппроксимировали их по взвешивающей функции, которая может быть легко вычислена для произвольной точки фазы. Как описано ниже, это может выполнить путем взвешивания уровней изображения, таким образом, чтобы аппроксимировать взвешивание, которое может быть получено при помощи ECR способа.
1
2
3
ЭКГ пациента схематически показана под ссылочной позицией 502 для примерно трех (3) циклов сердечных сокращений; представляющая интерес 504 случайно выбранная точка фазы в каждом из циклов обозначена как 504, 504 и 504.
2
На этапе 602 генерируют взвешивающую функцию 506 сердечной деятельности. Взвешивающая функция 506 сердечной деятельности выбирает проекционные данные, которые соответствуют относительно времени желаемой точке 504 фазы. В одном из вариантов осуществления генерируют взвешивающую функцию cos. Такая взвешивающая функция назначает относительно большие веса к лучам, которые относительно времени находятся рядом с точкой фазы и обращаются в ноль на границах окна синхронизации.
f
l
На этапе 604 генерируют взвешивание 508 освещения. Окно освещения для данных вокселей 510 представляет собой функцию геометрии системы наблюдения и характеризуется угловыми положениями θ и θ, которые определяют соответствующие первый и последний проекционные углы, под которыми освещаются воксели. Окно освещения может быть численно определено для каждого вокселя объема реконструкции. В одном из вариантов осуществления генерируют трапецеидальное взвешивание. Угловой диапазон, определяемый окном освещения, определяет количество избыточных данных, доступных для каждого вокселя в объеме реконструкции. Полученное множественное покрытие в пределах окна освещения используется при взвешивании сердечной деятельности.
На этапе 606 взвешивающие функции сердечной деятельности 506 и освещения 508 перемножают для генерации объединенной взвешивающей функции 510.
На этапе 608 объединенную взвешивающую функцию 510 нормализуют для генерации нормализованной взвешивающей функции 512 для каждого вокселя или для блоков или групп вокселей. Более конкретно, взвешивающие функции нормализуют относительно всех различных π партнеров. В контексте ECR способа, нормализованная взвешивающая функция 512 обычно применяется во время обратной проекции. Однако, как описано ниже, нормализованная взвешивающая функция служит в качестве референсной функции, которая используется при вычислении взвешивания уровней изображения.
На этапе 610 определяют уровни 308 изображения, сгенерированные только из проекций, которые освещают воксели.
2
На этапе 612 вычисляют взвешивание для определенных уровней изображения. Более конкретно, вычисляют веса уровней изображения для минимизации разницы между зависящим от проекции весовым профилем каждого вокселя и нормализированной взвешивающей функцией 512 для этого вокселя. Это может быть выполнено при помощи неитеративной замкнутой формулы для вычисления разницы, например, путем минимизации I нормы гильбертого пространства, зависящей от проекции разницы.
На этапе 614 отрицательные значения весов обрезаются до нуля.
На этапе 616 веса уровней изображения нормализируют. Более конкретно, веса изображения нормализируют таким образом, чтобы сумма весов для уровней, сдвинутых в потоке получения данных путем умножений на π, была равна единичному значению или, иными словами, общему значению. Выполнение этой нормализации и выбора углового смещения между последовательными уровнями изображения, которое хорошо разделяется по времени, равное половине поворота гентри, гарантирует, что сумма весов, полученных для всех проекций, используемых для реконструкции вокселя в конечном изображении и принадлежащая одному и тому же углу между 0 и π, также будет равна единичному значению. Как известно из ECR способа, вычисление весов изображения может быть выполнено быстрее путем вычисления весов для областей или блоков вокселей, чем при вычислении отдельных весов для каждого вокселя.
1
2
3
Также возможны и другие варианты. Как изображено на Фиг.4, точки 504 фазы показаны как представляющие собой одну и ту же точку в цикле 502 сердечных сокращений. Однако качество реконструированного изображения может быть чувствительным к изменениям сердечного ритма. Следовательно, точные местоположения точек 504, 504 и 504 фазы в соответствующих циклах сердечных сокращений могут быть модифицированы индивидуально перед вычислением временной взвешивающей функции 506.
1
2
3
Таким образом, пользователь получает возможность подгонять одну или несколько точек 504, 504 и 504 фазы относительно цикла сердечных сокращений до требуемой величины. Взвешивающие функции обновляют соответствующим образом и генерируют и интерактивно отображают обновленные изображения. Поскольку взвешивающая функция по существу может корректироваться непрерывно, пользователь получает дополнительную возможность уменьшения эффектов пертурбаций в цикле сердечных сокращений, что к тому же облегчается интерактивной природой процесса. Очевидно, что такая конструкция облегчает интерактивную оптимизацию качества изображения, позволяя пользователю просматривать результаты модификации в интерактивном режиме и по желанию применять дальнейшие модификации.
Хотя приведенное выше описание было изложено относительно референсных функций, относящихся к аппроксимации, которые могут быть получены при помощи ECR способа, также могут использоваться другие необходимые референсные взвешивающие функции.
Изобретение было раскрыто в отношении предпочтительных вариантов осуществления. Конечно, приведенное выше описание может подсказать специалисту другие модификации и изменения. Следует иметь в виду то, что настоящее изобретение следует трактовать как включающее в себя все такие модификации и изменения в той мере, в какой они предусмотрены объемом прилагаемой формулы изобретения.