Способ получения субтракционного ангиографического изображения

Номер патента: 24855

Опубликовано: 31.10.2016

Автор: Афанасенко Арсений Сергеевич

Есть еще 2 страницы.

Смотреть все страницы или скачать PDF файл.

Формула / Реферат

1. Способ получения ангиографического изображения, заключающийся в том, что посредством рентгенологического оборудования выполняют съемку пациента с получением предконтрастной серии из N цифровых изображений, вводят рентгеноконтрастное вещество в исследуемые сосуды пациента, выполняют рентгенографическую съемку с получением постконтрастной серии из M цифровых изображений, производят пространственное совмещение цифровых изображений предконтрастной и постконтрастной серий, вычитают цифровые изображения предконтрастной серии из цифровых изображений постконтрастной серии и полученное ангиографическое цифровое изображение передают на устройство вывода, отличающийся тем, что на этапе совмещения цифровых изображений выполняют поиск характерных деталей на каждом изображении, по смещению названных деталей определяют вероятные смещения органов пациента, осуществляют сегментацию изображения предконтрастной серии на области, в которых характерные детали смещены согласованно, для каждой области рассчитывают афинное преобразование и выполняют соответствующие афинные преобразования для названных областей на цифровых изображениях предконтрастной серии, на каждом из полученных цифровых изображений органы и ткани, затеняющие друг друга, разделяют по меньшей мере на два перекрывающихся изображения - слои, соответствующие органам и тканям, находящимся на различной глубине; совмещение изображений проводят независимо для каждого слоя и результаты затем суммируют.

2. Способ по п.1, отличающийся тем, что в качестве характерных деталей выбирают такие элементы костных структур и органов пациента, положение которых на цифровых изображениях на протяжении съемки однозначно определяется независимо от изменения их положения в пространстве.

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

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

5. Способ по любому из пп.1-4, отличающийся тем, что после получения предконтрастной серии из N цифровых изображений формируют единственное опорное предконтрастное цифровое изображение путем взвешенного суммирования цифровых изображений, первое из которых запоминают без изменений, а для каждого последующего перед суммированием с текущим опорным изображением выполняют пространственное совмещение с ним.

6. Способ по любому из пп.1-4, отличающийся тем, что в процессе получения постконтрастной серии цифровых изображений выбирают опорное изображение из постконтрастной серии и повторяют для него поиск характерных деталей, оценку смещений и сегментацию.

Рисунок 1

Текст

Смотреть все

СПОСОБ ПОЛУЧЕНИЯ СУБТРАКЦИОННОГО АНГИОГРАФИЧЕСКОГО ИЗОБРАЖЕНИЯ Изобретение относится к области цифровой обработки рентгенограмм и может быть использовано в цифровой субтракционной ангиографии для компенсации влияния непроизвольного движения пациента и естественного движения его внутренних органов на изображение сосудистой системы. Техническим результатом заявляемого изобретения является повышение диагностической ценности субтракционного ангиографического изображения за счет устранения артефактов изображения, вызванных подвижностью анатомических структур. Технический результат достигается за счет того, что на этапе совмещения цифровых изображений выполняют поиск характерных деталей на каждом изображении, по смещению названных деталей определяют вероятные смещения органов пациента, осуществляют сегментацию изображения предконтрастной серии на области, в которых характерные детали смещены согласованно, для каждой области рассчитывают афинное преобразование и выполняют соответствующие афинные преобразования для всех областей на цифровых изображениях предконтрастной серии. Афанасенко Арсений Сергеевич (RU) Михайлова Г.Л. (RU)(71)(73) Заявитель и патентовладелец: ЗАКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО "ИМПУЛЬС" (RU) Область техники Изобретение относится к области цифровой обработки рентгенограмм и может быть использовано в цифровой субтракционной ангиографии для компенсации влияния непроизвольного движения пациента и естественного движения его внутренних органов на изображение сосудистой системы. Предшествующий уровень техники Ангиография является методом визуализации кровеносных сосудов, основанным на введении контрастного вещества в исследуемую часть сосудистой системы пациента. На ангиографических снимках,помимо сосудов, присутствуют изображения близлежащих органов и тканей, на фоне которых сосуды,заполненные контрастным веществом, во многих случаях слабо различимы. В субтракционной ангиографии для уменьшения видимости близлежащих органов формируют предконтрастное изображение рентгенографическое изображение исследуемой части организма пациента, сформированное до начала ввода контрастного вещества в сосудистую систему. В дальнейшем выполняется съемка серии изображений сосудов по мере заполнения их рентгеноконтрастным веществом - постконтрастных изображений. Под термином постконтрастное изображение понимаются все изображения, зарегистрированные не ранее начала ввода контрастного вещества и отражающие фазы заполнения сосудов контрастным веществом, продвижения контрастного вещества по сосудам и постепенного растворения контрастного вещества в крови. Предконтрастное изображение вычитается из каждого постконтрастного изображения, благодаря чему компенсируются те детали изображения, которые не относятся к сосудистой системе - кости,мягкие ткани и т. д. Это значительно улучшает различимость сосудов и облегчает проведение лечебнодиагностических процедур. Поскольку получение предконтрастных и постконтрастных изображений происходят в различные моменты времени, любые изменения пространственного положения органов и тканей, присутствующих на изображении, приводят к появлению артефактов. Особенно отчетливо такие артефакты проявляются на контрастных границах подвижных объектов. В некоторых случаях контрастность артефактов превышает контрастность сосудов, что приводит к маскированию существенных фрагментов сосудистой системы и делает существенные детали неразличимыми. Известен способ совмещения изображений применительно к цифровой субтракционной ангиографии по патенту US7409078. В соответствии с этим способом ангиографическое изображение формируют следующим образом: получают предконтрастное изображение перед началом ввода контрастного вещества в сосудистую систему пациента; вводят контрастное вещество в сосудистую систему пациента; получают постконтрастное изображение; на предконтрастное изображение накладывают сетку с ячейками правильной геометрической формы; выполняют пространственное совмещение постконтрастного изображения с предконтрастным по контрольным точкам, размещенным в узлах сетки; после совмещения вычитают одно изображение из другого и отображают результат на дисплее; проверяют устойчивость свойств контрольных точек, и если она оказывается недостаточной, сетку разбивают на ячейки меньшего или большего размера. Недостатком способа является неполный учет характеристик обрабатываемого изображения, что приводит к субоптимальной коррекции артефактов движения. Известен способ и устройство для цифровой субтракционной ангиографии по патенту US5848121. В соответствии с этим способом ангиографическое изображение получают следующим образом: получают предконтрастное изображение перед началом ввода контрастного вещества в сосудистую систему пациента; вводят контрастное вещество в сосудистую систему пациента; получают постконтрастное изображение; находят множество характерных деталей на предконтрастном изображении и множество характерных деталей на постконтрастном изображении; определяют соответствие между найденными характерными деталями в двух изображениях; формируют модель локально-адаптивного афинного преобразования на основе координат характерных точек; применяют полученное преобразование для деформации предконтрастного изображения; формируют субтракционное изображение путем вычитания логарифма деформированного предконтрастного изображения из логарифма постконтрастного изображения. Недостатком способа является то, что ошибочное определение соответствия характерных точек предконтрастного и постконтрастного изображений приводит к появлению неустранимых дефектов на результирующем изображении. Наиболее близким к заявляемому является способ совмещения, вычитания и отображения рентгеновских изображений грудной клетки по заявке US2010/0266188. Способ состоит из следующих опе-1 024855 раций: получают предконтрастное изображение перед началом ввода контрастного вещества в сосудистую систему пациента; вводят контрастное вещество в сосудистую систему пациента; получают постконтрастное изображение; таким образом, получают не менее двух рентгеновских снимков, содержащих изображения одних и тех же органов; проводят предварительную обработку изображений, включающую нормировку по размерам, яркости и глубине цвета, а также выделение зон интереса; выполняют грубое совмещение изображений, включая компенсацию сдвига и поворота; выполняют точное совмещение изображений, используя корреляционное совмещение мелких деталей либо методы вычисления оптического потока; вычитают совмещенные изображения друг из друга таким образом, который обеспечивает наилучшую видимость необходимых деталей и наилучшее подавление посторонних деталей. Общими признаками с заявляемым способом являются: 1. Оценка афинных преобразований, связывающих совмещаемые изображения по ограниченному множеству контрольных областей. 2. Возможность компенсации сдвига, поворота, а также локальных деформаций исследуемых органов на изображении. 3. Многомасштабный анализ изображений при оценке движения. 4. Разностная обработка совмещенных изображений с целью подчеркивания определенных деталей(сосудов в заявляемом способе или легочных микрообразований в прототипе). Основным недостатком способа по заявке US2010/0266188 является ограниченная точность оценки движения при определенных характеристиках исследуемого объекта, что существенно влияет на качество получаемого ангиографического изображения. Это вызвано последовательной схемой компенсации движения на различных масштабах, от более грубого к более точному. Как показывает пространственно-частотный анализ большого числа изображений, встречаемых в клинической практике, многие изображения содержат существенные детали, проявляющиеся только при понижении пространственного разрешения. В таких случаях указанный способ дает точные результаты лишь до определенного момента, а в дальнейшем накапливает ошибку. Кроме того, по известному способу функция преобразования координат является гладкой в пределах области интереса, что негативно сказывается на точности компенсации движения в случае, когда в области интереса находится не менее двух движущих фрагментов, границы которых четко выражены. Раскрытие изобретения Техническим результатом заявляемого изобретения является повышение диагностической ценности субтракционного ангиографического изображения за счет устранения артефактов подвижности анатомических структур. Технический результат в способе получения субтракционного ангиографического изображения, заключающемся в том, что посредством рентгенологического оборудования выполняют съемку пациента с получением предконтрастной серии из N цифровых изображений, вводят рентгеноконтрастное вещество в исследуемые сосуды пациента, выполняют рентгенографическую съемку с получением постконтрастной серии из M цифровых изображений, производят пространственное совмещение цифровых изображений предконтрастной и постконтрастной серий, вычитают цифровые изображения предконтрастной серии из цифровых изображений постконтрастной серии и полученное ангиографическое цифровое изображение передают на устройство вывода достигается тем, что на этапе совмещения цифровых изображений выполняют поиск характерных деталей на каждом изображении, по смещению названных деталей определяют вероятные смещения органов пациента, осуществляют сегментацию изображения предконтрастной серии на области, в которых характерные детали смещены согласованно, для каждой области рассчитывают афинное преобразование и выполняют соответствующие афинные преобразования для всех областей на цифровых изображениях предконтрастной серии. Возможен вариант осуществления настоящего изобретения, где в качестве характерных деталей выбирают такие элементы костных структур и органов пациента, положение которых на цифровых изображениях на протяжении съемки однозначно определяется независимо от изменения их положения в пространстве. Возможен вариант осуществления настоящего изобретения, где при сегментации цифрового изображения предконтрастной серии на области, в которых характерные детали смещены согласованно,приписывают каждой характерной детали конечное множество возможных смещений относительно текущего изображения постконтрастной серии и выделяют области на изображении, в пределах которых характерные детали смещаются согласованно. Возможен вариант осуществления настоящего изобретения, где в изображениях предконтрастной серии выделяют неподвижные части рентгенологического оборудования и на этапе совмещения цифровых изображений этим областям присваивают нулевое смещение. Возможен вариант осуществления настоящего изобретения, согласно которому после получения предконтрастной серии из N цифровых изображений формируют единственное опорное предконтрастное цифровое изображение путем взвешенного суммирования цифровых изображений, первое из которых запоминают без изменений, а для каждого последующего перед суммированием с текущим опорным изображением выполняют пространственное совмещение с ним. Возможен вариант осуществления настоящего изобретения, согласно которому в процессе получения постконтрастной серии цифровых изображений выбирают опорное изображение также из постконтрастной серии и повторяют для него поиск характерных деталей, оценку смещений и сегментацию. Возможен вариант осуществления настоящего изобретения, согласно которому на каждом из полученных цифровых изображений органы и ткани, затеняющие друг друга, разделяют по меньшей мере на два перекрывающихся изображения - слои, соответствующие органам и тканям, находящимся на различной глубине; совмещение изображений проводят независимо для каждого слоя и результаты затем суммируют. Подробное описание изобретения Заявляемое техническое решение, возможность его технической реализации и достижение указанного технического результата поясняется фиг. 1-4. Фиг. 1. Схема аппарата. Фиг. 2. Состав изображения постконтрастной серии. Фиг. 3. Сегментация подвижных органов. Фиг. 4. Результаты обработки. Получение цифровых рентгенограмм осуществляют, например, с помощью устройства, показанного на фиг. 1. Оно содержит рентгеновскую трубку 1, которая испускает поток рентгеновского излучения 2. Поток рентгеновского излучения 2 проходит через тело пациента 3, находящегося на столе 4, и поступает на приемник 5. Приемник 5 осуществляет преобразование потока рентгеновского излучения в цифровое изображение. В одном из возможных вариантов реализации приемник содержит сцинтилляционный экран (на чертеже не показан), преобразующий рентгеновское излучение в видимый свет, и матрицу фоточувствительных элементов (на чертеже не показаны). Рентгеновская трубка 1 и приемник 5 неподвижно закреплены на штативе 6, который выполнен с возможностью перемещения относительно стола 4 в различных плоскостях. Согласно заявляемому способу выполняют съемку пациента с получением предконтрастной серии из N цифровых изображений (где N1). После ввода рентгеноконтрастного вещества в исследуемые сосуды пациента получают постконтрастную серию из M цифровых изображений (где M1). Для каждого изображения постконтрастной серии выбирают соответствующее опорное изображение, которое может принадлежать предконтрастной серии, принадлежать постконтрастной серии или быть результатом совместной обработки нескольких предшествующих изображений из обеих серий. В соответствии с различными реализациями изобретения возможны следующие варианты получения и обработки изображений предконтрастной и постконтрастной серий. 1. При неподвижном штативе получают предконтрастную серию из N цифровых изображений (гдеN1). Наиболее качественное изображение, выбранное вручную (или единственное предконтрастное изображение в случае N = 1), используют в качестве опорного при обработке всех изображений постконтрастной серии. 2. При неподвижном штативе получают предконтрастную серию из не менее чем двух цифровых изображений. Наиболее качественное изображение, выбранное вручную, запоминают в качестве опорного. Для всех последующих изображений предконтрастной серии выполняют пространственное совмещение с опорным изображением. Совмещенные изображения суммируют, результат суммирования запоминают в качестве обновленного опорного изображения. В результате на опорном изображении ослабляются шумы. При обработке всех изображений постконтрастной серии используется одно опорное изображение. 3. В случае, когда ангиографическое исследование требует перемещения штатива, например, при исследовании протяженных сосудов пациента, изображения предконтрастной серии получают при различных положениях штатива. Траектории перемещения штатива при съемке предконтрастной и постконтрастной серий совпадают. Каждое изображение предконтрастной серии является опорным при обработке изображения постконтрастной серии, полученного при таком же положении штатива. 4. В случаях, когда при получении постконтрастной серии произошло непредусмотренное перемещение штатива, существенное движение пациента, или по другим причинам изменилось содержание изображения, с этого момента в качестве опорного изображения выбирают вместо изображения предконтрастной серии одно из изображений постконтрастной серии, обработанное с целью уменьшения заметности сосудов, заполненных контрастным веществом. На фиг. 2 показан состав цифрового ангиографического изображения постконтрастной серии. Прямоугольное рабочее поле 7 приемника частично закрыто деталями 8 рентгеновского аппарата, неподвижными во время исследования. Изображение пациента содержит заполненные контрастным веществом кровеносные сосуды 9, а также органы и ткани 10, не являющиеся объектами ангиографического исследования. Для совмещения изображения постконтрастной серии с опорным изображением определяют такое афинное преобразование, при котором достигается наилучшее совпадение этих двух изображений. За счет предварительного выделения неподвижных областей ограничивают область, для которой требуется совмещение изображений. Это необходимо в тех случаях, когда часть рабочего поля 7 приемника закрыта створками коллиматора или в составе приемника используется электронно-оптический преобразователь круглого сечения. Границы неподвижных областей имеют априорно известную форму прямых или гладких кривых линий, а также высокий контраст, что позволяет использовать для их поиска известные алгоритмы выделения геометрических примитивов (например, преобразование Хафа) в сочетании с методами регрессионного анализа. Неподвижность данных областей позволяет выделить их единожды при обработке каждого опорного изображения. При совмещении изображений сначала выбирают конечное множество контрольных областей на опорном изображении. Контрольные области размещают таким образом, чтобы в них попадали характерные детали изображения. Характерная деталь должна обладать свойствами, позволяющими достоверно выделять ее на изображениях данного объекта независимо от яркостных и афинных преобразований, а также при наличии шумов. Для выделения контрольных областей выполняют следующие действия: 1. Формируют многомасштабное представление изображения, применяя пирамидальное разложение Гаусса или аналогичное преобразование. Число масштабов задают вручную или согласовывают с формой спектра пространственных частот изображения. 2. На каждом масштабе к изображению применяют математический оператор подчеркивания характерных деталей (далее, для краткости -"детектор"). В одном из вариантов реализации изобретения данный оператор основан на вычислении ковариационной матрицы текстуры в каждой точке изображения. Откликом детектора является произведение собственных значений данной матрицы. Локальные максимумы отклика детектора соответствуют перепадам яркости, изгибам и пересечениям границ объектов. 3. Ограничивают множество характерных деталей, выбирая детали с максимальным значением отклика детектора, найденные на различных масштабах. Для корректного сравнения значений отклика детектора на различных масштабах вводят поправочный множитель, монотонно возрастающий при переходе от более крупного масштаба к более мелкому. 4. В окрестности каждой из найденных характерных деталей формируют контрольную область,размер которой соответствует масштабу, на котором получен максимальный отклик детектора. Каждой области может быть дополнительно приписана пространственная ориентация, соответствующая главному собственному вектору ковариационной матрицы текстуры. Возможен упрощенный вариант реализации изобретения, при котором поиск характерных деталей не осуществляется. Вместо этого на изображение накладывается геометрическая сетка априорно заданной структуры (например, прямоугольная) и множество контрольных областей располагается в узлах сетки. Результатом анализа опорного изображения является совокупность контрольных областей, для каждой из которых сохраняют координаты центра, размеры и, в зависимости от реализации изобретения,дополнительные атрибуты. Для каждого цифрового изображения постконтрастной серии, а также для изображений предконтрастной серии, не являющихся опорными, осуществляют процедуру совмещения с опорным изображением, состоящую из следующих действий: 1. Осуществляют выделение контрольных областей на изображении тем же способом, что при обработке опорного изображения. 2. Выполняют сравнение атрибутов контрольных областей текущего изображения с сохраненными атрибутами контрольных областей опорного изображения. На основе корреляции определяют взаимнооднозначные соответствия между контрольными областями. В результате каждой контрольной области приписывают вектор смещения относительно опорного изображения. В общем случае, множество контрольных областей, выделенных на опорном изображении, может не полностью совпадать с множеством контрольных областей текущего изображения. Возникающие при этом ошибки в определении смещений устраняют статистическими методами при дальнейшей обработке. 3. Выделяют из множества контрольных областей группы, перемещение которых происходит согласованно и с достаточной точностью описывается выбранной моделью афинного преобразования. В одном из возможных вариантов реализации изобретения такой моделью является аффинное преобразование, подразумевающее различные комбинации сдвига, поворота, наклона и изменения масштаба изображения. Для примера на фиг. 3 показана сегментация (разбиение) латеральной проекции головы пациента на три области с различной подвижностью: череп 11, нижнюю челюсть 12 и плечи 13. Центрам контрольных областей 14 приписаны векторы смещения, в том числе ложные векторы 15. В результате кластерного анализа смещений приписывают каждой группе контрольных областей общую модель афинного преобразования. 4. Определяют границы сегментов изображения, содержащих выделенные группы контрольных об-4 024855 ластей, используя информацию о градиенте яркости изображения. 5. Преобразуют опорное изображение, применяя к выделенным сегментам изображения рассчитанные афинные преобразования. На границах сегментов возможно применение интерполяции для ослабления краевых эффектов. После того как выполнено совмещение изображений, осуществляют поэлементное арифметическое деление яркости текущего изображения на яркости преобразованного опорного изображения. Эта операция соответствует логарифмическому вычитанию. Результат является искомым субтракционным изображением. На фиг. 4 представлено сравнение ангиографического изображения 16, полученного без применения компенсации движения, и ангиографического изображения 17, полученного заявляемым способом. В первом случае на изображения сосудов 18, заполненных контрастным веществом, наложены артефакты движения 19, не являющиеся изображениями каких-либо существующих органов или тканей. Во втором случае артефакты 20 выражены существенно слабее. Благодаря сегментации поля смещений, полученного после сопоставления контрольных областей опорного изображения и текущего изображения, достигают меньшей остаточной видимости органов,ухудшающих различимость исследуемых сосудов, на ангиографическом изображении, что повышает диагностическую ценность изображений, получаемых заявляемым способом. Лучший вариант осуществления изобретения Рентгенографические изображения получают с помощью аппарата, у которого неподвижные части не проецируются на рабочее поле преемника. Выполняют рентгенографическую съемку с получением предконтрастной серии из 2-3 цифровых изображений. На первом цифровом изображении предконтрастной серии осуществляют поиск характерных деталей. Для этого над ним выполняют следующие операции. 1. Формируют многомасштабное представление изображения, состоящее из 4 уровней. Первым уровнем является само изображение. Далее для формирования каждого следующего уровня применяют сглаживающий цифровой фильтр с апертурой 33 элемента и прореживают изображение вдвое. 2. На каждом уровне многомасштабного изображения выделяют характерные детали. Для этого сначала вычисляют первые частные производные яркости по горизонтали и по вертикали Затем в каждом пикселе изображения вычисляют ковариационную матрицу текстуры где g - импульсная характеристика сглаживающего фильтра Гаусса с апертурой 55 элементов. Затем вычисляют отклик оператора подчеркивания углов, взвешенный с учетом масштаба где n[0,3] - номер текущего уровня многомасштабного представления,det C - определитель ковариационной матрицы текстуры,trC - след ковариационной матрицы текстуры. Затем выбирают такие точки на изображении, для которых значение h является локальным максимумом в пределах окрестности радиуса r. Величину r назначают равной 2% от минимального из линейных размеров изображения (с округлением до ближайшего целого числа). 3. Объединяют результаты выделения характерных деталей на всех уровнях многомасштабного представления и формируют начальное множество характерных деталей. При этом каждой характерной детали приписывают координаты и значение локального максимума h, а также номер уровня многомасштабного представления, на котором данная деталь была выделена. Исключают из множества те детали,вблизи которых на других уровнях располагаются детали с большим значением h (радиус окрестности для исключения повторяющихся деталей назначают равным 1% от минимального из линейных размеров изображения). Полученное множество сортируют по убыванию значения h и оставляют в итоговом множестве не более 100 первых деталей. Каждой детали ставят в соответствие контрольную область, размер которой тем меньше, чем крупнее масштаб, на котором выделенная характерная деталь. Вводят рентгеноконтрастное вещество в исследуемые сосуды пациента и выполняют рентгенографическую съемку с получением постконтрастной серии из M цифровых изображений (где M1). Для группы изображений, включая не менее одного изображения предконтрастной серии и все изображения постконтрастной серии, выполняют совмещение изображений, состоящее из следующих операций: 1. Формируют многомасштабное представление изображения, состоящее из 4 уровней. Первым уровнем является само изображение. Далее для формирования каждого следующего уровня применяют сглаживающий цифровой фильтр с апертурой 33 элемента и прореживают отсчеты полученного изображения в два раза. 2. Для каждой контрольной области корреляционным методом определяют смещение и поворот. Поиск начинают на том уровне многомасштабного представления, которому соответствует размер характерной детали, и продолжают на более крупном масштабе, последовательно уменьшая диапазон поиска и увеличивая точность. Поиск заключается в максимизации коэффициента взаимной корреляции между фрагментами опорного и текущего изображений где TR - фрагмент опорного изображения R, соответствующий контрольной области,IR - фрагмент текущего изображения I, соответствующий контрольной области,warp1(I,dx,dy,) - преобразование изображения I, состоящее из сдвига на вектор (dx,dy) и поворота на угол ,cov - ковариация отсчетов яркости текущего и опорного изображений в пределах фрагмента,D - дисперсия. 3. Осуществляют разбиение множества контрольных областей на K сегментов таким образом, чтобы для каждого сегмента смещения контрольных областей с наименьшей погрешностью описывались аффинным преобразованием. Аффинное преобразование координат имеет вид где x, y - исходные координаты точки,x', y' - преобразованные координаты,а 1 а 6 - коэффициенты преобразования. Погрешность для каждого сегмента определяют как где k[1,K] - номер сегмента,i - номер контрольной области,Pk - подмножество контрольных областей, относящихся к k-му сегменту. Задача сегментации состоит в минимизации суммарной погрешности В процессе минимизации выражения (7) определяются непересекающиеся подмножества Pk, образующие полное множество контрольных областей, а также определяются соответствующие матрицы аффинного преобразования Ak. Число сегментов K априорно неизвестно, но для улучшения сходимости процедуры сегментации его ограничивают значением 5. В ходе итерационного процесса оптимизации сегменты с близкими матрицами преобразования объединяют, благодаря чему итоговое число сегментов,в зависимости от характера подвижности контрольных областей, составляет от 1 до 5. Выполняют обработку текущего изображения с учетом движения. Если текущее изображение относится к предконтрастной серии, к нему применяют инвертированное преобразование, вычисленное на предыдущем шаге. Опорное изображение обновляют путем взвешенного арифметического суммирования с модифицированным текущим изображением: где R' - опорное изображение после обновления,R - опорное изображение до обновления,t - номер изображения в обрабатываемой последовательности,It - текущее изображение,warp2(It, A1, A2, , Ak) - совместное преобразование, в результате которого каждый сегмент изо-6 024855 бражения It подвергается аффинному преобразованию, заданному матрицей Ak. Если текущее изображение относится к постконтрастной серии, для чего выполняют субтракцию. К опорному изображению применяют преобразование, вычисленное на предыдущем шаге. Из логарифма текущего изображения вычитают взвешенный логарифм модифицированного опорного изображения: где s - коэффициент субтракции в диапазоне от 0 до 1. Полученное изображение F отображают на устройстве вывода. Регулируя s, рентгенолог может делать заметными на устройстве вывода не только исследуемые сосуды, но и окружающие органы, служащие анатомическими ориентирами. Таким образом, диагностическая ценность полученного ангиографического изображения повышается за счет устранения артефактов изображения, вызванных подвижностью анатомических структур. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ получения ангиографического изображения, заключающийся в том, что посредством рентгенологического оборудования выполняют съемку пациента с получением предконтрастной серии изN цифровых изображений, вводят рентгеноконтрастное вещество в исследуемые сосуды пациента, выполняют рентгенографическую съемку с получением постконтрастной серии из M цифровых изображений, производят пространственное совмещение цифровых изображений предконтрастной и постконтрастной серий, вычитают цифровые изображения предконтрастной серии из цифровых изображений постконтрастной серии и полученное ангиографическое цифровое изображение передают на устройство вывода, отличающийся тем, что на этапе совмещения цифровых изображений выполняют поиск характерных деталей на каждом изображении, по смещению названных деталей определяют вероятные смещения органов пациента, осуществляют сегментацию изображения предконтрастной серии на области, в которых характерные детали смещены согласованно, для каждой области рассчитывают афинное преобразование и выполняют соответствующие афинные преобразования для названных областей на цифровых изображениях предконтрастной серии, на каждом из полученных цифровых изображений органы и ткани, затеняющие друг друга, разделяют по меньшей мере на два перекрывающихся изображения слои, соответствующие органам и тканям, находящимся на различной глубине; совмещение изображений проводят независимо для каждого слоя и результаты затем суммируют. 2. Способ по п.1, отличающийся тем, что в качестве характерных деталей выбирают такие элементы костных структур и органов пациента, положение которых на цифровых изображениях на протяжении съемки однозначно определяется независимо от изменения их положения в пространстве. 3. Способ по п.1, отличающийся тем, что при сегментации цифрового изображения предконтрастной серии на области с согласованным смещением приписывают каждой характерной детали конечное множество возможных смещений относительно текущего изображения постконтрастной серии и выделяют области на изображении, в пределах которых характерные детали смещаются согласованно. 4. Способ по п.1, отличающийся тем, что в изображениях предконтрастной серии выделяют неподвижные части рентгенологического оборудования и на этапе совмещения цифровых изображений этим областям присваивают нулевое смещение. 5. Способ по любому из пп.1-4, отличающийся тем, что после получения предконтрастной серии изN цифровых изображений формируют единственное опорное предконтрастное цифровое изображение путем взвешенного суммирования цифровых изображений, первое из которых запоминают без изменений, а для каждого последующего перед суммированием с текущим опорным изображением выполняют пространственное совмещение с ним. 6. Способ по любому из пп.1-4, отличающийся тем, что в процессе получения постконтрастной серии цифровых изображений выбирают опорное изображение из постконтрастной серии и повторяют для него поиск характерных деталей, оценку смещений и сегментацию.

МПК / Метки

МПК: G06T 7/00, G06K 9/46, A61B 6/02

Метки: ангиографического, способ, изображения, получения, субтракционного

Код ссылки

<a href="http://easpatents.com/10-24855-sposob-polucheniya-subtrakcionnogo-angiograficheskogo-izobrazheniya.html" rel="bookmark" title="База патентов Евразийского Союза">Способ получения субтракционного ангиографического изображения</a>

Похожие патенты