RU2542096C2 - Система для представления информации о вентиляции легких - Google Patents

Система для представления информации о вентиляции легких Download PDF

Info

Publication number
RU2542096C2
RU2542096C2 RU2011135966/14A RU2011135966A RU2542096C2 RU 2542096 C2 RU2542096 C2 RU 2542096C2 RU 2011135966/14 A RU2011135966/14 A RU 2011135966/14A RU 2011135966 A RU2011135966 A RU 2011135966A RU 2542096 C2 RU2542096 C2 RU 2542096C2
Authority
RU
Russia
Prior art keywords
ventilation
images
image
phase
processor
Prior art date
Application number
RU2011135966/14A
Other languages
English (en)
Other versions
RU2011135966A (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 RU2011135966A publication Critical patent/RU2011135966A/ru
Application granted granted Critical
Publication of RU2542096C2 publication Critical patent/RU2542096C2/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/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30061Lung

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

Группа изобретений относится к медицине, а именно к системам и способам отображения информации о вентиляции легких. Система содержит устройство ввода и процессор. Устройство ввода предназначено для получения множества КТ-изображений легкого, при этом каждое КТ-изображение соответствует одной фазе из, по меньшей мере, двух разных фаз в дыхательном цикле. Процессор выполнен с возможностью сравнения КТ-изображений, соответствующих разным фазам в дыхательном цикле, для определения поля векторов деформаций для каждой фазы, формирования для каждой фазы, изображения вентиляции на основании соответствующего поля векторов деформаций, пространственного совмещения изображений вентиляции и формирования, для, по меньшей мере, одной общей позиции на каждом из совмещенных изображений вентиляции, функции изменения во времени значения вентиляции для упомянутой общей позиции, при этом каждое значение вентиляции в функции изменения во времени значения вентиляции основано на векторных полях деформаций, соответствующих совмещенным изображениям вентиляции. Использование изобретения обеспечивает повышение точности оценки локальной вентиляции легких. 3 н и 8 з.п. ф-лы, 10 ил.

Description

ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ
Настоящее изобретение относится к системе для отображения информации о вентиляции легких, при этом система содержит устройство ввода для получения множества компьютерных томографических изображений (КТ-изображений) легкого, причем каждое КТ-изображение соответствует одной фазе, по меньшей мере, двух разных фаз дыхательного цикла, и процессор, сконфигурированный с возможностью:
- сравнения КТ-изображений, соответствующих последовательным фазам, для определения векторного поля деформаций для каждой фазы,
- формирования, для каждой фазы, изображения вентиляции на основании соответствующего векторного поля деформаций, и
- пространственного совмещения изображений вентиляции,
- и дополнительного содержания дисплея для отображения информации о вентиляции легких.
Настоящее изобретение дополнительно относится к способу отображения информации о вентиляции легких и к компьютерному программному продукту.
УРОВЕНЬ ТЕХНИКИ ИЗОБРЕТЕНИЯ
В японской заявке на патент JP 2005-028121 предложена система для измерения локальной вентиляции легких. Упомянутая система использует трехмерное (3-мерное) компьютерное томографическое сканирование (КТ-сканирование) легких в моменты вдоха и выдоха. Затем изображение вдоха и изображение выдоха совмещают, и получают векторное поле смещений. Локальную вентиляцию легких вычисляют из векторного поля смещений. Интегрирование локальных значений выполняют для получения суммарного объема вентиляции. На фиг. 3 в заявке на патент JP 2005-028121 представлено изображение, на котором визуализируется локальная вентиляция легких. На данном изображении указаны зоны расширения, зоны сокращения и зоны, имеющие постоянный объем.
Недостаток системы в соответствии с заявкой на патент JP 2005-028121 заключается в том, что выдаваемое изображение обеспечивает информацию только о локальном изменении объема для полного дыхательного цикла. Более подробная информация может помочь врачу, терапевту или исследователю лучше оценить физическое состояние обследуемых легких.
ЦЕЛЬ ИЗОБРЕТЕНИЯ
Целью настоящего изобретения является повышение точности оценки локальной вентиляции легких.
СУЩНОСТЬ ИЗОБРЕТЕНИЯ
В соответствии с первым аспектом изобретения, данная цель достигается путем обеспечения системы для отображения информации о вентиляции легких, при этом, система содержит устройство ввода и процессор. Устройство ввода обеспечено для получения множества КТ-изображений легкого, причем, каждое КТ-изображение соответствует одной фазе из, по меньшей мере, двух разных фаз в дыхательном цикле. Процессор сконфигурирован с возможностью: сравнения КТ-изображений, соответствующих разным фазам в дыхательном цикле для определения векторного поля деформаций для каждой фазы, формирования для каждой фазы изображения вентиляции на основании соответствующего векторного поля деформаций, пространственного совмещения изображений вентиляции и формирования для, по меньшей мере, одной общей позиции в каждом из совмещенных изображений вентиляции, функции изменения во времени значения вентиляции для упомянутой общей позиции, при этом каждое значение вентиляции в функции основано на векторных полях деформаций, соответствующих совмещенным изображениям вентиляции. Для отображения, по меньшей мере, одного из изображений вентиляции или КТ-изображений и диаграммы сформированной функции может быть обеспечен дисплей.
Путем вычисления локальной вентиляции, врачу может представляться степень сокращения или расширения для каждого перехода фаз. Объем вентиляции можно представлять для любой выбранной пользователем позиции в легком и для любого момента в течение дыхательного цикла, а не как общий объем полного дыхательного цикла. Данная система дает возможность отображать релевантную информацию во всех пространственных измерениях (2 измерениях или 3 измерениях) и одновременно отображать информацию во временном измерении. Представляемые на дисплее изображения вентиляции показывают распределение вентиляции легких для изображаемой области, в то время как диаграмма или диаграммы показывают изменение во времени локальной вентиляции в, по меньшей мере, одной характерной позиции или области. Позиции или области могут выбираться автоматически или путем взаимодействия с пользователем. Например, можно выбрать область опухоли или опорную область со здоровой тканью. При использовании системы в соответствии с изобретением, можно выполнять не только сравнение состояний максимального вдоха и максимального выдоха, но можно сравнивать любые две фазы, что обеспечивает измерение во времени локальной вентиляции легких. Следует отметить, что, даже если вентиляция, вычисленная для перехода от максимального вдоха к максимальному выходу, является одинаковой для двух разных позиций или двух пациентов, вентиляция от фазы к фазе может различаться. Информация из всех изображений вентиляции от фазы к фазе обеспечивает более надежное средство для ранней диагностики рака.
В варианте осуществления системы в соответствии с изобретением для сформированной функции подбирают модель дыхательного цикла. Подобранную модель можно отображать вместе со сформированной функцией, чтобы визуализировать различия между теоретической моделью и измеренной совокупностью параметров вентиляции. Кроме того, параметры, использованные для подборки модели дыхательного цикла соответственно сформированной функции, могут содержать полезную информацию, например амплитуду вентиляции или фазовый сдвиг вентиляции. Если данные параметры определяются для разных позиций в легком, их можно визуализировать в виде атласа легкого, показывающего значения параметров в соответствующих позициях, например, с использованием цветокодирования.
В дополнительном варианте осуществления системы в соответствии с изобретением, совокупность параметров вентиляции определяют несколько раз, например, в течение курса лечения. Затем несколько сформированных функций можно использовать для анализа трендов. Анализ трендов может обеспечивать, в результате, отображение последовательности совокупностей параметров вентиляции. В предпочтительном варианте вычисляют и отображают тренды важных параметров, например амплитуды вентиляции и/или сдвиг фаз вентиляции.
Возможно также определение совокупности параметров вентиляции для нескольких пациентов и применение полученных данных для создания статистической модели функциональных данных легкого и представления статистических данных в виде функционального атласа легкого.
В предпочтительном варианте обеспечено средство пользовательского ввода, которое связано с процессором и дает пользователю возможность выбирать, по меньшей мере, одну общую позицию в отображаемом изображении вентиляции. При использовании упомянутого средства ввода, терапевт или обследующий врач получают возможность более детально изучить локальную вентиляцию легких в конкретной, представляющей интерес позиции.
В соответствии со вторым аспектом изобретения, предлагается способ отображения информации о вентиляции легких, при этом способ содержит этап получения множества КТ-изображений легкого, причем каждое КТ-изображение соответствует одной фазе из, по меньшей мере, двух разных фаз дыхательного цикла, этап сравнения КТ-изображений, соответствующих разным фазам в дыхательном цикле, для определения векторного поля деформаций для каждой фазы, этап формирования для каждой фазы изображения вентиляции на основании соответствующего векторного поля деформаций, этап пространственного совмещения изображений вентиляции, этап формирования для, по меньшей мере, одной общей позиции на каждом из совмещенных изображений вентиляции, диаграммы изменения во времени значения вентиляции для упомянутой общей позиции, причем каждое значение вентиляции на диаграмме основано на векторных полях деформаций, соответствующих совмещенным изображениям вентиляции, и этап отображения, по меньшей мере, одного из изображений вентиляции или КТ-изображений и диаграммы.
В соответствии с третьим аспектом изобретения предлагается компьютерный программный продукт для обеспечения выполнения процессором вышеописанного способа.
Специалистам в данной области техники будет ясно, что, по меньшей мере, два вышеупомянутых варианта осуществления, исполнения и/или аспекта изобретения можно объединять любым, представляющимся полезным, способом.
На основе настоящего описания, специалистом в данной области техники могут быть созданы модификации и варианты системы, способа, устройства для получения изображения, рабочей станции и/или компьютерного программного продукта, которые соответствуют описанным модификациям и вариантам системы или способа.
Приведенные и другие аспекты изобретения станут понятны из приведенного далее объяснения со ссылками на нижеописанные варианты осуществления.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
На чертежах:
Фиг. 1 - схема системы в соответствии с изобретением,
Фиг. 2 - блок-схема последовательности операций способа в соответствии с изобретением,
Фиг. 3 - последовательность изображений компьютерных томографических изображений (КТ-изображений) пары легких,
Фиг. 4 - изображение последовательности векторного поля деформаций, полученного по КТ-изображениям, представленным на фиг.3,
Фиг. 5 - последовательность изображений вентиляции легких, полученных из векторных полей деформаций, представленных на фиг.4,
Фиг. 6 - примерный дисплей системы в соответствии с изобретением,
Фиг. 7 - диаграмма сформированной функции вентиляции вместе с подобранной моделью дыхательного цикла,
Фиг. 8 - изображение примерного атласа легкого,
Фиг. 9 - КТ-изображение легкого, на котором выделены характерные зоны, и
Фиг. 10 - график изменения во времени параметра вентиляции в течение десятинедельного периода.
ПОДРОБНОЕ ОПИСАНИЕ ИЗОБРЕТЕНИЯ
На фиг. 1 схематически изображена система в соответствии с изобретением. Система содержит процессор 15 с устройством 12 ввода для получения КТ-изображений легкого. КТ-изображения получают компьютерным томографическим сканером (КТ-сканером) 11, который может быть частью системы или соединен с ней. В изображения, получаемые КТ-сканером 11, вводится временная метка. Во время компьютерного томографического сканирования (КТ-сканирования), признаки дыхания пациента получают с помощью маркерного блока, помещенного, например, на верхнюю часть брюшной полости. Изображения проекций ретроспективно разделяют, на основании дыхательных фаз, например, на десять групп 3-мерных данных КТ-изображений (т.е. от фазы 0% до фазы 90%, с интервалом 10%). Временные метки вводят посредством, например, внутренних часов 13 КТ-сканера 11. Процессор 15 может записывать полученные изображения в запоминающем средстве 16, например на жестком диске. Изображения могут обрабатываться немедленно при получении или могут вызываться из запоминающего средства 16 в более поздний момент времени для обработки. Обработка приводит к формированию диаграмм, показывающих изменения во времени локальной вентиляции легких в некоторой локализации исследуемого легкого. После обработки КТ-изображений, дисплей 17 будет представлять изображение вентиляции легких и, по меньшей мере, одну из сформированных диаграмм. Обработка КТ-изображений описана ниже со ссылкой на фигуру 2. Для представления пользователю возможности управления обработкой и/или выбора режимов отображения и опций дисплея, можно обеспечить средство 18 пользовательского ввода.
На фиг. 2 представлена блок-схема последовательности операций способа в соответствии с изобретением. На этапе 22 ввода, КТ-изображения 71 принимаются в устройстве 12 ввода системы. КТ-изображения 71 отображают, по меньшей мере, часть легкого. КТ-изображения 71 получают в разные моменты дыхательного цикла пациента. В предпочтительном варианте, получают, приблизительно, десять изображений 71 на один дыхательный цикл, что имеет следствием разделение дыхательного цикла на десять разных фаз, изменяющихся от максимального выдоха до максимального вдоха и содержащих между ними промежуточные фазы. Посредством получения изображений на разных глубинах, можно получить 3-мерное изображение 71 для каждой фазы.
На этапе 23 сравнения, изображения 71 последовательных фаз сравнивают с использованием алгоритмов распознавания изображений и/или совмещения изображений. Различия между изображениями 71 двух последовательных фаз сохраняются, например, в виде векторного поля деформаций. Следует отметить, что различия между изображениями 71 двух последовательных фаз могут проистекать либо из дыхательного движения, либо из перемещений объекта обследования как целого. Поскольку перемещения пациента как целого являются более или менее одинаковыми для всех частей легкого, то точно совмещенные различия между двумя изображениями 71 можно компенсировать для данных перемещений. Вместо сравнения изображений 71 последовательных фаз, можно также сравнивать каждую фазу, например, с одной опорной фазой.
На этапе 24 определения вентиляции векторные поля деформаций служат для вычисления изображений 72 вентиляции. Вентиляция характеризует приток воздуха. Более высокие значения вентиляции соответствуют более высокому измеренному притоку воздуха, например, в л/мин. Вентиляцию можно вычислять, например, посредством вычисления отклонения векторных полей деформаций или определителя Якоби векторных полей деформаций, или посредством взятия значения Хаунсфилда из выбранной фазы и значения Хаунсфилда из другой деформированной фазы, которое соответствует вычисленному векторному полю деформаций, с последующим вычислением относительного изменения. Для каждой фазы формируется изображение 72 вентиляции. Примерные изображения 72 вентиляции представлены на фиг. 5. Изображения 72 вентиляции показывают, для каждой фазы и для каждой локализации в наблюдаемом легком, имеет ли место локальное расширение или сокращение ткани легкого. Смещение вентиляции выбирают так, чтобы значение 0 соответствовало отсутствию притока воздуха.
На этапе 25 совмещения, изображения 72 вентиляции совмещают. Одно из изображений 72 используют как опорное изображение. Например, состояние максимального вдоха может использовать как опорную фазу. Другие изображения 72 вентиляции смещают так, чтобы неподвижные опорные точки в изображенной зоне находились в одной и той же позиции в каждом из совмещенных изображений 73. Опорные точки могут быть указаны пользователем на всех изображениях. В альтернативном варианте пользователь указывает одну опорную точку на одном изображении, и процессор 15 автоматически производит поиск соответствующих точек на изображениях, соответствующих разным фазам. Кроме того, одна опорная точка, подлежащая отысканию, может быть найдена процессором с использованием поискового алгоритма и некоторых предварительно заданных критериев поиска. В результате, каждое изображение 72 вентиляции задается в системе координат опорной фазы, тогда как, перед совмещением, i-тое изображение вентиляции было связано с системой координат, принадлежащей i-той фазе. В альтернативном варианте совмещение изображений можно выполнять перед формированием векторных полей деформаций. В таком случае, КТ-изображения 71 совмещаются, и совмещение изображений 72 вентиляции больше не требуется. В обоих случаях конечным результатом будет набор из совмещенных изображений 73 вентиляции.
На этапе 26 формирования диаграммы формируется, по меньшей мере, одна диаграмма для представления изменения во времени значения вентиляции для конкретной позиции в обследуемом легком. С данной целью вентиляцию легких в разных фазах можно анализировать для, по меньшей мере, одной стандартной позиции. В альтернативном варианте пользователь может выбрать, по меньшей мере, одну позицию, для которой пользователь имеет намерение просмотреть изменение вентиляции легких во времени.
На фиг. 3 приведена последовательность КТ-изображений 71 пары легких. На данной фигуре показано десять изображений, где каждое изображение представляет отличающуюся фазу дыхательного цикла. Несмотря на то что на КТ-изображениях 71, возможно, видно еще не очень четко, в данном случае, фаза 0% соответствует состоянию максимального вдоха, при этом состояния 10%, 20%, 30% и 40% представляют фазы во время выдоха. Состояние 50% представляет переход от выдоха к вдоху. В состояния от 60% до 90%, легкие расширяются по причине вдоха. На максимальном вдохе дыхательный цикл завершается, и легкие опять находятся в состоянии 0%. Следует отметить, что, в принципе, любой момент дыхательного цикла можно задать, как состояние 0%. Подходящей альтернативой состоянию максимального вдоха может быть состояние максимального выдоха. В данном примере существует около 5 состояний выдоха и 4 состояний вдоха. Для других пациентов или для других измерений одного и того же пациента, отношение фаз вдоха к фазам выдоха может быть другим.
На фиг. 4 приведена последовательность векторных полей деформаций, полученная по КТ-изображениям 71, показанным на фиг. 3. В данном примере отображена нижняя - верхняя составляющая. Для каждой фазы определяется критерий измерения смещения в нижнем - верхнем направлении. В предпочтительном варианте, смещение измеряют в 3 направлениях (x, y, z), для получения возможности выполнения 3-мерного анализа вентиляции легких. Смещение определяют с использованием алгоритма распознавания образов. В фазах 10%-50%, наибольшее смещение происходит в верхнем направлении. Более светлые зоны изображения представляют большие смещения, чем более темные зоны. Упомянутые пять фаз представляют фазу выдоха дыхательного цикла. Выдох вызывается, в основном, расслаблением дыхательной мускулатуры и перемещением вверх диафрагмы, что приводит к нажиму вверх на легочную ткань и выталкиванию воздуха из легких. В фазах 60%-100%/0%, вдох приводит к смещению вниз легочной ткани. В данном случае, более темные зоны представляют большие смещения. Движение диафрагмы в нижнем направлении вызывает всасывание воздуха в легкие. Применение цветокодирования может быть полезно для обозначения величины и направления смещения.
На фиг. 5 приведена последовательность изображений 72 вентиляции легких, полученных из векторных полей деформаций, показанных на фиг. 4. Данные изображения 72 показывают, для каждой позиции в исследуемом легком, значение расширения/сокращения. Например, значение 0 соответствует сохранению объема, а значение меньше (больше) 0 указывает на сокращение (расширение), что можно отображать разными цветами.
На фиг. 6 изображен примерный дисплей 17 системы в соответствии с изобретением. На дисплее 17 показано одно из КТ-изображений 71. КТ-изображение 71 представляет два легких. Показанное КТ-изображение представляет одну из фаз дыхательного цикла. Чтобы пользователь мог выбирать изображение 71, представляющее другую фазу, может быть обеспечен пользовательский интерфейс. В альтернативном варианте, одновременно могут представляться, по меньшей мере, два изображения 71. Вместо КТ-изображения 71, возможно также представление изображения 72 вентиляции. Для четырех позиций 62 на КТ-изображении 71 отображается соответствующая диаграмма 74. Диаграмма 74 показывает изменение во времени значения вентиляции для соответствующей позиции 62. Вентиляционное смещение является таким, что во время выдоха значение вентиляции меньше 1, и во время вдоха значение вентиляции больше 1. Значение вентиляции характеризует объем легких. Щелчком мыши на позициях на КТ-изображении, пользователь может вызвать диаграммы 74, показывающие значения вентиляции для соответствующих позиций. По желанию, отображаемые диаграммы 74 могут показывать указатель 61 фазы текущего отображаемого КТ-изображения 71 таким образом, чтобы пользователь видел, какая точка на диаграмме соответствует отображаемому изображению 71. Путем манипуляции данным указателем 61, пользователь может выбирать КТ-изображения других фаз для отображения, вместо выбора КТ-изображений для отображения в дополнение к текущему отображаемому изображению.
На фиг. 7 представлена диаграмма 74 сформированной функции 81 вентиляции, вместе с подобранной моделью 82 дыхательного цикла. Модель 82 дыхательного цикла является математическим представлением вентиляции легких в виде функции времени. Например, модель 82 дыхательного цикла можно задать в виде
V ( t ) = V 0 + b cos 2 n ( c 1 ( t φ ) + c 2 )
Figure 00000001
где:
- c 1
Figure 00000002
задано продолжительностью дыхательного цикла,
- c 2
Figure 00000003
означает начальный сдвиг фазы вентиляции, при t = φ
Figure 00000004
,
- V 0
Figure 00000005
означает вентиляционное смещение,
- b
Figure 00000006
означает амплитуду вентиляции,
- φ
Figure 00000007
означает момент начала отсчета времени, и
- характерное значение n
Figure 00000008
равно 1.
Для подбора под измеренный параметр могут быть пригодны различные модели 82. Выбор математической модели 82 может зависеть, например, от измеренного параметра и/или требуемой точности подбора. Значения параметров подобранной модели 82 могут вычисляться и отображаться для выбранных позиций в зоне легкого. Упомянутые параметры подобранной модели могут содержать информацию, полезную для лица, работающего с системой. Например, амплитуда вентиляции и сдвиг фазы вентиляции могут содержать ценную информацию для оценки физического состояния соответствующей зоны легких.
Для визуализации значений параметра подобранной модели 82 во множестве позиций в легком, можно вычислять и отображать атлас 90 легких. Например, цветокодированное изображение легких может показывать амплитуды вентиляции или сдвиг фазы вентиляции в клинически релевантных позициях на изображении легких. На фиг. 8 показан примерный атлас 90 легких, где каждый пиксель соответствует позиции в легком в некоторый момент дыхательного цикла. Параметр является амплитудой вентиляции. Цвет каждого пикселя указывает значение параметра. Подобные атласы легких можно формировать также для других параметров и областей легких.
Атлас легких может показывать информацию об изменении параметра вентиляции легких в течение периода времени в выбранных позициях в легком. На подобном атласе легких можно показать развитие состояния вентиляции, например, во время лечения или вследствие развивающегося заболевания. В альтернативном варианте, атлас легких можно формировать с использованием параметров вентиляции легких нескольких пациентов. Подобный атлас легких может обеспечивать статистическую информацию о функционировании легких определенных групп пациентов. Например, атлас легких можно обеспечивать для здоровых людей, курильщиков, некурящих, людей с раком легких и т.п., при этом отображаемое значение параметра может быть математическим ожиданием или среднеквадратичным отклонением параметра.
На фиг. 9 показано КТ-изображение 71 легкого, на котором выделены изображения характерных зон 91, 92, 93. На фиг. 10 показано изменение во времени параметра вентиляции в выделенных зонах 91, 92, 93 фиг. 9 за десятинедельный период. Характерные зоны 91, 92, 93 могут быть выбраны пользователем. Например, пользовательский интерфейс системы может давать пользователю возможность очерчивать выбранные зоны указательным устройством, например мышью или джойстиком. В дополнение к зонам, подлежащим контролю, пользователь может также выбирать, по меньшей мере, один параметр для контроля и/или шкалу времени контролируемой диаграммы. Зоны 91, 92, 93 могут выбираться также автоматически с использованием методов распознавания изображений. Например, система может распознавать и выбирать зону 91 опухоли. Затем вокруг зоны 91 опухоли выбирается зона 92 окрестности опухоли. По меньшей мере, одну опорную область 93 находят в другом легком или в том же легком, на расстоянии от зоны опухоли. В предпочтительном варианте, опорная область 93 содержит только здоровую ткань легкого.
На фиг. 10 представлен параметр вентиляции легких, амплитуда вентиляции, для выбранных областей 91, 92, 93. Значение параметра легкого для каждой области может быть средним значением из всех значений в области, максимальным значением в области и т.п. Следует отметить, что график может содержать параметры вентиляции для областей, которые не видны на КТ-изображении на фиг. 9, например, по той причине, что данные области расположены на отличающейся глубине. График на фиг. 10 показывает, что вентиляция легких в зоне 91 опухоли снижается в течение десятинедельного лечения, в то время как вентиляция легких в зоне 92 окрестности опухоли снижается лишь незначительно, и вентиляция легких в опорных областях 93 остается, по существу, постоянной.
Следует понимать, что изобретение относится также к компьютерным программам, в частности, к компьютерным программам на или в носителе, предназначенным для применения изобретения на практике. Программа может быть в форме исходного кода, объектного кода, источника промежуточного кода и объектного кода, например, в частично компилированной форме, или в любой другой форме, пригодной для использования в реализации способа в соответствии с изобретением. Кроме того, следует понимать, что упомянутая программа может иметь много разных архитектурных проектов. Например, управляющая программа, реализующая функции способа или системы в соответствии с изобретением, может быть подразделена на, по меньшей мере, одну подпрограмму. Специалисту будут очевидны много разных способов распределения функций между упомянутыми подпрограммами. Подпрограммы могут храниться вместе в одном исполняемом файле, в форме независимой программы. Данный исполняемый файл может содержать исполняемые компьютером команды, например команды процессора и/или команды программы-интерпретатора (например, команды программы-интерпретатора Java). В альтернативном варианте, одна, несколько или все подпрограммы могут храниться в, по меньшей мере, одном внешнем библиотечном файле и могут подключаться к основной программе как статически, так и динамически, например, во время выполнения. Основная программа содержит, по меньшей мере, одно обращение к, по меньшей мере, одной из подпрограмм. Кроме того, подпрограммы могут содержать обращения к функциям между собой. Вариант осуществления, относящийся к компьютерному программному продукту, содержит исполняемые компьютером команды, соответствующие каждому из этапов обработки данных, по меньшей мере, одного из вышеописанных способов. Упомянутые команды могут подразделяться на подпрограммы и/или храниться в, по меньшей мере, файле, который может подключаться статически или динамически. Другой вариант осуществления, относящийся к компьютерному программному продукту, содержит исполняемые компьютером команды, соответствующие каждому из средств из, по меньшей мере, одной/ого из вышеописанных систем и/или продуктов. Упомянутые команды могут быть подразделены на подпрограммы и/или храниться в, по меньшей мере, одном файле, который может подключаться статически или динамически.
Носитель компьютерной программы может быть любым объектом или устройством, способным содержать программу. Например, носитель может содержать среду для хранения, например устройство ROM (постоянное запоминающее устройство, например, CD ROM (ROM на компакт-диске), или полупроводниковое устройство ROM, или магнитный носитель записи, например, гибкий диск или жесткий диск. Кроме того, носитель может быть передающим носителем, например, электрическим или оптическим сигналом, который может передаваться по электрическому или оптическому кабелю или по радио, или другими средствами. Когда программа осуществлена в виде данного сигнала, то носитель может быть образован упомянутым кабелем или другим устройством или средством. В альтернативном варианте, носитель может быть интегральной схемой, в которую заложена программа, при этом микросхема предназначена для выполнения или для применения при выполнении соответствующего способа.
Следует отметить, что вышеупомянутые варианты осуществления иллюстрируют, а не ограничивают изобретение, и что специалисты в данной области техники смогут разработать множество альтернативных вариантов осуществления, не выходящих за пределы объема притязаний приложенной формулы изобретения. В формуле изобретения никакие условные обозначения, расположенные в скобках, не следует толковать как ограничивающие формулу изобретения. Применение глагола «содержать» и его спряжений не исключает присутствия элементов или этапов, отличающихся от элементов или этапов, заявленных в формуле изобретения. Элемент в единственном числе не исключает присутствия множества данных элементов. Изобретение можно реализовать с помощью аппаратного обеспечения, содержащего несколько отдельных элементов, и с помощью соответственно запрограммированного компьютера. В пункте формулы изобретения, касающемся устройства и перечисляющем несколько средств, некоторые из упомянутых средств могут быть осуществлены посредством одного и того же аппаратного изделия. Очевидное обстоятельство, что некоторые средства упомянуты во взаимно различающихся зависимых пунктах формулы изобретения, не означает, что в подходящих случаях нельзя использовать сочетание упомянутых средств.

Claims (11)

1. Система для отображения информации о вентиляции легких, при этом система содержит:
- устройство (12) ввода для получения множества КТ-изображений (71) легкого, причем каждое КТ-изображение (71) соответствует одной фазе из, по меньшей мере, двух разных фаз в дыхательном цикле,
- процессор (15), сконфигурированный с возможностью:
- сравнения КТ-изображений (71), соответствующих разным фазам в дыхательном цикле для определения векторного поля деформаций для каждой фазы,
- формирования для каждой фазы изображения (72) вентиляции на основании соответствующего векторного поля деформаций,
- пространственного совмещения изображений (72) вентиляции,
- формирования для, по меньшей мере, одной общей позиции (62) в каждом из совмещенных изображений (72) вентиляции, функции (81) изменения во времени значения вентиляции для упомянутой общей позиции (62), причем каждое значение вентиляции в функции (81) основано на векторных полях деформаций, соответствующих совмещенным изображениям (73) вентиляции, и
- подбора модели дыхательного цикла для функции (81).
2. Система по п. 1, в которой процессор (15) дополнительно сконфигурирован с возможностью получения амплитуды локальной вентиляции и/или сдвига фазы локальной вентиляции из подобранной модели дыхательного цикла.
3. Система по п. 1, в которой процессор (15) дополнительно сконфигурирован с возможностью формирования функции (81) несколько раз, например, в течение курса лечения, и выполнения анализа трендов на основании множества сформированных функций (81).
4. Система по п. 1, в которой процессор (15) дополнительно сконфигурирован с возможностью формирования функции (81) для множества объектов обследования и выполнения статистического анализа на основании множества сформированных функций (81).
5. Система по п. 1, дополнительно содержащая дисплей (17) для отображения, по меньшей мере, одного из изображений вентиляции или КТ-изображений (72, 71) и диаграммы (74) функции (81).
6. Система по п. 5, дополнительно содержащая средство (18) пользовательского ввода, связанное с процессором (15), для предоставления пользователю возможности выбора, по меньшей мере, одной общей позиции (62) на отображаемом изображении вентиляции или КТ-изображении (72, 71).
7. Система по п. 6, в которой процессор (15) сконфигурирован для обеспечения пользователя возможностью применения средства (18) пользовательского ввода для выбора одного из изображений вентиляции или КТ-изображений (72, 71), которое следует отобразить.
8. Система по п. 6, в которой процессор (15) сконфигурирован для обеспечения пользователя возможностью выбора опорной точки совмещения на отображаемом изображении (73) вентиляции и пространственного совмещения изображений (72) вентиляции относительно опорной точки совмещения.
9. Система по п. 1, дополнительно содержащая КТ-сканер (11) для получения КТ-изображений (71).
10. Способ отображения информации о вентиляции легких, при этом способ содержит этапы, на которых:
- (22) получают множество КТ-изображений (71) легкого, причем каждое КТ-изображение (71) соответствует одной фазе из, по меньшей мере, двух разных фаз дыхательного цикла,
- (23) сравнивают КТ-изображения (71), соответствующие разным фазам в дыхательном цикле, для определения векторного поля деформаций для каждой фазы,
- (24) формируют для каждой фазы, изображение (72) вентиляции на основании соответствующего векторного поля деформаций,
- (25) осуществляют пространственное совмещение изображений (72) вентиляции,
- (26) формируют для, по меньшей мере, одной общей позиции (62) на каждом из совмещенных изображений (73) вентиляции, функцию (81) изменения во времени значения вентиляции для упомянутой общей позиции (62), причем каждое значение вентиляции в функции (81) основано на векторных полях деформаций, соответствующих совмещенным изображениям (73) вентиляции,
- подбирают модель дыхательного цикла для функции (81).
11. Машиночитаемый носитель, содержащий исполняемые процессором команды, предназначенные предписывать процессору выполнять способ по п. 10.
RU2011135966/14A 2009-01-30 2010-01-25 Система для представления информации о вентиляции легких RU2542096C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP09151764.9 2009-01-30
EP09151764 2009-01-30
PCT/IB2010/050305 WO2010086776A1 (en) 2009-01-30 2010-01-25 System for providing lung ventilation information

Publications (2)

Publication Number Publication Date
RU2011135966A RU2011135966A (ru) 2013-03-10
RU2542096C2 true RU2542096C2 (ru) 2015-02-20

Family

ID=42077897

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2011135966/14A RU2542096C2 (ru) 2009-01-30 2010-01-25 Система для представления информации о вентиляции легких

Country Status (6)

Country Link
US (1) US8798343B2 (ru)
EP (1) EP2391271B1 (ru)
CN (1) CN102300503B (ru)
BR (1) BRPI1005358A2 (ru)
RU (1) RU2542096C2 (ru)
WO (1) WO2010086776A1 (ru)

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX2014002153A (es) 2011-08-30 2014-05-01 Koninkl Philips Nv Integracion de entradas de usuario y correcion de campo de vector de deformacion en la dinamica de trabajo del registro deformable de imagenes.
JP5844187B2 (ja) * 2012-03-23 2016-01-13 富士フイルム株式会社 画像解析装置および方法並びにプログラム
EP2829231B1 (en) * 2012-03-23 2018-08-01 Konica Minolta, Inc. Image-generating apparatus
JP2015527108A (ja) * 2012-06-27 2015-09-17 コーニンクレッカ フィリップス エヌ ヴェ 画質駆動型非剛体画像レジストレーション
JP2015521880A (ja) * 2012-06-27 2015-08-03 コーニンクレッカ フィリップス エヌ ヴェ 運動パラメータ推定
US20160019680A1 (en) * 2013-03-29 2016-01-21 Koninklijke Philips N.V. Image registration
US9943703B2 (en) 2014-07-28 2018-04-17 The University Of Maryland, Baltimore System and method for irradiation therapy using voxel based functional measurements of organs at risk
US9962086B2 (en) * 2015-03-31 2018-05-08 Toshiba Medical Systems Corporation Medical image data processing apparatus and method for determining the presence of an abnormality
JP6565422B2 (ja) * 2015-07-24 2019-08-28 富士通株式会社 画像処理プログラム、画像処理装置及び画像処理方法
CN108369642A (zh) * 2015-12-18 2018-08-03 加利福尼亚大学董事会 根据头部计算机断层摄影解释和量化急症特征
US20170329927A1 (en) * 2016-05-11 2017-11-16 InFluidS LLC System and method for analyzing airway-pulmonary response using computational fluid dynamics to diagnose and monitoring potential health anomalies
JP6686733B2 (ja) * 2016-06-23 2020-04-22 コニカミノルタ株式会社 動態解析システム
JP2018064848A (ja) * 2016-10-21 2018-04-26 コニカミノルタ株式会社 動態解析システム
US10549117B2 (en) 2017-07-21 2020-02-04 Varian Medical Systems, Inc Geometric aspects of radiation therapy planning and treatment
US10609806B2 (en) 2017-07-21 2020-03-31 Varian Medical Systems Particle Therapy Gmbh Energy modulation of a cyclotron beam
US10245448B2 (en) 2017-07-21 2019-04-02 Varian Medical Systems Particle Therapy Gmbh Particle beam monitoring systems and methods
US10183179B1 (en) 2017-07-21 2019-01-22 Varian Medical Systems, Inc. Triggered treatment systems and methods
US10843011B2 (en) 2017-07-21 2020-11-24 Varian Medical Systems, Inc. Particle beam gun control systems and methods
US10092774B1 (en) 2017-07-21 2018-10-09 Varian Medical Systems International, AG Dose aspects of radiation therapy planning and treatment
US11712579B2 (en) 2017-07-21 2023-08-01 Varian Medical Systems, Inc. Range compensators for radiation therapy
US11590364B2 (en) 2017-07-21 2023-02-28 Varian Medical Systems International Ag Material inserts for radiation therapy
US12048526B2 (en) 2017-12-15 2024-07-30 Respiratory Motion, Inc. Devices and methods of calculating and displaying continuously monitored tidal breathing flow-volume loops (TBFVL) obtained by non-invasive impedance-based respiratory volume monitoring
EP3545845B1 (en) * 2018-03-29 2020-09-09 Medizinische Hochschule Hannover Method for processing computed tomography imaging data of a suspect`s respiratory system
JP7345282B2 (ja) 2019-06-03 2023-09-15 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
CN111292309B (zh) * 2020-02-13 2023-08-18 中国人民解放军国防科学技术大学 一种肺组织异化程度判断方法及装置
US20230107737A1 (en) * 2021-10-05 2023-04-06 Koninklijke Philips N.V. Inter-and extrapolation of chest image and mechanical ventilation settings into a time lapse series for progression monitoring and outcome prediction during long term mechanical ventilation

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005028121A (ja) * 2003-06-17 2005-02-03 Iyo Gazo Kenkyusho:Kk 呼吸気ct画像による換気分布計測方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5790690A (en) * 1995-04-25 1998-08-04 Arch Development Corporation Computer-aided method for automated image feature analysis and diagnosis of medical images
AU7586698A (en) * 1997-05-23 1998-12-11 William Beaumont Hospital Method and apparatus for delivering radiation therapy during suspended ventilation
US6501981B1 (en) * 1999-03-16 2002-12-31 Accuray, Inc. Apparatus and method for compensating for respiratory and patient motions during treatment
JP4294881B2 (ja) * 2000-05-12 2009-07-15 富士フイルム株式会社 画像の位置合わせ方法および装置
US6915151B2 (en) * 2001-02-08 2005-07-05 Trustees Of The University Of Pennsylvania Quantitative pulmonary imaging
US7106891B2 (en) * 2001-10-15 2006-09-12 Insightful Corporation System and method for determining convergence of image set registration
JP4493408B2 (ja) * 2003-06-06 2010-06-30 富士フイルム株式会社 画像読影支援方法及び装置並びにプログラム
US7359535B2 (en) * 2003-06-20 2008-04-15 Ge Medical Systems Global Technology Company, Llc Systems and methods for retrospective internal gating
DE102004017478B4 (de) * 2004-04-08 2012-01-19 Siemens Ag Vorrichtung für die Gewinnung von Strukturdaten eines sich bewegenden Objekts
EP1744671B1 (en) * 2004-04-26 2013-10-23 David F. Yankelevitz Medical imaging system for accurate measurement evaluation of changes in a target lesion
US8447380B2 (en) 2004-05-17 2013-05-21 Siemens Aktiengesellschaft Method for diagnosis of functional lung illnesses
US7903849B2 (en) * 2005-04-15 2011-03-08 Kabushiki Kaisha Toshiba Medical image processing apparatus
US7991449B2 (en) * 2005-04-26 2011-08-02 Mayo Foundation For Medical Education And Research Imaging elastic properties of the lung with magnetic resonance elastography
US7668357B2 (en) * 2005-10-17 2010-02-23 Stanford University Method and system for using computed tomography to test pulmonary function
CN101512602B (zh) * 2006-08-29 2012-12-05 皇家飞利浦电子股份有限公司 胸部ct成像中心脏运动伪影的减少
US7835493B2 (en) * 2007-08-06 2010-11-16 Stanford University Method and system for four dimensional intensity modulated radiation therapy for motion compensated treatments
US8189885B2 (en) * 2008-02-15 2012-05-29 The University Of Iowa Research Foundation Apparatus and method for computing regional statistical distribution over a mean anatomic space
US8538111B2 (en) * 2008-04-21 2013-09-17 University Of South Florida Method and apparatus for pulmonary ventilation imaging using local volume changes
US8755635B2 (en) * 2008-08-11 2014-06-17 Siemens Aktiengesellschaft Method and system for data dependent multi phase visualization

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005028121A (ja) * 2003-06-17 2005-02-03 Iyo Gazo Kenkyusho:Kk 呼吸気ct画像による換気分布計測方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
П.М. Котляров "Лучевые методы в диагностике заболеваний органов дыхания" Русский медицинский журнал. ПУЛЬМОНОЛОГИЯ. ДИАГНОСТИКА. ОНКОЛОГИЯ. 03 марта 2001 г, N 5, с.197-201 *

Also Published As

Publication number Publication date
RU2011135966A (ru) 2013-03-10
EP2391271A1 (en) 2011-12-07
US20110286652A1 (en) 2011-11-24
EP2391271B1 (en) 2015-11-04
US8798343B2 (en) 2014-08-05
CN102300503B (zh) 2014-06-18
BRPI1005358A2 (pt) 2019-09-24
CN102300503A (zh) 2011-12-28
WO2010086776A1 (en) 2010-08-05

Similar Documents

Publication Publication Date Title
RU2542096C2 (ru) Система для представления информации о вентиляции легких
CN110167438B (zh) 用于将体内电读数重建成解剖结构的系统和方法
JP5534607B2 (ja) 医用画像処理装置及び、そのプログラム
US10952634B2 (en) Electrical impedance tomography system
JP5694323B2 (ja) オブジェクトデータの生成
EP2925216B1 (en) Stenosis therapy planning
EP2867861B1 (en) Motion parameter estimation
US20100208957A1 (en) Respiratory Motion Compensated Cardiac Wall Motion Determination System
US20210259574A1 (en) Systems, methods, and computer-readable media of estimating thoracic cavity movement during respiration
US8487933B2 (en) System and method for multi-segment center point trajectory mapping
EP3468668B1 (en) Soft tissue tracking using physiologic volume rendering
JP5415245B2 (ja) 医用画像表示装置および方法並びにプログラム
US20160180520A1 (en) Quantitative method for 3-d joint characterization
RU2565521C2 (ru) Обработка набора данных изображения
EP2038846B1 (en) Model-based determination of the contraction status of a periodically contracting object
JP7173747B2 (ja) 医用画像処理装置
JP2021083479A (ja) 生体情報表示装置、生体情報表示方法及び表示プログラム
JP2024500828A (ja) 末梢灌流測定
KR20140120234A (ko) 혈관의 직경 측정 방법
US10441192B2 (en) Impedance shift detection
JP7547950B2 (ja) 動態解析装置及びプログラム
Cohen et al. Detection of recurrent nerve paralysis: development of a computer aided diagnosis system
US20220084206A1 (en) Information processing apparatus, information processing method, and non-transitory storage medium