RU2401441C2 - Реконструкция в позитронной эмиссионной томографии в режиме времяпролетного списка с использованием функции отклика детектора - Google Patents
Реконструкция в позитронной эмиссионной томографии в режиме времяпролетного списка с использованием функции отклика детектора Download PDFInfo
- Publication number
- RU2401441C2 RU2401441C2 RU2008110176/28A RU2008110176A RU2401441C2 RU 2401441 C2 RU2401441 C2 RU 2401441C2 RU 2008110176/28 A RU2008110176/28 A RU 2008110176/28A RU 2008110176 A RU2008110176 A RU 2008110176A RU 2401441 C2 RU2401441 C2 RU 2401441C2
- Authority
- RU
- Russia
- Prior art keywords
- matching pair
- probability
- function
- flight
- time
- Prior art date
Links
- 238000005316 response function Methods 0.000 title claims abstract description 20
- 238000002600 positron emission tomography Methods 0.000 title description 19
- 238000000034 method Methods 0.000 claims abstract description 16
- 230000005855 radiation Effects 0.000 claims description 23
- 238000003384 imaging method Methods 0.000 claims description 15
- 230000005251 gamma ray Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims 2
- 238000012800 visualization Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 229920013655 poly(bisphenol-A sulfone) Polymers 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012879 PET imaging Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine (AREA)
- Image Analysis (AREA)
Abstract
Изобретение относится к усовершенствованному алгоритму реконструкции для процессов визуализации. РЕТ-система содержит усовершенствованный алгоритм реконструкции изображения на основе усовершенствованного моделирования функции рассеяния точки. Данные о времени пролета в PET используют для получения средней точки излучения и вероятностной функции времени пролета. Затем упомянутую информацию используют для моделирования функции рассеяния точки. Вероятностная функция времени пролета и функция отклика детектора применяются для определения объема вероятностей для данной линии срабатывания, которая затем применяется при реконструкции изображения. Технический результат - повышение разрешающей способности и общего качества изображения. 3 н. и 17 з.п. ф-лы, 4 ил.
Description
В области медицинской визуализации все шире используют позитронную эмиссионную томографию (PET). В визуализирующей РЕТ-системе в объект, подлежащий визуализации, вводят радионуклид. Радионуклид испускает позитроны. Когда испускаемый позитрон сталкивается с электроном, происходит акт, далее событие аннигиляции, в котором позитрон и электрон уничтожаются. В результате события аннигиляции порождаются два гамма-кванта, разлетающихся, по существу, под углом 180 градусов один от другого. Визуализирующие РЕТ-системы используют детекторы, расположенные напротив один другого, для детектирования гамма-квантов, испускаемых объектом. Обычно применяют кольцо из детекторов для регистрации гамма-квантов, приходящих под каждым углом.
Сбор большого числа событий формирует информацию, необходимую для оценки объекта методами томографической реконструкции. Упомянутые события, при привязке к детекторным элементам, формируют линии срабатывания (LOR), которые можно обрабатывать гистограммным методом в зависимости от их геометрических характеристик для определения проекций или синограмм, подлежащих реконструкции. События можно также вводить в изображение по отдельности.
Следовательно, фундаментальным элементом сбора данных и реконструкции изображения является LOR: линия, пересекающая апертуру для пациента в системе. Можно получать дополнительную информацию о местоположении события. Во-первых, известно, что, в ходе снятия отсчетов и реконструкции способность системы к реконструкции (определению координат) точки не является пространственно-инвариантной по полю видения; а именно, выше в центре и постепенно снижается к периферии. Для описания упомянутой характеристики обычно используют функцию рассеяния точки (PSF). Разработаны инструментальные средства для использования PSF в процессе реконструкции. Во-вторых, время пролета (TOF) или разность по времени между прилетами гамма-квантов к каждому детектору, участвующему в регистрации пары, можно использовать для ограничения по LOR, на которых событие могло произойти с большей вероятностью.
По своему характеру PSF является трехмерной и, обычно, изменяется от изотропной с размером 4-6 мм в центре сканера до сильно асимметрично эллипсоидальной, с размером, изменяющимся от 6 до 8 мм, на периферии. Информацию, касающуюся TOF, можно преобразовать в пространственную информацию (посредством преобразования скорости пролета двух гамма-квантов). Обычно, информация, касающаяся TOF, будет приводить к определению местонахождения события с точностью от 10-12 см до нескольких сантиметров с использованием очень быстродействующего сканера.
По существу, желательно создание способа и средства для использования всей имеющейся информации, касающейся LOR, во время реконструкции.
Настоящее изобретение относится к усовершенствованному алгоритму реконструкции для процессов визуализации. Визуализирующая система, например РЕТ-система, объединяет одновременно моделирование функции рассеяния точки детектора и знание TOF для повышения общего качества реконструированного изображения.
В одном варианте осуществления времяпролетные данные используют для определения средней точки излучения и вероятностной функции времени пролета. Средняя точка излучения и вероятностная функция времени пролета, в свою очередь, служат для определения функции отклика детектора. Затем вероятностную функцию времени пролета и функцию отклика детектора используют для определения объема вероятностей, например, эллипсоида. Элементы трехмерного изображения в объеме вероятностей применяются в алгоритме реконструкции для данной линии срабатывания. Указанный процесс повышает общее качество реконструированного изображения без необходимости в существенном увеличении времени обработки.
На прилагаемых чертежах, которые включены в настоящее описание и являются его частью, представлены варианты осуществления изобретения, которые, в сочетании с вышеприведенным общим описанием изобретения и нижеприведенным подробным описанием, служат для пояснения принципов настоящего изобретения. Специалисту в данной области техники должно быть ясно, что упомянутые наглядные варианты осуществления не предназначены для ограничения изобретения, а просто дают примеры, содержащие принципы изобретения.
Фигура 1 - изображение визуализирующей PET-системы.
Фигура 2 - изображение модели DRF (функции отклика детектора) для трехмерного (3-мерного) кольцевого РЕТ-сканера.
Фигура 3 - изображение данных с временной меткой, обеспечивающей среднюю точку излучения и (функцию вероятности времени пролета) probTOF.
Фигура 4 - изображение эллипсоида вероятностей для точки возникновения пары совпадающих гамма-квантов по модификации DRF, основанной на функции probTOF.
РЕТ-система, предлагаемая в настоящей заявке, обеспечивает способ и средство для точного моделирования PSF, чтобы PSF можно было использовать для обеспечения усовершенствованной реконструкции изображений в течение допустимого времени обработки. РЕТ-система моделирует PSF в сочетании с данными, собранными в режиме списка, о временах пролета (TOF) в PET. РЕТ-система фактически обеспечивает более точную модель PSF, а также ограничивает объем вычислений ограничением числа элементов трехмерного изображения, которые могут вносить вклад в точку возникновения зарегистрированных гамма-квантов.
На фигуре 1 представлена РЕТ-система 8, содержащая множество колец 10 детекторов, область 12 визуализации, объект 13, подлежащий визуализации, и консоль или корпус 14. Объект 13, подлежащий визуализации, внутри которого содержится радионуклид, испускающий позитроны, помещают в область 12 визуализации в корпусе 14. Когда позитроны испускаются из радионуклида, они сталкиваются с электронами, и при этом происходит событие 16 аннигиляции. Событие 16 аннигиляции порождает два гамма-кванта 17, 18, разлетающихся в противоположных направлениях вдоль линии 31 отклика. Затем гамма-кванты 17, 18 детектируются кольцами 10 детекторов, и сигнал передается в процессорный блок 54. РЕТ-система 8 может содержать любое число колец 10 детекторов, из которых пять (r1, r2, r3, r4, r5) показаны на фигуре 1.
Процессорный блок 54 использует информацию из детекторов 10 вместе с алгоритмами реконструкции для формирования изображения объекта. Алгоритм реконструкции включает в себя моделирование PSF для определения точки возникновения каждого зарегистрированного гамма-кванта. В двумерной PET отклик детектора является двумерной функцией Х и ε, как показано на фигуре 2, где Х является радиальной координатой линии срабатывания (LOR), и ε, которое равно tA/tB, является относительной тангенциальной координатой места излучения гамма-кванта вдоль LOR. PSF в центре поперечного поля видения имеют более узкий и центрированный профиль, чем PSF, которые находятся на периферии поперечного поля видения. PSF более или менее асимметрична относительно Х для периферических LOR. PSF является более узкой при ε=1 (точка излучения находится в середине LOR) и расширяется с уменьшением ε (точка излучения находится ближе к одному из детекторов). В этом отношении, неточное моделирование DRF создает проблемы качества изображения, которые становятся более серьезными на периферии поперечного поля видения.
Что касается трехмерной PET, PSF является четырехмерной функцией X, ε, Z' и β, как показано на фигуре 2. Как показано, Z' является осью, перпендикулярной к LOR в осевом направлении, β означает угол наклона LOR к оси, и плоскость X-Z' ортогональна к LOR.
Так как РЕТ-системы создают гамма-кванты, которые всегда разлетаются под углом 180 градусов один к другому, то для облегчения определения точки возникновения вдоль LOR можно применять времяпролетные (TOF) алгоритмы. Например, как показано на фигуре 3, поле 68 видения PET показано с LOR 70, найденной по паре совпадающих гамма-квантов. Исходя из разности времен регистрации совпадающей пары и известного диаметра поля 68 видения, можно определить среднюю точку 72 излучения вдоль LOR 70. Однако точку излучения невозможно определить абсолютно точно. Вместо этого точку излучения можно выразить в виде функции вероятности (probTOF) 75 для всех точек вдоль LOR 70. Как показано на фигуре 3, probTOF 75 напоминает функцию Гаусса с максимальной вероятностью в средней точке 72 излучения.
Данные от детекторов 10 могут подаваться в процессорный блок 54 и могут содержать временную метку для каждого зарегистрированного гамма-кванта. Временные метки могут служить для определения средней точки 72 излучения и probTOF 75. Среднюю точку 72 излучения и probTOF 75 можно использовать для более точного моделирования DRF. На фигуре 4 представлен алгоритм реконструкции, который использует среднюю точку 72 излучения и probTOF 75 для модификации DRF. Оси координат обозначены позицией 76. Регистрация совпадения дает LOR 70, и временные метки служат для определения средней точки 72 излучения и probTOF 75, при этом probTOF является функцией L. DRF 78 можно вычислить по средней точке 72 излучения, при этом DRF является функцией X, Z' при L=1. Для данной LOR 70 вносить вклад в LOR 70 могут только элементы трехмерного изображения внутри эллипсоида 80. Границы эллипсоида 80 вероятностей определяются функциями probTOF 75 и DRF 78. В частности, границы эллипсоида 80 вероятностей определяются ненулевой областью probTOF 75 в направлении LOR и ненулевой областью DRF 78 в плоскости X-Z', L=1. В альтернативном варианте объем 80 вероятностей может быть ограничен заданным участком probTOF, заданным участком DRF или заданным участком probTOF и заданным участком DRF. При изменении заданного(ных) участка(ков) изменяется форма объема вероятностей, и можно уменьшить число элементов трехмерного изображения в объеме вероятностей. Например, probTOF и/или DRF могут быть вычислены при одном или двух стандартных отклонениях от средней точки излучения.
Поскольку PSF может быть функцией относительной тангенциальной координаты ε, контур эллипсоида для каждого фиксированного L=1 может быть разным. Весовой множитель точки пересечения каждого элемента трехмерного изображения с LOR является простым произведением probTOF(L) и PFS в плоскости (X, Z', L=1). Поэтому в передней проекции, для данной LOR, все элементы трехмерного изображения в эллипсоиде 80 вероятностей должны проецироваться с весовым множителем их probTOF и DRF. В обратной проекции, для данной LOR, скорректированная LOR получает соответствующий весовой множитель и распределяется каждому элементу трехмерного изображения в пределах цилиндра для эллипсоида вероятностей. Следует отметить, что объем 80 вероятностей может иметь любую трехмерную форму, определяемую функциями probTOF и DRF, и, следовательно, предполагается, что объем настоящего изобретения не должен ограничиваться конкретной геометрической формой объема вероятностей.
Следовательно, с использованием данных TOF, полученных в режиме списка, для каждой измеренной LOR координата излучения совпадающих гамма-квантов известна и находится в пределах временного разрешения сканера. PSF моделируется для ограниченного диапазона элементов трехмерного изображения объекта вдоль LOR, исходя из средней точки 72 излучения и probTOF 75. Уменьшение объема необходимых вычислений значительно сокращает время реконструкции изображения для алгоритмов, которые точно моделируют DRP. Точное моделирование PSF допускает формирование изображений с намного более высокими разрешающей способностью и общим качеством.
Выше изобретение описано со ссылкой на, по меньшей мере, один предпочтительный вариант осуществления. Безусловно, специалистам, прочитавшим и изучившим настоящее описание, станут очевидными модификации и варианты. Например, предполагается, что настоящее изобретение не ограничено РЕТ-системами, но может включать в себя визуализирующие системы, использующие разные способы, например комбинированные системы РЕТ/СТ (компьютерной томографии) или систему PET/MR (ядерного магнитного резонанса). Предполагается охват всех подобных модификаций и изменений в той степени, насколько они находятся в пределах объема прилагаемой формулы изобретения или ее эквивалентов.
Claims (20)
1. Визуализирующая РЕТ-система, содержащая
(a) по меньшей мере, один детектор для регистрации гамма-квантов, излучаемых из подлежащего визуализации объекта; и
(b) процессорный блок, соединенный с упомянутым, по меньшей мере, одним детектором, при этом упомянутый процессорный блок содержит
i) средство для присвоения временных меток каждому событию регистрации гамма-кванта;
ii) средство для определения совпадающих пар по упомянутым временным меткам;
iii) средство для определения средней точки излучения и вероятностной функции времени пролета для каждой совпадающей пары;
iv) средство для определения функции рассеяния точки детектора по средней точке излучения и вероятностной функции времени пролета для каждой совпадающей пары;
v) средство для определения объема вероятностей по вероятностной функции времени пролета и функции отклика детектора для каждой совпадающей пары; и
vi) средство для использования объема вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения.
(a) по меньшей мере, один детектор для регистрации гамма-квантов, излучаемых из подлежащего визуализации объекта; и
(b) процессорный блок, соединенный с упомянутым, по меньшей мере, одним детектором, при этом упомянутый процессорный блок содержит
i) средство для присвоения временных меток каждому событию регистрации гамма-кванта;
ii) средство для определения совпадающих пар по упомянутым временным меткам;
iii) средство для определения средней точки излучения и вероятностной функции времени пролета для каждой совпадающей пары;
iv) средство для определения функции рассеяния точки детектора по средней точке излучения и вероятностной функции времени пролета для каждой совпадающей пары;
v) средство для определения объема вероятностей по вероятностной функции времени пролета и функции отклика детектора для каждой совпадающей пары; и
vi) средство для использования объема вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения.
2. Визуализирующая РЕТ-система по п.1, в которой упомянутый объем вероятностей является эллипсоидом.
3. Визуализирующая РЕТ-система по п.1, в которой средство для использования объема вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения содержит применение весового множителя точки пересечения вероятностной функции времени пролета и функции отклика детектора для каждой совпадающей пары.
4. Визуализирующая РЕТ-система по п.1, в которой объем вероятностей для каждой совпадающей пары частично ограничен ненулевой областью вероятностной функции времени пролета для каждой совпадающей пары.
5. Визуализирующая РЕТ-система по п.1, в которой объем вероятностей для каждой совпадающей пары частично ограничен ненулевой областью функции отклика детектора для каждой совпадающей пары.
6. Визуализирующая РЕТ-система по п.1, в которой объем вероятностей для каждой совпадающей пары ограничен ненулевой областью вероятностной функции времени пролета и ненулевой областью функции отклика детектора для каждой совпадающей пары.
7. Визуализирующая РЕТ-система по п.1, в которой объем вероятностей для каждой совпадающей пары ограничен заданным участком вероятностной функции времени пролета для каждой совпадающей пары.
8. Визуализирующая РЕТ-система по п.1, в которой объем вероятностей для каждой совпадающей пары ограничен заданным участком функции отклика детектора для каждой совпадающей пары.
9. Визуализирующая РЕТ-система по п.1, в которой объем вероятностей для каждой совпадающей пары ограничен заданным участком вероятностной функции времени пролета и заданным участком функции отклика детектора для каждой совпадающей пары.
10. Визуализирующая РЕТ-система по п.9, в которой, по меньшей мере, один из упомянутого заданного участка вероятностной функции времени пролета и упомянутого заданного участка функции отклика детектора для каждой совпадающей пары, приблизительно равен двум стандартным отклонениям от упомянутой средней точки излучения.
11. Визуализирующая РЕТ-система по п.1, дополнительно содержащая визуализирующую систему, использующую разные способы.
12. Способ формирования изображения в PET, содержащий следующие этапы:
(a) регистрируют излучения гамма-квантов из подлежащего визуализации объекта;
(b) присваивают временные метки каждому из упомянутых зарегистрированных гамма-квантов;
(c) определяют совпадающие пары гамма-квантов;
(d) определяют среднюю точку излучения и вероятностную функцию времени пролета для каждой совпадающей пары;
(e) определяют функцию отклика детектора на основании средней точки излучения и вероятностной функции времени пролета для каждой совпадающей пары;
(f) определяют объем вероятностей на основании вероятностной функции времени пролета и функцию отклика детектора для каждой совпадающей пары; и
(g) применяют объем вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения.
(a) регистрируют излучения гамма-квантов из подлежащего визуализации объекта;
(b) присваивают временные метки каждому из упомянутых зарегистрированных гамма-квантов;
(c) определяют совпадающие пары гамма-квантов;
(d) определяют среднюю точку излучения и вероятностную функцию времени пролета для каждой совпадающей пары;
(e) определяют функцию отклика детектора на основании средней точки излучения и вероятностной функции времени пролета для каждой совпадающей пары;
(f) определяют объем вероятностей на основании вероятностной функции времени пролета и функцию отклика детектора для каждой совпадающей пары; и
(g) применяют объем вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения.
13. Способ по п.12, в котором упомянутый объем вероятностей для каждой совпадающей пары является эллипсоидом.
14. Способ по п.12, в котором объем вероятностей для каждой совпадающей пары ограничен заданным участком вероятностной функции времени пролета для каждой совпадающей пары.
15. Способ по п.12, в котором объем вероятностей для каждой совпадающей пары ограничен заданным участком функции отклика детектора для каждой совпадающей пары.
16. Способ по п.12, в котором объем вероятностей для каждой совпадающей пары ограничен заданным участком вероятностной функции времени пролета и заданным участком функции отклика детектора для каждой совпадающей пары.
17. Способ по п.16, в котором, по меньшей мере, один из упомянутого заданного участка вероятностной функции времени пролета и упомянутого заданного участка функции отклика детектора для каждой совпадающей пары приблизительно равен двум стандартным отклонениям от упомянутой средней точки излучения.
18. Процессор для реконструкции изображения, содержащий
(a) средство для присвоения временных меток событию регистрации гамма-кванта;
(b) средство для определения совпадающих пар по упомянутым временным меткам;
(c) средство для определения средней точки излучения и вероятностной функции времени пролета для каждой совпадающей пары;
(d) средство для определения функции отклика детектора по средней точке излучения и вероятностной функции времени пролета для каждой совпадающей пары;
(e) средство для определения объема вероятностей по вероятностной функции времени пролета и функции отклика детектора для каждой совпадающей пары; и
(f) средство для использования объема вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения.
(a) средство для присвоения временных меток событию регистрации гамма-кванта;
(b) средство для определения совпадающих пар по упомянутым временным меткам;
(c) средство для определения средней точки излучения и вероятностной функции времени пролета для каждой совпадающей пары;
(d) средство для определения функции отклика детектора по средней точке излучения и вероятностной функции времени пролета для каждой совпадающей пары;
(e) средство для определения объема вероятностей по вероятностной функции времени пролета и функции отклика детектора для каждой совпадающей пары; и
(f) средство для использования объема вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения.
19. Процессор для реконструкции изображения по п.18, в котором упомянутый объем вероятностей для каждой совпадающей пары является эллипсоидом.
20. Процессор для реконструкции изображения по п.18, в котором средство для использования объема вероятностей для каждой совпадающей пары в алгоритме реконструкции для формирования изображения содержит применение весового множителя точки пересечения вероятностной функции времени пролета и функции отклика детектора для каждой совпадающей пары.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US59593005P | 2005-08-18 | 2005-08-18 | |
| US60/595,930 | 2005-08-18 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2008110176A RU2008110176A (ru) | 2009-09-27 |
| RU2401441C2 true RU2401441C2 (ru) | 2010-10-10 |
Family
ID=37496419
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2008110176/28A RU2401441C2 (ru) | 2005-08-18 | 2006-07-27 | Реконструкция в позитронной эмиссионной томографии в режиме времяпролетного списка с использованием функции отклика детектора |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US7365335B2 (ru) |
| EP (1) | EP1922566B1 (ru) |
| JP (1) | JP5069683B2 (ru) |
| CN (1) | CN101243331B (ru) |
| RU (1) | RU2401441C2 (ru) |
| WO (1) | WO2007020543A1 (ru) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2597074C2 (ru) * | 2011-06-16 | 2016-09-10 | Конинклейке Филипс Н.В. | Усовершенствование пространственной выборки для сбора данных пэт в виде списка, использующее спланированное перемещение стола/гентри |
| RU2717563C1 (ru) * | 2019-07-10 | 2020-03-24 | Акционерное общество "Научно-исследовательский институт технической физики и автоматизации" (АО "НИИТФА") | Способ определения функции рассеяния точки системы рентгеновской визуализации |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8314796B2 (en) * | 2006-02-24 | 2012-11-20 | The Board Of Trustees Of The Leland Stanford Junior University | Method of reconstructing a tomographic image using a graphics processing unit |
| DE102007023657B4 (de) * | 2007-05-22 | 2014-03-20 | Siemens Aktiengesellschaft | Verfahren zur Datenaufnahme bei einer funktionellen Gehirnuntersuchung mit einem kombinierten Magnetresonanz-PET-Gerät |
| GB2459870B (en) * | 2008-05-07 | 2012-06-20 | Petrra Ltd | Positron emission detection and imaging |
| US8326075B2 (en) | 2008-09-11 | 2012-12-04 | Google Inc. | System and method for video encoding using adaptive loop filter |
| US8767908B2 (en) | 2008-10-21 | 2014-07-01 | University Of Southern California | Exact and approximate rebinning of time-of-flight PET positron emission tomography data |
| JP2011002306A (ja) * | 2009-06-18 | 2011-01-06 | Shimadzu Corp | Ｐｅｔ装置の逐次近似画像再構成法 |
| JP2011252855A (ja) * | 2010-06-03 | 2011-12-15 | Toshiba Corp | 核医学イメージング装置 |
| EP2685284A4 (en) * | 2011-03-10 | 2014-11-05 | Shimadzu Corp | POSITRON TOMOGRAPHY DEVICE AND RECONSTRUCTION METHOD |
| US8780996B2 (en) | 2011-04-07 | 2014-07-15 | Google, Inc. | System and method for encoding and decoding video data |
| US9332952B2 (en) * | 2011-10-06 | 2016-05-10 | Koninklijke Philips N.V. | Data-driven optimization of event acceptance/rejection logic |
| GB2487017A (en) * | 2012-03-30 | 2012-07-04 | Petrra Ltd | Positron emission detection and imaging |
| CN103006252B (zh) * | 2012-12-10 | 2014-08-27 | 沈阳东软医疗系统有限公司 | 一种pet系统点扩散函数的测量方法及其测量装置 |
| EP3125765A4 (en) | 2014-03-04 | 2018-04-18 | The Trustees of Columbia University in the City of New York | Regularization of images |
| CN104408763B (zh) * | 2014-10-29 | 2017-08-11 | 沈阳东软医疗系统有限公司 | 一种图像重建方法及装置 |
| CN109480892B (zh) * | 2018-12-29 | 2022-11-04 | 上海联影医疗科技股份有限公司 | 一种图像的生成方法 |
| US20230260171A1 (en) * | 2022-02-17 | 2023-08-17 | Canon Medical Systems Corporation | Event property-dependent point spread function modeling and image reconstruction for pet |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6173083A (ja) * | 1984-09-17 | 1986-04-15 | クレイトン フアウンデイシヨン フオ− リサ−チ | 時経過陽電子放射カメラによる放射線検出方法及びその実施に使用するカメラ |
| CA1303256C (en) | 1988-12-14 | 1992-06-09 | Royal Institution For The Advancement Of Learning (The) | Scintillation crystals for positron emission tomography having a non reflecting band |
| US5323007A (en) * | 1992-02-07 | 1994-06-21 | Univ. Of Chicago Development Corp. Argonne National Laboratories | Method of recovering tomographic signal elements in a projection profile or image by solving linear equations |
| US6448560B1 (en) | 1996-02-02 | 2002-09-10 | Tumay O. Tumer | Method and apparatus for gamma ray detection |
| WO2001087140A2 (en) | 2000-05-16 | 2001-11-22 | Crosetto Dario B | Method and apparatus for anatomical and functional medical imaging |
| US6484051B1 (en) | 2001-05-15 | 2002-11-19 | James Daniel | Coincident multiple compton scatter nuclear medical imager |
| US7205546B1 (en) * | 2003-09-12 | 2007-04-17 | Kennedy Stephen W | Spherical PET mammography scanner |
| US7173247B2 (en) * | 2003-09-24 | 2007-02-06 | Radiation Monitoring Devices, Inc. | Lu1-xI3:Cex—a scintillator for gamma ray spectroscopy and time-of-flight PET |
| US7518114B2 (en) | 2004-11-10 | 2009-04-14 | General Electric Company | Method and system for imaging a patient |
| US7057178B1 (en) | 2004-11-15 | 2006-06-06 | General Electric Company | Method and system for imaging using a filter for Time-of-Flight PET |
| US7129495B2 (en) | 2004-11-15 | 2006-10-31 | General Electric Company | Method and apparatus for timing calibration in a PET scanner |
| US7301144B2 (en) | 2004-12-29 | 2007-11-27 | General Electric Company | Method and system for calibrating a positron emission tomography system |
| US7227149B2 (en) | 2004-12-30 | 2007-06-05 | General Electric Company | Method and system for positron emission tomography image reconstruction |
| US7129496B2 (en) | 2005-01-21 | 2006-10-31 | General Electric Company | Method and system for scattered coincidence estimation in a time-of-flight positron emission tomography system |
-
2006
- 2006-07-27 EP EP06780231A patent/EP1922566B1/en active Active
- 2006-07-27 JP JP2008526572A patent/JP5069683B2/ja not_active Expired - Fee Related
- 2006-07-27 CN CN200680030108.9A patent/CN101243331B/zh active Active
- 2006-07-27 WO PCT/IB2006/052583 patent/WO2007020543A1/en active Application Filing
- 2006-07-27 RU RU2008110176/28A patent/RU2401441C2/ru active
- 2006-08-15 US US11/464,534 patent/US7365335B2/en active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2597074C2 (ru) * | 2011-06-16 | 2016-09-10 | Конинклейке Филипс Н.В. | Усовершенствование пространственной выборки для сбора данных пэт в виде списка, использующее спланированное перемещение стола/гентри |
| RU2717563C1 (ru) * | 2019-07-10 | 2020-03-24 | Акционерное общество "Научно-исследовательский институт технической физики и автоматизации" (АО "НИИТФА") | Способ определения функции рассеяния точки системы рентгеновской визуализации |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2008110176A (ru) | 2009-09-27 |
| US20070040123A1 (en) | 2007-02-22 |
| WO2007020543A1 (en) | 2007-02-22 |
| CN101243331A (zh) | 2008-08-13 |
| JP2009505093A (ja) | 2009-02-05 |
| US7365335B2 (en) | 2008-04-29 |
| EP1922566A1 (en) | 2008-05-21 |
| JP5069683B2 (ja) | 2012-11-07 |
| EP1922566B1 (en) | 2013-02-13 |
| CN101243331B (zh) | 2012-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2401441C2 (ru) | Реконструкция в позитронной эмиссионной томографии в режиме времяпролетного списка с использованием функции отклика детектора | |
| JP5378787B2 (ja) | 粗い角度及びスライスリビニングによる三次元飛行時間ｐｅｔ | |
| JP2009505093A5 (ru) | ||
| EP3224801B1 (en) | Multi-modality imaging system and method | |
| CN106466188B (zh) | 用于发射断层显像定量的系统和方法 | |
| RU2620862C2 (ru) | Быстрая оценка рассеяния при реконструкции посредством позитронно-эмиссионной томографии | |
| US7417231B2 (en) | Fourier re-binning of time-of-flight positron emission tomography data | |
| US7402807B2 (en) | Method for reducing an electronic time coincidence window in positron emission tomography | |
| Markiewicz et al. | NiftyPET: a high-throughput software platform for high quantitative accuracy and precision PET imaging and analysis | |
| JP4208284B2 (ja) | 核像形成方法及び装置 | |
| RU2427037C2 (ru) | Способ и система для реконструкции рет-изображения с использованием суррогатного изображения | |
| JP2000180550A (ja) | Ｍｌ―ｅｍ画像再構成法及び医用画像形成装置 | |
| EP2984631B1 (en) | Method for modeling and accounting for cascade gammas in images | |
| CN101903908A (zh) | 基于特征的2d/3d图像配准 | |
| US10215864B2 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
| CN109887048B (zh) | Pet散射校正方法、图像重建方法、装置及电子设备 | |
| EP2080166A2 (en) | Dirty isotope pet reconstruction | |
| CN106415317B (zh) | 单光子发射计算机化断层摄影术中的多个发射能量 | |
| KR101766294B1 (ko) | 부호화 구경을 이용한 방사선원의 분포에 대한 3차원 영상 처리 장치 및 방법 | |
| JP2010204755A (ja) | 画像処理装置、画像再構成システム、画像処理方法およびプログラム | |
| US11898972B2 (en) | Imaging method using jointly a PET reconstruction and a compton reconstruction, preferably in 3D compton | |
| JP2003294842A (ja) | 放射線投影データの補正方法 | |
| US20200085397A1 (en) | Attenuation coefficient image estimation method, attenuation coefficient image estimation program, and positron ct apparatus equipped with the same | |
| US7890282B2 (en) | Estimation of crystal efficiency with axially compressed sinogram | |
| McLennan et al. | SORTEO: Monte Carlo-based simulator with list-mode capabilities |