RU2575943C2 - Получение данных компьютерной томографии (ст) - Google Patents

Получение данных компьютерной томографии (ст) Download PDF

Info

Publication number
RU2575943C2
RU2575943C2 RU2013137507/28A RU2013137507A RU2575943C2 RU 2575943 C2 RU2575943 C2 RU 2575943C2 RU 2013137507/28 A RU2013137507/28 A RU 2013137507/28A RU 2013137507 A RU2013137507 A RU 2013137507A RU 2575943 C2 RU2575943 C2 RU 2575943C2
Authority
RU
Russia
Prior art keywords
integration
integration period
gantry
timing
signal
Prior art date
Application number
RU2013137507/28A
Other languages
English (en)
Other versions
RU2013137507A (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 RU2013137507A publication Critical patent/RU2013137507A/ru
Application granted granted Critical
Publication of RU2575943C2 publication Critical patent/RU2575943C2/ru

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5294Devices using data or image processing specially adapted for radiation diagnosis involving using additional data, e.g. patient information, image labeling, acquisition parameters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/545Control of apparatus or devices for radiation diagnosis involving automatic set-up of acquisition parameters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/547Control of apparatus or devices for radiation diagnosis involving tracking of position of the device or parts of the device

Landscapes

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

Abstract

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

Description

ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ
Приведенное ниже относится, в целом, к получению данных системы визуализации, и описано с конкретным применением для компьютерной томографии (CT).
УРОВЕНЬ ТЕХНИКИ
Сканер компьютерной томографии (CT) включает в себя поворотную часть, поддерживаемую с возможностью поворота неподвижной частью. Поворотная часть поддерживает рентгеновскую трубку, которая испускает излучение, которое проходит через область исследования и, в ходе этого, через часть объекта или субъекта, и матрицу детекторов, которая обнаруживает излучение, проходящее через область исследования, и генерирует данные проекции, показывающие обнаруженное излучение. Реконструктор восстанавливает данные проекции и генерирует объемные данные изображения, показывающие часть объекта или субъекта в области исследования. На основании данных изображения может быть сгенерировано одно или более изображений.
Матрица детекторов включает в себя пиксели детекторов, которые преобразуют обнаруженные фотоны рентгеновского излучения в электрические сигналы, указывающие на них. Для каждого оборота поворотного гентри пиксели детекторов обнаруживают и преобразуют фотоны рентгеновского излучения для множества периодов интегрирования, каждый из которых соответствует различному диапазону угловых положений. Временная продолжительность периода интегрирования зависит от скорости поворота поворотного гентри и числа периодов интегрирования для каждого оборота сканирования. При использовании набора интегрирующих детекторов, в начале каждого периода интегрирования интеграторы для пикселей детекторов перезапускаются, и затем интеграторы принимают и интегрируют электрические сигналы за периоды интегрирования. Интегрированные сигналы формируют данные проекции для этого периода интегрирования.
Реконструктор обрабатывает данные проекции для заданного периода интегрирования как соответствующие конкретному диапазону угловых положений. Однако если временная привязка периода интегрирования относительно углового положения является неточной, то фактические данные проекции соответствуют другому диапазону угловых положений, и это может привести к снижению качества изображения для данных восстановленного изображения. К сожалению, имеется естественная погрешность положения поворотного гентри, в которой преобладающим нарушением скорости поворотного гентри является дисбаланс, который является функцией угла поворотного гентри, и величина угловой ошибки вследствие дисбаланса является функцией, обратно пропорциональной кубу скорости поворота.
Устройство высокоточного определения положения поворотного гентри большого диаметра использовалось для определения фактического положения поворотного гентри, которое использовалось для синхронизации временной привязки периода интегрирования и диапазона угловых положений поворотного гентри. Такое устройство включало в себя кодировщик и датчик положения высокого разрешения, вспомогательную электронику, источник питания, и контактное кольцо для передачи данных между поворотной и неподвижной частями сканера. К сожалению, такие устройства обычно являются сложными и дорогими и увеличивают общую стоимость сканера.
СУЩНОСТЬ ИЗОБРЕТЕНИЯ
Аспекты настоящей заявки решают указанные выше и другие вопросы.
В соответствии с одним из аспектов, система визуализации включает в себя неподвижный гентри и поворотный гентри, поддерживаемый с возможностью вращения неподвижным гентри и выполненный с возможностью вращения вокруг области исследования. Система также содержит источник излучения, находящийся на поворотном гентри, который выполнен с возможностью испускания излучения, которое проходит через область исследования, и матрицу детекторов, удерживаемую поворотным гентри и расположенную напротив источника излучения через область исследования, которая обнаруживает излучение, которое проходит через область исследования. Матрица детекторов обнаруживает излучение для множества периодов интегрирования во время вращения поворотного гентри, при этом множество периодов интегрирования соответствует различным диапазонам угловых положений, и матрица детекторов генерирует сигнал, указывающий на обнаруженное излучение, соответственно, для множества периодов интегрирования. Система также содержит контроллер периодов интегрирования, который генерирует сигнал временной привязки периода интегрирования, который включает в себя временную привязку для начала каждого периода интегрирования для оборота поворотного гентри, на основании по меньшей мере временной продолжительности предыдущего оборота поворотного гентри вокруг области исследования, при этом сигнал временной привязки интегрирования используется для инициации множества периодов интегрирования.
В соответствии с другим аспектом, способ включает в себя генерацию сигнала временной привязки периода интегрирования, включающего в себя временную привязку множества периодов интегрирования получения данных для оборота поворотного гентри системы визуализации, при этом сигнал временной привязки периода интегрирования генерируется на основании по меньшей мере временной продолжительности предыдущего оборота поворотного гентри.
В соответствии с другим аспектом, способ включает в себя корректировку предсказанной временной привязки получения данных периода интегрирования по отношению к поворотному гентри системы визуализации из-за дисбаланса поворотного гентри.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Изобретение может принимать форму различных компонентов и компоновок компонентов, а также различных этапов и расположения этапов. Чертежи предназначены только для целей иллюстрации предпочтительных вариантов осуществления и не должны рассматриваться как ограничивающие изобретение.
Фигура 1 схематически иллюстрирует систему визуализации в связи с контроллером периодов интегрирования.
Фигура 2 схематически иллюстрирует типовой контроллер периодов интегрирования.
На фигуре 3 показан примерный график временной привязки IP и ошибки IP как функции скорости поворота поворотного гентри.
Фигура 4 схематически иллюстрирует примерную систему визуализации, использующую один датчик и одну метку для способствования определению сигнала временной привязки IP.
На фигуре 5 показана примерная временная диаграмма датчика с фигуры 4.
Фигура 6 схематически иллюстрирует примерную систему визуализации, использующую множество датчиков и множество меток для способствования определению сигнала временной привязки IP и его корректировке.
На фигуре 7 показана примерная диаграмма временной привязки для датчиков с фигуры 6.
На фигуре 8 показана примерная кривая ошибок, представляющая ошибку IP как функцию углового положения.
Фигура 9 иллюстрирует примерный способ для определения сигнала временной привязки IP.
Фигура 10 иллюстрирует примерный способ для определения скорректированного с учетом дисбаланса гентри сигнала временной привязки IP.
ПОДРОБНОЕ ОПИСАНИЕ ВАРИАНТОВ ОСУЩЕСТВЛЕНИЯ
Фигура 1 схематически иллюстрирует систему визуализации, такую как сканер 100 компьютерной томографии (CT). Сканер 100 включает в себя неподвижный гентри 102 и поворотный гентри 104, который удерживается с возможностью поворота неподвижным гентри 102. Поворотный гентри 104 поворачивается вокруг области исследования 106 вокруг продольной или z-оси 108 один или более раз в течение одного или более циклов получения данных. Опора 114 пациента, такая как кушетка, поддерживает пациента в области исследования 106.
Источник 110 излучения, такой как рентгеновская трубка, поддерживается поворотным гентри 104 и поворачивается вместе с ней вокруг области исследования 106. Источник 110 излучения испускает излучение, которое коллимируется коллиматором источника с целью выдачи, как правило, веерообразного, клиновидного или конического пучка излучения, который проходит через область исследования 106. Матрица 112 чувствительных к излучению детекторов включает в себя одно- или двумерную матрицу пикселей детекторов, которые, соответственно, обнаруживают излучение, которое проходит через область исследования 106, и генерируют электрические сигналы (например, ток или напряжение), указывающие на обнаруженное излучение.
Интегратор 116, который включает в себя один или более субинтеграторов, таких как субинтегратор для каждого пикселя детектора матрицы детекторов, интегрирует электрические сигналы от множества пикселей детекторов, соответственно, для множества периодов интегрирования (видов, проекций и т.д.), при этом каждый период интегрирования соответствует различному диапазону угловых положений при обороте поворотного гентри 104, для одного или более оборотов поворотного гентри 104 во время получения данных, и выдает данные проекции.
Контроллер 118 периодов интегрирования (IP) управляет интегратором 116. В качестве примера, проиллюстрированный контроллер IP 118 инициирует перезапуск интегратора 116 в начале каждого периода интегрирования наряду с передачей произвольного интегрированного сигнала для предыдущего периода интегрирования.
Как подробнее описано ниже, контроллер IP 118 определяет, когда перезапустить интегратор 116, и выдает сигнал, который перезапускает интегратор 116 на основании информации о предыдущем обороте поворотного гентри 104, полученной до (например, при предварительном сканировании) и/или во время сканирования. В проиллюстрированном варианте осуществления, определитель 120 информации оборота гентри определяет такую информацию. В одном из случаев определенная информация оборота гентри позволяет осуществить синхронизацию временной привязки периодов интегрирования (IP) (то есть, интегратор перезапускается) с диапазоном угловых положений поворотного гентри 104 с корректированием и/или без корректирования сигналов временной привязки IP из-за дисбаланса поворотного гентри.
Следует понимать, что такая синхронизация может быть достигнута без использования дополнительного устройства определения положения поворотного гентри. Таким образом, контроллер IP 118 может способствовать сохранению или снижению общей стоимости сканера, например, посредством избегания включения высокоточного устройства определения положения поворотного гентри в систему визуализации 100, при этом достигая предварительно заданного качества изображения (например, заданного пространственного разрешения и/или избегания внесения артефактов, вызываемого неточностями временной привязки) относительно конфигурации, в которой контроллер IP 118 не используется. Естественно, также может быть использовано устройство определения положения поворотного гентри.
Реконструктор 122 восстанавливает данные проекции и генерирует объемные данные изображения, показывающие область исследования 106. Вычислительная система общего назначения служит в качестве консоли 124 оператора, и включает в себя устройство вывода, такое как дисплей, и устройство ввода, такое как клавиатура, мышь и/или т.п. Программное обеспечение, имеющееся в консоли 124, позволяет оператору управлять функционированием системы 100, например, позволяя оператору начать сканирование и т.д.
Несмотря на то, что интегратор 116 и контроллер IP 118 показаны как блоки, отдельные от матрицы детекторов 112, следует понимать, что один или более элементов из интегратора 116 или контроллера IP 118 могут являться частью матрицы детекторов 112. Кроме того, определитель 120 информации оборота гентри может быть частью поворотного и/или стационарного гентри 104 и 102.
Более того, следует понимать, что по меньшей мере субчасть контроллера IP 118 может быть реализована через один или более процессоров, выполняющих читаемые компьютером инструкции, закодированные на машиночитаемом носителе информации.
Дополнительно или альтернативно, одна или более читаемых компьютером инструкций могут переноситься сигналом, несущей волной и т.п.
Фигура 2 иллюстрирует пример контроллера IP 118 в связи с определителем 120 информации оборота гентри, консолью 124 и интегратором 116.
Определитель 120 информации оборота гентри определяет информацию временной привязки о вращении поворотного гентри 104. Обычно соответствующий определитель 120 информации оборота гентри может содержать один или более датчиков, прикрепленных к поворотному гентри 104 (или неподвижному гентри 102) и одну или более меток, прикрепленных к неподвижному гентри 102 (или поворотному гентри 104), при этом датчик(-и) и метка(-и) расположены таким образом, что между ними имеет место относительное перемещение, когда поворотный гентри 104 поворачивается, и датчик обнаруживает метку, проходящую в процессе поворота через датчик, и генерирует сигнал, указывающий на обнаруженную метку. Таймер и т.п. может быть использован для определения объема времени между обнаружением меток для одного или более датчиков. Неограничивающий пример реализации подробнее описан ниже.
Предсказатель 202 временной привязки IP предсказывает сигнал временной привязки IP для конкретного сканирования. В проиллюстрированном варианте осуществления предсказатель 202 временной привязки IP предсказывает сигнал временной привязки IP для оборота на основании информации оборота гентри для предыдущего оборота и информации, указывающей скорость поворота поворотного гентри 104, и количества периодов интегрирования (IP) для представляющего интерес протокола сканирования. В качестве примера, в одном из случаев предсказатель 202 временной привязки IP делит объем времени для предыдущего оборота на скорость поворота и количество IP для вычисления сигнала временной привязки IP, указывающего на множество IP, находящихся на равном угловом расстоянии друг от друга. Подходящая временная продолжительность или ширина IP включает, но не ограничивается, ширину временного интервала в диапазоне от около пятидесяти микросекунд (50 мкс) до около пятисот микросекунд (500 мкс). Естественно, в настоящем описании предусматриваются и другие значения ширины временного интервала IP.
В проиллюстрированном варианте осуществления предсказатель 202 временной привязки IP может получать скорость поворота гентри и количество IP с консоли 124. Если сканер 100 включает в себя протоколы с одними и теми же скоростью поворота гентри и количеством IP, то скорость поворота гентри и количество IP может дополнительно или альтернативно храниться в памяти контроллера IP 118 и/или быть доступными ему. В другом варианте осуществления контроллер IP 118 дополнительно или альтернативно хранит и/или может иметь доступ к памяти, в которой хранится соответствие (например, справочная таблица (LUT)) между конкретными протоколами сканирования и скоростью поворота гентри и количеством IP для каждого протокола сканирования, и использует соответствие для определения поворота гентри и количества IP для идентифицированного протокола, например, в информации из консоли 124 и/или другого компонента.
Определитель 204 ошибки временной привязки IP генерирует сигнал ошибки IP, показывающий колебания скорости поворота поворотного гентри 104, например, вследствие дисбаланса поворотного гентри 104 и/или по другой причине, которые могут привести к неточности временной привязки между предсказанным сигналом временной привязки IP и фактическим диапазоном угловых положений поворотного гентри 104. В проиллюстрированном варианте осуществления определитель 204 ошибки временной привязки IP генерирует сигнал ошибки на основании информации от определителя 120 информации оборота гентри и/или предсказателя 202 временной привязки IP. В одном из случаев определитель 204 ошибки временной привязки IP определяет разницу между предсказанной временной привязкой IP и измеренными фактическими IP для подмножества угловых положений, определяет значения ошибок на основании разностей, и использует значения ошибок для генерации корректировки для предсказанного сигнала IP. Типовая реализация определителя 204 ошибки временной привязки IP подробнее описана ниже.
Корректор 206 временной привязки IP корректирует предсказанный сигнал временной привязки IP от предсказателя 202 временной привязки IP на основании сигнала ошибки временной привязки IP от определителя 204 ошибки временной привязки IP. В одном неограничивающем варианте осуществления корректор 206 временной привязки IP может просто суммировать или объединять другим способом сигнал временной привязки IP и сигнал ошибки временной привязки IP.
Генератор 208 импульсов IP генерирует импульсы перезапуска интегратора на основании скорректированного сигнала временной привязки IP, который перезапускает интегратор 116 для каждого периода интегрирования для каждого поворота поворотного гентри 104 во время сканирования. В проиллюстрированном варианте осуществления генератор 208 импульсов IP генерирует поток последовательных импульсов, который инициирует передачу данных от интегратора 116 или перезапускает его. На выходе интегратора 116 выдаются данные проекции, которые восстанавливаются реконструктором 122, как обсуждалось выше.
Также предусматриваются вариации.
Если ошибка IP находится в пределах предварительно определенного диапазона (например, пятьдесят (50) угловых секунд или менее), то определитель 204 ошибки временной привязки IP и корректор 206 временной привязки IP может быть исключен или может не использоваться. В этом случае предсказанный сигнал временной привязки от предсказателя 202 временной привязки IP выдается и используется генератором 208 импульсов IP для выдачи импульсов перезапуска интегратора.
На фигуре 3 показаны примеры кривой 300 временной привязки IP и кривой 302 ошибок IP как функции скорости поворота поворотного гентри 104, в которых y-ось 304 представляет собой ширину временного интервала IP, и x-ось 306 представляет собой скорость поворота поворотного гентри. Как показано на фигуре 3, временная привязка IP и ошибка IP снижается с возрастанием скорости поворота поворотного гентри. Таким образом, определитель 204 ошибки временной привязки IP и корректор 206 временной привязки IP могут использоваться, не включаться или не использоваться при более высоких скоростях поворота, и могут использоваться при более низких скоростях поворота для корректирования ошибки IP.
В другой вариации значения ошибок IP для сканера определяются заранее, например, во время калибровки или другой процедуры, и хранятся в таблице, такой как справочная таблица (LUT), со сканером. В данном случае LUT может использоваться корректором 206 временной привязки IP для корректировки сигнала временной привязки IP от предсказателя 202 временной привязки IP.
В другой вариации, адаптивный алгоритм вычисления ошибки может использоваться для нестационарных ошибок. Это было бы полезным, например, если в системе возникает непредсказанное нарушение, такое как меняющийся дисбаланс или поправки крутящего момента двигателя во время получения.
В другой вариации, стандартный кодировщик двигателя может быть использован для осуществления движения поворотного гентри с постоянной скоростью, что может скомпенсировать или снизить крутящий момент дисбаланса, что может привести к снижению ошибки временной привязки IP.
В другой вариации, один цикл на амплитуду оборота и фазу тока мотора может быть использован для измерения неточности IP и улучшения временной привязки IP аналогично и, возможно, совместно с дополнительными датчиками положения, как описано выше. Крутящий момент может быть использован после выполнения получения данных для улучшения временной привязки.
Фигуры 4 и 5 иллюстрируют неограничивающий пример, в котором определитель 204 ошибки временной привязки IP и корректор 206 временной привязки IP опускаются или не используются.
Обратимся сначала к фигуре 4; датчик 402 прикреплен к поворотному гентри 104, и метка 404 прикреплена к неподвижному гентри 102. Датчик 402 и метка 404 размещены друг относительно друга таким образом, что датчик 402 проходит в процессе поворота через метку 404 при каждом обороте поворотного гентри 104. Датчик 402 сконфигурирован для обнаружения метки 404 при повороте датчика 402 с прохождением метки 404, и для генерации выходного сигнала, указывающего на эту метку.
На фигуре 5 показан пример выхода датчика 402. На фигуре 5 x-ось 502 представляет собой время, и y-ось 504 представляет собой выход датчика 402. Выходной сигнал 506 представляет собой выход датчика 402. В отношении выходного сигнала 506, первый выходной импульс 508 соответствует прохождению датчика 402 через метку 404 для оборота, и второй выходной импульс 510 соответствует прохождению датчика 402 через метку 404 для следующего оборота.
Продолжим рассмотрение фигур 4 и 5; определитель 120 информации оборота гентри включает в себя таймер 406. В данном примере таймер 406 перезапускается и начинает отсчет времени оборота на нарастающем фронте импульса 508. На нарастающем фронте следующего импульса 510, таймер 406 выдает сигнал 516 времени, который показывает время оборота, предсказателю 202 временной привязки IP, и таймер 406 перезапускается и начинает отсчет времени следующего оборота.
Предсказатель 202 временной привязки IP определяет сигнал временной привязки IP для следующего оборота, как описано выше, например, посредством деления сигнала 516 времени на количество IP и скорость поворота поворотного гентри 104 для следующего оборота, которые могут быть получены с консоли 124 или другим способом. Сигнал временной привязки IP выдается генератору 208 импульсов IP, который использует сигнал временной привязки IP для генерации импульсов перезапуска интегратора 116 для оборота.
Неограничивающая реализация может включать в себя счетчик, который инициируется выходным сигналом датчика 402 и который определяет временную продолжительность оборота, триггер-защелку, которая также инициируется выходом датчика 402 и в которой задерживается сигнал временной привязки IP (который рассчитывается путем деления временной продолжительности на произведение количества IP и скорости поворота), и таймер, который выдает поток последовательных импульсов инициации на основании сигнала временной привязки IP. Сумматор может быть подсоединен до таймера с целью добавления сигнала корректировки ошибки IP к сигналу временной привязки IP.
Фигуры 6 и 7 иллюстрируют неограничивающий пример, в котором определитель 204 ошибки временной привязки IP и корректор 206 временной привязки IP используются для корректировки предсказанного сигнала временной привязки IP. Вспомним, что при более низких скоростях поворота может иметься необходимость корректировки предсказанного сигнала временной привязки IP относительно, например, вариаций в угловом положении вследствие дисбаланса гентри и/или по другим причинам.
Обратимся сначала к фигуре 6; по меньшей мере первый датчик 602 и второй датчик 604 прикреплены к поворотному гентри 104, и по меньшей мере первая метка 606 и вторая метка 608 прикреплены к неподвижному гентри 102 (или поворотному гентри 104). По меньшей мере две метки 606 и 608 имеют угловое смещение друг относительно друга на угловое расстояние (α) 610. В проиллюстрированном варианте осуществления α составляет около 90 градусов. По меньшей мере два датчика 602 и 604 имеют угловое смещение друг относительно друга на угловое расстояние (β) 612. В проиллюстрированном варианте осуществления β составляет около ста восьмидесяти градусов.
Следует понимать, что проиллюстрированный вариант осуществления не является ограничивающим. В других вариантах осуществления может присутствовать больше или меньше датчиков и/или меток, и расположение датчиков и/или меток друг относительно друга и поворотном и/или неподвижном гентри может быть различным.
На фигуре 7 показан примерный выходной сигнал датчиков 602 и 604. На фигуре 7 x-ось 702 представляет время, и y-ось 704 представляет выходной сигнал по меньшей мере датчиков 602 и 604. Первый выходной сигнал 706 представляет выход датчика 602 (или датчика 604), и второй выходной сигнал 708 представляет выход датчика 604 (или датчика 602).
В отношении выходного сигнала 706, первый выходной импульс 710 соответствует прохождению первого датчика 602 через первую метку 606, и второй выходной сигнал 712 соответствует прохождению первого датчика 602 через вторую метку 608. В отношении выходного сигнала 708, первый выходной импульс 714 соответствует прохождению второго датчика 604 через первую метку 606, и второй выходной сигнал 716 соответствует прохождению второго датчика 604 через вторую метку 608. Импульс 711 соответствует прохождению первого датчика 602 через первую метку 606 для последующего оборота.
В данном примере, и аналогично примеру с фигур 4-6, таймер 406 используется для определения временной продолжительности 718 оборота на основании разницы во времени между импульсом 710 и импульсом 711. Таймер 406 также может быть использован для определения временного расстояния 720 между метками 606 и 608 на основании разницы во времени между импульсами 710 и импульсом 712, и временного расстояния 722 между датчиками 602 и 604 на основании разницы между импульсами 714 и импульсом 716.
В проиллюстрированном варианте осуществления, метка 606 присваивается опорному угловому положению ноль (0). Угловое расстояние (α) между меткой 606 и меткой 608 может быть определено как отношение временного расстояния 720 между метками 606 и 608 и временного расстояния 718 одного оборота. Угловое расстояние (β) между датчиками 602 и 604 может быть определено как отношение временного расстояния 722 между датчиками 602 и 604 при прохождении метки 608 и временного расстояния 718 одного оборота. Угловое расстояние между датчиком 602 при прохождении метки 606 и датчиком 604 при прохождении метки 608 может быть определено через суммирование угловых расстояний α и β (α+β). Моменты времени α и β можно назвать откалиброванными моментами времени измерений, и они могут быть определены путем измерения времени между метками во время состояния минимального нарушения, которое обычно имеет место при максимальной скорости ротора.
Определитель 204 ошибки временной привязки IP (фигура 2) может определить ошибку IP для предсказанного сигнала временной привязки IP на основании предсказанной временной привязки IP и измеренной временной привязки IP 718-722 для угловых положений 0, α, β и α+β, как показано в уравнении 1:
Уравнение 1:
Figure 00000001
где Θ представляет собой угловые положения 0, α, β и α+β, T^predicted представляет собой предсказанную временную привязку IP для угловых положений Θ, T^measured представляет собой измеренную временную привязку для угловых положений Θ, I является неизвестным и представляет собой инерциальную амплитуду Mo/Jω3, где Mo представляет собой статический дисбаланс, J представляет собой инерцию ротора вокруг оси поворота z, и ω представляет собой скорость поворота, Φ является неизвестным и представляет собой фазу дисбаланса, и
Figure 00000002
представляет собой среднее T^ideal - T^measured.
С использованием тригонометрического тождества Isin(Θ+Φ)=IcosΦsinΘ+IsinΦcosΘ, и подстановки γ=IcosΦ, σ=IsinΦ, S=sinΘ, C=cosΘ и Y=T^predicted-T^measured-
Figure 00000002
, уравнение 1 может быть переписано как показано в уравнении 2:
Уравнение 2:
Figure 00000003
где γ и σ являются неизвестными. Умножение уравнения 2 на транспонированную S (то есть, St) и транспонированную С (то есть, Ct), соответственно, дает уравнения 3 и 4:
Уравнение 3:
Figure 00000004
и
Уравнение 4:
Figure 00000005
или, два уравнения с двумя неизвестными, что можно решить как систему уравнений. В матричной форме уравнения 3 и 4 могут быть представлены, как показано в уравнениях 5 и 6:
Уравнение 5:
Figure 00000006
и
Уравнение 6:
Figure 00000007
которые можно решить относительно γ и σ.
С полученными значениями γ и σ неизвестные I и Φ из уравнения 1 могут быть определены, соответственно, из уравнений 7 и 8:
Уравнение 7:
Figure 00000008
и
Уравнение 8:
Figure 00000009
Выше, I представляет собой инерциальную амплитуду Mo/Jω3, где Mo представляет собой статический дисбаланс, J представляет собой инерцию ротора вокруг оси поворота z, и ω представляет собой скорость поворота, и может быть использована для определения ошибки IP, как показано в уравнении 9:
Уравнение 9:
Figure 00000010
На фигуре 8 показана кривая 802 ошибок, представляющая ошибку IP как функцию углового положения, при этом y-ось 804 представляет величину ошибки (или отклонение от фактического углового положения), и x-ось 806 представляет угловое положение. Точки 808, 810, 812 и 814 соответствуют угловым положениям θ (или 0, α, β и α+β). Кривая 802 ошибок подгоняется по точкам 808, 810, 812 и 814. В точке 808, исходном угловом положении IP для оборота, ошибка углового положения IP равна нулю (0). Ошибка в заданном угловом диапазоне может быть определена по отклонению в ошибке относительно значения в точке 808.
Определитель 204 ошибки временной привязки IP выдает кривую ошибок временной привязки IP корректору 206 временной привязки IP, который использует кривую ошибок временной привязки IP для корректировки предсказанного сигнала временной привязки IP, предсказанного предсказателем 202 временной привязки IP. В качестве примера, в одном неограничивающем случае, корректор 206 временной привязки IP прибавляет кривую ошибок временной привязки IP к сигналу временной привязки IP, в результате чего происходит сдвиг временной привязки IP таким образом, что IP синхронизируются с угловым положением поворотного гентри.
Следует понимать, что использование дополнительных датчиков и/или меток может улучшить оценку ошибки дисбаланса. Например, использование дополнительных датчиков и/или меток может снизить шум. В зависимости от количества датчиков и/или меток, контроллер IP 118 может функционировать, по существу, аналогично кодировщику с высоким разрешением.
Фигура 9 иллюстрирует примерный способ для определения сигнала временной привязки IP.
Следует понимать, что порядок действий в способах, описанных в настоящем описании, не является ограничивающим. Таким образом, в настоящем описании предусматриваются и другие порядки. Кроме того, одно или более действий может быть пропущено и/или может быть включено одно или более дополнительных действий.
На этапе 902 определяется временная продолжительность оборота поворотного гентри 104. Как описано в настоящем описании, этого можно достичь посредством датчика, метки и таймера.
На этапе 904 получают количество IP для одного или более последовательных оборотов.
На этапе 906 получают скорость поворота поворотного гентри для одного или более последовательных оборотов.
На этапе 908 сигнал временной привязки IP для одного или более последовательных оборотов определяется путем деления времени оборота на произведение количества IP и скорости поворота поворотного гентри.
На этапе 910 сигнал временной привязки IP используется во время получения данных для определения IP во время по меньшей мере одного или более последовательных оборотов.
Фигура 10 иллюстрирует примерный способ для определения скорректированного сигнала временной привязки IP.
Следует понимать, что порядок действий в способах, описанных в настоящем описании, не является ограничивающим. Таким образом, в настоящем описании предусматриваются и другие порядки. Кроме того, одно или более действий может быть пропущено и/или может быть включено одно или более дополнительных действий.
На этапе 1002 определяется первая временная продолжительность оборота поворотного гентри 104, в соответствии с описанным в настоящем описании.
На этапе 1004 определяется по меньшей мере вторая временная продолжительность от начального углового положения оборота до по меньшей мере первого опорного углового положения, которое находится на расстоянии, меньшем полного оборота (2π) от начального углового положения. Как описано в настоящем описании, этого можно достичь через использование датчика, метки и таймера.
На этапе 1006 предсказывается временная привязка IP для опорного углового положения.
На этапе 1008 измеряется временная привязка IP для опорного углового положения.
На этапе 1010 определяется ошибка временной привязки IP на основании различия между предсказанным и измеренным угловыми положениями.
На этапе 1012 генерируется сигнал ошибки IP путем подгонки значений ошибок в начальном угловом положении и по меньшей мере первом опорном угловом положении, кривой, охватывающей оборот.
На этапе 1014 сигнал временной привязки IP для одного или более последовательных оборотов определяется на основании первой временной продолжительности, количества IP и скорости поворота для последующего оборота.
На этапе 1016 сигнал временной привязки IP корректируется на основании сигнала ошибки IP. Как описано в настоящем описании, этого можно достичь путем суммирования сигналов.
На этапе 1018 скорректированный сигнал временной привязки IP используется во время получения данных для определения IP во время по меньшей мере одного или более последовательных оборотов.
Описанное выше может быть реализовано посредством машиночитаемых инструкций, которые при выполнении процессором(-ами) компьютера вызывают исполнение процессором(-ами) описанных действий. В этом случае инструкции хранятся на машиночитаемом носителе информации, ассоциированном с соответствующим компьютером или доступным с него.
Изобретение в настоящей заявке было описано со ссылкой на различные варианты осуществления. Модификации и изменения могут быть понятны другим после прочтения представленного описания. Предполагается, что изобретение ограничивается как включающее в себя все такие модификации и изменения в той степени, в которой они находятся в пределах объема прилагаемых пунктов формулы изобретения и их эквивалентов.

Claims (15)

1. Система визуализации (100), содержащая:
неподвижный гентри (102);
поворотный гентри (104), поддерживаемый с возможностью вращения неподвижным гентри и выполненный с возможностью вращения вокруг области исследования;
источник (110) излучения, удерживаемый поворотным гентри, выполненный с возможностью испускания излучения, которое проходит через область исследования;
матрицу (112) детекторов, удерживаемую поворотным гентри и расположенную напротив источника излучения через область исследования, которая обнаруживает излучение, которое проходит через область исследования,
при этом матрица детекторов обнаруживает излучение в течение множества периодов интегрирования во время вращения поворотного гентри, при этом множество периодов интегрирования соответствует различным диапазонам угловых положений, и при этом матрица детекторов генерирует сигнал, указывающий на обнаруженное излучение, соответственно, для упомянутого множества периодов интегрирования; и
контроллер (118) периодов интегрирования, который генерирует сигнал временной привязки периода интегрирования, который включает в себя временную привязку для начала каждого периода интегрирования для оборота поворотного гентри на основании по меньшей мере временной продолжительности предыдущего оборота поворотного гентри вокруг области исследования, при этом сигнал временной привязки интегрирования используется для инициации множества периодов интегрирования.
2. Система визуализации по п. 1, в которой контроллер периодов интегрирования определяет сигнал временной привязки периода интегрирования путем деления временной продолжительности предыдущего оборота на произведение числа периодов интегрирования и скорости поворота поворотного гентри для последующего оборота.
3. Система визуализации по п. 2, дополнительно содержащая:
по меньшей мере одну метку (404), прикрепленную к неподвижному гентри;
по меньшей мере один датчик (402), прикрепленный к поворотному гентри, при этом датчик обнаруживает метку в ответ на поворачивание посредством метки и генерирует выходной сигнал, указывающий это; и
таймер, который определяет временной интервал между первым сигналом выхода датчика, представляющим первый момент времени, когда датчик поворачивается, посредством метки, и вторым сигналом выхода датчика, представляющим следующий момент времени, когда датчик поворачивается, посредством метки, при этом временной интервал представляет собой временной интервал предыдущего оборота.
4. Система визуализации по п. 1, в которой контроллер периодов интегрирования генерирует сигнал ошибки периода интегрирования и корректирует сигнал временной привязки периода интегрирования на основании сигнала ошибки интегрирования.
5. Система визуализации по п. 4, в которой контроллер периодов интегрирования корректирует сигнал временной привязки периода интегрирования посредством суммирования сигнала временной привязки периода интегрирования и сигнала ошибки интегрирования.
6. Система визуализации по любому из пп. 4 или 5, в которой сигнал ошибки периода интегрирования указывает на колебания скорости поворотного гентри во время вращения.
7. Система визуализации по любому из пп. 4 или 5, в которой сигнал ошибки периода интегрирования указывает на дисбаланс поворотного гентри.
8. Система визуализации по любому из пп. 4 или 5, дополнительно содержащая:
определитель (204) ошибки временной привязки периода интегрирования, который определяет ошибки временной привязки для двух или более угловых положений поворотного гентри на основании различия между предсказанной временной привязкой периода интегрирования для двух или более угловых положений поворотного гентри и измеренной временной привязкой периода интегрирования для двух или более угловых положений поворотного гентри.
9. Система визуализации по п. 6, дополнительно содержащая:
две или более меток (606, 608), прикрепленных к неподвижному гентри и смещенные по углу друг от друга на первый угол;
два или более датчиков (602, 604), прикрепленных к поворотному гентри и смещенные по углу друг от друга на первый угол, при этом два или более датчиков, соответственно, обнаруживают метку в ответ на прохождение метки при повороте и генерируют выходной сигнал, указывающий это; и
таймер, который определяет временные интервалы между двумя метками и/или двумя датчиками, при этом одно или более угловых положений определяются путем деления соответствующих временных интервалов на временную продолжительность оборота.
10. Система визуализации по п. 8, в которой определитель ошибки временной привязки периода интегрирования определяет сигнал ошибки периода интегрирования путем подгонки ошибок временной привязки для двух или более угловых положений поворотного гентри к кривой, охватывающей предварительно определенный угловой диапазон оборота.
11. Способ визуализации, содержащий генерацию сигнала временной привязки периода интегрирования, причем способ содержит:
генерацию сигнала временной привязки периода интегрирования, включающего в себя временную привязку множества периодов интегрирования получения данных для оборота поворотного гентри системы (100) визуализации, при этом сигнал временной привязки периода интегрирования генерируется на основании по меньшей мере временной продолжительности предыдущего оборота поворотного гентри.
12. Способ по п. 11, дополнительно содержащий:
генерацию сигнала временной привязки периода интегрирования посредством деления временной продолжительности на количество периодов интегрирования и скорость поворота поворотного гентри.
13. Способ по любому из пп. 11 или 12, дополнительно содержащий:
корректировку сигнала временной привязки периода интегрирования с использованием сигнала ошибки положения интегрирования.
14. Способ по п. 13, дополнительно содержащий:
корректировку сигнала временной привязки периода интегрирования посредством суммирования сигнала временной привязки периода интегрирования и сигнала ошибки положения интегрирования.
15. Способ по п. 13, дополнительно содержащий:
предсказание временной привязки периода интегрирования для множества угловых местоположений;
измерение временной привязки периода интегрирования для множества угловых местоположений; и
определение ошибок определения положения периода интегрирования для множества угловых местоположений; и
генерацию сигнала ошибки положения интегрирования на основании определенных ошибок определения положения периодов интегрирования для множества угловых местоположений.
RU2013137507/28A 2011-01-10 2012-01-06 Получение данных компьютерной томографии (ст) RU2575943C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161431081P 2011-01-10 2011-01-10
US61/431,081 2011-01-10
PCT/IB2012/050079 WO2012095773A2 (en) 2011-01-10 2012-01-06 Computed tomography (ct) data acquisition

Publications (2)

Publication Number Publication Date
RU2013137507A RU2013137507A (ru) 2015-02-20
RU2575943C2 true RU2575943C2 (ru) 2016-02-27

Family

ID=45529158

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2013137507/28A RU2575943C2 (ru) 2011-01-10 2012-01-06 Получение данных компьютерной томографии (ст)

Country Status (6)

Country Link
US (1) US9173627B2 (ru)
EP (1) EP2663236B1 (ru)
JP (1) JP6071902B2 (ru)
CN (1) CN103501700B (ru)
RU (1) RU2575943C2 (ru)
WO (1) WO2012095773A2 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2715835C1 (ru) * 2019-06-19 2020-03-03 федеральное государственное бюджетное образовательное учреждение высшего образования "Донской государственный технический университет" (ДГТУ) Способ снижения структурной погрешности традиционного цифрового датчика физической величины в аналого-цифровой системе автоматического управления или контроля

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2899648A1 (en) * 2012-09-20 2015-07-29 Kabushiki Kaisha Toshiba Data processing device, data management system, data processing method, and program
CN104042231A (zh) * 2013-03-15 2014-09-17 Ge医疗系统环球技术有限公司 用于ct设备的机架的定位系统以及包括其的ct设备
CN103767725B (zh) 2013-12-31 2016-06-29 沈阳东软医疗系统有限公司 一种用于平衡ct机架的方法和装置
JP6564330B2 (ja) * 2016-01-22 2019-08-21 株式会社日立製作所 フォトンカウンティングct装置
DE102016224940B3 (de) * 2016-12-14 2017-08-24 Bruker Axs Gmbh Röntgenvorrichtung mit motorischer Drehmomentkompensation am Detektorkreis des Goniometers und Verfahren zum Betrieb einer Röntgenvorrichtung
WO2018182616A1 (en) * 2017-03-30 2018-10-04 Analogic Corporation System and method for sampling frequency adjustment for radiation imaging system
CN113995431B (zh) * 2021-10-27 2022-12-02 赛诺威盛科技(北京)股份有限公司 Ct扫描系统、方法、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5175754A (en) * 1992-05-11 1992-12-29 General Electric Company Gantry position reference for tomographic scanners
US6169778B1 (en) * 1998-03-24 2001-01-02 Siemens Aktiengesellschaft Computed tomography device
US6671345B2 (en) * 2000-11-14 2003-12-30 Koninklijke Philips Electronics N.V. Data acquisition for computed tomography
RU2384866C2 (ru) * 2005-04-22 2010-03-20 Конинклейке Филипс Электроникс Н.В. Сканер позитронно-эмиссионной томографии и магнитно-резонансной визуализации со способностью определения времени полета

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR9408119A (pt) * 1993-12-06 1997-08-05 Analogic Corp Aparelho para medir um ou mais parâmetros cinemáticos geométricos e posicionais de um dispositivo de rotação aparelho para medir qualquer deslocamento lateral de um centro geométrico de um dispositivo de rotação aparelho para medir substancialmente a posição angular instantânea de um dispositivo de rotação para um ponto de referência fixado do dispositivo aparelho para medir substancialmente a velocidade angular instantânea de um dispositivo de rotação à medida que um dispositivo de rotação gira em torno de um eixo de rotação aparelho incluindo um objeto de rotação montado nos dispositivos de suporte para girar em torno de um eixo rotacional sistema de tomografia de raio-x método de determinação de parâmetros associados com o movimento de um disco circular montado nos dispositivos de suporte para rotação em um plano de rotação em torno de um eixo rotacional
US6574301B1 (en) 2002-02-11 2003-06-03 Ge Medical Systems Global Technology Company, Llc CT data acquisition system trigger jitter filter
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
US7813473B2 (en) * 2002-07-23 2010-10-12 General Electric Company Method and apparatus for generating temporally interpolated projections
JP2006525060A (ja) 2003-05-06 2006-11-09 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 回転エンコーダによるcアームシステムにおいて画像収集を開始する方法及び装置
US7327822B2 (en) * 2005-07-20 2008-02-05 Purdue Research Foundation Methods, apparatus, and software for reconstructing an image
US7826587B1 (en) * 2009-09-11 2010-11-02 General Electric Company System and method of fast kVp switching for dual energy CT

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5175754A (en) * 1992-05-11 1992-12-29 General Electric Company Gantry position reference for tomographic scanners
US6169778B1 (en) * 1998-03-24 2001-01-02 Siemens Aktiengesellschaft Computed tomography device
US6671345B2 (en) * 2000-11-14 2003-12-30 Koninklijke Philips Electronics N.V. Data acquisition for computed tomography
RU2384866C2 (ru) * 2005-04-22 2010-03-20 Конинклейке Филипс Электроникс Н.В. Сканер позитронно-эмиссионной томографии и магнитно-резонансной визуализации со способностью определения времени полета

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2715835C1 (ru) * 2019-06-19 2020-03-03 федеральное государственное бюджетное образовательное учреждение высшего образования "Донской государственный технический университет" (ДГТУ) Способ снижения структурной погрешности традиционного цифрового датчика физической величины в аналого-цифровой системе автоматического управления или контроля

Also Published As

Publication number Publication date
US9173627B2 (en) 2015-11-03
WO2012095773A3 (en) 2012-10-26
CN103501700A (zh) 2014-01-08
CN103501700B (zh) 2016-05-11
WO2012095773A2 (en) 2012-07-19
US20130287165A1 (en) 2013-10-31
JP6071902B2 (ja) 2017-02-01
JP2014501596A (ja) 2014-01-23
RU2013137507A (ru) 2015-02-20
EP2663236B1 (en) 2019-04-24
EP2663236A2 (en) 2013-11-20

Similar Documents

Publication Publication Date Title
RU2575943C2 (ru) Получение данных компьютерной томографии (ст)
US7471760B2 (en) Device for obtaining structure data of a moving object
JP5572221B2 (ja) X線画像内のマーカ位置からの患者運動ベクトルの抽出
US9155516B2 (en) Apparatus and method for count loss calibration of photon-counting detectors in spectral computed tomography imaging
US7801264B2 (en) Method for calibrating a dual -spectral computed tomography (CT) system
US20030117136A1 (en) Method and apparatus for measuring and correcting motion effects using navigator echoes
WO2005099577A1 (en) Dynamic dose control for computed tomography
CN103284740A (zh) Ct机的扫描控制方法、扫描控制装置及ct机
JP2016198507A (ja) X線コンピュータ断層撮影装置及び医用画像処理装置
CN104068885A (zh) Ct机的扫描触发控制方法、装置及ct机
US20150036905A1 (en) Three-dimensional x-ray imaging
US8023713B2 (en) Method and system for reducing artifact due to time delay in data acquisition system in computer tomography
EP3513383A1 (en) Motion compensated reconstruction for helical computer tomography
JP4316017B2 (ja) X線ct装置
US12109052B2 (en) Radiographic imaging device and radiographic treatment device
EP4201336A1 (en) Method for operating a ct imaging system
US9345445B2 (en) Registration of 4th-generation detectors relative to 3rd-generation CT system coordinates using 4th-generation detector shadow pattern
JP2007216043A (ja) X線ct装置
CN102428500B (zh) 免内插扇形到平行射束重分装
US20080267344A1 (en) Signal Processing Unit for Producing Images
CN102855603B (zh) 一种自适应插值方法和装置
US9216005B2 (en) Synchronized view data acquisition between stationary spectral detectors and rotating energy integrating detectors for spectral computer-aided tomography
JP2007236777A (ja) X線ct装置及びx線ct装置の制御方法
JP3696155B2 (ja) 回転角度検出装置およびx線コンピュータ断層撮影装置
JP2014068717A (ja) X線ct装置

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20180107