RU2605519C2 - Двухпроходная коррекция металлического артефакта с компенсацией движения для изображений срезов компьютерной томографии - Google Patents

Двухпроходная коррекция металлического артефакта с компенсацией движения для изображений срезов компьютерной томографии Download PDF

Info

Publication number
RU2605519C2
RU2605519C2 RU2014121096/28A RU2014121096A RU2605519C2 RU 2605519 C2 RU2605519 C2 RU 2605519C2 RU 2014121096/28 A RU2014121096/28 A RU 2014121096/28A RU 2014121096 A RU2014121096 A RU 2014121096A RU 2605519 C2 RU2605519 C2 RU 2605519C2
Authority
RU
Russia
Prior art keywords
image
images
projection
slice
artifact
Prior art date
Application number
RU2014121096/28A
Other languages
English (en)
Other versions
RU2014121096A (ru
Inventor
Михаэль ГРАСС
Хольгер ШМИТТ
Дирк СХАФЕР
Original Assignee
Конинклейке Филипс Н.В.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Конинклейке Филипс Н.В. filed Critical Конинклейке Филипс Н.В.
Publication of RU2014121096A publication Critical patent/RU2014121096A/ru
Application granted granted Critical
Publication of RU2605519C2 publication Critical patent/RU2605519C2/ru

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/12Arrangements for detecting or locating foreign bodies
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A61B6/5241Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5264Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/541Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Optics & Photonics (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Data Mining & Analysis (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

Изобретение относится к компьютерной томографии (CT), в частности к коррекции изображений, получаемых с помощью CT. Способ включает в себя этапы, на которых принимают изображение среза и множество изображений СТ-проекций. Выполняют сегментацию в срезе области артефакта. Используют сегментированную таким образом область артефакта для идентификации в каждом из множества изображений СТ-проекций отпечатка, соответствующего области артефакта. Сопоставляют в каждом из множества изображений СТ-проекций идентифицированный отпечаток в изображении СТ-проекции с соответствующим отпечатком в этом изображении СТ-проекции. Причем соответствующий отпечаток представляет положение элемента во время получения этого изображения СТ-проекции. Формируют по пользовательскому запросу новое изображение среза путем реконструирования на основании множества соответствующих отпечатков нового изображения среза, компенсируя движение элемента во время получения изображений СТ-проекций. Повторяют предыдущие этапы для заново реконструированного изображения среза с компенсированным движением вместо предыдущего изображения среза. Выполняют интерполяцию для каждого из множества изображений СТ-проекций и используют соответствующий отпечаток в соответствующем изображении СТ-проекции в качестве ограничения для интерполяции, чтобы таким образом получить множество интерполированных изображений СТ-проекций. Используют множество интерполированных изображений СТ-проекций для реконструкции изображения корректированной версии изображения среза. Выводят корректированную версию изображения среза. 4 н. и 8 з.п. ф-лы, 5 ил.

Description

Область техники, к которой относится изобретение
Настоящее изобретение относится к устройству для коррекции изображения поперечного среза для артефакта в изображении, к способу коррекции изображения поперечного среза для артефакта в изображении, к медицинской системе обработки изображений для коррекции изображения поперечного среза для артефакта в изображении, к компьютерному продукту и к машиночитаемому носителю.
Уровень техники
Все более возрастающее количество пациентов, которые подвергаются исследованию при помощи компьютерной томографии (СТ) сердца, имеют в своем сердце металлические имплантаты, такие как искусственные клапаны сердца, бивентрикулярные кардиостимуляторы, стенты, устройства Амплатцера или подобные им. В перспективе компьютерно-томографическое сканирование сердца со стробированием дает хорошие результирующие изображения для мягкой ткани, поскольку сканирующие проходы выполняются во время фазы покоя сердца, чтобы таким образом минимизировать вызванные движением артефакты изображения. К сожалению, небольшое, но все же заметное движение металлического имплантата внутри окна стробирования может привести к артефактам изображения. Система компьютерной томографии сердца раскрыта в заявке WO 2007/060572, принадлежащей заявителю настоящей заявки. В US 2010/0183214 описан способ уменьшения артефактов в изображениях рентгеновской компьютерной томографии, вызванных материалами с высоким значением коэффициента ослабления. Этот способ включает в себя объединение проекций, полученных под одинаковыми углами обзора, для формирования набора данных плоскости проекции, из которого формируется переформатированная проекция.
Раскрытие изобретения
Задача настоящего изобретения состоит в создании альтернативного устройства и способа коррекции в компьютерно-томографических (СТ) изображениях срезов для артефактов.
Задача настоящего изобретения решается объектами независимых пунктов формулы изобретения, причем дополнительные варианты осуществления включены в зависимые пункты формулы изобретения.
Следует отметить, что описанные ниже аспекты настоящего изобретения одинаково применимы к способу коррекции изображения среза для артефакта в изображении и к медицинской системе обработки изображений для коррекции изображения поперечного среза для артефакта в изображении, к компьютерному программному элементу и к машиночитаемому носителю.
Изобретение соответствует прилагаемой формуле изобретения.
В соответствии с первым аспектом предложено устройство для коррекции изображения поперечного среза для артефакта в изображении. Изображение среза основано на изображениях СТ-проекций объекта, полученных ранее в разное время и под различными углами в то время, как часть объекта с большим декрементом ослабления находится в движении, вызывая таким образом артефакт в изображении. Устройство включает в себя блок ввода, обрабатывающий блок и блок вывода. Блок ввода выполнен с возможностью приема требующего коррекции изображения среза, а также выполнен с возможностью осуществления доступа к изображениям СТ-проекций. Обрабатывающий блок включает в себя блок сегментирования («сегментатор»), блок идентификации («идентификатор»), блок сопоставления («сопоставитель»), блок интерполяции («интерполятор»), а также блок реконструкции изображения среза (реконструктор изображения).
Сегментатор выполнен с возможностью сегментирования области артефакта в срезе.
Идентификатор выполнен с возможностью использования сегментированной области артефакта для идентифицирования в каждом из множества изображений СТ-проекций отпечатка, соответствующего области артефакта.
Сопоставитель выполнен с возможностью сопоставления в каждом из множества изображений СТ-проекций идентифицированного отпечатка в изображении компьютерно-томографической проекции с соответствующим отпечатком в этом изображении СТ-проекции. Соответствующий отпечаток представляет собой положение части во время получения этого изображения СТ-проекции.
Интерполятор выполнен с возможностью формирования для каждого из множества изображений СТ-проекций интерполированного изображения СТ-проекции путем использования соответствующего отпечатка в изображении СТ-проекции в качестве ограничения для интерполяции в этом изображении СТ-проекции, чтобы таким образом получить множество интерполированных изображений компьютерно-томографических проекций.
Реконструктор изображения выполнен с возможностью использования множества интерполированных изображений СТ-проекций для реконструкции корректированной версии изображения среза.
Блок вывода выполнен с возможностью вывода корректированной версии изображения среза.
Корректированное изображение среза предварительно реконструируется по изображениям СТ-проекций, полученных при сканировании, выполненном СТ-сканирующей системой.
Часть объекта с большим декрементом ослабления имеет более высокий коэффициент ослабления используемого сканером СТ-излучения (например, рентгеновского), чем средний коэффициент ослабления остальной части объекта. Примером является металлический элемент, встроенный в органическую мягкую ткань, такую как сердце, проходящее через свои циклы сердечных сокращений, и тем самым перемещающаяся относительно направляющей проволоки катетера, или искусственный клапан сердца, находящийся в сердце во время сканирования.
Однако это является только примером, и настоящее изобретение может использоваться для любого другого объекта, имеющего элемент с большим декрементом ослабления, который изменяет свое положение или ориентацию во время получения компьютерно-томографического изображения.
Изображения СТ-проекций могут рассматриваться в качестве «исходных изображений» для «целевого» изображения, образуемого срезом. Устройство работает преимущественно на уровне исходного изображения, а не на уровне изображения среза.
Артефакт изображения в предварительно реконструированном изображении среза является результатом используемого алгоритма реконструкции и движения элемента с большим декрементом ослабления во время получения изображений СТ-проекций. Поскольку устройство работает преимущественно на уровне исходного изображения, то есть на изображениях СТ-проекций, может быть обеспечена более совершенная коррекция, поскольку устройство работает на том уровне, где искажения и размытие являются обычно менее выраженными по сравнению с искажениями и размытием в изображении среза. Вместо того, чтобы рассматривать металлические артефакты («металлические тени») в изображении среза как простую неприятность, устройство использует информацию, обеспеченную самим артефактом изображения, для того чтобы получить корректированные изображения в высоком качестве. Артефакт может рассматриваться как визуальное представление «реконструированного среднего положения» металлического объекта, то есть «изображения» элемента, усредненного во времени по различным положениям во время получения изображений СТ-проекций. Область артефакта определяет «полосу», образуемую различными положениями, которые элемент принимал во время получения изображений.
Каждое из изображений проекций получается сканером компьютерного томографа в конкретный момент времени и под конкретным углом проекции.
Каждый из соответствующих отпечатков в любом отдельном изображении СТ-проекций представляет собой проекцию части в «истинном» мгновенном положении в конкретный момент времени. В противоположность этому, «положение в изображении» металлического элемента в изображении среза, представленное сегментированной областью артефакта, является видимой «областью металлического объекта». Идентификация отпечатка области артефакта на уровне изображения проекции производится в одном варианте осуществления прямой проекцией, которая является операцией, обратной по отношению к обратной проекции, используемой в предыдущем первом проходе реконструкции изображения среза по изображениям СТ-проекций.
Соответствующие отпечатки, кодирующие истинные положения, находятся сопоставителем на уровне изображения проекции путем использования формы отпечатка прямоспроектированной области артефакта, реконструированной из уровня изображения среза.
Область артефакта включает в себя области изображения среза, представляющие металлический элемент, но может также включать в себя области воксела, не представляющие надлежащим образом металлический элемент, но являющиеся чистыми артефактами, такими как полосы, пятна и искажения, вызванные алгоритмом реконструкции и движением металлического элемента во время получения изображения СТ-проекции. В соответствии с одним вариантом осуществления сегментатор конфигурирован таким образом, чтобы ограничить сегментацию областями воксела, представляющими металлический элемент надлежащим образом.
В одном варианте осуществления сегментация может включать в себя подобласти, которые являются несвязанными, то есть для любой пары точек внутри сегментации путь, соединяющий эту пару, может в части своей длины не принадлежать сегментированной области.
В еще одном варианте осуществления сегментатор осуществляет отображение только «связанных областей», то есть любая пара точек внутри сегментации может быть соединена путем, который по всей своей длине принадлежит сегментированной области.
В соответствии с одним вариантом осуществления область артефакта в изображении среза или в новом изображении среза представляет собой среднее во времени положение части во время ее движения при получении изображений СТ-проекций, причем область артефакта является результатом предыдущей реконструкции из изображений СТ-проекций изображения среза или нового изображения среза.
В одном варианте осуществления производится процесс регистрации на изображениях проекции и оценивается смещение отпечатка артефакта или металлической тени в каждом из изображений проекции по отношению к прямо спроектированному усредненному отпечатку металлической тени. В процессе регистрации геометрическая форма прямо спроектированной металлической тени используется для того, чтобы упорядочить регистрацию. «Истинные» положения движущейся части, закодированные соответствующими отпечатками, найденными сопоставителем, затем используются во втором проходе реконструкции, чтобы таким образом получить корректированную версию изображения среза.
В соответствии с одним вариантом осуществления множество интерполированных изображений СТ-проекций формируется путем замены в каждом из множества изображений СТ-проекций информации изображения внутри соответствующего отпечатка в этом изображении СТ-проекции на интерполированную информацию изображения. Интерполированная информация изображения либо получается путем использования усредненной информации изображения (пиксельных значений) снаружи отпечатка, чтобы таким образом существенно исключить путем усреднения отпечаток проволоки и получить корректированное изображение среза, на котором проволока не видна. Затем используется пиксельная информация для окружающей ткани, и выбранное изображение будет показано пользователю так, как будто оно изначально не содержало металлического элемента. В другом варианте осуществления достигается противоположное, а именно ищется улучшенное представление с корректированным контрастом металлического элемента, внедренного в окружающую ткань. Корректированный контраст может быть достигнут путем замены интерполированных пикселей более высокими или более низкими пиксельными значениями, чем пиксельные значения снаружи каждого из соответствующих отпечатков. В этом случае осуществляются манипуляции с информацией изображения внутри соответствующих отпечатков, и она заменяется корректированными значениями интенсивности для лучшего контраста. При реконструкции используется шкала Хаунсфилда.
В соответствии с одним вариантом осуществления по пользовательскому запросу и до интерполяции и реконструкции корректированного изображения формируется новое изображение среза для возврата в устройство в качестве входа. Реконструктор выполнен с возможностью реконструировать новое изображение среза на основании множества соответствующих отпечатков, компенсируя таким образом движение элемента во время получения изображений СТ-проекций. Блок ввода затем получает заново реконструированное изображение среза с компенсированным движением, и затем запускается обрабатывающий блок, который обрабатывает заново реконструированное изображение среза вместо предыдущего изображения среза.
Другими словами, этап предварительной обработки для предварительной обработки изображений срезов реализуется в соответствии с одним вариантом осуществления для получения новых изображений срезов, которые затем используются для возврата в устройство в качестве входа и обрабатываются вместо предыдущего изображения среза. Таким образом может быть получено корректированное контрастное изображение области артефакта с меньшим размытием.
Устройство затем использует это новое корректированное контрастное изображение среза в качестве входного и обрабатывает его для того, чтобы снова получить улучшенную версию изображения среза с корректированным контрастом.
В одном варианте осуществления новое изображение среза выводится в качестве окончательного результата, обходя таким образом интерполятор.
Реконструктор изображения в устройстве может быть выделенным реконструктором изображения, однако в других вариантах осуществления конструктор изображения, ранее использованный для получения изображения среза, также может использоваться для выполнения реконструкции. В этом случае устройство содержит подходящие интерфейсные средства и средства программного интерфейса приложения (API) для сопряжения с собственной системой реконструктора изображения. В этом случае могут быть необходимы соответствующее перепрограммирование и конфигурирование собственного реконструктора изображения.
В общем, устройство в одном варианте осуществления обеспечивает двухпроходную схему уменьшения металлического артефакта с компенсацией движения.
В одном варианте осуществления устройство применяется к перспективно стробируемому получению изображений сердца. Данное устройство уменьшает артефакты в сканированных изображениях сердца, когда металлические объекты остаются в движении внутри окна стробирования. Данное устройство также может быть применено в инвазивных компьютерно-томографических процедурах для коррекции артефактов, образующихся из-за перемещения игл или катетеров. Следует понимать, что настоящее изобретение не ограничивается компьютерно-томографической визуализацией сердца или металлическими объектами, но настоящее изобретение может быть применено к любому анатомическому или неанатомическому объекту, включающему в себя особенности с большим декрементом ослабления, где элемент движется во время получения компьютерной томограммы.
Определения
Термин «артефакт изображения» или «металлическая тень» включает в себя пятна, полосы и искажения, видимые наблюдателю изображения среза.
Термин «отпечатки» означает пиксельные области, представляющие проекции либо части, либо прямо спроектированного (функционально) изображения артефакта.
Термин «сегментация» означает область вокселов или пикселей, которые вместе определяют объект в плоскости изображения. О каждом пикселе или вокселе в плоскости изображения затем можно сказать, находится ли он внутри объекта или снаружи объекта. Сегментация также определяет внешний контур, границу и, следовательно, «форму» объекта.
Термины «движение» или «изменение положения» металлического объекта включают в себя любое их строгое преобразование, то есть вращение и/или перенос объекта, при условии, что изменение положения является различимым в изображении проекции. Проекция движения элемента может быть выражена как плоское строгое преобразование на плоскости изображения проекции отпечатка элемента.
Краткое описание чертежей
Примерные варианты осуществления настоящего изобретения будут описаны далее со ссылками на прилагаемые чертежи, в которых:
Фиг. 1 показывает вид в перспективе компьютерно-томографической сканирующей системы;
Фиг. 2 показывает более подробный вид компьютерно-томографической сканирующей системы, изображенной на Фиг.1, в боковой проекции;
Фиг. 3 показывает работу реконструктора изображения среза компьютерного томографа;
Фиг. 4 показывает устройство для коррекции артефакта в изображении поперечного среза в соответствии с одним вариантом осуществления настоящего изобретения;
Фиг. 5 показывает блок-схему способа коррекции артефакта в изображениях поперечного среза в соответствии с одним вариантом осуществления настоящего изобретения.
Осуществление изобретения
На Фиг. 1 показана примерная компьютерно-томографическая сканирующая система («третьего поколения»). Система 100 включает в себя корпус 102, имеющий канал. Гентри 120 расположен внутри канала в виде жесткой цилиндрической структуры с возможностью вращения. Гентри 120 включает в себя в жесткой и противоположной взаимосвязи рентгеновскую трубку 130 и массив 115 детекторов. Диагностический стол 110 проходит в гентри 120 через центр цилиндра гентри 120. Интересующий объект, такой как человеческое сердце 135 пациента, лежащего на столе 110, позиционируется в центральной точке гентри 120. Сердце 135 содержит металлический объект, такой как направляющая проволока 140 катетера. Из-за позиционирования сердца 135 и, следовательно, проволоки 140 в центре или в «фокусе» гентри 120 проволока 140 во время сканирования всегда располагается между рентгеновской трубкой 130 и массивом 115 детекторов, в то время как гентри 120 вращается вокруг проволоки 140. Из-за жесткого расположения в гентри 120 рентгеновская трубка 130 и детектор 115 всегда находятся в противоположном пространственном отношении друг к другу. Компьютерно-томографическая сканирующая система 100 позволяет получать изображения поперечных срезов сердца 135 и проволоки 140 в плоскости 125 среза.
Фиг. 2 представляет собой вид сканера 100 сбоку в разрезе и в виде сверху на плоскость 125 среза. Компьютерно-томографическая сканирующая система 100 дополнительно включает в себя рабочую станцию 240 для того, чтобы управлять работой сканера 100, систему 215 сбора данных (DAS), реконструктор 220 изображения и базу 225 данных для хранения полученных изображений среза. Изображения среза могут быть получены из базы 225 данных и выведены рабочей станцией 240 для представления на экране 230.
По мере того как гентри 120 вращается вокруг сердца 135, получают серию изображений проекции при различных углах проекции α («альфа»), причем каждое изображение получают в разное время во время вращения гентри по мере того, как гентри 120 вращается в угловых интервалах от угла альфа до угла β («бета»).
В операции получения изображений рентгеновский источник 2 30 проецирует веерообразный рентгеновский пучок на детектор 115 на противоположной стороне цилиндра гентри 120. Веерообразный пучок рентгеновских лучей формируется из набора отдельных остронаправленных рентгеновских лучей 210а-с. Остронаправленные лучи 210a-b испускаются из рентгеновского источника 130 с начальной интенсивностью. Остронаправленные лучи 210а-с ослабляются по мере того, как они проходят через сердце 135, и ослабленные таким образом остронаправленные лучи выходят из сердца 135 с его противоположной стороны и затем попадают на массив 115 детекторов. Массив 115 детекторов образован детекторными элементами. Примерный детекторный элемент 115а показан на Фиг. 3.
Каждый детекторный элемент 115а формирует электрический ток обратно пропорционально ослаблению, которое претерпевает остронаправленный рентгеновский луч, попадающий на этот детекторный элемент 115а. Электрический ток затем считывается и принимается системой 215 сбора данных и преобразуется с помощью подходящего аналого-цифрового преобразования в численное цифровое значение. Набор цифровых значений (″необработанные данные детектора″) для каждого детекторного элемента 115а затем записывается в матричной структуре, формирующей цифровое изображение 310а,b проекции, по одному изображению 310а или 310b проекции для каждого угла проекции альфа или бета. Позиция строки и столбца i, j («пиксель») в матрице соответствует позиции строки и столбца i, j соответствующего детекторного элемента 115а в массиве 115 детекторов. Для каждого изображения 310а,b проекции регистрируются его время получения и угол проекции и последовательность изображений 310а,b проекции затем либо сохраняется в базе 225 данных и/или передается в реконструктор 220 изображения.
Реконструктор 220 изображения использует известные алгоритмы реконструирования изображения, основанные на обратном преобразовании Радона, для того чтобы сформировать изображение 320 среза, показывающее поперечное сечение в плоскости 125 изображения поперечного сечения сердца 135 с проволокой 140. Другими словами, изображение 320 среза является результатом обратного преобразования Радона (или одной из его многих числовых реализаций), примененного к последовательности изображений 310а,b для каждого угла проекции альфа и бета. Реконструированное таким образом изображение 320 среза затем может быть сохранено в базе 225 данных.
Поскольку проволока 140 катетера имеет более высокое значение коэффициента ослабления рентгеновского излучения, чем окружающая мягкая ткань, остронаправленный рентгеновский луч 210b, проходящий через проволоку 140, ослабляется больше, чем остронаправленные рентгеновские лучи 210а,с, проходящие не через проволоку 140, а через окружающую мягкую ткань сердца. По мере того как сердце проходит через свой цикл сердечных сокращений, находящейся в нем проволоке 140 придается движение.
В компьютерно-томографической системе, показанной на Фиг.2, используется схема визуализации сердца с перспективным стробированием. Для лучшего отображения сердца в конкретной фазе цикла сердечных сокращений используется подходящее электрокардиооборудование 205. Электродные прокладки, присоединенные к груди пациента вокруг области сердца, принимают электрические сигналы, указывающие на текущую фазу цикла сердечных сокращений. Когда обнаруживается электрический сигнал, соответствующий требуемой фазе цикла сердечных сокращений (например, фазе систолы или диастолы), соответствующая управляющая команда подается от электрокардиооборудования 205 на рабочую станцию 240, которая, в свою очередь, управляет включением/выключением рентгеновской трубки 130 для получения «снимка» изображения проекции точно в нужное время, когда сердце находится в требуемой фазе цикла сердечных сокращений. Таким образом, изображения 310а,b проекции получают для каждого угла проекции альфа в требуемой фазе цикла сердечных сокращений, и тяжелые артефакты изображения, вызванные движением сердца, могут быть по меньшей мере частично устранены, поскольку рентгеновский источник 130 выключается, когда сердце не находится в целевом цикле или фазе сердечных сокращений (в покое). Небольшие движения проволоки 140, тем не менее, остаются при перемещении гентри 120 от одного угла проекции к другому из-за задержек включения/выключения рентгеновской трубки 130 и других несовершенств процедуры перспективного стробирования.
Далее со ссылкой на фиг. 3 более подробно объясняется работа компьютерно-томографического сканера 100 по мере того, как проволока 140 изменяет свое положение из-за двигательной активности сердца. Левая часть Фиг. 3 показывает положение проволоки 140 (для ясности иллюстрации окружающая ткань сердца не показана) в то время, когда изображение 310а проекции получают с углом проекции альфа. Пропорциональный электрический сигнал, формируемый элементом 115а после возбуждения падающим на него остронаправленным рентгеновским лучом 210b, направляется затем в систему 215 сбора данных, где изображение 310а проекции для угла альфа составляется путем регистрации каждого преобразованного в цифру аналогового сигнала для каждого из детекторных элементов. Поскольку проволока 140 имеет более высокий коэффициент ослабления, чем окружающая ее ткань сердца, на цифровом изображении 310а проекции формируется или «отпечатывается» отпечаток 315а проволоки.
За то время, которое необходимо для того, чтобы гентри 120 повернулся к следующему углу получения данных бета, положение проволоки 140 изменяется, как вверху справа на фиг. 3. По мере того как гентри 120 поворачивается от одного угла проекции альфа к следующему углу проекции бета, проволока 140 принимает различные положения в различные моменты времени получения данных.
Вследствие этого отпечаток 315b во втором изображении 310b проекции отличается от отпечатка 315а в предыдущем изображении 310а проекции, полученном с углом проекции альфа. Полученное таким образом множество изображений 310а,b проекции, каждое из которых показывает отпечаток 315а,b или пиксельное значение истинного положения отпечатка проволоки 140 в соответствующий получению изображений проекции момент времени и с соответствующим углом проекции, затем подаются в реконструктор 220 изображения. Реконструктор 220 изображения с помощью обратного преобразования Радона преобразует одно за другим последовательность изображений 310а,b проекции в изображение 320 поперечного среза, представляющее собой поперечное сечение сердца 135 в плоскости 125 изображения среза. Поскольку проволока 140 изменяет свое положение во время получения данных, реконструированная область изображения, представляющая проволоку 140, выглядит как артефакт 330 размытого изображения, окружающий упомянутую область изображения. Истинная граница изображения проволоки 140 на изображении может быть неразличимой даже для натренированного медицинского глаза. Артефакт формируется из «металлической тени», размазанной поперек и вокруг области изображения, представляющей реконструированную проволоку на изображении среза. Артефакт 330 схематично показан на Фиг. 3. В результате вызванных движением размытия и искажения артефакт 330, соответствующий проволоке 140, выглядит намного более толстым, чем фактическая толщина отпечатков 315а,b проволоки на соответствующих изображениях 310а,b проекции.
Фиг.4 показывает один вариант осуществления устройства для коррекции изображения среза. Это устройство включает в себя блок ввода или средство 415 входного интерфейса и блок вывода или средство 470 выходного интерфейса. Устройство 400 дополнительно включает в себя обрабатывающий блок 400, включающий в себя сегментатор 405, идентификатор 410, сопоставитель 415 и интерполятор 420. Также имеется реконструктор 430 изображения, который может принимать выходные сигналы от сопоставителя 415 или интерполятора 420. Средство 450 интерфейса обеспечивает взаимодействие с системой 215 сбора данных или базой 225 данных приема подлежащего коррекции содержащего артефакт изображения 330. Средство 450 интерфейса также выполнено с возможностью осуществления доступа к множеству изображений 310а,b проекции, которые ранее использовались реконструктором 220 изображения компьютерного томографа для того, чтобы реконструировать изображение 330 среза, содержащее артефакт 330.
Компоненты 450, 405, 410, 415, 420, 430 и 470 устройства выполняются как программные процедуры на обрабатывающем блоке 400. Однако в альтернативных вариантах осуществления также рассматривается распределенная архитектура устройства, в которой все компоненты соединены подходящей коммуникационной сетью. Компоненты также могут быть организованы как выделенные программируемые логические матрицы (FPGA) или аппаратно соединенные автономные микросхемы. Компоненты могут быть запрограммированы на подходящей платформе научных вычислений, такой как Matlab® или Simulink®, и могут затем быть транслированы в, например, процедуры на языке С++ или С, сохраняемые в библиотеке и связываемые при вызове обрабатывающим блоком 400.
Вообще говоря, устройство 400 принимает на средстве 450 интерфейса изображение 320 среза, имеющее артефакт 330 изображения, и изображения 310а,b проекции, использованные ранее для реконструирования среза 320. Изображение 320 среза обрабатывается на основании изображений 310а,b СТ-проекций и изображение среза с корректированным артефактом выводится на средства 470 вывода. Корректированное изображение либо не содержит представления металлической проволоки (оно удаляется, включая артефакт, вызванный движением проволоки 140 во время получения данных компьютерной томографии), либо представление металлической проволоки улучшается, поскольку она показывается с более высоким контрастом и с уменьшенными артефактами. В одном варианте осуществления устройство также выполнено с возможностью итеративно улучшать по пользовательскому запросу представление артефакта, показываемого в изображении 330, путем выполнения этапа предварительной обработки, причем новое изображение среза показывает артефакт менее размытым, чем в первоначальном срезе 320. Изображение среза с улучшенным артефактом затем может быть рассмотрено пользователем и, в случае его одобрения, оно становится новым предварительно обработанным изображением среза с артефактом, которое затем возвращается обратно в сегментатор 405 и обрабатывается обрабатывающим блоком 4 00 для того, чтобы получить корректированное изображение среза.
Работа
Изображение 320 среза, имеющее артефакт 330, принимается средством 450 ввода и передается в сегментатор 405, и устанавливается местоположение в хранилище 225 изображений 310а,b проекции для обеспечения доступа обрабатывающему блоку 400 к упомянутым изображениям 310а,b проекции.
Сегментатор 405 сегментирует артефакт 330 изображения путем использования известных методов сегментирования, таких как сегментация по порогу интенсивности значения серого. Каждый воксел, чье значение серого превышает предопределенное и конфигурируемое значение порога серого, помечается и рассматривается как являющийся частью области 330 артефакта изображения и/или области изображения, представляющей металлический элемент 140. Пороговая сегментация может использоваться, поскольку металлическая проволока 140 имеет более высокий коэффициент ослабления рентгеновского излучения, чем окружающий органический материал ткани сердца и область 330 артефакта изображения, включающая в себя область, представляющую проволоку 140, будет казаться более темной (или более светлой в зависимости от цветового кодирования, используемого для изображения 320 среза), чем остальная часть изображения. Поскольку предполагается, что проволока 140 движется во время получения изображений 310а,b проекции, используемых для предыдущей реконструкции изображения 330, область артефакта можно рассматривать как графическое представление усредненного положения проволоки 140 во время получения изображений. Другими словами, артефакт можно считать результатом суперпозиции истинных положений проволоки 140 в любой отдельно взятый момент времени получения соответствующих изображений 310а,b проекции. Структура данных сегментации, кодирующая положения в срезе этих вокселов, которые были установлены сегментатором 405 как принадлежащие области 330 артефакта, затем направляется в идентификатор 410.
Идентификатор 410 может быть реализован как прямой проектор, то есть как численная реализация, основанная на преобразовании Радона. Идентификатор 410 тогда использует структуру данных сегментации для того, чтобы прямо спроектировать на каждое из изображений 310а,b проекции вокселы, составляющие сегментированную область 330 артефакта. Область пикселей в каждом из изображений 310а,b проекции, формирующая изображение, в математическом смысле, прямой проекции, примененной к сегментированной области артефакта, идентифицирует отпечаток области артефакта в каждом из изображений 310а,b проекции. Если F является функцией прямой проекции, a S является набором вокселов в плоскости среза, образующим область артефакта, тогда отпечаток область 330 артефакта является изображением F(S) функции F. Идентификация может быть реализована битовыми масками, ограничивающими действие прямой проекции только вокселами внутри сегментированной области 330 артефакта. Таким образом может быть обеспечено соответствие области артефакта отпечаткам в каждом из изображений 310а,b проекции, хотя область артефакта сама по себе существует только в срезе 320 как результат предыдущей реконструкции изображения реконструктором 220 компьютерно-томографического изображения. Идентифицированные таким образом пиксельные области в каждом из изображений 310а,b СТ-проекций составляют отпечаток, один для каждого изображения 310а,b проекции, области артефакта 330 на срезе. Отпечаток 330а артефакта 330 примерно показан для изображения 310а проекции на фиг.4. Изображения 310а,b СТ-проекций и подходящая структура данных, помечающая пиксели, составляющие соответствующие идентифицированные отпечатки области артефакта, затем передаются сопоставителю 415.
В каждом изображении 310а,b проекции сопоставитель 415 соотносит отпечаток области артефакта (330а для изображения 310а) с соответствующим истинным отпечатком 315а,b проволоки в этом изображении 310а,b проекции. Поскольку область артефакта в изображении среза рассматривается как представление усредненного положения и усредненной формы истинных положений и формы (в представлении проекции) проволоки 140, можно ожидать, что общая форма отпечатка области артефакта будет подобна истинным отпечаткам 315а,b в каждом из изображений проекции 310а,b.
Сопоставитель 415 может быть выполнен с возможностью реализовать алгоритм блочного сопоставления для сканирования каждого из изображений 310а,b проекции с тем, чтобы соотнести идентифицированную область артефакта с отпечатком 315а,b проволоки в соответствующем изображении 310а,b проекции. Сопоставитель 415 использует выбираемую меру схожести для сопоставления в каждой проекции 310а,b формы отпечатка артефакта в этом изображении проекции с формой соответствующего отпечатка 315а или b проволоки, существующего в соответствующем изображении 310а или b проекции, соответственно. Например, в отношении проволоки 140 ожидается, что ее отпечатки артефакта (например, 330а), идентифицированные в каждом из изображений 310а,b проекции будут иметь форму криволинейной полосы, но из-за размытости полоса будет более широкой, чем аналогичные истинные отпечатки 315а,b проволоки в форме полосы.
Положение отпечатка артефакта относительно положения истинного отпечатка 315а,b проволоки в плоскости соответствующего изображения 310а,b проекции может быть описано вектором и матрицей, представляющими строгое преобразование на плоскости. Это преобразование, будучи примененным к истинному отпечатку 315а проволоки, дает существенное совпадение (с точностью до меры подобия) отпечатка 315а проволоки с отпечатком 330а области артефакта проволоки. То же самое применяется к каждому из других изображений 310а,b проекции, к каждому из других отпечатков 315b и соответствующих отпечатков 330b области артефакта. Операция сопоставления схематично показана на Фиг.4 криволинейной двунаправленной стрелкой в подробном виде изображения 310 проекции в правой части Фиг 4. Координатное описание в подходящей структуре данных сопоставленных таким образом истинных отпечатков 315а,b проволоки затем передается в интерполятор 420.
Интерполятор 420 затем заменяет пиксельные значения, составляющие истинные отпечатки 315а,b проволоки, интерполированными значениями. Используя пиксельную область, определенную отпечатками 315а,b, в качестве ограничений для интерполяции, устройство 400 может быть сконфигурировано пользователем так, чтобы полностью удалить путем усреднения металлический проволочный элемент в корректированном срезе или сформировать корректированный срез, показывающий металлическую проволоку с более высоким контрастом по сравнению с первоначальным срезом. Если металлический элемент должен быть полностью удален путем усреднения, интерполятор 420 использует усредненную пиксельную информацию снаружи каждого из истинных отпечатков проволоки, чтобы таким образом заменить пиксельные значения внутри отпечатка проволоки средними пиксельными значениями.
В соответствии с одним вариантом осуществления интерполяция производится путем усреднения по линиям, охватывающим с двух сторон отпечаток 315а,b проволоки в плоскости 310а,b изображения проекции. Первоначально пиксельные значения на линии к любой стороне охватываемой с двух сторон области 315а или 315b отпечатка, но вне области 315а или 315b отпечатка, усредняются и используются для того, чтобы заменить пиксельные значения на линии внутри области 315а или 315b отпечатка. Продолжая аналогичным образом вдоль каждой из линий путем использования скользящих средних значений вдоль каждой из линий, пиксельные значения внутри области 315а или 315b отпечатка последовательно заполняются от границы области отпечатка к ее центру, также используя ранее интерполированные пиксельные значения в области 315а или 315b отпечатка, соответственно.
Если представление металлической проволоки в окончательном корректированном изображении среза должно быть сохранено, интерполятор 420 сбрасывает пиксельные значения внутри отпечатка к такому уровню, чтобы улучшить интенсивность/контраст по сравнению с внешней областью в каждом из изображений проекции и корректирует физические эффекты, такие как шумовое голодание (недостаток соотношения сигнал/шум) и увеличение жесткости излучения. Интерполированные таким образом изображения проекции затем передаются в реконструктор 430 изображения.
Реконструктор 430 реализует операцию, обратную операции прямого проектора 410. Реконструктор 430 изображения реконструирует во «втором проходе» (предыдущая реконструкция в компьютерно-томографической сканирующей системе 100 первоначального среза 320 образует «первый проход») путем обратного проектирования вдоль прямых линий значений из интерполированных изображений проекции на плоскость изображения, чтобы таким образом получить корректированное изображение среза, которое затем выводится средствами 470 вывода и передается для хранения в базу 225 данных. Корректированное изображение может быть затем обработано подходящим средством графического рендеринга для вывода на экран 230.
В одном варианте осуществления устройство выполнено с возможностью обеспечения процедуры для улучшения представления самой области 330 артефакта путем формирования нового изображения среза, которое затем обрабатывается, как было описано ранее, вместо первоначального изображения среза. Другими словами, полученное изображение 330 среза сначала предварительно обрабатывается для получения нового изображения среза, показывающего область артефакта с корректированным по сравнению с первоначальным изображением среза контрастом, и этот новый срез затем возвращается обратно в устройство к сегментатору 205 и обрабатывается в последовательности, описанной выше в отношении первоначального среза. Формирование улучшенного изображения среза таким образом может итеративно повторяться до тех пор, пока не будет получено удовлетворительное изображение среза, показывающее артефакт в соответствии с ожидаемым стандартом с более высокой резкостью.
Если запрашивается новое изображение среза, сопоставитель 415 направляет изображения 310а,b проекции с размеченными соответствующими отпечатками 315а,b проволоки в реконструктор 430 вместо того, чтобы передать изображения 310а,b проекции в интерполятор 420, как было описано ранее. Реконструктор 430 затем использует размеченные истинные отпечатки 315а,b проволоки, найденные сопоставителем 415, для компенсации движения проволоки 140 во время получения изображений 310а,b проекции. Реконструктор 430 осуществляет компенсацию путем регистрации изображений 310а,b проекции вдоль общей системы координат и расположения изображений 310а,b СТ-проекций в соответствии с их временем получения.
Реконструктор 430 затем прослеживает очевидное движение отпечатков 315а,b от самого раннего изображения проекции к самому последнему изображению проекции путем использования ранее зарегистрированных строгих преобразований, описывающих очевидное движение отпечатков проволоки относительно предыдущего отпечатка 330 области артефакта. Прослеженное таким образом движение затем используется в преобразовании Радона, реализованном в реконструкторе 430, при реконструкции улучшенного изображения среза из изображений 310а,b проекции.
Вместо того чтобы использовать прямые линии для реконструкции на основании преобразования Радона, реконструктор 430 использует линии, изогнутые в соответствии с отслеженным движением, для того, чтобы отобразить смещение в соответствии с упомянутым движением в обратной проекции, выполняя таким образом второй проход компенсации движения. В другом варианте осуществления реконструктор 430 осуществляет обратную проекцию вдоль прямых линий проекции, но сдвигает каждое из изображений проекции в соответствии с движением, чтобы таким образом аналогично выполнить второй проход компенсации движения. В одном варианте осуществления строгие преобразования, ранее установленные сопоставителем 415, используются для операций изгиба линии и сдвига плоскости изображения. Реконструированное таким образом изображение среза будет затем показывать область артефакта более резко. Пользователь может затем рассмотреть новый срез после его вывода на экран 230 и решить, удовлетворяет ли его выведенная на экран область артефакта, и если не удовлетворяет, то новый срез затем возвращается обратно в сегментатор 405, обрабатывается как было описано ранее, и передается сопоставителем 415 в реконструктор изображения 430, чтобы таким образом итеративно формировать дополнительно улучшенные изображения среза. Как только пользователь удовлетворен улучшенным изображением среза, изображение среза обрабатывается как было описано ранее и передается сопоставителем 415 интерполятору 420 и затем выводится в качестве окончательного корректированного изображения.
На фиг.5 показана блок-схема способа коррекции артефакта 330 изображения в изображениях поперечных срезов.
На этапе S505 принимают срез 320 и изображения 310а,b проекции или файловые ссылки на них.
На этапе S510 изображение среза сегментируют для получения сегментированной области артефакта.
На этапе S520 используют сегментированную область артефакта для идентификации в каждом из множества изображений проекции отпечаток, соответствующий области артефакта.
На этапе S530 для каждого из множества изображений проекции идентифицированный отпечаток в соответствующих изображениях СТ-проекций сопоставляют с соответствующим отпечатком в этом изображении СТ-проекции. Соответствующий отпечаток представляет собой положение металлического элемента во время получения соответствующего изображения СТ-проекции.
На этапе S540 для каждого из множества изображений СТ- проекций информацию об изображении интерполируют с использованием области изображения, заключенной в соответствующем отпечатке, в качестве ограничения для интерполяции в этом изображении СТ-проекции. Интерполяция обеспечивает множество интерполированных изображений СТ-проекций.
На этапе S550 множество интерполированных таким образом компьютерно-томографических изображений затем реконструируют на втором проходе для получения корректированной версии изображения среза.
На этапе S560 реконструированный таким образом срез затем выводят в качестве корректированной версии изображения среза.
В одном варианте осуществления способ дополнительно содержит возможный этап S570, выполняемый перед вышеописанными этапами S540 и S550. На этапе S570 для реконструкции нового изображения среза используют соответствующие отпечатки в изображениях проекции для того, чтобы таким образом компенсировать движение элемента во время получения изображений СТ-проекций.
В одном варианте осуществления новый срез, реконструированный на этапе S570, напрямую передается на этап S560 в качестве окончательного результата.
Следует понимать, что в вышеприведенном описании проволока используется в качестве примера элемента с высоким коэффициентом ослабления в среде со сравнительно низким коэффициентом ослабления. Описанное устройство может также использоваться для любого другого элемента, металлического или неметаллического, погруженного в объект любого вида или автономного, который движется во время выполнения компьютерной томографии.
В соответствии с одним вариантом осуществления реконструктор отслеживает соответствующие отпечатки во множестве изображений проекции для того, чтобы отследить движение объекта во время получения изображений проекции. Отслеженное движение затем учитывается при реконструкции новых срезов.
В другом примерном варианте осуществления настоящего изобретения предложена компьютерная программа или компьютерный программный элемент, который отличается тем, что он выполнен с возможностью выполнения на подходящей системе этапов способа в соответствии с одним из предыдущих вариантов осуществления.
Компьютерный программный элемент, следовательно, может быть сохранен на компьютерном блоке, который также может быть частью варианта осуществления настоящего изобретения. Этот вычислительный блок может быть выполнен с возможностью исполнения или инициирования исполнения этапов способа, описанного выше. Кроме того, он может быть выполнен с возможностью оперировать компонентами вышеописанного устройства. Вычислительный блок может быть выполнен с возможностью автоматической работы и/или выполнения команд пользователя. Компьютерная программа может быть загружена в рабочую память процессора. Процессор, таким образом, может быть оборудован для выполнения способа по настоящему изобретению.
Этот примерный вариант осуществления настоящего изобретения покрывает как компьютерную программу, которая с самого начала использует настоящее изобретение, так и компьютерную программу, которая посредством обновления превращает существующую программу в программу, которая использует настоящее изобретение.
Далее, компьютерный программный элемент может обеспечивать все необходимые этапы для выполнения процедуры примерного варианта осуществления описанного выше способа.
В соответствии с дополнительным примерным вариантом осуществления настоящего изобретения предложен машиночитаемый носитель, такой как CD-ROM, причем на машиночитаемом носителе сохранен компьютерный программный элемент, который описан выше.
Компьютерная программа может быть сохранена и/или распространена на подходящем носителе, таком как оптический носитель данных или твердотельный носитель, поставляемый совместно или как часть других аппаратных средств, но может также распространяться в других формах, например через Интернет или другие проводные или беспроводные телекоммуникационные системы.
Однако компьютерная программа может также быть представлена по сети, такой как всемирная паутина, и может быть загружена в рабочую память процессора из такой сети. В соответствии с дополнительным примерным вариантом осуществления настоящего изобретения предложен носитель для обеспечения возможности загрузки компьютерного программного элемента, предназначенного для выполнения способа в соответствии с одним из ранее описанных вариантов осуществления настоящего изобретения.
Следует отметить, что варианты осуществления настоящего изобретения описываются в отношении различных сущностей. В частности, некоторые варианты осуществления описываются в отношении способа, тогда как другие варианты осуществления описываются в отношении устройства. Однако специалисту в данной области техники будет понятно, что если не указано иное, в дополнение к любой комбинации функций, принадлежащих одной сущности, также любая комбинация функций, относящихся к другим сущностям, рассматривается как раскрытая в этой заявке. Однако все функции могут сочетаться, обеспечивая синергетические эффекты, которые являются более чем простой суммой функций.
В то время как настоящее изобретение было проиллюстрировано и подробно описано на чертежах и вышеприведенном описании, такие чертежи и описание следует считать иллюстративными или примерными, но не ограничивающими. Настоящее изобретение не ограничивается раскрытыми вариантами осуществления. Другие вариации раскрытых вариантов осуществления могут быть поняты и осуществлены специалистами в данной области техники при осуществлении заявленного изобретения на основании изучения чертежей, описания и зависимых пунктов формулы изобретения.
В формуле изобретения слово «содержащий» не исключает других элементов или этапов, а формы единственного числа элементов не исключают множественного числа. Один процессор или другой блок может выполнять функции нескольких блоков, упомянутых в формуле изобретения. Тот факт, что некоторые средства упоминаются в различных зависимых пунктах формулы изобретения, не указывает на то, что сочетание этих средств не может использоваться с достижением преимущества. Любые ссылочные позиции в формуле изобретения не должны рассматриваться как ограничивающие объем.

Claims (12)

1. Устройство для коррекции изображения поперечного среза для артефакта в изображении, причем указанное изображение среза основано на изображениях (310а,b) компьютерно-томографических (СТ) проекций объекта, полученных в разное время и когда элемент (140) с высоким значением коэффициента ослабления объекта (135) находится в движении, причем артефакт вызван упомянутым движением, причем устройство содержит:
блок (450) ввода;
обрабатывающий блок (400);
блок (470) вывода;
причем
блок ввода выполнен с возможностью приема изображения (320) среза и осуществления доступа к множеству изображений СТ-проекций;
обрабатывающий блок содержит:
сегментатор (405), выполненный с возможностью сегментирования области (330) артефакта в срезе;
идентификатор (410), выполненный с возможностью использования сегментированной таким образом области артефакта для идентификации в каждом из множества изображений СТ-проекций отпечатка (330а,b), соответствующего области артефакта;
сопоставитель (415), выполненный с возможностью сопоставления в каждом из множества изображений СТ-проекций идентифицированного отпечатка (330а,b) в изображении СТ-проекции с соответствующим отпечатком (315а,b) в этом изображении СТ-проекции, причем соответствующий отпечаток представляет собой положение элемента во время получения этого изображения СТ-проекции;
реконструктор (430) изображения, выполненный с возможностью реконструкции по пользовательскому запросу нового изображения среза на основании множества соответствующих отпечатков, компенсируя тем самым движение элемента во время получения изображений СТ-проекций, причем сегментатор, идентификатор и сопоставитель затем обрабатывают упомянутое новое изображение среза вместо изображения (320) среза; и
интерполятор (420), выполненный с возможностью формирования для каждого из множества изображений СТ-проекций интерполированного изображения СТ-проекции путем использования соответствующего отпечатка в изображении СТ-проекции в качестве ограничения для интерполяции в этом изображении СТ-проекции, чтобы таким образом получить множество интерполированных изображений СТ-проекций;
реконструктор (430) изображения, выполненный с возможностью использования множества интерполированных изображений СТ-проекций для реконструкции корректированной версии изображения среза;
блок вывода, выполненный с возможностью вывода корректированной версии изображения среза.
2. Устройство по п. 1, в котором множество интерполированных изображений СТ-проекций формируется путем замены в каждом из множества изображений СТ-проекций информации изображения внутри соответствующего отпечатка в этом изображении СТ-проекции информацией интерполированного изображения.
3. Устройство по п. 1, в котором реконструктор выполнен с возможностью отслеживания соответствующих отпечатков во множестве изображений СТ-проекций, чтобы таким образом отслеживать движение элемента, и выполнен с возможностью использования отслеженного движения для реконструкции нового среза, осуществляя таким образом компенсацию движения.
4. Устройство по любому из пп. 1-3, в котором идентификатор реализован как прямой проектор для прямой проекции области артефакта на каждое из множества изображений СТ-проекций, идентифицируя таким образом отпечаток области артефакта; и
сопоставитель реализован путем использования алгоритма регистрации изображения.
5. Устройство по п. 1, в котором область артефакта в изображении среза или в новом изображении среза представляет собой усредненное во времени положение элемента при его движении во время получения изображений СТ-проекций, причем область артефакта является результатом предыдущей реконструкции из изображений СТ-проекций изображения среза или нового изображения среза.
6. Устройство по п. 1, в котором объект является сердцем человека или животного, и изображения СТ-проекций получают посредством СТ-сканера сердца с перспективным стробированием, и смещение элемента происходит внутри окна стробирования при получении изображений.
7. Устройство по п. 1, в котором элемент является металлическим и является, в частности, любым из искусственного клапана, бивентрикулярного кардиостимулятора, стента, направляющей проволоки для введения стента и устройства Амплатцера.
8. Способ коррекции изображения поперечного среза для артефакта в изображении, причем изображение среза основано на изображениях компьютерно-томографических (СТ) проекций объекта, полученных в разное время и когда элемент с высоким значением коэффициента ослабления объекта находится в движении, причем артефакт вызван упомянутым движением, причем способ содержит этап обработки, включающий в себя этапы, на которых:
принимают (S505) изображение среза и множество изображений СТ-проекций;
выполняют (S510) сегментацию в срезе области артефакта;
используют сегментированную таким образом область артефакта для идентификации (S520) в каждом из множества изображений СТ-проекций отпечатка (330а,b), соответствующего области артефакта;
сопоставляют (S530) в каждом из множества изображений СТ-проекций идентифицированный отпечаток (330а,b) в изображении СТ-проекции с соответствующим отпечатком (315а,b) в этом изображении СТ-проекции, причем соответствующий отпечаток представляет положение элемента во время получения этого изображения СТ-проекции;
формируют по пользовательскому запросу новое изображение среза путем реконструирования (S570) на основании множества соответствующих отпечатков нового изображения среза, компенсируя таким образом движение элемента во время получения изображений СТ-проекций,
повторяют предыдущие этапы S505-S530 для заново реконструированного изображения среза с компенсированным движением вместо предыдущего изображения среза;
выполняют (S540) интерполяцию для каждого из множества изображений СТ-проекций и используют соответствующий отпечаток в соответствующем изображении СТ-проекции в качестве ограничения для интерполяции, чтобы таким образом получить множество интерполированных изображений СТ-проекций;
используют множество интерполированных изображений СТ-проекций для реконструкции (S550) изображения корректированной версии изображения среза; и
выводят (S560) корректированную версию изображения среза.
9. Способ по п. 8, в котором множество интерполированных изображений СТ-проекций формируют путем замены в каждом из множества изображений СТ-проекций информацию изображения внутри соответствующего отпечатка в этом изображении СТ-проекции на информацию интерполированного изображения.
10. Способ по п. 8, в котором этап реконструирования (S570), применяемый к новому изображению среза, включает в себя этап, на котором отслеживают соответствующие отпечатки во множестве изображений СТ-проекций, чтобы таким образом отследить движение элемента, и используют отслеживаемое движение на этапе реконструирования (S570) для нового среза, осуществляя таким образом компенсацию движения.
11. Медицинская система обработки изображений для коррекции изображения поперечного среза для артефакта в изображении, причем изображение среза основано на изображениях компьютерно-томографических (СТ) проекций объекта, полученных в разное время и под разными углами относительно объекта, когда элемент объекта с высоким значением коэффициента ослабления находится в движении, причем артефакт вызван упомянутым движением, причем система содержит:
устройство по любому из пп. 1-7;
систему базы данных, содержащую изображение среза и изображения СТ-проекций.
12. Машиночитаемый носитель, на котором сохранен компьютерный программный элемент для управления устройством по одному из пп. 1-7, который при его исполнении обрабатывающим блоком выполнен с возможностью выполнения этапов способа по пп. 8-10.
RU2014121096/28A 2011-10-24 2012-10-22 Двухпроходная коррекция металлического артефакта с компенсацией движения для изображений срезов компьютерной томографии RU2605519C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161550509P 2011-10-24 2011-10-24
US61/550,509 2011-10-24
PCT/IB2012/055799 WO2013061239A2 (en) 2011-10-24 2012-10-22 Motion compensated second pass metal artifact correction for ct slice images

Publications (2)

Publication Number Publication Date
RU2014121096A RU2014121096A (ru) 2015-12-10
RU2605519C2 true RU2605519C2 (ru) 2016-12-20

Family

ID=47278919

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2014121096/28A RU2605519C2 (ru) 2011-10-24 2012-10-22 Двухпроходная коррекция металлического артефакта с компенсацией движения для изображений срезов компьютерной томографии

Country Status (8)

Country Link
US (1) US9275454B2 (ru)
EP (1) EP2747655B1 (ru)
JP (1) JP6144268B2 (ru)
CN (1) CN103889329B (ru)
BR (1) BR112014009434A2 (ru)
IN (1) IN2014CN02779A (ru)
RU (1) RU2605519C2 (ru)
WO (1) WO2013061239A2 (ru)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10679104B2 (en) 2012-10-26 2020-06-09 Varian Medical Systems, Inc. Forward projection cylinder volume method
US9672600B2 (en) * 2012-11-19 2017-06-06 Carestream Health, Inc. Clavicle suppression in radiographic images
DE102013202313A1 (de) * 2013-02-13 2014-08-14 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Korrektur von Bewegungsartefakten bei einem computertomographischen Bild
JP6886606B2 (ja) * 2015-03-03 2021-06-16 株式会社ニコン 測定処理装置、x線検査装置、測定処理方法、測定処理プログラム、制御装置、および構造物の製造方法
JP6554353B2 (ja) * 2015-07-27 2019-07-31 株式会社日立製作所 X線撮像装置及び透視撮影制御方法
CN105374014B (zh) * 2015-08-25 2018-01-12 上海联影医疗科技有限公司 图像校正方法及装置、医学图像生成方法及装置
GB2542764A (en) * 2015-09-23 2017-04-05 Pathxl Ltd Image processing method and apparatus for normalisation and artefact correction
CN108369745B (zh) 2015-12-15 2023-11-07 皇家飞利浦有限公司 条纹伪影预测
JP6780948B2 (ja) * 2016-03-28 2020-11-04 キヤノンメディカルシステムズ株式会社 X線ct装置
EP3500171B1 (en) 2016-08-22 2020-06-03 Koninklijke Philips N.V. Model regularized motion compensated medical image reconstruction
KR101874950B1 (ko) * 2016-11-23 2018-07-05 연세대학교 산학협력단 재귀적 액티브 컨투어를 이용한 3차원 콘빔 ct 영상의 세그먼테이션 방법 및 ct 영상 장치
KR101870890B1 (ko) * 2016-11-23 2018-06-25 연세대학교 산학협력단 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치
US10255694B2 (en) * 2016-12-21 2019-04-09 Uih America, Inc. Methods and systems for emission computed tomography image reconstruction
US10360724B2 (en) 2016-12-21 2019-07-23 Uih America, Inc. Methods and systems for emission computed tomography image reconstruction
US10966670B2 (en) * 2017-08-07 2021-04-06 Carestream Health, Inc. Imaging system and method for dual-energy and computed tomography
CN109431528B (zh) * 2017-11-08 2022-05-31 上海联影医疗科技股份有限公司 校正投影图像的系统和方法
WO2019090541A1 (en) * 2017-11-08 2019-05-16 Shenzhen United Imaging Healthcare Co., Ltd. Systems and methods for correcting projection images in computed tomography image reconstruction
KR102489406B1 (ko) * 2018-02-27 2023-01-18 기초과학연구원 전산화 단층촬영 영상의 빔 경화 인공물 보정 방법 및 그 방법을 이용한 전산화 단층촬영 장치
KR20190103816A (ko) * 2018-02-28 2019-09-05 삼성전자주식회사 컴퓨터 단층 촬영 영상을 보정하는 방법 및 장치
CN112368715B (zh) * 2018-05-15 2024-08-23 蒙纳士大学 用于磁共振成像的运动校正的方法和系统
CN110660064A (zh) * 2018-06-29 2020-01-07 通用电气公司 心脏ct图像的处理方法和装置,非暂态计算机可读存储介质
JP7199643B2 (ja) 2019-05-10 2023-01-06 日本精工株式会社 車両用操向装置
JP7317639B2 (ja) * 2019-09-05 2023-07-31 富士フイルムヘルスケア株式会社 放射線画像処理システム及び画像処理方法
CN112233158B (zh) * 2020-10-14 2022-02-15 俐玛精密测量技术(苏州)有限公司 微纳米ct投影图像的二次投影配准方法
CN117830456B (zh) * 2024-03-04 2024-05-28 中国科学技术大学 用于校正图像金属伪影的方法、装置及电子设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6721387B1 (en) * 2001-06-13 2004-04-13 Analogic Corporation Method of and system for reducing metal artifacts in images generated by x-ray scanning devices
US20080247624A1 (en) * 2007-04-04 2008-10-09 Siemens Aktiengesellschaft Method for determining a three-dimensional reconstruction of an examination object
US20100183214A1 (en) * 2009-01-20 2010-07-22 Mccollough Cynthia H System and Method for Highly Attenuating Material Artifact Reduction in X-Ray Computed Tomography
US20110007956A1 (en) * 2009-07-07 2011-01-13 Esther Meyer Sinogram processing to reduce metal artifacts in computed tomography

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0819533A (ja) * 1994-07-05 1996-01-23 Hitachi Medical Corp X線ct装置
US6263096B1 (en) * 1999-06-23 2001-07-17 The Board Of Trustees Of The University Of Illinois Multilevel domain decomposition method for fast reprojection of images
US8275091B2 (en) * 2002-07-23 2012-09-25 Rapiscan Systems, Inc. Compact mobile cargo scanning system
JP4584550B2 (ja) * 2003-06-05 2010-11-24 株式会社日立メディコ X線計測装置
US20090074278A1 (en) * 2004-10-12 2009-03-19 Universite Laval Method and apparatus for metal artifact reduction in computed tomography
CN100361632C (zh) * 2005-03-22 2008-01-16 东软飞利浦医疗设备系统有限责任公司 一种自动去除黑心伪影的x-射线计算机层析成像机
DE102005037367B3 (de) * 2005-08-08 2007-04-05 Siemens Ag Verfahren für eine Röntgeneinrichtung
DE102005049862A1 (de) * 2005-10-18 2007-04-26 Siemens Ag Verfahren und Vorrichtung zur Bewegungskorrektur bei der Bildgebung des Herzens
AU2006331312A1 (en) * 2005-11-17 2007-07-05 Brain Research Institute Pty Ltd Apparatus and method for detection and monitoring of electrical activity and motion in the presence of a magnetic field
US8184883B2 (en) 2005-11-24 2012-05-22 Koninklijke Philips Electronics N.V. Motion compensated CT reconstruction of high contrast objects
US7746976B2 (en) * 2005-12-30 2010-06-29 Carestream Health, Inc. Bone mineral density assessment using mammography system
JP2010501856A (ja) * 2006-08-22 2010-01-21 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 動きアーティファクト画像についてのアーティファクト補正
WO2008084352A1 (en) * 2007-01-04 2008-07-17 Koninklijke Philips Electronics N. V. Apparatus, method and computer program for producing a corrected image of a region of interest from acquired projection data
KR100923098B1 (ko) 2008-01-15 2009-10-22 (주)이우테크놀로지 엑스 레이 씨티 촬영 영상의 메탈 아티팩트를 제거하는방법
JP5303154B2 (ja) * 2008-02-20 2013-10-02 株式会社日立メディコ X線ct装置
JP5243160B2 (ja) * 2008-09-16 2013-07-24 株式会社日立メディコ X線ct装置
CN102307535B (zh) 2009-01-05 2014-09-24 皇家飞利浦电子股份有限公司 用于电磁跟踪系统的动态金属畸变补偿的系统和方法
US8503750B2 (en) * 2009-10-06 2013-08-06 General Electric Company Method and apparatus for reduction of metal artifacts in CT images

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6721387B1 (en) * 2001-06-13 2004-04-13 Analogic Corporation Method of and system for reducing metal artifacts in images generated by x-ray scanning devices
US20080247624A1 (en) * 2007-04-04 2008-10-09 Siemens Aktiengesellschaft Method for determining a three-dimensional reconstruction of an examination object
US20100183214A1 (en) * 2009-01-20 2010-07-22 Mccollough Cynthia H System and Method for Highly Attenuating Material Artifact Reduction in X-Ray Computed Tomography
US20110007956A1 (en) * 2009-07-07 2011-01-13 Esther Meyer Sinogram processing to reduce metal artifacts in computed tomography

Also Published As

Publication number Publication date
WO2013061239A2 (en) 2013-05-02
EP2747655A2 (en) 2014-07-02
BR112014009434A2 (pt) 2017-04-11
CN103889329A (zh) 2014-06-25
JP2014530709A (ja) 2014-11-20
IN2014CN02779A (ru) 2015-07-03
WO2013061239A3 (en) 2013-11-07
US9275454B2 (en) 2016-03-01
EP2747655B1 (en) 2019-02-06
RU2014121096A (ru) 2015-12-10
JP6144268B2 (ja) 2017-06-07
US20140270450A1 (en) 2014-09-18
CN103889329B (zh) 2016-10-26

Similar Documents

Publication Publication Date Title
RU2605519C2 (ru) Двухпроходная коррекция металлического артефакта с компенсацией движения для изображений срезов компьютерной томографии
US20210056688A1 (en) Using deep learning to reduce metal artifacts
Meilinger et al. Metal artifact reduction in cone beam computed tomography using forward projected reconstruction information
US6810102B2 (en) Methods and apparatus for truncation compensation
US9934597B2 (en) Metal artifacts reduction in cone beam reconstruction
US7378660B2 (en) Computer program, method, and system for hybrid CT attenuation correction
CN111915696B (zh) 三维图像数据辅助的低剂量扫描数据重建方法及电子介质
US20110007956A1 (en) Sinogram processing to reduce metal artifacts in computed tomography
US11580678B2 (en) Systems and methods for interpolation with resolution preservation
CN117522747A (zh) 一种用于ct图像的金属伪影校正方法与系统
Kim et al. Cardiac motion correction for helical CT scan with an ordinary pitch
CN113226184B (zh) 在x射线牙科体积断层扫描中金属伪影减少的方法
US20220071578A1 (en) Improved method of acquiring a radiographic scan of a region-of-interest in a metal containing object
CN115564858A (zh) 一种金属伪影的校正方法、系统和x射线图像重建方法
US8526757B2 (en) Imaging system and imaging method for imaging a region of interest
CN116457823A (zh) 锥形束计算机断层摄影中伪影的减少
KR20180003400A (ko) 의료 영상 처리 장치 및 방법
Agrawal et al. Metal artifact reduction in cone-beam extremity images using gated convolutions
CN112258596A (zh) 图像生成方法、装置、控制台设备及ct系统
El Naqa et al. Automated breathing motion tracking for 4D computed tomography
EP3776478B1 (en) Tomographic x-ray image reconstruction
US11666299B2 (en) Controlling a medical X-ray device
CN116635888A (zh) 图像分割方法及装置、图像引导系统、放射治疗系统
CN117252943A (zh) 一种基于概率扩散模型(ddmp)的金属伪影去除方法
CN111415394A (zh) 骨硬化伪影校正方法、装置、计算机设备和可读存储介质

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20171023