RU2582055C2 - Совмещение данных изображения для динамической перфузионной компьютерной томографии - Google Patents
Совмещение данных изображения для динамической перфузионной компьютерной томографии Download PDFInfo
- Publication number
- RU2582055C2 RU2582055C2 RU2012152513/08A RU2012152513A RU2582055C2 RU 2582055 C2 RU2582055 C2 RU 2582055C2 RU 2012152513/08 A RU2012152513/08 A RU 2012152513/08A RU 2012152513 A RU2012152513 A RU 2012152513A RU 2582055 C2 RU2582055 C2 RU 2582055C2
- Authority
- RU
- Russia
- Prior art keywords
- image data
- contrast
- dimensional
- dimensional image
- similarity
- Prior art date
Links
- 230000010412 perfusion Effects 0.000 title claims description 22
- 238000002591 computed tomography Methods 0.000 title description 11
- 238000000034 method Methods 0.000 claims abstract description 28
- 239000002872 contrast media Substances 0.000 claims description 15
- 230000000747 cardiac effect Effects 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 2
- 238000011478 gradient descent method Methods 0.000 abstract description 2
- 239000000126 substance Substances 0.000 abstract 1
- 230000005855 radiation Effects 0.000 description 12
- 238000003384 imaging method Methods 0.000 description 11
- 210000001519 tissue Anatomy 0.000 description 8
- 238000013459 approach Methods 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 5
- 230000002107 myocardial effect Effects 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 210000003484 anatomy Anatomy 0.000 description 4
- 238000011084 recovery Methods 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 239000012634 fragment Substances 0.000 description 3
- 230000000241 respiratory effect Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000002386 leaching Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 208000001778 Coronary Occlusion Diseases 0.000 description 1
- 206010011086 Coronary artery occlusion Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 229940039231 contrast media Drugs 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 210000002751 lymph Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/507—Apparatus 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 determination of haemodynamic parameters, e.g. perfusion CT
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10076—4D tomography; Time-sequential 3D tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Dentistry (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
Изобретение относится к совмещению данных изображения. Техническим результатом является повышение точности определения движения подвижной структуры. Способ совмещения набора данных 4-мерного изображения с контрастированием содержит этап, на котором: совмещают данные изображения, соответствующие множеству разных периодов времени, с опорным изображением из одного из периодов времени, данные изображения совмещаются с использованием итерационного алгоритма совмещения на основе градиентов и показателя подобия, алгоритм совмещения основывается на алгоритме оптимизации методом стохастического градиентного спуска, и показатель подобия включает в себя показатель подобия, основывающийся на нормированной кросс-корреляции с нулевым средним. 2 н. и 9 з.п. ф-лы, 7 ил.
Description
Нижеследующее в целом относится к совмещению данных изображения и описывается с конкретным применением к компьютерной томографии (КТ); однако в этом документе также предполагаются другие применения по визуализации.
Компьютерный томограф (КТ) обычно включает в себя неподвижный гентри и вращающийся гентри. Вращающийся гентри поддерживается с возможностью вращения неподвижным гентри и перемещает рентгеновскую трубку. Детектор располагается напротив рентгеновской трубки, поперек области обследования. Вращающийся гентри вращается вокруг области обследования относительно оси z, и излучение, выпущенное рентгеновской трубкой, пересекает область обследования и анатомию пациента в ней и обнаруживается детектором. Детектор формирует проекционные данные, указывающие облученную анатомию пациента. Блок восстановления восстанавливает проекционные данные и сформированные восстановленные данные изображения, указывающие облученную анатомию пациента.
КТ, а также другие устройства получения изображений могут использоваться для визуализации перфузии. Как правило, визуализация перфузии является подходом к визуализации для сбора информации, соответствующей протеканию жидкости (например, крови, лимфы и т.д.) через анатомическую ткань, и определению количества информации для упрощения идентификации состояния здоровья этой ткани или другой ткани. Например, визуализация перфузии миокарда может использоваться для измерения снабжения кислородом мышечной ткани в сердце. По существу, сердце пациента может изображаться в трех измерениях (3D) до и во время впрыскивания рентгеноконтрастного вещества, и отличия в восстановленных изображениях могут приписываться контрастному веществу, которое промывается в ткани. Мышечную ткань затем можно классифицировать в соответствии со степенью поглощения контрастного вещества, и можно обнаружить патологии, например инфаркт миокарда или коронарную окклюзию.
Проспективная ЭКГ-синхронизированная КТ использовалась ранее для визуализации перфузии миокарда. По этой методике, хотя некоторая или конкретная фаза сердца (например, систола, диастола и т.д.) выбирается для проспективной синхронизации, состояния движения сердца в разных периодах времени сердечного цикла могут отличаться. Это обычно обусловлено тем, что синхронизация с ЭКГ опирается на электрический сигнал, а не на фактическое механическое движение, и что изменение длины цикла сердца не учитывает проспективную синхронизацию. В результате периоды времени динамических измерений могут не выравниваться по вокселям, и такое несовпадение может вносить артефакты (искажения) в кривые временного поглощения контрастного вещества, что может приводить к ошибочным результатам анализа.
Аспекты настоящей заявки решают вышеуказанные и другие вопросы.
В соответствии с одним аспектом способ включает в себя совмещение набора данных 4-мерного изображения с контрастированием. Способ включает в себя совмещение данных изображения, соответствующих множеству разных периодов времени, с опорным изображением из одного из периодов времени. Набор данных 4-мерного изображения с контрастированием включает в себя данные изображения одного и того же исследуемого объема, полученные в разные периоды времени с изменяющимся контрастированием, причем исследуемый объем включает в себя подвижную структуру, и разные периоды времени соответствуют заранее установленной исследуемой фазе движения в разных циклах движения подвижной структуры.
В соответствии с другим аспектом система включает в себя компонент деформирования, который деформирует данные изображения в наборе данных 4-мерного изображения, чтобы совместить данные изображения с данными опорного изображения из набора данных 4-мерного изображения на основе алгоритма градиентного спуска. Данные изображения и данные опорного изображения соответствуют одному и тому же исследуемому объему и разным периодам времени. Система дополнительно включает в себя определитель подобия, который определяет значение, указывающее подобие между данными деформированного изображения и данными опорного изображения.
В соответствии с другим аспектом машиночитаемый носитель информации, кодированный командами, которые при выполнении процессором компьютера побуждают процессор: эластично совместить данные 4-мерного проспективного ЭКГ-синхронизированного изображения с контрастированием на основе алгоритма оптимизации методом стохастического градиентного спуска с адаптивным предсказанием величины шага в сочетании с показателем подобия, основывающийся на нормированной кросс-корреляции с нулевым средним.
Изобретение может принимать форму различных компонентов и схем компонентов, и различных этапов и схем этапов. Чертежи предназначаются только для целей иллюстрации предпочтительных вариантов осуществления и не должны толковаться как ограничивающие изобретение.
Фиг. 1 иллюстрирует систему визуализации применительно к компоненту совмещения данных изображения.
Фиг. 2 иллюстрирует примерный компонент совмещения.
Фиг. 3 иллюстрирует примерную пространственную фильтрацию данных 4-мерного изображения.
Фиг. 4 иллюстрирует примерную временную фильтрацию данных 4-мерного изображения.
Фиг. 5 и 6 иллюстрируют примерное совмещение.
Фиг. 7 иллюстрирует примерный способ.
Фиг. 1 иллюстрирует систему 100 визуализации, например компьютерный томограф (КТ). Система 100 визуализации включает в себя неподвижный гентри 102 и вращающийся гентри 104, который поддерживается с возможностью вращения неподвижным гентри 102. Вращающийся гентри 104 вращается вокруг области 106 обследования относительно продольной оси или оси z. В проиллюстрированном варианте осуществления вращающийся гентри 104 конфигурируется для достаточно быстрого вращения для сканирования бьющегося сердца человека, например, для сканирования всего сердца в течение одного сердечного цикла. В одном примере проиллюстрированный вращающийся гентри 104 конфигурируется для поворота по меньшей мере за двести семьдесят миллисекунд (270 мс). Вращающийся гентри 104 также может поворачиваться за меньшее и большее время обращения.
Источник 108 излучения, например рентгеновская трубка, поддерживается вращающимся гентри 104 и вращается вместе с вращающимся гентри 104, и испускает излучение. Коллиматор источника коллимирует излучение, чтобы образовать конический, веерный, клиновидный или иным образом оформленный пучок излучения, который пересекает область 106 обследования. Чувствительная к излучению детекторная матрица 110, расположенная напротив источника 108, обнаруживает излучение, которое пересекает область 106 обследования, и формирует указывающие его проекционные данные. В проиллюстрированном варианте осуществления чувствительная к излучению детекторная матрица 110 достаточно большая, чтобы обнаружить излучение, пересекающее все сердце пациента за один раз. Для выполнения этого в одном варианте осуществления детекторная матрица 110 содержит по меньшей мере двести пятьдесят шесть (256) слоев или рядов детекторов. В другом варианте осуществления детекторная матрица 110 может иметь больше или меньше слоев или рядов.
Контрольный прибор или датчик движения 114 конфигурируется для восприятия состояния движения по меньшей мере заранее установленного исследуемого объема подвижного объекта или предмета в области 106 обследования и формирования указывающего его сигнала. Датчик 114 движения может включать в себя датчик сердечного, и/или дыхательного, и/или другого движения. Универсальная вычислительная система служит в качестве пульта 116 оператора, который включает в себя человекочитаемые устройства вывода, например дисплей и/или принтер, и устройства ввода, например клавиатуру и/или мышь. Программное обеспечение, постоянно находящееся в пульте 116, позволяет оператору управлять работой системы 100, например, позволяя оператору выбирать протокол сканирования (например, синхронизированный с ЭКГ или иначе протокол с контрастированием, который синхронизирует сканирование на основе сигнала от контрольного прибора 114 движения), инициировать сканирование и т.д. Впрыскиватель 112 конфигурируется для впрыскивания контрастного вещества (веществ), например, для процедуры визуализации с контрастированием. Проиллюстрированный впрыскиватель 112 управляется пультом 116. В другом примере контрастное вещество вводится вручную.
Блок 118 восстановления восстанавливает проекционные данные и формирует объемные данные изображения, указывающие область 106 обследования. Блок 118 восстановления может применять алгоритмы восстановления, которые охватывают различные угловые диапазоны, например, от ста восьмидесяти градусов (180°) плюс угол веерного пучка до трехсот шестидесяти градусов (360°) плюс угол веерного пучка, и/или другие диапазоны. Набор четырехмерных (4D) данных можно восстановить на основе данных, полученных в разные периоды времени. В одном примере каждый период времени может соответствовать разному циклу движения (например, сердечного, дыхательного и т.д.) или субфрагменту разных циклов движения. Субфрагмент может соответствовать заранее установленной и/или выбранной (например, сердечной, дыхательной и т.д.) исследуемой фазе (фазам) в каждом цикле. При визуализации сердца фаза может быть систолой, диастолой или другой фазой. В случае данных изображения с контрастированием данные изображения, соответствующие разному циклу движения/периодам времени, могут иметь разное контрастирование, представляющее поглощение и вымывание контрастного вещества в ткани во время сканирования.
Опора 124, например кушетка, поддерживает объект или предмет в области 106 обследования. Опора 124 может использоваться для различного расположения объекта или предмета в области 106 обследования относительно осей х, у и/или z до, во время и/или после сканирования. В одном примере опора 124 используется для расположения исследуемого объема, например сердца или другой ткани, в области 106 обследования, и затем весь исследуемый объем (или его существенная часть) сканируется с опорой 124 в том же положении. Примером такого сканирования является сканирование сердца, например проспективное ЭКГ-синхронизированное сканирование сердца с контрастированием. Для такого сканирования коллиматор источника коллимирует выпущенное излучение, чтобы создать пучок излучения с подходящим углом раствора пучка (направления х/у и/или z) для охвата исследуемого объема, детекторная матрица имеет подходящее количество слоев или рядов детекторов по оси z, чтобы обнаружить излучение, пересекающее исследуемый объем, и сердце сканируется во время выбранной фазы (фаз) каждого сердечного цикла или для некоторого выбранного подмножества сердечных циклов, например менее каждого из них.
Компонент 120 совмещения совмещает данные 4-мерного изображения или их заранее установленный субфрагмент. Как отмечалось выше, данные 4-мерного изображения могут включать в себя данные изображения, соответствующие разным циклам движения, и соответственно периодам времени. При исследованиях с контрастированием данные изображения, соответствующей одной и той же фазе движения в разных циклах движения и периодах времени, могут иметь разное контрастирование, отражающее поглощение и вымывание контрастного вещества в течение процедуры визуализации и времени. Как подробнее описано ниже, в одном примере алгоритм совмещения включает в себя итерационное эластичное совмещение на основе градиента с показателем подобия, который упрощает определение того, когда данные изображения совмещаются соответствующим образом. С помощью данных с контрастированием такой алгоритм может учитывать изменения в концентрации контрастного вещества в разные периоды времени и/или несогласованные состояния движения. Такой алгоритм может уменьшить артефакты из-за геометрического несовпадения, посредством этого предоставляя большую количественную точность для различных процедур визуализации, например визуализации КТ перфузии миокарда.
Анализатор 122 может использоваться для анализа совмещенных данных изображения. При перфузионных сканированиях анализатор 122 может использоваться для формирования перфузионных карт и/или перфузионных параметров, таких как время максимального накопления, площадь под кривой (то есть объем крови, циркулирующей через ткань), локальный максимум интенсивности или усиление, средний возрастающий наклон, максимальный подъем, время до максимального подъема и т.д. Такая информация может формироваться по вокселям, по исследуемым областям или на другой основе, например, по временной оси. Сформированные перфузионные карты и/или перфузионные параметры могут представляться по-разному, например, в виде перфузионных карт с цветной маркировкой, графиков и т.д., сохраняться в запоминающем устройстве, записываться на пленку, передаваться для обработки и т.д.
Нужно учитывать, что компонент 120 совмещения и/или анализатор 122 могут быть частью системы 100 (как показано) или быть удаленными от нее, например, в вычислительной системе, такой как рабочая станция или т.п. В любом случае один или несколько процессоров могут выполнять машиночитаемые команды, кодированные и/или воплощенные на локальном или удаленном машиночитаемом носителе информации, например запоминающем устройстве, чтобы реализовать компонент 120 совмещения и/или анализатор 122.
Фиг. 2 иллюстрирует пример компонента 120 совмещения.
Проиллюстрированный компонент 120 совмещения включает в себя фильтр 202, который применяется для фильтрации данных 4-мерного изображения, которые включают в себя данные изображения, полученные в разные моменты времени, или данные изображения, соответствующие разным периодам времени. Проиллюстрированный фильтр 202 конфигурируется для пространственной фильтрации данных изображения, и затем временной фильтрации данных изображения. Фиг. 3 показывает пример пространственной фильтрации данных 4-мерного изображения, соответствующих проспективному ЭКГ-синхронизированному перфузионному сканированию сердца с контрастированием, по пространственным координатам (х, у, и z), а фиг. 4 показывает пример временной фильтрации данных 4-мерного изображения по временным координатам (время). Возвращаясь к фиг. 2, такая многомерная фильтрация изображения может увеличить отношение сигнал-шум (SNR) у данных изображения наряду с сохранением пространственной и временной четкости данных изображения. Увеличение SNR (то есть уменьшение шума изображения) может привести к более точному количественному анализу перфузии миокарда. В другом варианте осуществления фильтр 202 отсутствует, или используется другой фильтр. Пример подходящего фильтра обсуждается в статье Н. Bruder et al., "Spatio-temporal filtration of dynamic CT data using diffusion filter", in Proc. of SPIE Med. Imag. Conf., vol. 7258, 2009.
Идентификатор 204 опорного изображения идентифицирует данные опорного изображения в наборе 4-мерных данных. По отношению к данным 4-мерного изображения с контрастированием идентификатор 204 опорного изображения может выполнять это путем идентификации периода времени или данных изображения, имеющих наибольшее (максимальное) контрастирование. Это может выполняться автоматически, например, из кривых максимального увеличения или иным образом. Как правило, данные изображения, соответствующие разным сканированным сердечным циклам, будут иметь разные величины контрастирования, когда контрастное вещество поглощается и вымывается, и данные изображения с наибольшим поглощением контрастного вещества соответствуют времени, когда обе камеры сердца полностью заполнены контрастным веществом и, скорее всего, покажут большую четкость на границе относительно других данных изображения. Этот период времени можно идентифицировать автоматически с помощью компонента 120 совмещения, например путем отыскания периода времени с максимальной энергией изображения. В другом примере идентификатор 204 опорного изображения идентифицирует данные опорного изображения на основе сигнала, указывающего выбранные пользователем данные изображения, и/или иным образом.
Компонент 206 деформирования деформирует или преобразует данные 4-мерного изображения на основе данных опорного изображения. В проиллюстрированном варианте осуществления компонент 206 деформирования применяет эластичный или нежесткий алгоритм. В другом варианте осуществления дополнительно или в качестве альтернативы применяется жесткий алгоритм. Компонент 206 деформирования может деформировать данные изображения один или несколько раз, например, путем деформирования ранее деформированных данных изображения. Например, когда определяется, что данные деформированного изображения и опорного изображения не достаточно похожи структурно, компонент 206 деформирования может деформировать данные деформированного изображения. Определитель 208 подобия определяет значение, указывающее подобие между данными деформированного изображения и данными опорного изображения. Оптимизатор 210 облегчает минимизацию критерия подобия.
Компонент 212 принятия решений определяет, когда достигнуто подходящее совмещение. Если оно не достигнуто, то компонент 120 совмещения деформирует данные деформированного изображения, определитель 208 подобия определяет другое значение, и компонент 212 принятия решений на основе нового значения определяет, достигнуто ли подходящее совмещение. Эти этапы могут повторяться, пока совмещение определяется как приемлемое. В одном примере компонент 212 принятия решений определяет значение разности между этим значением и значением, определенным для предыдущего деформирования данных изображения, и сравнивает значение разности с заранее установленной пороговой величиной. В другом примере компонент 212 принятия решений сравнивает значение с заранее установленным (относительным или абсолютным) пороговым значением. В еще одном примере компонент 212 принятия решений определяет, когда выполнено заранее установленное количество итераций над данными изображения. Другие критерии могут использовать дополнительно или в качестве альтернативы для определения, когда достигается подходящее совмещение.
Фиг. 5 и 6 иллюстрируют примерное совмещение с помощью компонента 120 совмещения. В этом примере компонент 120 совмещения определяет поле деформации, которое удовлетворяет Уравнению 1.
Уравнение 1
fw(x)=ft(g(x))=fr,
где fw представляет данные деформированного изображения, ft представляет данные изображения, совмещаемые с данными опорного изображения, g представляет поле деформации, a fr представляет данные опорного изображения. Компонент 120 совмещения может использовать все или подмножество значений интенсивности вокселей у совмещаемых изображений.
На фиг. 5 данные 502 изображения совмещаются с данными 504 опорного изображения. В этом примере компонент 206 деформирования включает в себя интерполятор, который деформирует данные 502 изображения, формируя деформированные данные 508 изображения. Подходящий алгоритм интерполяции включает в себя кубический би-сплайн или другой алгоритм. Компонент 208 подобия определяет значение подобия между данными 508 деформированного изображения и данными 504 опорного изображения. В проиллюстрированном варианте осуществления критерий подобия включает в себя нормированную кросс-корреляцию с нулевым средним между изображениями fr и ft. Оптимизатор 210 пытается найти наилучшую деформацию g, которая минимизирует критерий подобия. Подходящий оптимизатор включает в себя адаптивную оптимизацию методом стохастического градиентного спуска. Этот критерий хорошо подходит для обработки меняющейся концентрации контрастного вещества в данных 4-мерного изображения. Соответствие между дискретными и непрерывными версиями изображений и полем g 510 деформации может устанавливаться с использованием кубических би-сплайнов. Как отмечалось в этом документе, совмещение может включать в себя одну или несколько итераций 512. Как показано на фиг. 6, данные изображения из одного или нескольких периодов времени 600 совмещаются на основе и/или с данными 504 опорного изображения.
В проиллюстрированном варианте осуществления алгоритм адаптивной оптимизации методом стохастического градиентного спуска применяется вместе с адаптивным предсказанием величины шага. Подходящий алгоритм оптимизации с предсказанием величины шага обсуждается в S. Klein, J. P. W. Pluim, M. Staring and M.
A. Viergever, "Adaptive stochastic gradient descent optimisation for image registration", International Journal of Computer Vision, 81(3), pp. 227-239, 2009. Такой алгоритм может уменьшить время вычисления в расчете на итерацию, не влияя на скорость сходимости, окончательную точность или надежность. Этот алгоритм аналогичен детерминированному градиентному спуску с отличием в том, что градиент критерия подобия заменяется приближением. Стохастическое приближение производной критерия подобия может определяться с использованием нового, случайно выбранного подмножества вокселей в каждой итерации процесса оптимизации. Таким образом, избегают сдвига в ошибке приближения.
Скорость и точность совмещения могут зависеть от качества градиентного приближения, полученного путем случайной выборки. Подмножество случайных местоположений вокселей может определяться с использованием либо равномерной выборки, либо неравномерной выборки, что может улучшить градиентные приближения. Границы изображения могут сильно влиять на оценки совмещения на основе интенсивности. По существу, выборочное распределение, которое подчеркивает границы изображения, может использоваться для улучшения градиентных приближений. В одном примере двоичная маска используется для случайной выборки только вокселей, которые имеют наибольшие величины градиента изображения. Другой подходящий подход для случайной выборки описывается в Bhagalia R, Fessler JA, Kim В, "Accelerated nonrigid intensity-based image registration using importance sampling", IEEE Trans Med Imaging, 2009 Aug; 28 (8): 1208-16.
Подход с переменной разрешающей способностью может применяться для повышения надежности и эффективности алгоритма совмещения. С помощью такого алгоритма сначала может выполняться совмещение с использованием грубых субдискретизированных изображений и определения поля деформации с уменьшенным количеством степеней свободы. Потом результаты могут распространяться на следующий более мелкий уровень. Эта итерационная процедура может поочередно расширять сетки контрольных точек би-сплайна у изображений и у поля деформации, пока не достигнут самый мелкий уровень. Подходящий подход с переменной разрешающей способностью обсуждается в М. Unser, А. Aldroubi, and М. Eden, "The L2-polynomial spline pyramid", IEEE Trans. Pattern Anal. Mach. Intell. 15, 364-379, 1993. Этот подход использует 3-мерные операторы уменьшения/расширения би-сплайна с коэффициентом два (2) для построения пирамиды или набор постепенно уменьшенных версий исходных изображений и поля деформации, который оптимален.
В случаях, в которых существует высокая степень свободы при совмещении, в сущности, может иметь место некорректная задача, которая могла бы привести к нереальному свертыванию полей деформации при отсутствии подходящих ограничений. По существу, различные ограничения могут применяться при оценивании поля g деформации. Например, в одном примере одно ограничение может состоять в том, что поле g деформации является обратимым, когда движение органа и ткани человека является обратимым. Один подходящий подход для 3-мерных совмещений можно найти в S. Y. Chun and J. A. Fessler, "A simple regularizer for b-spline nonrigid image registration that encourages local invertibility", IEEE J. Sel. Top. Signal Process. 3, 159-169, 2009. Штрафная функция, раскрытая в работе Fessler, поощряет положительные якобианы путем ограничения отличий двух соседних коэффициентов деформации в направлении х, у, z. Путем ограничения только отличий вместо коэффициентов даже большие деформации g с градиентами в тех пределах включаются в пространство поиска решений.
Фиг. 7 иллюстрирует способ. Для этого примера компонент 120 совмещения совмещает данные 4-мерного изображения, полученные посредством проспективного ЭКГ-синхронизированного перфузионного сканирования сердца с контрастированием. Такое сканирование синхронизируется на основе заранее установленной фазы (фаз) сердца, например систолической фазы, диастолической фазы и/или другой сердечной фазы в каждом сердечном цикле или сердечном сокращении (то есть в каждом интервале R-R). Фаза (фазы) может выбираться оператором сканера, идентифицироваться по выбранному протоколу визуализации и/или определяться иным образом. Кроме того, сканирование может синхронизироваться так, что фаза (фазы) сердца сканируется в каждом цикле сердца, каждом втором цикле сердца или некотором другом количестве циклов сердца меньше каждого цикла сердца. Сканирование может инициироваться до, во время или вскоре после введения контрастного вещества.
На этапе 702 получаются данные 4-мерного перфузионного изображения. На этапе 704 к набору 4-мерных данных применяется динамическая диффузионная фильтрация, которая описана применительно к фиг. 2 или иным образом. На этапе 706 фильтрованное изображение совмещается, как описано в этом документе. Например, в одном примере это включает в себя эластичное совмещение данных изображения с использованием данных опорного изображения, используя способ оптимизации методом стохастического градиентного спуска вместе с адаптивным предсказанием величины шага в сочетании с показателем подобия, основывающимся на нормированной кросс-корреляции с нулевым средним. На этапе 708 обрабатываются результирующие совмещенные данные изображения. Это включает в себя определение различных перфузионных карт и/или параметров. Используя этот способ, артефакты вследствие геометрическое несовпадения можно устранить или уменьшить при совмещении данных изображения с контрастированием у подвижного объекта или анатомической структуры. По существу, перфузионные применения, такие как визуализация КТ перфузии миокарда, могут предоставлять данные, хорошо подходящие для относительно высокой количественной точности.
Вышеописанные действия могут быть реализованы посредством машиночитаемых команд, которые при выполнении процессором (процессорами) компьютера побуждают процессор (процессоры) выполнять описанные в этом документе действия. В таком случае команды сохраняются на машиночитаемом носителе информации, например в запоминающем устройстве, ассоциированном и/или доступном иным образом для соответствующего компьютера.
Изобретение описано со ссылкой на предпочтительные варианты осуществления. После прочтения и осмысления предшествующего подробного описания могут возникнуть его модификации и изменения. Подразумевается, что изобретение должно создается как включающее в себя все такие модификации и изменения в той мере, в которой они входят в объем прилагаемой формулы изобретения либо ее эквивалентов.
Claims (11)
1. Способ совмещения набора данных 4-мерного изображения с контрастированием, полученного перфузионным сканированием с контрастированием посредством контрастного вещества, в котором набор данных 4-мерного изображения с контрастированием включает в себя данные изображения одного и того же исследуемого объема, полученные в разные периоды времени с изменяющимся контрастированием, при этом исследуемый объем включает в себя подвижную структуру, и разные периоды времени соответствуют заранее установленной исследуемой фазе движения в разных циклах движения подвижной структуры, причем способ содержит этап, на котором:
совмещают данные изображения, соответствующие множеству разных периодов времени, с опорным изображением из одного из периодов времени, при этом данные изображения совмещаются с использованием итерационного алгоритма совмещения на основе градиентов и показателя подобия, причем алгоритм совмещения основывается на алгоритме оптимизации методом стохастического градиентного спуска, и показатель подобия включает в себя показатель подобия, основывающийся на нормированной кросс-корреляции с нулевым средним.
совмещают данные изображения, соответствующие множеству разных периодов времени, с опорным изображением из одного из периодов времени, при этом данные изображения совмещаются с использованием итерационного алгоритма совмещения на основе градиентов и показателя подобия, причем алгоритм совмещения основывается на алгоритме оптимизации методом стохастического градиентного спуска, и показатель подобия включает в себя показатель подобия, основывающийся на нормированной кросс-корреляции с нулевым средним.
2. Способ по п. 1, в котором показатель подобия нормализованной взаимной корреляции с нулевым средним обрабатывает меняющуюся концентрацию контрастного вещества в наборе данных 4-мерного изображения с контрастированием.
3. Способ по п. 1 или 2, в котором совмещение данных изображения содержит:
деформацию данных изображения для конкретного периода времени, чтобы привести в соответствие с данными опорного изображения;
определение значения, указывающего подобие между данными деформированного изображения и данными опорного изображения; и
совмещение данных изображения с данными опорного изображения, когда значение удовлетворяет заранее установленным критериям, и деформацию данных деформированного изображения, когда значение не удовлетворяет заранее установленным критериям.
деформацию данных изображения для конкретного периода времени, чтобы привести в соответствие с данными опорного изображения;
определение значения, указывающего подобие между данными деформированного изображения и данными опорного изображения; и
совмещение данных изображения с данными опорного изображения, когда значение удовлетворяет заранее установленным критериям, и деформацию данных деформированного изображения, когда значение не удовлетворяет заранее установленным критериям.
4. Способ по п. 3, в котором заранее установленные критерии включают в себя пороговое значение.
5. Способ по п. 3, в котором заранее установленные критерии включают в себя максимальное количество итераций.
6. Способ по п. 1 или 2, дополнительно содержащий этап, на котором фильтруют данные 4-мерного изображения до совмещения данных 4-мерного изображения, причем фильтрация включает в себя пространственную фильтрацию набора данных 4-мерного изображения и временную фильтрацию данных 4-мерного изображения.
7. Способ по п. 1 или 2, дополнительно содержащий этап, на котором формируют по меньшей мере одно из параметрической карты и/или перфузионного параметра на основе набора совмещенных данных 4-мерного изображения.
8. Способ по п. 1 или 2, в котором набор данных 4-мерного изображения соответствует проспективному ЭКГ-синхронизированному перфузионному сканированию с контрастированием.
9. Способ по п. 1 или 2, в котором данные опорного изображения соответствуют данным изображения, имеющим заранее установленное целевое контрастирование.
10. Способ по п. 1 или 2, в котором набор данных 4-мерного изображения включает в себя данные изображения сердца и соответствует конкретной заранее установленной исследуемой сердечной фазе.
11. Машиночитаемый носитель информации, кодированный исполняемыми компьютером командами, которые при их исполнении процессором компьютера предписывают процессору:
эластично совмещать данные 4-мерного проспективного ЭКГ-синхронизированного изображения с контрастированием на основе алгоритма оптимизации методом стохастического градиентного спуска с адаптивным предсказанием величины шага в сочетании с показателем подобия, основывающимся на нормированной кросс-корреляции с нулевым средним.
эластично совмещать данные 4-мерного проспективного ЭКГ-синхронизированного изображения с контрастированием на основе алгоритма оптимизации методом стохастического градиентного спуска с адаптивным предсказанием величины шага в сочетании с показателем подобия, основывающимся на нормированной кросс-корреляции с нулевым средним.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US33187110P | 2010-05-06 | 2010-05-06 | |
| US61/331,871 | 2010-05-06 | ||
| PCT/IB2011/051579 WO2011138694A1 (en) | 2010-05-06 | 2011-04-12 | Image data registration for dynamic perfusion ct |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2012152513A RU2012152513A (ru) | 2014-06-20 |
| RU2582055C2 true RU2582055C2 (ru) | 2016-04-20 |
Family
ID=44310075
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2012152513/08A RU2582055C2 (ru) | 2010-05-06 | 2011-04-12 | Совмещение данных изображения для динамической перфузионной компьютерной томографии |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US9002089B2 (ru) |
| EP (1) | EP2567359B1 (ru) |
| JP (1) | JP5833637B2 (ru) |
| CN (1) | CN102884550B (ru) |
| RU (1) | RU2582055C2 (ru) |
| WO (1) | WO2011138694A1 (ru) |
Families Citing this family (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8643642B2 (en) | 2009-08-17 | 2014-02-04 | Mistretta Medical, Llc | System and method of time-resolved, three-dimensional angiography |
| US8654119B2 (en) | 2009-08-17 | 2014-02-18 | Mistretta Medical, Llc | System and method for four dimensional angiography and fluoroscopy |
| WO2011091300A2 (en) | 2010-01-24 | 2011-07-28 | Mistretta Medical, Llc | System and method for implementation of 4d time-energy subtraction computed tomography |
| US8655041B2 (en) * | 2010-09-23 | 2014-02-18 | Siemens Aktiengesellschaft | Automatic detection of contrast injection |
| US8768031B2 (en) * | 2010-10-01 | 2014-07-01 | Mistretta Medical, Llc | Time resolved digital subtraction angiography perfusion measurement method, apparatus and system |
| CN103562959B (zh) * | 2011-03-02 | 2017-05-24 | 皇家飞利浦有限公司 | 用于提供支持对血管介入程序中的介入设备的准确引导的图像表示的医疗成像***和方法 |
| WO2012174263A2 (en) | 2011-06-15 | 2012-12-20 | Mistretta Medical, Llc | System and method for four dimensional angiography and fluoroscopy |
| US9031288B2 (en) * | 2012-04-18 | 2015-05-12 | International Business Machines Corporation | Unique cardiovascular measurements for human identification |
| CN104684482B (zh) * | 2012-10-01 | 2019-02-12 | 东芝医疗***株式会社 | 图像处理装置以及x射线ct装置 |
| RU2015120325A (ru) * | 2012-10-31 | 2016-12-20 | Конинклейке Филипс Н.В. | Визуализация перфузии |
| WO2014199995A1 (ja) * | 2013-06-11 | 2014-12-18 | 株式会社 東芝 | Ｘ線コンピュータ断層撮影装置、およびスキャン開始タイミング決定方法 |
| KR101604812B1 (ko) * | 2014-01-15 | 2016-03-18 | 삼성전자주식회사 | 의료 영상 처리 장치 및 그에 따른 의료 영상 처리 방법 |
| JP6576690B2 (ja) | 2014-05-29 | 2019-09-18 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
| US20170325755A1 (en) * | 2014-12-17 | 2017-11-16 | Koninklijke Philips N.V. | Perfusion imaging |
| DE102015208929B3 (de) | 2015-05-13 | 2016-06-09 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Verfahren zur 2D-3D-Registrierung, Recheneinrichtung und Computerprogramm |
| CN107735814A (zh) | 2015-06-30 | 2018-02-23 | 皇家飞利浦有限公司 | 用于处理计算机断层摄影成像数据的设备和方法 |
| US10223331B2 (en) * | 2016-04-08 | 2019-03-05 | Goodrich Corporation | Warp models for registering multi-spectral imagery |
| US11996129B2 (en) | 2016-06-10 | 2024-05-28 | Cornell University | Semiconductor circuits and devices based on low-energy consumption semiconductor structures exhibiting multi-valued magnetoelectric spin hall effect |
| CN107884164B (zh) * | 2017-09-25 | 2019-07-26 | 广西电网有限责任公司电力科学研究院 | 一种ncc-sar算法的断路器弹簧性能检测方法 |
| US20210081698A1 (en) * | 2018-02-09 | 2021-03-18 | Nano Techgalaxy, Inc. D/B/A Galaxy.Ai | Systems and methods for physical object analysis |
| JP7164964B2 (ja) * | 2018-04-04 | 2022-11-02 | キヤノン株式会社 | 情報処理装置、放射線撮影装置、放射線撮影システム、情報処理方法及びプログラム |
| JP7270453B2 (ja) * | 2019-04-26 | 2023-05-10 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
| RU2732895C1 (ru) * | 2019-05-27 | 2020-09-24 | Общество с ограниченной ответственностью "ПЛАТФОРМА ТРЕТЬЕ МНЕНИЕ" | Метод для выделения и классификации типов клеток крови с помощью глубоких сверточных нейронных сетей |
| WO2021253294A1 (zh) * | 2020-06-17 | 2021-12-23 | 深圳高性能医疗器械国家研究院有限公司 | 一种ct图像的分离重建方法及应用 |
| WO2022165526A1 (en) * | 2021-01-28 | 2022-08-04 | Cornell University | Resolution-enhanced optical coherent tomography |
| CN116309922B (zh) * | 2023-05-22 | 2023-08-11 | 杭州脉流科技有限公司 | Ct灌注影像的去伪影方法、装置、设备和存储介质 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2000118170A (ru) * | 2000-07-12 | 2002-06-27 | Научно-исследовательский институт неврологии Российской АМН | Способ совмещения изображений биологического объекта |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2989325B2 (ja) * | 1991-06-12 | 1999-12-13 | 日本放送協会 | 動きベクトル信頼性の判定装置 |
| EP1658588A1 (en) * | 2003-08-21 | 2006-05-24 | Philips Intellectual Property & Standards GmbH | Device and method for combined display of angiograms and current x-ray images |
| US7672540B2 (en) * | 2005-07-13 | 2010-03-02 | Siemens Medical Solutions USA, Inc, | Nonrigid registration of cardiac perfusion MR images using adaptive local template matching |
| US20090169080A1 (en) | 2005-08-09 | 2009-07-02 | Koninklijke Philips Electronics, N.V. | System and method for spatially enhancing structures in noisy images with blind de-convolution |
| DE102006020864A1 (de) * | 2006-05-04 | 2007-11-08 | Siemens Ag | Verfahren zur Ermittlung und Anzeige von wenigstens einer Information über ein Zielvolumen |
| US7693349B2 (en) * | 2006-08-15 | 2010-04-06 | General Electric Company | Systems and methods for interactive image registration |
| DE102006051919B4 (de) | 2006-10-31 | 2016-09-01 | Siemens Healthcare Gmbh | Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes des Herzens eines Patienten |
| US8098911B2 (en) | 2006-12-05 | 2012-01-17 | Siemens Aktiengesellschaft | Method and system for registration of contrast-enhanced images with volume-preserving constraint |
| DE102007028270B4 (de) * | 2007-06-15 | 2013-11-14 | Siemens Aktiengesellschaft | Verfahren zur Segmentierung von Bilddaten zur Erkennung einer Leber |
| DE102008016892A1 (de) | 2008-04-02 | 2009-10-15 | Siemens Aktiengesellschaft | Betriebsverfahren für eine Bildgebungsanlage zur zeitaufgelösten Abbildung eines sich iterativ bewegenden Untersuchungsobjekts |
| US8208709B2 (en) * | 2008-04-17 | 2012-06-26 | The Ohio State University Research Foundation | System and method for improved real-time cine imaging |
| US8755635B2 (en) * | 2008-08-11 | 2014-06-17 | Siemens Aktiengesellschaft | Method and system for data dependent multi phase visualization |
| JP5551960B2 (ja) * | 2009-09-30 | 2014-07-16 | 富士フイルム株式会社 | 診断支援システム、診断支援プログラムおよび診断支援方法 |
-
2011
- 2011-04-12 JP JP2013508584A patent/JP5833637B2/ja not_active Expired - Fee Related
- 2011-04-12 US US13/643,167 patent/US9002089B2/en not_active Expired - Fee Related
- 2011-04-12 CN CN201180022575.8A patent/CN102884550B/zh not_active Expired - Fee Related
- 2011-04-12 RU RU2012152513/08A patent/RU2582055C2/ru not_active IP Right Cessation
- 2011-04-12 WO PCT/IB2011/051579 patent/WO2011138694A1/en active Application Filing
- 2011-04-12 EP EP11725511.7A patent/EP2567359B1/en not_active Not-in-force
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2000118170A (ru) * | 2000-07-12 | 2002-06-27 | Научно-исследовательский институт неврологии Российской АМН | Способ совмещения изображений биологического объекта |
Non-Patent Citations (1)
| Title |
|---|
| Isola et al, "Fully automatic nonrigid registration-based local motion estimation for motion-corrected iterative cardiac CT reconstruction", Medical Physics, vol. 37, N3, 03.2010. Breeuwer M. et al, "Automatic QUantative Analysis of Cardiac MR Perfusion Images", Proceedings of Spie, The International Society for Optical Ebgineering Spie, USA, vol. 4322, 19.02.2001. Bruder H. et al, "Spatio-temporal filtration of dynamic CT data using diffusion filters", Proceedings of Spie, vol. 7258, 01.01.2009. Stanz K. M. et al, "In-vivo regional myocardial perfusion measurements in a porcine model by ECG-gated multi-slice computed tomography", Proceedings of Spie, The International Society for Optical Engineering Spie, USA, vol. 5031, 01.01.2003. * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20130039559A1 (en) | 2013-02-14 |
| EP2567359A1 (en) | 2013-03-13 |
| US9002089B2 (en) | 2015-04-07 |
| WO2011138694A1 (en) | 2011-11-10 |
| RU2012152513A (ru) | 2014-06-20 |
| JP2013525056A (ja) | 2013-06-20 |
| EP2567359B1 (en) | 2014-10-29 |
| JP5833637B2 (ja) | 2015-12-16 |
| CN102884550B (zh) | 2016-03-16 |
| CN102884550A (zh) | 2013-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2582055C2 (ru) | Совмещение данных изображения для динамической перфузионной компьютерной томографии | |
| US10748293B2 (en) | Tomography apparatus and method for reconstructing tomography image thereof | |
| CN109035355B (zh) | 用于pet图像重建的***和方法 | |
| US7782998B2 (en) | Method and apparatus for correcting motion in image reconstruction | |
| US7382852B2 (en) | Method and apparatus for correcting motion in image reconstruction | |
| Isola et al. | Fully automatic nonrigid registration‐based local motion estimation for motion‐corrected iterative cardiac CT reconstruction | |
| CN104424659B (zh) | 图像数据的基于图像的运动补偿 | |
| US20080267455A1 (en) | Method for Movement Compensation of Image Data | |
| US20110044559A1 (en) | Image artifact reduction | |
| US9595079B2 (en) | High resolution volumetric image generation based on super resolution post-processing and point spread function | |
| Lossau et al. | Motion estimation and correction in cardiac CT angiography images using convolutional neural networks | |
| US10034648B2 (en) | System and method for motion artifacts reduction | |
| US11940578B2 (en) | Super resolution in positron emission tomography imaging using ultrafast ultrasound imaging | |
| WO2007072279A1 (en) | Iterative image reconstruction of a moving object from projection data | |
| WO2013130086A1 (en) | Integrated image registration and motion estimation for medical imaging applications | |
| RU2596982C2 (ru) | Проверка и определение местоположения сердца в эмиссионных изображениях | |
| Van Stevendaal et al. | A motion‐compensated scheme for helical cone‐beam reconstruction in cardiac CT angiography | |
| KR101783964B1 (ko) | 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법 | |
| Pack et al. | Four‐dimensional computed tomography of the left ventricle, Part I: motion artifact reduction | |
| Parker et al. | Respiratory motion correction in gated cardiac SPECT using quaternion‐based, rigid‐body registration | |
| Perez-Liva et al. | Ultrafast ultrasound imaging for super-resolution preclinical cardiac PET | |
| JP2003190147A (ja) | Ｘ線ｃｔ装置及び画像処理装置 | |
| US20230079025A1 (en) | Information processing apparatus, information processing method, and non-transitory computer readable medium | |
| Hsieh | Local motion detection, characterization, and quantification for x-ray CT | |
| Asif | Enhanced CNN for Image Enhancement Improved Approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20170413 |