RU2542096C2 - Система для представления информации о вентиляции легких - Google Patents
Система для представления информации о вентиляции легких Download PDFInfo
- 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
Links
- 238000009423 ventilation Methods 0.000 title claims abstract description 128
- 210000004072 lung Anatomy 0.000 title claims abstract description 83
- 230000000241 respiratory effect Effects 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000010586 diagram Methods 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims 1
- 239000003814 drug Substances 0.000 abstract 2
- 230000000694 effects Effects 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 206010028980 Neoplasm Diseases 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000008602 contraction Effects 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000003434 inspiratory effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 201000005202 lung cancer Diseases 0.000 description 1
- 208000020816 lung neoplasm Diseases 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 210000003019 respiratory muscle Anatomy 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 210000004887 upper abdominal cavity Anatomy 0.000 description 1
Images
Classifications
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices 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/5235—Devices 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
-
- 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/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- 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/10081—Computed x-ray tomography [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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
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 дыхательного цикла можно задать в виде
где:
Для подбора под измеренный параметр могут быть пригодны различные модели 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).
- устройство (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).
- (22) получают множество КТ-изображений (71) легкого, причем каждое КТ-изображение (71) соответствует одной фазе из, по меньшей мере, двух разных фаз дыхательного цикла,
- (23) сравнивают КТ-изображения (71), соответствующие разным фазам в дыхательном цикле, для определения векторного поля деформаций для каждой фазы,
- (24) формируют для каждой фазы, изображение (72) вентиляции на основании соответствующего векторного поля деформаций,
- (25) осуществляют пространственное совмещение изображений (72) вентиляции,
- (26) формируют для, по меньшей мере, одной общей позиции (62) на каждом из совмещенных изображений (73) вентиляции, функцию (81) изменения во времени значения вентиляции для упомянутой общей позиции (62), причем каждое значение вентиляции в функции (81) основано на векторных полях деформаций, соответствующих совмещенным изображениям (73) вентиляции,
- подбирают модель дыхательного цикла для функции (81).
11. Машиночитаемый носитель, содержащий исполняемые процессором команды, предназначенные предписывать процессору выполнять способ по п. 10.
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)
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)
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)
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 |
-
2010
- 2010-01-25 WO PCT/IB2010/050305 patent/WO2010086776A1/en active Application Filing
- 2010-01-25 US US13/146,885 patent/US8798343B2/en active Active
- 2010-01-25 RU RU2011135966/14A patent/RU2542096C2/ru active
- 2010-01-25 CN CN201080005802.1A patent/CN102300503B/zh active Active
- 2010-01-25 EP EP10703511.5A patent/EP2391271B1/en active Active
- 2010-01-25 BR BRPI1005358A patent/BRPI1005358A2/pt not_active IP Right Cessation
Patent Citations (1)
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)
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 |