EA016862B1 - Способ и устройство для оценки начального момента времени электронного импульса, предназначенного для использования в позитронно-эмиссионной томографии - Google Patents

Способ и устройство для оценки начального момента времени электронного импульса, предназначенного для использования в позитронно-эмиссионной томографии Download PDF

Info

Publication number
EA016862B1
EA016862B1 EA201070554A EA201070554A EA016862B1 EA 016862 B1 EA016862 B1 EA 016862B1 EA 201070554 A EA201070554 A EA 201070554A EA 201070554 A EA201070554 A EA 201070554A EA 016862 B1 EA016862 B1 EA 016862B1
Authority
EA
Eurasian Patent Office
Prior art keywords
pulse signal
analog
digital
time
curve
Prior art date
Application number
EA201070554A
Other languages
English (en)
Other versions
EA201070554A1 (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 EA201070554A1 publication Critical patent/EA201070554A1/ru
Publication of EA016862B1 publication Critical patent/EA016862B1/ru

Links

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
    • 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)
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • 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/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647—Processing of scintigraphic data
    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography

Landscapes

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

Abstract

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

Description

Перекрестная ссылка на родственную заявку
Данная заявка претендует на приоритет предварительной патентной заявки США №60985083, поданной 2 ноября 2007 г., содержание которой целиком включено в настоящий документ по ссылке.
Уровень техники
Возможность создания изображений внутренних областей живого организма без инвазивного хирургического вмешательства стала главным достижением в медицине за последнюю сотню лет. Технологии создания изображений, такие как рентгеновская компьютерная томография (СТ) и магниторезонансная интроскопия (ММ), дали врачам и ученым возможность наблюдать изображение анатомических структур внутри тела с высоким разрешением. Хотя это и привело к совершенствованию диагностики и лечения болезней, большой набор болезней вызывает изменения в анатомической структуре только на последних стадиях заболевания или не вызывает их вовсе. Это привело к становлению направления в области создания медицинских изображений, согласно которому фиксируется определенная метаболическая активность внутри живого тела. Позитронно-эмиссионная томография (РЕТ) относится именно к этому направлению в области создания медицинских изображений.
Позитронно-эмиссионная томография
РЕТ представляет собой способ создания медицинских изображений, в котором используют преимущества, которые дает анализ результатов радиоактивного распада для измерения определенной метаболической активности внутри живых организмов. Системы создания изображений типа РЕТ содержат три основных компонента, схематически показанных на фиг. 1: радиоактивный индикатор, который вводится в объект, подлежащий сканированию; сканер, предназначенный для обнаружения местоположения радиоактивного индикатора (косвенным образом, как обсуждается ниже); и систему обработки для создания томографических изображений.
Первым шагом является создание и введение радиоактивного индикатора 90, содержащего радиоактивный изотоп и метаболически активную молекулу. Индикатор 90 вводится в тело 91, подлежащее сканированию. По окончании времени, необходимого для концентрирования индикатора 90 в определенных тканях, тело 91 позиционируют надлежащим образом внутри сканера 92. Событием радиоактивного распада для индикаторов, используемых в исследованиях РЕТ, является позитронная эмиссия. Эмитированный позитрон проходит короткое расстояние в ткани тела, пока не вступит во взаимодействие с электроном. Позитронно-электронное взаимодействие во время события аннигиляции создает два антипараллельных фотона с энергией 511 килоэлектрон-вольт (кэВ). Сканер 92 приспособлен для обнаружения, по меньшей мере, некоторых фотонов, возникающих из события аннигиляции.
Сканер 92 (вторая компонента системы РЕТ) включает в себя кольцо из датчиков, которые обнаруживают фотоны с энергией 511 кэВ, и электронную аппаратуру предварительной обработки, которая обрабатывает сигналы, созданные датчиками. Датчики, как правило, содержат сцинтилляционные кристаллы или сцинтилляторы 93 и фотоэлектронные умножители (РМТ), кремниевые фотоумножители (δίΜΡ) или лавинные фотодиоды (ΑΡΌ) 94. Сцинтилляционный кристалл 93 преобразует фотоны высокой энергии 511 кэВ в множество фотонов низкой энергии, как правило, фотонов видимого света. РМТ, δίΜΡ или ΆΡΌ 94 обнаруживают фотоны видимого света и создают соответствующий электрический импульс. Импульсы РМТ обрабатываются электронной аппаратурой предварительной обработки для определения параметров или характеристик импульса (то есть энергии, временных параметров). Для удобства здесь ссылки на РМТ следует понимать так, что они включают в себя любой механизм или устройство для обнаружения фотонов высокой энергии, например фотонов с энергией 511 кэВ, и создания в ответ фотонов низкой энергии, таких как фотоны видимого света.
Наконец, данные посылают в главный компьютер 95, который выполняет реконструкцию томографического изображения, преобразуя данные в трехмерное (3-Ό) изображение.
Радиофармацевтика
Для синтезирования индикатора 90 к метаболически активной молекуле прикрепляют радиоактивный изотоп с малым периодом полураспада. Малый период полураспада сокращает экспонирование объекта ионизирующим излучением, но обычно при этом необходимо, чтобы индикатор 90 находился рядом со сканером. Наиболее часто используемым индикатором является фтор-18 фтордеоксиглюкоза ([Е18]ЕЭС) - аналог глюкозы, имеющий период полураспада 110 мин. |Е-18|ЕЭС аналогичен глюкозе, которая фосфорелирована клетками, использующими глюкозу, но не подвергнута гликолизу. Таким образом, радиоактивная часть молекулы захватывается в ткани. Клетки, которые потребляют много глюкозы, такие как раковые клетки и клетки мозга, аккумулируют со временем больше |Е-18|ЕЭС. чем другие ткани.
После того как прошло время, достаточное для получения исследуемой тканью достаточного количества индикатора 90, используют сканер 92 для обнаружения событий радиоактивного распада, то есть путем обнаружения фотонов с энергией 511 кэВ. Излученный позитрон, как правило, проходит несколько миллиметров ткани, прежде чем аннигилирует с электроном, создавая два фотона с энергией 511 кэВ, которые направлены под углом 180°±0,23° друг относительно друга.
Фотонная сцинтилляция
Фотон с энергией 511 кэВ обладает существенной энергией и может пройти сквозь множество ма
- 1 016862 териалов, включая ткань тела. Хотя это, как правило, позволяет фотону проходить через тело насквозь, фотоны с высокой энергией трудно обнаружить. Обнаружение фотонов является задачей сцинтиллятора 93. Сцинтиллятор 93 поглощает фотоны с высокой энергией и излучает фотоны с низкой энергией, как правило, фотоны видимого света. Сцинтиллятор 93 может быть выполнен из различных материалов, в том числе пластмассы, органических и неорганических кристаллов и органических жидкостей. Каждый тип сцинтиллятора имеет разную плотность, коэффициент преломления, временные характеристики и длину волны максимального излучения.
В общем случае плотность кристалла сцинтиллятора определяет, насколько хорошо материал задерживает фотоны с высокой энергией. Коэффициент преломления кристалла сцинтиллятора и длина волны излученного света предопределяют, насколько легко может быть уловлен свет от кристалла. Также необходимо, чтобы длина волны излученного света соответствовала устройству, которое будет преобразовывать свет в электрический импульс (например, РМТ), чтобы оптимизировать эффективность системы. Временные характеристики сцинтиллятора определяют, сколько времени понадобится для того, чтобы видимый свет достиг своего максимума (время нарастания), и длительность спада (время спада). Времена нарастания и спада являются важными характеристиками, поскольку чем больше сумма этих двух времен, тем меньше количество событий, которое детектор сможет обработать за заданный период, и, следовательно, тем дольше будет сканирование для получения одного и того же количества отсчетов. Также, чем больше указанные временные характеристики, тем больше вероятность того, что два события перекроются (наложатся друг на друга) и данные будут потеряны.
Примером современного материала для сцинтиллятора является Ьи2§1о5(Се), или Ь8О, который представляет собой неорганический кристалл. Ь8О, как сообщается, имеет постоянную нарастания 30 пкс и постоянную спада 40 нс. Эти значения могут незначительно изменяться из-за вариаций геометрии кристалла и электронной аппаратуры, с которой они соединены. Ь8О является современным материалом для сцинтилляторов, который отличается малым временем отклика и значительной светоотдачей.
Фотоэлектронные умножители
Сцинтиллятор 93 прикреплен к электронным устройствам, которые преобразуют фотоны видимого света от сцинтиллятора 93 в электронные импульсы. Наиболее часто используют два типа устройств: РМТ и ΆΡΌ. РМТ представляет собой вакуумную трубку с фотокатодом, несколькими динодами и анодом, имеющим высокий коэффициент усиления, позволяющий обнаруживать свет с очень низкими уровнями энергии. ΆΡΌ является полупроводниковой версией РМТ. Другой технологией, исследуемой в настоящее время с целью использования в сканерах ΡΕΤ, являются фотоумножители δίΡΜ. Фотоумножители δίΡΜ содержат матрицу полупроводниковых фотодиодов, которые работают в режиме Гейгера, так что при взаимодействии фотона и создании несущего сигнала возникает короткий импульс тока. В примерном δίΡΜ матрица фотодиодов содержит примерно 103 диода на мм2. Все диоды подсоединены к общей кремниевой подложке, так что выход матрицы равен сумме выходных сигналов всех диодов. Следовательно, выходной сигнал может принимать значение от минимального, когда включается один фотодиод, до максимального, когда включаются все фотодиоды. Это позволяет получить линейный выход, даже если устройства выполнены из цифровых приборов.
В примерной системе используется РМТ с двенадцатью каналами: шесть в направлении 'х' и шесть в направлении 'у', как показано на фиг. 2. Отдельные каналы позволяют более точно определять место события. Например, если событие обнаружено в левом верхнем углу РМТ, то тогда каналы Υ1 и XI будут иметь большой сигнал, причем по мере последовательного увеличения номера канала сигнал канала будет постепенно уменьшаться. В каналах Υ6 и Х6 сигналы будут фактически отсутствовать.
При обнаружении достаточного количества совпадающих событий может начаться реконструкция изображения. Обнаруженные события по существу разделяются на параллельные линии отклика (интерпретированная траектория фотонной пары), которые можно использовать для создания трехмерного изображения с применением компьютерной томографии.
Хотя технологии ΡΕΤ, ΜΚΙ и СТ являются общеизвестными способами создания изображений в медицине, информация, получаемая от этих различных подходов, сильно отличается. Технологии ΜΚ1 и СТ предоставляют анатомическую или структурную информацию. То есть они создают картину внутри тела. Это большое достижение с точки зрения таких проблем, как переломы костей, разрывы связок или что-либо подобное, что представляет собой нарушение структуры. Однако технологии ΜΚΙ и СТ не показывают метаболическую активность. Это является сферой томографии ΡΕΤ. Использование метаболически активных индикаторов означает, что изображения, созданные с помощью технологии ΡΕΤ, предоставляют функциональную или биохимическую информацию.
В настоящее время томография ΡΕΤ чаще всего применяется в онкологии (исследование опухолей). Некоторые канцерогенные ткани поглощают больше глюкозы, чем нормальная ткань. [Ε-18]ΕΌΟ по своим свойствам достаточно близко к глюкозе, которую канцерогенные клетки легко поглощают, и, следовательно, они будут иметь высокую радиоактивность по сравнению с фоновой тканью во время сканирования. Сканирование ΡΕΤ позволяет обнаружить некоторые виды опухоли, прежде чем те станут достаточно большими, чтобы их можно было обнаружить при сканировании согласно технологии ΜΚΙ. Информация, полученная при сканировании ΡΕΤ, также очень полезна для контроля за ходом лечения, так
- 2 016862 как можно отслеживать количество введенного индикатора по ходу лечения. Если сканирование указывает на низкую активность в той же канцерогенной ткани после терапии, это указывает на то, что терапия дает результат.
Томография РЕТ также полезна в неврологии (исследования нервной системы) и кардиологии (исследования сердца). Интересным приложением неврологии является ранняя диагностика болезни Паркинсона. Были разработаны индикаторы, которые концентрируются в клетках мозга, производящих допамин, который является нейропередатчиком. У пациентов с болезнью Паркинсона уменьшается количество нейронов, производящих допамин. Поэтому сканирование пациента с болезнью Паркинсона показывает меньшую активность, чем у здорового пациента. Это может обеспечить раннюю диагностику, поскольку множество других ранних признаков болезни Паркинсона схожи с другими заболеваниями.
Существует потребность в продолжении работ, направленных на снижение стоимости и повышение эффективности и точности систем РЕТ.
Сущность изобретения
Данный раздел предусмотрен в качестве введения к выбору концепций в упрощенной форме, которые дополнительно описываются ниже в подробном описании изобретения. Данный раздел не претендует на идентификацию ключевых признаков заявленного предмета изобретения и не предполагает возможность его использования при определении объема заявленного предмета изобретения.
Раскрыт способ для оценки начального момента времени электронных импульсов, таких как импульсы, созданные при обнаружении фотонов высокой энергии в позитронно-эмиссионной томографии, где полезно иметь очень точную информацию о начальном моменте времени. Обеспечены детекторы для обнаружения внешнего события, такого как попадание фотонов на детектор, и создания электронного аналогового импульсного сигнала. Выбирают параметризованную идеальную форму кривой для представления аналоговых импульсных сигналов, созданных детекторами. После приема аналогового импульсного сигнала его оцифровывают с помощью аналого-цифрового преобразователя (АЦП) для создания цифрового импульсного сигнала. Затем амплитуду цифрового импульсного сигнала нормализуют на основе вычисленной площади цифрового импульсного сигнала. Например, дискретные значения, которые содержит цифровой импульсный сигнал, можно масштабировать, используя отношение вычисленной площади к площади нормализованной кривой. Используя по меньшей мере одну точку из нормализованного импульсного сигнала, кривую на основе параметризованной идеальной формы кривой используют для оценки начального момента времени принятого аналогового импульсного сигнала. Следует иметь в виду, что в общем случае этот начальный момент времени является промежуточной точкой в цифровом импульсном сигнале. Затем может быть записана временная метка для аналогового импульсного сигнала.
В одном варианте параметризованная идеальная форма кривой содержит первый экспоненциальный участок, имеющий заранее определенную постоянную времени нарастания, и второй экспоненциальный участок, имеющий заранее определенную постоянную времени спада.
В одном варианте изобретения аналоговый импульсный сигнал создается кремниевым фотоумножителем.
В одном варианте изобретения для задания кривой используют только первую точку цифрового импульсного сигнала.
В одном варианте изобретения аналоговый импульсный сигнал фильтруют фильтром нижних частот до оцифровки сигнала.
В одном варианте изобретения используют обратный поиск, где в одной или нескольких поисковых таблицах хранится оценка начального момента времени для набора нормализованных точечных амплитуд сигнала.
В одном варианте изобретения детекторы содержат сцинтилляционный кристалл, связанный с одним из: фотоэлектронным умножителем, лавинным фотодиодом или кремниевым фотоумножителем.
Описание чертежей
Вышеуказанные аспекты и многие сопутствующие преимущества этого изобретения станут более очевидными и понятными со ссылками на последующее подробное описание вместе с сопроводительными чертежами, на которых фиг. 1 - внешний вид, показывающий систему со сканером РЕТ согласно настоящему изобретению;
фиг. 2 - выходные каналы примерной двенадцатиканальной системы РМТ с шестью каналами в направлении 'х' и шестью каналами в ортогональном направлении 'у';
фиг. 3 - блок-схема, показывающая архитектуру электронной аппаратуры предварительной обработки для одного варианта сканера РЕТ с высоким разрешением согласно настоящему изобретению;
фиг. 4 - блок-схема, показывающая архитектуру электронной аппаратуры предварительной обработки для второго варианта сканера РЕТ с высоким разрешением согласно настоящему изобретению;
фиг. 5 - распределение разности временных меток для импульсов при частоте дискретизации 70 МГц при использовании способа по настоящему изобретению;
фиг. 6 - корреляция между отфильтрованными импульсами и амплитудой импульсов, с указанием наилучшей линейной оценки корреляции;
- 3 016862 фиг. 7 - график стандартного отклонения точек отфильтрованного импульса, который дискретизируется с помощью АЦП с частотой 70 МГц, где линия показывает отфильтрованный импульс (инвертированный) со ссылкой на положение каждой точки в импульсе; и фиг. 8 - блок-схема, показывающая архитектуру схемы датчика временных характеристик, реализованной в вентильной матрице, программируемой пользователем (ЕРОА), для системы по фиг. 4.
Подробное описание изобретения
Далее следует описание конкретных вариантов системы РЕТ согласно настоящему изобретению со ссылками на фигуры, где одинаковые ссылочные позиции указывают на одинаковые элементы. Обратимся снова к фиг. 1, где раскрыт сканер 92 РЕТ с высоким разрешением с детекторами, содержащими сцинтилляторы 93 и фотоэлектронные умножители РМТ 94. Данные датчиков фильтруются фильтром 96 нижних частот, оцифровываются аналого-цифровым преобразователем 97 и оцифрованные данные первоначально обрабатываются вентильными матрицами 98, программируемыми пользователем (ЕРОА).
Аналоговые импульсы, созданные фотоэлектронными умножителями РМТ 94, содержат информацию, используемую для создания изображения РЕТ. Аналоговые импульсы обрабатывают для получения начального момента времени, местоположения и общей энергии. Устройство для выполнения этой начальной обработки называют здесь аппаратурой предварительной обработки, которая включает в себя фильтры 96, АЦП 97 и матрицы ЕРОА 98. Аналоговый импульс, полученный из РМТ 94, фильтруется фильтром 96 нижних частот для удаления шума, а затем оцифровывается преобразователем АЦП 97 для обработки матрицей ЕРОА 98. Хотя некоторые современные АЦП могут выполнять дискретизацию с частотой вплоть до 400 миллионов отсчетов в секунду (М8Р8), в данном варианте изобретения выбран последовательный АЦП 97 с частотой дискретизации 70 М8Р8. Такой выбор значительно снижает сложность, стоимость и энергопотребление установки РЕТ. Другим аспектом является количество входов в матрицу ЕРОА 98. Быстродействующие АЦП имеют параллельный выход, для которого требуется 10-12 бит на канал, то есть от десятков до сотен каналов для матрицы ЕРОА. Таким образом, количество входов будет превышать возможности обработки современными матрицами ЕРОА. Поэтому в данной системе используют АЦП 97 с последовательным выходом, что ограничивает частоту дискретизации примерно до 100 М8Р8. Однако для систем, где требуется меньше АЦП на одну матрицу ЕРОА, для обеспечения лучшего временного разрешения можно использовать более быстродействующие АЦП.
После оцифровки данных аналогового импульса в матрице ЕРОА 98 могут быть выделены необходимые параметры импульса. Например, общую энергию импульса можно получить путем суммирования отсчетов значений импульсов и вычитания из базового уровня (значение выхода АЦП 97 без входного импульса).
Начальный момент времени импульса важен для определения совпадающих пар, то есть двух обнаруженных фотонов, которые возникают в результате одного события аннигиляции. Специалистам в данной области техники очевидно, что многие фотоны, созданные в результате событий аннигиляции, не обнаруживаются сканером 92. Например, созданный фотон может либо поглощаться, либо рассеиваться тканью тела, или может пройти по траектории, не пересекающейся со сцинтиллятором 93. Для создания изображения РЕТ требуется обнаружение обоих фотонов, получившихся в результате события аннигиляции, поскольку, как известно, траектория обнаруженных фотонов фактически представляет прямую линию между ними. Если сканер 92 обнаружил только один из двух эмитированных фотонов, то невозможно определить место появления указанного события. Если два обнаруженных фотона являются совпадающей парой, они должны быть обнаружены в течение определенного времени друг относительно друга, причем каждый из детекторов должен находиться в поле зрения другого детектора.
Был сконструирован сканер РЕТ с высоким разрешением для небольших животных, содержащий кольцо из восемнадцати детекторных кассет, где каждая кассета имеет четыре матрицы сцинтилляторов 93, прикрепленных к четырем фотоэлектронным умножителям РМТ 94. Каждая детекторная кассета подсоединена к специализированному набору электронной аппаратуры предварительной обработки.
На фиг. 3 показана блок-схема электронной аппаратуры 100 предварительной обработки для примерного сканера для небольших животных. Электронная аппаратура 100 предварительной обработки содержит несколько узлов (показано два), содержащих пару (показано две) устройств: микропроцессор 102 и ЕРОА 104. Каждый узел поддерживает два фотоэлектронных умножителя РМТ 106. Все узлы соединены в цепочку вместе с использованием стандартного высокоскоростного коммуникационного интерфейса 108, например, интерфейса 1394а стандарта 1ЕЕЕ (например, в реализации ЕйеХУйе®) для создания соединения с главным компьютером (не показан). Главный компьютер содержит программное обеспечение, которое собирает и обрабатывает данные от всех узлов. Эти данные обрабатываются фактически автономно для создания требуемого изображения.
Как показано на фиг. 3, в электронной аппаратуре предварительной обработки 100 имеется множество дискретных блоков. Фотоумножители 106 обнаруживают свет от сцинтилляционных кристаллов (на фиг. 3 не показаны) и создают соответствующие аналоговые импульсы. Каждый фотоэлектронный умножитель РМТ 106 выдает двенадцать сигналов: шесть для направления 'х' и шесть для направления 'у' (смотри фиг. 2). В данном варианте для сокращения входов/выходов двенадцать сигналов сокращены до четырех с помощью платы 110 суммирования. Специализированные интегральные схемы А81С 112 по
- 4 016862 лучают данные и реализуют алгоритм, определяющий время обнаруженного события. Сигналы оцифровываются преобразователями АЦП 114 и посылаются в матрицы ЕРСА 104. Может, но не обязательно, быть предусмотрен блок 109 совпадения для идентификации пар событий, которые могут совпадать с фильтрацией событий, которые очевидно бесполезны для дальнейшей обработки.
Вдобавок к матрицам ЕРСА 104 имеется микропроцессор 102 (в данном примере используется микроконтроллер КаЬЫ™1) и микросхема 108 физического уровня коммуникационного интерфейса, используемая для связи с главным компьютером. Матрица ЕРСА 104 и микропроцессор 102 играют ключевую роль в сборе данных. Микропроцессор 102 обеспечивает общее управление каждым отдельным узлом. Это включает в себя конфигурирование матриц ЕРСА 104, обмен данными с главным компьютером, инициализацию системы и взаимодействие с ЕРСА 104 для настройки системы. В то время как микропроцессоры 102 обеспечивают обработку данных для управления системой, матрицы ЕРСА 104 выстраивают основную часть тракта данных. При нормальной работе матрицы ЕРСА 104 решают две основные задачи: обработку импульсов и пакетирование данных для коммуникационного интерфейса 108. При обработке импульсов на первом шаге определяют, совпадает ли данное событие с другим событием по другую сторону от сканера 92. Если обнаружен случай совпадения, то интегрируется импульс для определения энергии и к данным прикладывается метка с грубым временным разрешением. Данные об энергии и времени (грубый отсчет из матрицы ЕРСА и точный отсчет из схемы А81С) посылают в главный компьютер через коммуникационный интерфейс 108.
Время от времени необходима настройка сканера 92 для установки коэффициентов усиления усилителей в схемах А81С 112. Это необходимо для того, чтобы сцинтилляционные кристаллы 93 имели разную светоотдачу и светоулавливающую способность, а фотоэлектронные умножители РМТ 106 имели разные характеристики усиления. Настройка нормализует разницу между датчиками и корректирует дрейф во времени. Для настройки сканера микропроцессор 102 реконфигурирует матрицы ЕРСА 104, используя алгоритм настройки, и инициализирует матрицы ЕРСА 104 для сохранения значений энергии для каждого импульса. К сожалению, краевые эффекты возле периферии матрицы 93 сцинтилляционных кристаллов и РМТ 106 могут вносить ошибки в алгоритм настройки, и поэтому они игнорируются. Для того чтобы отфильтровать событие, появляющееся на краях матрицы кристаллов, необходимо декодировать местоположение события. Четыре сигнала, поступающие из платы 110 суммирования, содержат достаточно информации для определения местоположения события в матрице 93 кристаллов.
Как только определено, что событие имело место в интересующей области матрицы 93 сцинтилляционных кристаллов, вычисляют общую энергию события и данные об этом помещают в разные столбцы гистограммы энергии. В результате создается гистограмма энергии с пиком, называемым фотопиком (поскольку он представляет энергию для того случая, когда вся энергия фотона 511 кэВ попадает в сцинтилляционный кристалл). Отсчеты, лежащие ниже фотопика, представляют рассеянные фотоны, когда только часть энергии фотона попадает в кристалл или фотон претерпевает комптоновское рассеяние в изображаемом объекте, прежде чем он достигнет кристалла.
Если система настроена, то фотопики должны быть выровнены для всех детекторов. Если пики имеют некоторые отклонения, микропроцессор изменяет коэффициенты усиления и вновь выполняет алгоритм настройки. Все это делается автоматически, так что, как только оператор даст команду машине выполнить самонастройку, микропроцессор 102 выдаст команду матрицам ЕРСА 104 на выполнение алгоритма настройки. После того как ЕРСА 104 просигнализировал о том, что процедура настройки завершена, микропроцессор 102 считывает упомянутую гистограмму из памяти и определяет местоположение фотопика. Если фотопик смещен, то микропроцессор 102 регулирует коэффициенты усиления в схеме А81С 112 и выполняет итерации до тех пор, пока фотопики не выровняются для всех датчиков.
Во втором варианте системы сканирования РЕТ, показанном на фиг. 4, фотоумножительные устройства 144 являются твердотельными умножителями 81РМ, которые имеют один выход на кристалл, что значительно увеличивает количество каналов, используемых на один датчик. На фиг. 4 представлена блок-схема, показывающая электронную аппаратуру 150 предварительной обработки для указанного второго варианта. Для этого варианта были выбраны матрица ЕРСА 148 типа 8ΐηΙίχ® II ЕР860 и последовательный АЦП 147 на 70 МГц, которые могут легко осуществлять связь друг с другом с использованием специализированной схемы параллельно-последовательного/последовательно-параллельного преобразования и контуров фазовой автоподстройки частоты. В этом варианте большинство функциональных возможностей электронной аппаратуры предварительной обработки реализуется матрицей ЕРСА 148, что обеспечивает более простую и компактную архитектуру. Фильтры 146 показаны в виде простых К.С-фильтров, хотя очевидно, что можно использовать другие фильтры, известные в данной области техники.
Имеются другие аспекты этих более современных матриц ЕРСА 148, которые делают их весьма подходящими для применения в качестве приложения для предварительной обработки в системе РЕТ. Например, микропроцессор 102 КаЬЬй ™ заменяют на встроенный процессор с программным ядром Νίοδ® II (152) в матрице ЕРСА 148, исключая связь на низкой скорости между ЕРСА 104 и микропроцессором 102 в архитектуре, показанной на фиг. 3. Как следует из сравнения фиг. 4 с фиг. 2, архитектура
- 5 016862
150 электронной аппаратуры предварительной обработки в сканере по второму варианту является более компактной благодаря интеграции отдельных частей в матрице ΡΡΘΆ 148, которая включает в себя логический блок 154 временных датчиков, логический блок 156 вычисления энергии, а также ядро 158 Е1ге^1ге® вдобавок к встроенному процессору 152.
Другим следствием большого количества каналов является неприемлемая стоимость и недопустимо большое пространство на плате, необходимое для схемы А81С, что определяет временные характеристики для каждого канала. Это привело к разработке алгоритма для выполнения полного определения временных характеристик внутри матрицы ΡΡΘΆ 148.
Важной целью создания высококачественных изображений ΡΕΤ является точное определение временных характеристик фотонов, взаимодействующих со сцинтилляционным кристаллом 93. Временное разрешение непосредственно коррелируется с количеством несовпадающих событий, которые воспринимаются как правильные события и, следовательно, добавляются к шуму конечного изображения. В описываемом втором варианте определение временных характеристик выполняется в матрице ΕΡΘΆ 148 с выборочными данными, что исключает необходимость иметь схемы Л81С на каждый канал.
Критерием качества измерения временных характеристик является распределение временных меток. Другими словами, для конфигурации с двумя детекторами и источником, расположенным точно по центру между ними (так что фотоны попадают на оба детектора одновременно), анализируют распределение разностей между временными метками для каждого детектора. В случае идеальной системы разность между временными метками будет нулевой. Однако шум в системе вызывает ошибки. Для имитации этого временные метки вычислялись для множества различных отсчетов одного и того же импульса. Временные метки для каждого импульса затем сравнивались с временными метками всех других импульсов для создания распределения 160 разностей временных меток, как показано на фиг. 5.
Для определения временных характеристик фотона, взаимодействующего с сцинтилляторами 93, используют схему 156 временного датчика. Схема 156 временного датчика присваивает временную метку конкретному признаку обнаруженного импульсного сигнала, полученного от РМТ. Например, этим признаком может быть начальный момент времени импульса, момент пикового значения импульса или момент, когда импульс пересекает заранее определенный уровень напряжения. Два традиционных способа измерения временных характеристик предусматривают использование дискриминатора переднего фронта и дискриминатора постоянной составляющей (СЕО). Передний фронт просто определяется по моменту пересечения импульсом некоторого фиксированного порогового напряжения. Для этого требуется аналоговая схема, которая обнаруживает указанное пересечение. Недостатком этого способа является то, что время достижения порогового значения зависит от амплитуды импульса. Этот эффект усугубляется, когда порог срабатывания установлен высоко.
На сегодняшний день измерение временных характеристик для систем ΡΕΤ выполняют с помощью аналоговых дискриминаторов СЕО, поскольку они не чувствительны к изменению амплитуды импульса. СЕО реализует схему для следующего уравнения:
где δ(ί) - входной сигнал. Вычисление уравнения выполняется путем разделения аналогового импульса на две копии и задержки одной копии на I). Другая копия инвертируется и ослабляется на постоянную дробную величину (как правило, примерно на 0,2). Наконец, две измененные копии складывают для создания импульса с переходом через ноль, который можно обнаружить и снабдить временной меткой. Переход через ноль появляется на постоянной части амплитуды импульса для импульсов одной и той же формы. Дискриминацию как по варианту СЕО, так и по переднему фронту, как правило, осуществляют в специализированных схемах А81С, причем необходимо иметь схему для преобразования сигнала переключения во временную метку. Дискриминаторы СЕО могут обеспечить субнаносекундное временное разрешение.
Способ измерения временных характеристик
Здесь раскрыт способ использования известных характеристик импульсов для вычисления начала импульса, обеспечивающий тем самым временное разрешение на уровне частоты субдискретизации. Например, для сцинтилляционных кристаллов Б8О время нарастания предопределено характеристикой фотоэлектронного умножителя РМТ, в то время как время спада зависит от сцинтилляционного кристалла. На основе этих положений начальный момент времени импульса можно определить путем фильтрации идеального импульса с получением дискретизированного импульса и путем использования идеального импульса для интерполяции начальной точки импульса.
Для проверки этого алгоритма определения временных характеристик на реальных данных был использован осциллограф с частотой 25 гигаотсчетов/с для дискретизации девятнадцати импульсов из фотоэлектронного умножителя РМТ, который был соединен с кристаллом Б8О. Для создания импульсов был использован источник (22\а) с энергий 511 кэВ. Затем данные из осциллографа были импортированы в программу МАТБАВ®.
Модельная кривая, обеспечивающая хорошее соответствие с данными импульсов, состоит из двух экспоненциальных кривых, например:
- 6 016862 [-η*Τ, -η·Τ, ехр г* -ехр г'
В качестве начального шага была выдвинута гипотеза о том, что если создан импульс с двумя экспоненциальными кривыми (одна для переднего фронта и одна для заднего фронта) и найдена амплитуда, временной сдвиг, спадающая экспоненциальная кривая и нарастающая экспоненциальная кривая, которые обеспечивают наилучшее совпадение с измеренным данными цифровых импульсов по методу наименьших квадратов, то можно определить идеальный импульс и использовать его для интерполяции начальной точки импульса. При использовании этого способа решения в лоб получено стандартное отклонение 1,0 нс при измерении временных характеристик с помощью АЦП на 70 МГц. Хотя это является хорошим временным разрешением, пространство поиска слишком велико для того, чтобы матрица РРОА могла выполнять вычисления в реальном времени. На основе метода решения в лоб было обнаружено, что время нарастания находится в диапазоне от 0,1 до 0,5 нс, время спада находится в диапазоне от 28 до 38 нс, а амплитуда находится в диапазоне от 0,082 до 0,185 В. Для перекрытия этих диапазонов при разумном временном шаге (примерно 40 пс) потребуется по меньшей мере 215000 операций вычисления по методу наименьших квадратов и сравнения для каждого импульса (11 шагов для времени спада, 5 шагов для времени нарастания, 11 шагов для амплитуды и 357 временных шагов).
Для разработки более эффективного алгоритма и способа на основе РРОЛ предположим сначала, что времена нарастания и спада (τκ, τΡ) импульсов РМТ/81РМ являются константами, и изменчивость импульсов определяется амплитудой импульса и белым шумом.
Например, в нашем тестовом устройстве времена нарастания и спада, которые дали наилучшее совпадение по методу наименьших квадратов для всех не отфильтрованных и не дискретизированных данных (то есть необработанные данные от осциллографа), составили 310 пс и 34,5 нс соответственно. При использовании метода решения в лоб с фиксированными значениями времени нарастания и времени спада стандартное отклонение при измерении временных характеристик увеличивается до 1,1 нс. Однако даже после исключения операций поиска постоянной времени для каждого случая все еще потребуется почти 4000 операций поиска.
Для дальнейшего упрощения способа исключаем расхождение в амплитуде для эталонного импульса (определенного уравнением с двумя экспонентами) и импульса входных данных с использованием прямой корреляции между площадью и амплитудой импульса. Между площадью импульса, порожденного событием, и амплитудой этого импульса существует значительная прямая корреляция 170, как показано на фиг. 6. Для нормализации амплитуды события по эталонному импульсу вычисляют отношение площади эталонного импульса к площади порожденного событием импульса. Затем порожденный событием импульс можно масштабировать с помощью этого отношения для выравнивания амплитуды эталонного и порожденного событием импульсов. Например, можно вычислить площадь оцифрованного порожденного событием импульса и сравнить ее с площадью идеального эталонного импульса. Затем отсчитанные точки порожденного событием импульса масштабируют или нормализуют согласно отношению этих двух площадей. Затем порожденный событием нормализованный цифровой импульс можно использовать для оценки начального момента порожденного событием импульса путем сравнения значений нормализованного порожденного событием импульса с эталонным импульсом.
Путем дискретизации каждого из девятнадцати импульсов с множеством различных начальных точек и корреляции площади, полученной для каждого отсчета, с известной амплитудой для данного полного импульса была определена функция для преобразования площади в амплитуду. При использовании этой оценки стандартное отклонение измерения временных характеристик ухудшилось до 1,2 нс.
При использовании этих двух аппроксимаций большинство размерностей поиска методом полного перебора было исключено с потерей всего лишь 20% временного разрешения. Однако этот алгоритм все еще требует 357 операций поиска для каждого возможного временного сдвига. При условии, что данные импульса согласуются с эталонной кривой, имеющей известные времена нарастания и спада, и что амплитуду вычисляют, исходя из площади импульса, поиск в лоб можно преобразовать в обратный просмотр.
Положим, например, что импульс имеет длительность порядка 2 х 10-7 с, и тогда при частоте дискретизации 70 МГц получится около тринадцати точек отсчета. Для каждого возможного входного напряжения предварительно вычисляют момент его появления в эталонном импульсе. Таким образом, каждое приходящее напряжение можно преобразовать во временной сдвиг с помощью простой операции запоминания. Это выполняется для каждого импульса, так что после просмотра любые или все тринадцать точек отсчета можно использовать для оценки момента времени, с которого начался импульс. Если эти тринадцать начальных моментов времени усреднить, то временное разрешение значительно ухудшается (до 2,84 нс).
После тщательной проверки результатов, полученных описываемым способом просмотра, стало очевидным, что некоторые точки отсчета дают гораздо лучшие результаты, чем другие. Это показано на фиг. 7, где изображено стандартное отклонение 180 для вычисленных начальных моментов времени для
- 7 016862 каждой из тринадцати точек отсчета. Из фиг. 7 следует, что стандартное отклонение 180 коррелирует со спадом отфильтрованного импульса 182 и расстоянием от начала импульса. Точки рядом с пиком (отсчеты 4 и 5) имеют малый уклон, и поэтому небольшое изменение напряжения приводит к большому временному сдвигу. Хвост импульса также имеет большое отклонение. Однако, если использовать только одну точку отсчета, то стандартное отклонение измерения временных характеристик составит 1,03 нс, что фактически совпадает с методом поиска в лоб.
На основе этой информации шаг обратного просмотра был изменен в том отношении, что использовалась только первая точка отсчета выше 0,05 В в порожденном событием импульсе. Однако при использовании более быстродействующих АЦП или импульсов с большей длительностью их нарастания можно усреднить или иным образом скоррелировать большее количество точек отсчета для получения лучшего конечного результата.
Таким образом, в данном алгоритме определения временных характеристик используется одна постоянная спада, одна постоянная нарастания, вычисляется амплитуда импульса, исходя из его площади и используется поиск по схеме напряжение-время для первого отсчета. Проведенное тестирование показало, что алгоритм дает стандартное отклонение при измерении временных характеристик всего лишь 1,03 нс. Распределение результатов для конечного алгоритма показано на фиг. 5 для АЦП с частотой 70 МГ ц. Хотя использование только первой точки отсчета дает очень хорошие результаты, специалистам в данной области техники должно быть очевидно, что в некоторых случаях оценку измерения можно улучшить, если использовать взвешенное среднее для более чем одной точки отсчета. Например, в другом варианте настоящего изобретения можно использовать первые две или три точки отсчета с экспериментально полученными весовыми коэффициентами для дальнейшего повышения согласованности измерения временных характеристик.
Реализованная архитектура алгоритма определения временных характеристик показана на фиг. 8. Фотоэлектронный умножитель РМТ 200, например фотоумножитель 8ίΡΜ, принимает импульсный сигнал от сцинтиллятора и создает выходной импульс (представленный в виде зашумленного сигнала 200'), который подвергается фильтрации, как правило, фильтром 202 нижних частот (представлен сглаженным сигналом 202'), а затем оцифровывается аналого-цифровым преобразователем 204 (представлено цифровым сигналом 204'). Цифровой сигнал 204' посылают в матрицу ΡΡΟΑ 206. ΡΡΟΑ 206 конфигурируется для вычисления площади обнаруженного сигнала, а для оценки амплитуды 208 сигнала используют корреляцию между площадью и амплитудой. Затем сигнал может быт нормализован (210) по площади эталонного импульса. Выбранная из цифрового сигнала 204' точка данных (или набор точек), например первая точка над заданным напряжением, используется затем при обратном просмотре 212 для нахождения эталонной кривой 214 импульса, которую используют для определения точного начального момента 216 для обнаруженного импульса.
Заметим, что хотя раскрытый алгоритм при использовании АЦП на 70 МГц может обеспечить несколько меньшее временное разрешение, чем аналоговый метод на основе СРИ, временное разрешение будет повышаться при совершенствовании технологий аналого-цифрового преобразования. Учитывая, что разрешение метода на основе СРИ не соответствует технологии (рабочие характеристики СРИ остаются фактически постоянными последние десять или более лет), настоящий алгоритм по прогнозам превзойдет метод на основе СРИ при использовании АЦП на 500 МГц (что имеется в настоящее время в параллельных АЦП и ожидается в ближайшем будущем в последовательных АЦП).
Следует иметь в виду, что даже в тех случаях, когда настоящий способ не совпадает со способами на основе СРИ по временному разрешению, для способов на основе СРИ потребуется заказной логический блок на каждый канал в виде нерегулируемых схем А81С. Предлагаемый полностью цифровой способ позволяет избежать этих затрат, что является весьма существенным в случае сканеров ΡΕΤ, которые могут содержать, например, 128 каналов на одну матрицу 6ΡΟΑ.
Короче говоря, томография ΡΕΤ является приложением, которое хорошо подходит для реализации на основе матриц 6ΡΟΑ. Матрицы 6ΡΟΑ являются идеальным средством для разработки таких алгоритмов, как цифровое определение временных характеристик, но они также обеспечивают большинство элементов, необходимых для усовершенствованной системы сбора и обработки данных для ΡΕΤ. В настоящей системе используются возможности реконфигурирования матриц 6ΡΟΑ для разработки алгоритма определения временных характеристик, которые под управлением микропроцессора могут регулировать коэффициенты усиления и устанавливать регистры для учета изменений в различных частях сканера и усложненный ввод/вывод для взаимодействия с быстродействующими последовательными АЦП, что позволяет обрабатывать данные во множестве каналов. Эти каналы также могут обрабатываться параллельно в реконфигурируемой архитектуре, что повышает скорость отсчета, которую может обрабатывать сканер. Увеличение производительности вычислений в современных матрицах 6ΡΟΑ по сравнению с предыдущим поколением позволяет определять временные характеристики непосредственно в матрице 6ΡΟΑ и исключить схемы А81С.
Здесь также раскрыт новый, чисто цифровой механизм измерения временных характеристик, который отличается более высоким временным разрешением, чем известные подходы, с учетом существую
- 8 016862 щих и будущих технологий аналого-цифрового преобразования. Множество функциональных особенностей современных матриц ЕРСА можно в полной мере использовать для поддержки законченной комплексной системы обработки сигналов в важной области электроники.
Хотя вышеописанный способ дает очень хорошие результаты измерения временных характеристик, предполагается, что измерение временных характеристик можно дополнительно усовершенствовать на основе использования альтернативного способа определения подходящего эталонного импульса при сохранении нормализации амплитуды и методики нахождения временных характеристик, описанной выше. В альтернативном способе матрицы РРСЛ конфигурируют для фиксации и запоминания множества появившихся импульсов и использования этих импульсов для формирования эталонного импульса. В частности, зафиксированные данные рассеивают путем дифференцирования амплитуд и временных сдвигов частоты субдискретизации (то есть начального момента импульса относительно интервала отсчета импульса). Для формирования составного эталонного импульса используют двухшаговый процесс. Сначала принятые импульсы нормализуют по амплитудам, как обсуждалось выше, и создают начальный эталонный импульс путем усреднения данных. Затем каждый отдельный импульс выравнивают по этому эталонному импульсу путем сдвига во времени, выполняя итерации, пока все кривые не окажутся максимально коррелированными. Затем выровненные импульсы используют для формирования конечного эталонного импульса. Затем конечный эталонный импульс можно использовать вместо кривой эталонного импульса, состоящей из двух экспонент, как обсуждалось выше.
Хотя раскрытый способ для точной цифровой оценки начального момента времени для появившихся аналоговых импульсных сигналов был разработан для позитронно-эмиссионной томографии, специалистам в данной области техники должно быть очевидно, что этот общий способ можно использовать для весьма точной оценки начального момента импульсных сигналов в других контекстах, и поэтому несомненно, что этот способ подойдет для использования в других приложениях, где требуется иметь информацию о начальном моменте времени для быстро протекающих событий.
Хотя здесь были показаны и описаны иллюстративные варианты изобретения, следует иметь в виду, что в них могут быть внесены различные изменения, не выходящие за рамки существа и объема изобретения.

Claims (27)

  1. ФОРМУЛА ИЗОБРЕТЕНИЯ
    1. Способ оценки начального момента времени электронного импульса, созданного в ответ на обнаруженное событие, причем способ содержит обнаружение детектором внешнего события и реагирование детектором на обнаруженное событие, создавая электронный аналоговый импульсный сигнал;
    выбор параметризованной идеальной формы кривой для представления аналоговых импульсных сигналов, созданных детектором;
    прием аналогового импульсного сигнала, созданного детектором;
    оцифровку принятого аналогового импульсного сигнала для создания цифрового импульсного сигнала, имеющего амплитуду;
    нормализацию амплитуды цифрового импульсного сигнала на основе вычисленной площади цифрового импульсного сигнала;
    использование нормализованной амплитуды по меньшей мере одной точки данных нормализованного цифрового импульсного сигнала в обратном просмотре для задания кривой исходя из параметризованной идеальной формы кривой для представления принятого аналогового импульсного сигнала;
    интерполяцию заданной кривой для оценки начального момента времени принятого аналогового импульсного сигнала и запись временной метки, указывающей оцененный начальный момент времени принятого аналогового импульсного сигнала.
  2. 2. Способ по п.1, в котором параметризованная идеальная форма кривой содержит первый экспоненциальный участок, имеющий заранее определенную постоянную времени нарастания, и второй экспоненциальный участок, имеющий заранее определенную постоянную времени спада.
  3. 3. Способ по п.1, в котором принятый аналоговый импульсный сигнал формируется кремниевым фотоумножителем.
  4. 4. Способ по п.1, в котором шаг использования по меньшей мере одной точки цифрового импульсного сигнала для задания кривой содержит использование только первой точки цифрового импульсного сигнала для задания кривой.
  5. 5. Способ по п.1, дополнительно содержащий шаг фильтрации принятого аналогового импульсного сигнала фильтром нижних частот.
  6. 6. Способ по п.1, в котором шаг использования нормализованного цифрового импульсного сигнала для задания кривой исходя из параметризованной идеальной формы кривой содержит обратный просмотр, где в одной или нескольких поисковых таблицах хранится оценка начального момента времени для набора нормализованных точечных амплитуд сигнала.
    - 9 016862
  7. 7. Способ по п.1, в котором детекторы содержат сцинтилляционный кристалл, связанный с одним из: фотоэлектронным умножителем, лавинным фотодиодом или кремниевым фотоумножителем.
  8. 8. Способ оценки начального момента времени для импульса, обнаруженного при позитронноэмиссионной томографии, причем способ содержит обнаружение детектором фотонов, имеющих энергию порядка 511 кэВ, и создание детектором в ответ аналогового импульсного сигнала;
    выбор параметризованной идеальной формы кривой для представления аналоговых импульсных сигналов, созданных детектором;
    прием аналогового импульсного сигнала, созданного детектором;
    оцифровку принятого аналогового импульсного сигнала для создания цифрового импульсного сигнала, имеющего амплитуду;
    нормализацию амплитуды цифрового импульсного сигнала на основе вычисленной площади цифрового импульсного сигнала;
    использование нормализованной амплитуды по меньшей мере одной точки данных нормализованного цифрового импульсного сигнала в обратном просмотре для задания кривой исходя из параметризованной идеальной формы кривой для представления принятого аналогового импульсного сигнала;
    интерполяцию заданной кривой для оценки начального момента времени принятого аналогового импульсного сигнала и запись временной метки, указывающей оцененный начальный момент времени принятого аналогового импульсного сигнала.
  9. 9. Способ по п.8, в котором параметризованная идеальная форма кривой содержит первый экспоненциальный участок, имеющий заранее определенную постоянную времени нарастания, и второй экспоненциальный участок, имеющий заранее определенную постоянную времени спада.
  10. 10. Способ по п.8, в котором принятый аналоговый импульсный сигнал создается кремниевым фотоумножителем.
  11. 11. Способ по п.8, в котором шаг использования по меньшей мере одной точки цифрового импульсного сигнала для задания кривой содержит использование только первой точки цифрового импульсного сигнала для задания кривой.
  12. 12. Способ по п.8, дополнительно содержащий шаг фильтрации принятого аналогового импульсного сигнала фильтром нижних частот.
  13. 13. Способ по п.8, в котором шаг использования нормализованного цифрового импульсного сигнала для задания кривой исходя из параметризованной идеальной формы кривой содержит обратный просмотр, где в одной или нескольких поисковых таблицах хранится оценка начального момента времени для набора нормализованных точечных амплитуд сигнала.
  14. 14. Способ по п.8, в котором детекторы содержат сцинтилляционный кристалл, связанный с одним из: фотоэлектронным умножителем, лавинным фотодиодом или кремниевым фотоумножителем.
  15. 15. Способ по п.8, дополнительно содержащий шаг предварительной обработки электронной аппаратурой, содержащей аналого-цифровой преобразователь для оцифровки принятого аналогового импульса и вентильную матрицу, программируемую пользователем, которая обрабатывает цифровой импульсный сигнал.
  16. 16. Способ по п.15, в котором разрешение оцененного начального момента времени принятого аналогового импульсного сигнала меньше интервала дискретизации, используемого для оцифровки принятого аналогового сигнала.
  17. 17. Способ для идентификации совпадающих пар в позитронно-эмиссионной томографии, причем способ содержит обнаружение множества фотонов и создание аналогового импульсного сигнала в ответ на каждый обнаруженный фотон с использованием множества детекторов, скомпонованных в форме кольца, так что некоторые из детекторов расположены в поле зрения друг друга;
    выбор параметризованной идеальной формы кривой для представления аналоговых импульсных сигналов, созданных детектором;
    оцифровку аналоговых импульсных сигналов для создания цифровых импульсных сигналов, имеющих амплитуду;
    нормализацию амплитуды каждого цифрового импульсного сигнала на основе вычисленной площади цифрового импульсного сигнала;
    использование нормализованной амплитуды по меньшей мере одной точки данных нормализованного цифрового импульсного сигнала в обратном просмотре для задания кривой исходя из параметризованной идеальной формы кривой для представления аналогового импульсного сигнала;
    интерполяция заданной кривой для оценки начального момента времени аналогового импульсного сигнала для каждого аналогового импульсного сигнала;
    запись временной метки, указывающей оцененный начальный момент времени аналоговых импульсных сигналов; и сравнение временных меток аналоговых импульсных сигналов от детекторов, расположенных в по
    - 10 016862 ле зрения друг друга, для идентификации совпадающих пар.
  18. 18. Способ по п.17, в котором параметризованная идеальная форма кривой содержит первый экспоненциальный участок, имеющий заранее определенную постоянную времени нарастания, и второй экспоненциальный участок, имеющий заранее определенную постоянную времени спада.
  19. 19. Способ по п.17, в котором шаг использования по меньшей мере одной точки цифрового импульсного сигнала для задания кривой содержит использование только первой точки цифрового импульсного сигнала для задания кривой.
  20. 20. Способ по п.17, дополнительно содержащий шаг фильтрации принятого аналогового импульсного сигнала фильтром нижних частот.
  21. 21. Способ по п.17, в котором шаг использования нормализованного цифрового импульсного сигнала для задания кривой исходя из параметризованной идеальной формы кривой содержит обратный просмотр, где в одной или нескольких поисковых таблицах хранится оценка начального момента времени для набора нормализованных точечных амплитуд сигнала.
  22. 22. Способ по п.17, в котором детекторы содержат сцинтилляционный кристалл, связанный с одним из: фотоэлектронным умножителем, лавинным фотодиодом или кремниевым фотоумножителем.
  23. 23. Способ по п.17, дополнительно содержащий шаг предварительной обработки электронной аппаратурой, содержащей аналого-цифровой преобразователь для оцифровки принятого аналогового импульса и вентильную матрицу, программируемую пользователем, которая обрабатывает цифровой импульсный сигнал.
  24. 24. Способ по п.23, в котором разрешение оцененного начального момента времени принятого аналогового импульсного сигнала меньше интервала дискретизации, используемого для оцифровки принятого аналогового сигнала.
  25. 25. Сканер для позитронно-эмиссионной томографии, содержащий множество детекторов, скомпонованных в форме кольца, причем каждый детектор содержит по меньшей мере один сцинтиллятор и по меньшей мере один фотоумножитель;
    систему электронной аппаратуры предварительной обработки, содержащую аналого-цифровой преобразователь, способный принимать аналоговые сигналы от детекторных фотоумножителей, причем системы электронной аппаратуры предварительной обработки включают в себя аналого-цифровые преобразователи, которые преобразуют принятые аналоговые сигналы в цифровые сигналы, и вентильные матрицы, программируемые пользователем, которые принимают цифровые сигналы и вычисляют начальные моменты времени аналоговых сигналов;
    где вентильные матрицы, программируемые пользователем, вычисляют начальные моменты времени аналоговых сигналов путем ί) нормализации амплитуды каждого цифрового импульсного сигнала на основе вычисленной площади цифрового импульсного сигнала; и) использования нормализованной амплитуды по меньшей мере одной точки данных нормализованного цифрового импульсного сигнала в обратном просмотре для задания кривой исходя из параметризованной идеальной формы кривой для представления аналогового импульсного сигнала; ш) интерполяции заданной кривой для оценки начального момента времени принятого аналогового импульсного сигнала для каждого аналогового импульсного сигнала.
  26. 26. Сканер для позитронно-эмиссионной томографии по п.25, в котором каждый из по меньшей мере одного фотоумножителя содержит одно из: электронный фотоумножитель, кремниевый фотоумножитель или лавинный фотодиод.
  27. 27. Сканер для позитронно-эмиссионной томографии по п.25, дополнительно содержащий фильтр, фильтрующий аналоговый сигнал перед его оцифровкой.
EA201070554A 2007-11-02 2008-11-03 Способ и устройство для оценки начального момента времени электронного импульса, предназначенного для использования в позитронно-эмиссионной томографии EA016862B1 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US98508307P 2007-11-02 2007-11-02
PCT/US2008/082273 WO2009059312A2 (en) 2007-11-02 2008-11-03 Data acquisition for positron emission tomography

Publications (2)

Publication Number Publication Date
EA201070554A1 EA201070554A1 (ru) 2010-12-30
EA016862B1 true EA016862B1 (ru) 2012-08-30

Family

ID=40591800

Family Applications (1)

Application Number Title Priority Date Filing Date
EA201070554A EA016862B1 (ru) 2007-11-02 2008-11-03 Способ и устройство для оценки начального момента времени электронного импульса, предназначенного для использования в позитронно-эмиссионной томографии

Country Status (9)

Country Link
US (2) US8003948B2 (ru)
EP (1) EP2208088B1 (ru)
JP (1) JP5667446B2 (ru)
KR (1) KR101583272B1 (ru)
CN (1) CN101903798B (ru)
AU (1) AU2008318371B2 (ru)
CA (1) CA2706388C (ru)
EA (1) EA016862B1 (ru)
WO (1) WO2009059312A2 (ru)

Families Citing this family (57)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006029475A1 (en) 2004-09-16 2006-03-23 Southern Innovation International Pty Ltd Method and apparatus for resolving individual signals in detector output data.
US8543356B2 (en) * 2008-01-24 2013-09-24 National Institute Of Standards And Technology Low cost multi-channel data acquisition system
JP2011516838A (ja) 2008-03-31 2011-05-26 サザン イノヴェーション インターナショナル プロプライアトリー リミテッド スクリーニング方法および装置
JP2011519415A (ja) 2008-03-31 2011-07-07 サザン イノヴェーション インターナショナル プロプライアトリー リミテッド 個別信号解像度を用いた放射線画像形成法
CN104020500A (zh) 2008-03-31 2014-09-03 南方创新国际股份有限公司 用于钻孔测井的方法和装置
AU2009328651A1 (en) 2008-12-18 2011-06-30 Southern Innovation International Pty Ltd Method and apparatus for resolving piled-up pulses by using a mathematical transform
US9040924B2 (en) 2009-10-27 2015-05-26 University Of Washington Through Its Center For Commercialization Optical-interface patterning for radiation detector crystals
KR101041020B1 (ko) * 2010-02-09 2011-06-13 서울대학교산학협력단 위치민감형 다중양극 광전자증배관의 이득 불균일 보정 시스템 및 방법
JP5458976B2 (ja) * 2010-03-10 2014-04-02 株式会社島津製作所 放射線信号処理装置、およびそれを備えた放射線検出器
GB201004922D0 (en) * 2010-03-24 2010-05-12 Sensl Technologies Ltd Silicon photomultiplier and readout method
JP5835990B2 (ja) * 2010-08-09 2015-12-24 株式会社東芝 核医学イメージング装置および核医学イメージング方法
US20120041695A1 (en) 2010-08-16 2012-02-16 Csi Technology, Inc. Integrated vibration measurement and analysis system
CN102772217B (zh) * 2011-05-13 2014-07-16 上海生物医学工程研究中心 一种针对pet符合***的测试方法及装置
CN102338880B (zh) * 2011-06-02 2012-12-26 中国科学技术大学 核脉冲幅度数字化方法及***
JP6012721B2 (ja) * 2011-06-14 2016-10-25 サザン イノヴェーション インターナショナル プロプライアトリー リミテッド 検出器出力データ中のパルスを特定するための方法および装置
KR101217867B1 (ko) * 2011-09-29 2013-01-02 한국원자력의학원 방사성분자체 및 양전자방출단층촬영기를 이용한 내부 장기의 움직임 계측 장치 및 그 방법
US9244179B2 (en) 2011-12-06 2016-01-26 The Board Of Trustees Of The Leland Stanford Junior University Time-based digitizer for PET photodetector
US9140805B2 (en) * 2012-10-22 2015-09-22 Kabushiki Kaisha Toshiba Apparatus and method for improving uniformity of performance in positron emission tomography
WO2014074666A1 (en) * 2012-11-07 2014-05-15 Massachusetts Institute Of Technology Inter-detector scatter enhanced emission tomography
KR101912715B1 (ko) 2012-11-20 2018-10-29 삼성전자주식회사 방사선이 방출된 위치의 분포를 추정하는 방법 및 장치
KR102026737B1 (ko) 2013-01-25 2019-09-30 삼성전자주식회사 영상 생성 장치 및 방법
KR101964891B1 (ko) * 2013-01-28 2019-08-07 삼성전자주식회사 실리콘 광증배관 디텍터 셀
CN103099639B (zh) * 2013-03-01 2014-12-24 江苏中惠医疗科技股份有限公司 Pet成像***的环形拓扑结构及其实现方法
US9442201B2 (en) * 2013-09-12 2016-09-13 Siemens Medical Solutions Usa, Inc. CMOS SPAD array with mixed timing pick-off for time-of-flight positron emission tomography
CN103607205B (zh) * 2013-11-12 2017-02-08 沈阳东软医疗***有限公司 一种信号处理方法、装置及设备
CN104656115B (zh) * 2013-11-19 2018-04-24 苏州瑞派宁科技有限公司 一种时间标记组合的方法与***
KR102347023B1 (ko) * 2013-12-11 2022-01-03 서던 이노베이션 인터내셔널 피티와이 엘티디 데이터에서의 신호들을 분해하기 위한 방법 및 장치
CN104688261B (zh) * 2015-03-19 2017-06-16 明峰医疗***股份有限公司 基于FPGA和ADC的全数字化SiPM PET数据采集***
US9606245B1 (en) 2015-03-24 2017-03-28 The Research Foundation For The State University Of New York Autonomous gamma, X-ray, and particle detector
GB2536650A (en) 2015-03-24 2016-09-28 Augmedics Ltd Method and system for combining video-based and optic-based augmented reality in a near eye display
AT517499B1 (de) * 2015-08-12 2018-02-15 Avl List Gmbh Verfahren und Vorrichtung zur Detektion von Signalpulsen
ITUB20154173A1 (it) 2015-10-01 2017-04-01 Datalogic IP Tech Srl Sensore optoelettronico e metodo di funzionamento di un sensore optoelettronico
WO2017091697A1 (en) * 2015-11-23 2017-06-01 Berr Stuart S Positron emission tomography systems and methods
US9939536B2 (en) 2016-02-19 2018-04-10 Sensi Technologies Ltd. Semiconductor photomultiplier with baseline restoration for a fast terminal signal output including output loads to correct an overshoot of an output signal (as amended)
US10690735B2 (en) 2016-04-26 2020-06-23 Aivitae LLC Wireless detection coil system
WO2017214766A1 (zh) * 2016-06-12 2017-12-21 上海联影医疗科技有限公司 正电子发射断层成像***及其图像重建方法
CN106372403A (zh) * 2016-08-29 2017-02-01 合肥菲力姆数码科技有限公司 一种基于无创成像的医用影像输出***
CN106443757B (zh) * 2016-09-23 2018-12-07 沈阳东软医疗***有限公司 一种获取湮灭事件发生时间的装置及探测器
EP3529628B1 (en) * 2016-10-24 2022-03-16 Xia LLC Ratio-reference measurement of the arrival time and/or amplitude of a digitized electronic pulse
US10416295B2 (en) 2016-10-24 2019-09-17 Xia Llc Interpolation measurement of the arrival time and/or amplitude of a digitized electronic pulse
US10817588B2 (en) 2016-10-24 2020-10-27 Xia Llc Ratio-reference measurement of the arrival time and/or amplitude of a digitized electronic pulse
US10673204B2 (en) 2017-03-07 2020-06-02 Sensl Technologies Ltd. Laser driver
CN107024711B (zh) * 2017-04-17 2019-02-26 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化信号的拟合方法
CN109816740B (zh) * 2017-11-18 2020-10-16 苏州瑞派宁科技有限公司 一种闪烁脉冲事件的符合处理方法
CN108152767A (zh) * 2017-11-30 2018-06-12 华东师范大学 一种基于fpga的磁共振信号实时处理方法
US11980507B2 (en) 2018-05-02 2024-05-14 Augmedics Ltd. Registration of a fiducial marker for an augmented reality system
CO2018007716A1 (es) 2018-07-25 2020-01-31 Univ Antonio Narino Receptor para señales ópticas de baja potencia con funcionamiento en condiciones de alta incidencia de luz de fondo y aplicación en comunicación con luz visible
US11766296B2 (en) 2018-11-26 2023-09-26 Augmedics Ltd. Tracking system for image-guided surgery
US11172906B2 (en) * 2019-01-18 2021-11-16 Parto Negar Persia (PNP) Company Normalization of a positron emission tomography scanner
US11061147B2 (en) 2019-03-01 2021-07-13 University Of Washington Accurate photon depth-of-interaction decoding and calibration of multiplexed detector modules
CN110226943B (zh) * 2019-07-05 2023-08-15 上海联影医疗科技股份有限公司 光子到达探测器的参数计算方法、装置和计算机设备
US11980506B2 (en) 2019-07-29 2024-05-14 Augmedics Ltd. Fiducial marker
US11382712B2 (en) 2019-12-22 2022-07-12 Augmedics Ltd. Mirroring in image guided surgery
CN111158039B (zh) * 2020-01-02 2022-01-04 苏州瑞派宁科技有限公司 信号采样、重建方法及装置
US11296713B1 (en) * 2020-08-10 2022-04-05 Pacific MicroCHIP Corp. Circuitry for event-driven data acquisition
CN113433579B (zh) * 2021-05-18 2023-01-20 中国工程物理研究院激光聚变研究中心 一种大灵敏面x射线光谱平响应二极管探测器
US11896445B2 (en) 2021-07-07 2024-02-13 Augmedics Ltd. Iliac pin and adapter

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005164334A (ja) * 2003-12-01 2005-06-23 Toshiba Corp 核医学診断装置
US7129497B2 (en) * 2004-05-11 2006-10-31 Scott David Wollenweber Method and system for normalization of a positron emission tomography system

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0173419A1 (en) * 1984-06-05 1986-03-05 CODMAN & SHURTLEFF INC. Transcutaneous nerve stimulation device using a common controller for pulse production and parameter display
JP2577386B2 (ja) * 1987-07-09 1997-01-29 株式会社東芝 放射線エネルギスペクトル測定装置
US5367209A (en) * 1992-01-07 1994-11-22 Hauck Scott A Field programmable gate array for synchronous and asynchronous operation
US5894565A (en) * 1996-05-20 1999-04-13 Atmel Corporation Field programmable gate array with distributed RAM and increased cell utilization
US5751000A (en) * 1997-01-08 1998-05-12 Smv America, Inc. Prefilter collimator for PET gamma camera
JP3103047B2 (ja) * 1997-09-01 2000-10-23 シュランベルジェ、ホールディング、リミテッド 核分光システムにおけるパルス波形調整及び弁別方法
US6493646B1 (en) * 2000-02-16 2002-12-10 Ge Medical Systems Global Technology Company, Llc High order primary decay correction for CT imaging system detectors
US7260044B2 (en) * 2002-09-06 2007-08-21 Ricoh Company, Ltd. Recording method for a phase-change optical recording medium
US7045802B2 (en) * 2003-11-26 2006-05-16 General Electric Company Method and apparatus for coincidence imaging digital triggering
EP1706759A1 (en) * 2004-01-13 2006-10-04 Koninklijke Philips Electronics N.V. Analog to digital conversion shift error correction
US7630528B2 (en) * 2004-03-02 2009-12-08 Koninklijke Philips Electronics N.V. Motion compensation
WO2006029475A1 (en) * 2004-09-16 2006-03-23 Southern Innovation International Pty Ltd Method and apparatus for resolving individual signals in detector output data.
JP4611106B2 (ja) * 2005-03-11 2011-01-12 住友重機械工業株式会社 放射線検出回路及び放射線検査装置
KR100715803B1 (ko) * 2005-04-15 2007-05-10 재단법인서울대학교산학협력재단 3개 이상의 결정층이 형성된 섬광검출기 및 이를 이용한양전자 방출 단층촬영장치
JP2008539414A (ja) * 2005-04-27 2008-11-13 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 心動態のモデル化におけるｅｃｇゲートの時間サンプリング
US7211799B2 (en) * 2005-09-14 2007-05-01 General Electric Company Method and system for calibrating a time of flight positron emission tomography system
US7626171B2 (en) * 2006-01-09 2009-12-01 Koninklijke Philips Electronics N.V. Method of constructing time-in-flight pet images
JP4844196B2 (ja) * 2006-03-30 2011-12-28 株式会社日立製作所 放射線検出回路
US7736319B2 (en) * 2007-01-19 2010-06-15 Cardiac Pacemakers, Inc. Ischemia detection using heart sound timing

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005164334A (ja) * 2003-12-01 2005-06-23 Toshiba Corp 核医学診断装置
US7129497B2 (en) * 2004-05-11 2006-10-31 Scott David Wollenweber Method and system for normalization of a positron emission tomography system

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ABU-AITA, Rami Rimon, "Count Rate and Spatial Resolution Performance of a 3-Dimensional Dedicated Positron Emission Tomography Scanner", Master's Thesis of University of Florida, 2004, see the whole document *
HASELMAN, Michael D. et al., "FPGA-Based Data Acquisition System for a Positron Emission Tomography Scanner", 16th ACM/SIGDA International Symposium on Field-Programmable Gate Arrays, Feb. 24-26, 2008, p. 264, see the whole document *
HASELMAN, Michael D. et al., "FPGA-Based Front-End Electronics for Positron Emission Tomography", 17th ACM/SIGDA International Symposium on Field-Programmable Gate Arrays, Feb. 22-24, 2009, pp. 93-102, see the whole document *
HASELMAN, Michael D. et al., "Simulation of Algorithms for Pulse Timing in FPGAs", 2007, IEEE Science Symposium Conf. Rec, Oct. 26 - Nov. 3, 2007, pp 3161-3165, see the whole document *
STREUN, M. et al., "Effects of Crosstalk and Gain Nonuniformity Using Multichannel PMTs in the Clearpet Scanner", Nucl. Inst., and Methods in Physics Research A, 2005, vol. 537, pp. 402-405, see the whole document *

Also Published As

Publication number Publication date
EP2208088B1 (en) 2019-05-08
JP2011503550A (ja) 2011-01-27
AU2008318371A1 (en) 2009-05-07
CN101903798B (zh) 2012-12-12
US20110301918A1 (en) 2011-12-08
EP2208088A4 (en) 2014-10-29
CN101903798A (zh) 2010-12-01
KR101583272B1 (ko) 2016-01-07
US20090224158A1 (en) 2009-09-10
EP2208088A2 (en) 2010-07-21
EA201070554A1 (ru) 2010-12-30
AU2008318371B2 (en) 2013-10-17
US8003948B2 (en) 2011-08-23
WO2009059312A2 (en) 2009-05-07
KR20100103490A (ko) 2010-09-27
CA2706388A1 (en) 2009-05-07
CA2706388C (en) 2016-12-06
WO2009059312A3 (en) 2009-06-18
US8309932B2 (en) 2012-11-13
JP5667446B2 (ja) 2015-02-12

Similar Documents

Publication Publication Date Title
EA016862B1 (ru) Способ и устройство для оценки начального момента времени электронного импульса, предназначенного для использования в позитронно-эмиссионной томографии
Vaquero et al. Positron emission tomography: current challenges and opportunities for technological advances in clinical and preclinical imaging systems
RU2381525C2 (ru) Система позитрон-эмиссионной томографии
US8716669B2 (en) Line of response estimation for high-resolution PET detector
RU2582887C2 (ru) Калибровка рет с изменяющимися интервалами совпадений
CN110772274A (zh) 一种多核素的伽马成像***与方法
US8340377B2 (en) Method for energy calculation and pileup determination for continuously sampled nuclear pulse processing
Wang et al. Characterization of the LBNL PEM camera
US8507842B2 (en) Method for identifying 3-D location of gamma interaction and flat panel gamma imaging head apparatus using the same
US9804206B2 (en) Method and a device for measuring parameters of an analog signal
JP2000131440A (ja) 放射線検出処理システム
US9804279B2 (en) Method for determining parameters of a reaction of a gamma quantum within a scintillator of a PET scanner
JP7019286B2 (ja) データ収集装置及びｘ線ｃｔ装置
Yamada et al. Development of a small animal PET scanner using DOI detectors
D'ascenzo et al. Recent advances in digital positron emission tomography
Bugalho et al. Experimental characterization of the Clear-PEM scanner spectrometric performance
Albuquerque et al. Characterization of the Clear-PEM breast imaging scanner performance
Trindade et al. Breast cancer imaging studies by Monte Carlo simulation with Clear-PEM
CN210697662U (zh) 一种契连柯夫事件与伽马事件符合成像装置
Ziegler Development of a small-animal PET system
Ay et al. Pet/Ct
Momsen Digital Waveform Acquisition and Adaptive Collimation Methods for Spect
Wang Calibration and characterization of the LBNL positron emission mammography camera
Tai PET/microPET imaging
Pani et al. A Dual-Modality Ultrasound-Scintigraphic Detector for Molecular Imaging in Medicine

Legal Events

Date Code Title Description
MM4A Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s)

Designated state(s): BY KZ KG MD TJ TM

MM4A Lapse of a eurasian patent due to non-payment of renewal fees within the time limit in the following designated state(s)

Designated state(s): AM AZ RU