EA016381B1 - Способ сведения и линеаризации амплитудных характеристик мультисенсорных детекторов - Google Patents
Способ сведения и линеаризации амплитудных характеристик мультисенсорных детекторов Download PDFInfo
- Publication number
- EA016381B1 EA016381B1 EA201100963A EA201100963A EA016381B1 EA 016381 B1 EA016381 B1 EA 016381B1 EA 201100963 A EA201100963 A EA 201100963A EA 201100963 A EA201100963 A EA 201100963A EA 016381 B1 EA016381 B1 EA 016381B1
- Authority
- EA
- Eurasian Patent Office
- Prior art keywords
- sensor
- sensors
- tabular
- detector
- response
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
- H04N25/60—Noise processing, e.g. detecting, correcting, reducing or removing noise
- H04N25/67—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response
- H04N25/671—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction
- H04N25/672—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction between adjacent sensors or output registers for reading a single image
Abstract
Изобретение относится к области рентгенотехники, а именно к цифровым детекторам рентгеновского изображения, и предназначено для сведения и линеаризации амплитудных характеристик отдельных сенсоров мультисенсорных детекторов. Технический результат - создание методики сведения и линеаризации амплитудных характеристик мультисенсорных детекторов в условиях неравномерной энергетической освещённости (не требует использования источника излучения с плоским рентгеновским (световым) полем). Способ основан на вычислении табличных функций преобразования интенсивности выходного сигнала сенсоров детектора. В результате применения данных табличных функций преобразования получают амплитудные характеристики сенсоров, которые в пределах точности измерений являются одинаковыми и линейными. Основная идея, используемая для вычисления табличных функций сведения, состоит в использовании того факта, что при наличии медленно изменяющегося по площади детектора неоднородного рентгеновского (светового) поля отклики любых двух смежных сенсоров, обладающих одинаковыми амплитудными характеристиками, должны возле стыка данных сенсоров иметь близкие значения.
Description
Область техники
Изобретение относится к области рентгенотехники, а именно к цифровым детекторам рентгеновского изображения, и предназначено для коррекции рентгеновского изображения. Более конкретно, данное изобретение предназначено для сведения и линеаризации амплитудных характеристик отдельных сенсоров мультисенсорных детекторов.
Предшествующий уровень техники
В настоящее время применение цифровых мультисенсорных детекторов в медицине, астрономии и других областях непрерывно расширяется |Но\те11 8.В. НапбЬоок οί ОСИ Айтопоту. СатЬпбде Ишует8Йу Рге55. 2000; Οίηο М. Ио18е, Иоще, Иоще; Мирошниченко С.И., Невгасимый А.А. Теория и техника многосенсорных цифровых рентгеновских приёмников. Биотехносфера, № 4/10/2010; Яне Б. Цифровая обработка изображений, М., Техносфера, 2007].
В мультисенсорных детекторах изображение может быть сформировано несколькими ПЗС (ССИ) или КМОП (СМО8) сенсорами. Каждый сенсор мультисенсорного детектора, в свою очередь, может состоять из более мелких сенсоров со своим собственным трактом преобразования сигнала. Амплитудные характеристики каждого сенсора такого детектора, определяющие выходные значения интенсивности сигнала формируемого изображения, в силу технологических и др. причин неизбежно будут различаться. Различия в амплитудных характеристиках приводят к общей неоднородности в отклике (неоднородности в выходном изображении) детектора, появляются разрывы непрерывности (т.н. швы) на стыках сенсоров в результирующем цифровом изображении. Поэтому на практике при использовании мультисенсорных детекторов возникает задача сведения амплитудных характеристик сенсоров, составляющих единый детектор. Для получения изображения, однородного по отклику и не содержащего швов на стыках сенсоров, необходимо применить процедуру калибровки, которая сводила бы все характеристики сенсоров к заданной. Кроме сведения для выполнения операции калибровки плоского поля должна быть решена задача линеаризации амлитудных характеристик.
Существуют несколько подходов к решению задачи сведения характеристик. В случае линейных откликов сенсоров стандартным подходом является т. н. двухточечная калибровка. Известны способы сведения и линеаризации их характеристик при существенно нелинейных характеристиках СМО8сенсоров. Так, в источнике [Ваблсоп 1тадшд Согр. ΑΝ08: Ро1упот1а1 Саш Соггесбоп Рот ВабЕуе Зепкога, (\у\у\у.габ-1соп.сот/рбГ/Ваб1соп ΑΝ08.ράί)] описан способ попиксельного сведения и линеаризации характеристик СМО8-сенсоров, при котором используя источник светового поля с равномерным распределением интенсивности по полю детектора, получают два калибровочных снимка на двух уровнях входного сигнала, причём первый уровень входного сигнала выбирают в два раза меньше следующего;
находят параметры квадратичной зависимости, моделирующей отклик сенсора, и, используя полученные калибровочные снимки, строят корректирующую функцию, линеаризующую и сводящую характеристики пикселей СМО8-сенсора.
В источнике |Τίμ С., 1отд Р. А Ртасйса1 №п-1теаг Оат Соттесйоп Ме11юб Рог Нщ11-ге5о1ибоп СМО8 Епащпд Эе1ес1ог5| описан другой способ сведения и линеаризации характеристик пикселей СМО8сенсоров, при котором используя световое поле с равномерным распределением интенсивности по полю детектора, получают три калибровочных снимка на трёх различных уровнях входного сигнала;
моделируют отклики пикселей сенсора кусочно-квадратичной гладкой зависимостью из трёх сегментов;
находят параметры модельной зависимости и, используя полученные калибровочные снимки, строят корректирующую функцию, линеаризующую и сводящую характеристики пикселей СМО8-сенсора.
Наиболее близким к заявляемому техническому решению является способ, описанный в источнике [Кобак. МиШр1е Ои1ри1 Зепкога 8еат§ Соггесйоп. Аррбсабоп №1е, 2009], в котором осуществляют посенсорную линеаризацию и сведение характеристик детектора.
В этом методе получают при возрастающих значениях интенсивности излучения с использованием источника светового поля с равномерным распределением интенсивности по полю детектора серию из Ν калибровочных снимков;
измеряют отклики сенсоров в виде табличных функций, описывающих зависимость выходного сигнала от значений входного, линеаризуют измеренные отклики сенсоров и сводят их к одному отклику, произвольно выбранному из этих линеаризованных откликов. Во всех перечисленных выше методах сведения и линеаризации амплитудных характеристик, включая метод [Кобак. МиШр1е Ои1ри1 Зепюге 8еат§ Соггесбоп. Аррбсабоп №1е, 2009], используют источник светового поля с равномерным распределением интенсивности по полю детектора. Однако в ряде случаев применение такого источника может быть неудобным либо невозможным в принципе. Неудобство связано с трудностью его создания. Невозможность использования источника светового поля с равномерным распределением интенсивности по полю детектора может быть обусловлена, например, наличием у разработчика детектора с уже интегрированным сцинтилляционным
- 1 016381 экраном, преобразующим рентгеновское излучение в свет. В последнем случае для сведения и линеаризации амплитудных характеристик сенсоров мультисенсорного детектора оказывается доступным только неплоское рентгеновское поле, приводящее к неравномерному распределению интенсивности по полю детектора (неравномерной энергетической (световой) освещённости).
Задачей настоящего изобретения является разработка способа сведения и линеаризации амплитудных характеристик детектора в условиях неравномерной энергетической освещённости.
Раскрытие изобретения
Техническим результатом заявляемого изобретения является создание способа сведения и линеаризации амплитудных характеристик мультисенсорных детекторов в условиях неравномерной энергетической освещённости (не требует использования источника излучения с равномерным распределением интенсивности по полю детектора).
Технический результат в способе сведения и линеаризации амплитудных характеристик мультисенсорных детекторов, заключающемся в том, что получают при возрастающих значениях интенсивности излучения серию из N калибровочных снимков с равномерным распределением интенсивности излучения по полю детектора, измеряют отклики сенсоров в виде табличных функций, описывающих зависимость выходного сигнала от значений входного, линеаризуют измеренные отклики сенсоров и сводят их к одному отклику, произвольно выбранному из этих линеаризованных откликов, достигается тем, что накапливают вдоль стыков смежных сенсоров М значений их откликов на полученных калибровочных снимках детектора методом вычисления скользящего среднего, получая для каждой пары смежных сенсоров ансамбль из ΜχΝ табличных функций преобразования, осуществляющих прямой перевод граничных значений сигнала смежных сенсоров; полученные табличные функции преобразования интерполируют на весь динамический диапазон интенсивностей выходного сигнала детектора, получая ансамбль из М табличных функций сведения смежных сенсоров; который затем усредняют, получая для каждого сенсора табличную функцию сведения его отклика к отклику смежного сенсора; выбирают отклик любого из сенсоров в качестве опорного и последовательно модифицируют табличные функции сведения смежных сенсоров, получая табличные функции сведения откликов сенсоров к отклику опорного сенсора, при этом каждую табличную функцию сведения для текущего сенсора строят с учётом табличной функции сведения предыдущего сенсора; вычисляют табличную функцию линеаризации, осуществляющую линеаризацию отклика опорного сенсора, и последовательно модифицируют табличные функции сведения с помощью найденной функции линеаризации, получая итоговые табличные функции сведения откликов сенсоров детектора к линеаризованному отклику выбранного опорного сенсора.
Возможен вариант осуществления изобретения, при котором вместо накопления М значений откликов сенсоров методом вычисления скользящего среднего вдоль стыков сенсоров применяют скользящую локальную линейную аппроксимацию сигнала плоскостью с экстраполяцией его М значений за пределы сенсора, что позволяет уменьшить влияние градиента входного рентгеновского (светового) поля на стыках сменных сенсоров.
Также возможен вариант, согласно которому осуществляют сведение и линеаризацию амплитудных характеристик сенсоров детектора, засвеченного входным световым полем с неравномерной интенсивностью.
Возможно разбиение по меньшей мере одного сенсора на части и выполнение вначале внутрисенсорного сведения с последующим междусенсорным сведением характеристик.
Главная отличительная особенность настоящего изобретения состоит в том, что сведение и линеаризация характеристик мультисенсорных детекторов не требует использования источника излучения с равномерным распределением интенсивности по полю детектора.
Подробное описание изобретения
Заявляемое техническое решение, возможность его технической реализации и достижение указанного технического результата поясняется фиг. 1-6.
На фиг. 1 показано устройство для получения рентгенограмм;
на фиг. 2 изображена общая схема оценки откликов сенсоров, используемая в настоящем изобретении для их сведения;
на фиг. 3 показан график, схематически изображающий сведение и линеаризацию характеристик сенсоров детектора, состоящего всего из двух смежных сенсоров;
на фиг. 4 приведены графики табличных функций преобразования сведения и линеаризации амплитудных характеристик;
на фиг. 5 показан исходный рентгеновский снимок размером 5000x4000 пикселей, полученный с использованием детектора, состоящего из 24 сенсоров;
на фиг. 6 продемонстрирована эффективность предлагаемой в данном изобретении методики сведения и линеаризации амплитудных характеристик для детектора, состоящего из 24 сенсоров.
Получение рентгенограмм осуществляют, например, с помощью устройства, показанного на фиг. 1. Оно содержит рентгеновскую трубку 1, которая испускает пучок рентгеновского излучения 2. Пучок рентгеновского излучения 2 поступает на детектор 3. Детектор 3 содержит сцинтилляционный экран (на
- 2 016381 чертеже не показан) и матрицу фоточувствительных элементов (на чертеже не показаны). Сцинтилляционный экран оптически связан с поверхностью активной области матрицы фоточувствительных элементов. Матрица фоточувствительных элементов (на чертеже не показана) состоит по меньшей мере из одного сенсора. Она может состоять из любого конечного числа сенсоров, например из двух.
Пучок рентгеновского излучения 2 попадает на детектор 3, сцинтилляционный экран преобразует его в видимый свет, который преобразуется сенсорами детектора в цифровую форму. Данная форма представляет собой цифровое изображение, состоящее из нескольких частей, количество которых соответствует числу сенсоров в детекторе. На фиг. 2 схематически показано цифровое изображение 4, сформированное детектором 3, состоящим из двух смежных сенсоров 5 и 6.
Согласно заявляемому способу сначала делают серию из N калибровочных снимков. Калибровочные снимки выполняют с равномерным распределением интенсивности излучения по полю детектора, без поглощающих объектов. Снимки генерируют при возрастающих значениях интенсивности излучения, с произвольным шагом по интенсивности от нулевой экспозиции (снимок считывания) до наступления сигнала насыщения с сенсоров детектора. В силу того что поле излучения в общем случае является неплоским, сенсоры детектора 3 уходят в насыщение неравномерно. Для уменьшения влияния шума в последующем этапе оценки откликов сенсоров для каждого значения интенсивности получают несколько снимков, которые затем усредняют. Одновременно с получением серии снимков замеряют с помощью дозиметра (на чертеже не показан) значения интенсивности излучения в произвольно выбранном месте детектора 3, располагая его, например, сверху над детектором 3 или рядом с ним.
Затем накапливают М значений откликов сенсоров методом вычисления скользящего среднего с помощью блока оценки сигнала 7 небольшого размера, например 30x30 пикселей, вдоль стыков смежных сенсоров. Блок оценки сигнала 7 в положении 8 соответствует отклику сенсора 5 под номером 1, положение 9 (фиг. 2) блока оценки сигнала 7 соответствует отклику сенсора 5 под номером М. Аналогично для сенсора 6 положение 10 блока оценки сигнала 7 соответствует первому отклику сенсора 6, а положение 11 блока оценки сигнала 7 соответствует отклику сенсора 6 под номером М. Для оценки сигнала применяют выборочное среднее или с целью повышения устойчивости к шуму медиану. На фиг. 3 приведены графики функций, где 12 - отклик сенсора 5 в одном из М положений блока оценки сигнала 7 на стыке смежных сенсоров 5 и 6, 13 - отклики сенсора 6 в одном из М положений блока оценки сигнала 7, 14 - желаемый общий линейный отклик соседних сенсоров 5 и 6. По горизонтальной оси графика (фиг. 3) показаны измеренные значения интенсивности (в относительных единицах, например дозах (О), нормированных на максимальное значение интенсивности), по вертикальной оси - вычисленное значение интенсивности (I) сигнала в выбранной области оценки сигнала.
Главная идея сбора значений откликов сенсоров вдоль их границ заключается в использовании того факта, что вблизи границы сенсоров интенсивности излучения практически одинаковы при плавно изменяющемся по площади детектора 3 входном рентгеновском (световом) поле, поэтому значения откликов откорректированных сенсоров 5 и 6 на стыках имеют близкие значения. Для реализации данной идеи на основе полученных М значений откликов сигналов сенсоров для N калибровочных снимков получают ансамбль из ΜχΝ табличных функций преобразования, осуществляющих прямой перевод граничных значений сигнала одного из смежных сенсоров в граничные значения сигнала другого сенсора. Например, каждой выходной интенсивности сигнала сенсора 6, вычисленной на его границе в калибровочном изображении, ставят в соответствие значение граничной интенсивности сенсора 5.
Данные табличные функции преобразования интерполируют на весь динамический диапазон интенсивностей выходного сигнала детектора, получая ансамбль из М табличных функций сведения.
Ансамбль из М табличных функций сведения усредняют с целью дальнейшего уменьшения ошибки в оценке сигнала, вызванной наличием шума на калибровочных снимках.
Таким образом получают для каждого сенсора табличную функцию сведения его отклика к отклику смежного сенсора. Для этого выбирают отклик любого из сенсоров детектора (например, отклик сенсора 5 на фиг. 2), который далее будем называть опорным, и последовательно, от опорного к остальным сенсорам, модифицируют табличные функции сведения так, чтобы осуществить сведение амплитудных характеристик сенсоров к амплитудной характеристике выбранного опорного сенсора 5.
В результате получают по калибровочным снимкам такие табличные функции сведения, которые обеспечивают сведение амплитудных характеристик сенсоров мультисенсорного детектора.
Помимо сведения амплитудных характеристик сенсоров 5 и 6 выполняют операцию калибровки плоского поля, амплитудную характеристику опорного сенсора 5 линеаризуют и преобразуют амплитудную характеристику сенсора 6 к линеаризованной амплитудной характеристике опорного сенсора 5. А именно, вычисляют табличную функцию линеаризации, осуществляющую линеаризацию отклика опорного сенсора. Затем последовательно модифицируют табличные функции сведения с помощью найденной функции линеаризации, получая итоговые табличные функции сведения откликов сенсора 6 детектора к линеаризованному отклику выбранного опорного сенсора 5.
На фиг. 4 показан график табличных функций сведения 15 отклика сенсора 5 в линейный отклик 14 фиг. 3, 16 - график табличной функции сведения отклика сенсора 6 в линейный отклик 14 фиг. 3, по го
- 3 016381 ризонтальной оси (фиг. 4) отложена входная интенсивность сигнала (I), по вертикальной изображена выходная интенсивность Т(1).
Описанная схема сведения и линеаризации амплитудных характеристик предполагает, что отклик каждого отдельного сенсора определяется, главным образом, одной характеристикой. Возможны различные комбинации сведения и линеаризации амплитудных характеристик мультисенсорных детекторов. Если характеристика сенсора обладает различной нелинейностью внутри отдельного сенсора, излагаемая здесь схема сведения и линеаризации допускает возможность разбиения сенсора на части и выполнение вначале внутрисенсорного сведения, с последующим междусенсорным сведением характеристик. С целью уменьшения влияния градиента, в особенности при использовании входного поля с большим градиентом, на стыках сенсоров вместо накопления М значений откликов сенсоров методом вычисления скользящего среднего вдоль стыков сенсоров на полученных калибровочных снимках применяют метод скользящей локальной линейной аппроксимации сигнала плоскостью с экстраполяцией его М значений за пределы сенсора.
Результат применения заявляемого в изобретении способа показан на фиг. 5-6. На фиг. 5 приведён снимок размером 5000x4000 пикселей, состоящий из 24 сенсоров, где 17 - изображение дозиметра, 18 и 19 - изображения, полученные с двух соседних сенсоров детектора 3. На фиг. 5 хорошо видны различия в отклике сенсоров, выражающиеся, в частности, в наличии стыков между изображениями 18 и 19, полученными со смежных сенсоров. На фиг. 6 приведён пример применения заявляемой в настоящем изобретении способа сведения и линеаризации амплитудных характеристик сенсоров: отклики сенсоров стали одинаковыми и линейными, границы между сенсорами незаметны.
Лучший вариант осуществления изобретения
Возможен следующий вариант реализации заявляемого способа.
Сначала получают серию из N калибровочных снимков с равномерным распределением интенсивности излучения по полю детектора (фиг. 2). Снимки генерируют при возрастающих значениях интенсивности излучения Н, с произвольным шагом по интенсивности излучения, от нулевой экспозиции (снимок считывания) до насыщения сигналов с сенсоров детектора. Интенсивность излучения измеряется дозиметром, расположенным в произвольном месте на поверхности детектора. Для уменьшения влияния шума для каждого значения интенсивности получают и усредняют несколько снимков. Затем с помощью блока оценки сигнала 7 накапливают вдоль стыков смежных сенсоров М значений их откликов на полученных калибровочных снимках детектора методом вычисления скользящего среднего. Рассмотрим для простоты ситуацию, при которой детектор состоит всего из двух сенсоров (фиг. 2). Оценку сигнала в данном случае осуществляют по следующим формулам:
= тейгаг{81 (х ± Δ, у ± Δ)) , = τηβάίαΐ^β2 (х ± Δ, .у ± Δ)), где - массив оценок сигнала первого сенсора - массив оценок сигнала второго сенсора , Δ радиус блока оценки сигнала, * = — количество оценок сигнала на стыке сенсоров, 7 = 1,Ν _ общее количество калибровочных снимков, (х, у) - координаты центра блока оценки сигнала. Вблизи стыка при использовании гладкого рентгеновского поля интенсивность излучения практически одинакова, поэтому при использовании сенсоров с одинаковыми амплитудными характеристиками должно выполняться приближённое равенство
(1)
В предлагаемом изобретении осуществляют построение табличных функций преобразования на основе использования приближённого равенства (1).
Для этого
1) для смежной пары сенсоров получают ансамбль из ΜχΝ табличных функций преобразования вида
2) интерполируют табличные функции преобразования Т на всём динамическом диапазоне I (на-
Т^/Шегро/а^Т^Т^к)
3) усредняют табличные функции сведения, т.е. получают одну табличную функцию сведения Т отклика сенсора 6 в сенсор 5
ϊ=1,М (2)
При усреднении в формуле (2) используют усреднение по ансамблю функций, заданных на одной и той же сетке значений зависимого аргумента (входной яркости е ^). Данная табличная функция сведения Т представляет собой искомую табличную функцию (ЬИТ, 1оок-ир 1аЫе), которая сводит различающиеся
- 4 016381 амплитудные характеристики сенсоров 5 и 6 к одной (к характеристике опорного сенсора 5) (фиг. 2).
Затем для выполнения процедуры линеаризации амплитудных характеристик сенсоров выбирают отклик любого из сенсоров, например, пусть это будет отклик К2^ сенсора 6, и применяют к нему найденную табличную функцию преобразования
1) рассматривая отклик как функцию измеренной интенсивности излучения (дозы), находят параметр наклона α желаемой линейной характеристики сенсоров детектора из следующей формулы [Кобак. МиШр1е ΟιιΙριιΙ 8епзог8 8еатз Соггесйоп. Аррйсайоп Νοίβ, 2009]:
N
Σι-/=1 к Л аР/
X
2) строят табличную функцию линеаризации, линеаризующую сенсор 6 с интерполяцией её на весь динамический диапазон детектора
Л = /п1е/ро/аё(Т11п)= ке/
3) итоговые табличные функции сведения откликов сенсоров 5 и 6 находят соответственно формулам
(сенсор 5 не изменялся).
Изложенная здесь схема последовательного сведения и линеаризации амплитудных характеристик сенсоров легко распространяется на случай более двух сенсоров в детекторе.
Литература
1. 11о\ме11 8.В. НапбЬоок ок ССО ЛзЛопошу. СатЬпбде Итуегзйу Рге88, 2000.
2. Кобак. Ми1йр1е ΟιιψιιΙ 8епзог8 8еатз Соггесйоп. Лррйеайоп №1е, 2009.
(^^^.кобак.сот/д1оЬа1/р1и§1П8/асгоЬа1/еп/Ьи81пе88/188/8иррог1бос8/МиЙ1р1е Ои1ри18еп8ог88еат8Соггесбоп.рб1).
3. Каб-1соп 1тадт§ Согр. ΑΝ08: Ро1упот1а1 Оат Соггесйоп ког КабЕуе 8епзог8 (^^^.габ1соп.сот/рбк/Каб1соп АЖ8.рбк).
4. Είρ С., .1огд Р. А Ргасйса1 ^п-йпеаг Оат Соггесйоп МеФоб ког Нщй-гезоМоп СМО8 1тадт§ 1)е1ес1ог8. (Ийр8://^^^.коег1есй.йаи8бе8Йоегеп8-о1бепЬигд.бе/бдтр2008/аЬ81гас1/Сао.рбк&гс1=1&д=Е1]1 Ргасйса1 №п-1теаг Оат Соггесйоп МеФоб ког Н1§й-ге8о1ийоп рбк).
5. Ото М. №18е, №1зе, №1зе (Ьир://^^^.а81горЬу8-а88181.еот/е0иеа1е/по18е/по18е.Ь1т).
6. Мирошниченко С.И., Невгасимый А.А. Теория и техника многосенсорных цифровых рентгеновских приёмников. Биотехносфера, № 4/10/2010.
7. Яне Б. Цифровая обработка изображений. М., Техносфера, 2007, 583 с.
Claims (4)
- ФОРМУЛА ИЗОБРЕТЕНИЯ1. Способ сведения амплитудных характеристик мультисенсорных детекторов, заключающийся в том, что получают при возрастающих значениях интенсивности излучения серию из Ν калибровочных снимков с равномерным распределением интенсивности излучения по полю детектора, измеряют отклики сенсоров в виде табличных функций, описывающих зависимость выходного сигнала от значений входного, линеаризуют измеренные отклики сенсоров и сводят их к одному отклику, произвольно выбранному из этих линеаризованных откликов, отличающийся тем, что накапливают вдоль стыков смежных сенсоров М значений их откликов на полученных калибровочных снимках детектора методом вычисления скользящего среднего, получая для каждой смежной пары сенсоров ансамбль из МxN табличных функций преобразования, осуществляющих прямой перевод граничных значений сигнала смежных сенсоров; полученные табличные функции преобразования интерполируют на весь динамический диапазон интенсивностей выходного сигнала детектора, получая ансамбль из М табличных функций сведения смежных сенсоров, который затем усредняют, получая для каждого сенсора табличную функцию сведения его отклика к отклику смежного сенсора; выбирают отклик любого из сенсоров в качестве опорного и последовательно модифицируют табличные функции сведения смежных сенсоров, получая табличные функции сведения откликов сенсоров к отклику опорного сенсора, при этом каждую табличную функцию сведения для текущего сенсора строят с учётом табличной функции сведения предыдущего сенсора; вычисляют табличную функцию линеаризации, осуществляющую линеаризацию отклика опорного сенсора, и последовательно модифицируют табличные функции сведения с помощью найденной функции линеаризации, получая итоговые табличные функции сведения откликов сенсоров детектора к линеаризованному отклику выбранного опорного сенсора.- 5 016381
- 2. Способ по п.1, отличающийся тем, что после получения серии калибровочных снимков М значений сигнала оценивают скользящей локальной линейной аппроксимацией сигнала плоскостью с экстраполяцией его М значений за пределы сенсора.
- 3. Способ по п.1, отличающийся тем, что осуществляют сведение и линеаризацию амплитудных характеристик сенсоров детектора, засвеченного входным световым полем с неравномерной интенсивностью.
- 4. Способ по любому из пп.1-3, отличающийся тем, что по меньшей мере один сенсор разбивают на части и выполняют внутрисенсорное сведение его характеристик с последующим междусенсорным сведением характеристик.
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EA201100963A EA016381B1 (ru) | 2011-07-18 | 2011-07-18 | Способ сведения и линеаризации амплитудных характеристик мультисенсорных детекторов |
| CN2012100890288A CN102891952A (zh) | 2011-07-18 | 2012-03-29 | 对多传感器检测器的增益特性进行拼接的方法 |
| JP2012091963A JP2013024861A (ja) | 2011-07-18 | 2012-04-13 | マルチセンサ検出器ゲイン特性のステッチング方法 |
| KR1020120041360A KR20130010429A (ko) | 2011-07-18 | 2012-04-20 | 멀티센서 검출기들의 이득 특성의 바느질 방식의 결합 및 선형화 방법 |
| EP12169306A EP2549744A2 (en) | 2011-07-18 | 2012-05-24 | Method of stitching and linearization of multisensor detectors gain characteristics |
| US13/550,451 US8606052B2 (en) | 2011-07-18 | 2012-07-16 | Method of stitching and linearization of multisensor detectors gain characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EA201100963A EA016381B1 (ru) | 2011-07-18 | 2011-07-18 | Способ сведения и линеаризации амплитудных характеристик мультисенсорных детекторов |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| EA201100963A1 EA201100963A1 (ru) | 2012-03-30 |
| EA016381B1 true EA016381B1 (ru) | 2012-04-30 |
Family
ID=45908256
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| EA201100963A EA016381B1 (ru) | 2011-07-18 | 2011-07-18 | Способ сведения и линеаризации амплитудных характеристик мультисенсорных детекторов |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US8606052B2 (ru) |
| EP (1) | EP2549744A2 (ru) |
| JP (1) | JP2013024861A (ru) |
| KR (1) | KR20130010429A (ru) |
| CN (1) | CN102891952A (ru) |
| EA (1) | EA016381B1 (ru) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6105903B2 (ja) * | 2012-11-09 | 2017-03-29 | キヤノン株式会社 | 画像処理装置、画像処理方法、放射線撮影システム及びプログラム |
| EP3025125B1 (de) | 2013-07-23 | 2018-03-21 | Balluff GmbH | Verfahren zur dynamischen linearisierung von sensorsignalen eines magnetband-längenmesssystems |
| JP6169626B2 (ja) * | 2014-03-10 | 2017-07-26 | 富士フイルム株式会社 | 放射線画像処理装置、方法およびプログラム |
| US9448134B2 (en) | 2014-06-23 | 2016-09-20 | Exxonmobil Upstream Research Company | Systems for detecting a chemical species and use thereof |
| US9471969B2 (en) | 2014-06-23 | 2016-10-18 | Exxonmobil Upstream Research Company | Methods for differential image quality enhancement for a multiple detector system, systems and use thereof |
| WO2015199914A1 (en) * | 2014-06-23 | 2015-12-30 | Exxonmobil Upstream Research Company | Methods for calibrating a multiple detector system |
| EP3158320B1 (en) | 2014-06-23 | 2018-07-25 | Exxonmobil Upstream Research Company | Methods and systems for detecting a chemical species |
| US10499863B2 (en) | 2014-11-17 | 2019-12-10 | Carestream Health, Inc. | Tiled digital radiography detectors for long-length imaging |
| US9820703B2 (en) | 2014-11-17 | 2017-11-21 | Carestream Health, Inc. | Tiled digital radiography detectors for long-length imaging |
| US10638986B2 (en) | 2014-11-17 | 2020-05-05 | Carestream Health, Inc. | Modular single shot digital radiography for long-length imaging |
| US9609197B1 (en) * | 2016-08-19 | 2017-03-28 | Intelligent Security Systems Corporation | Systems and methods for dewarping images |
| US9547883B1 (en) | 2016-08-19 | 2017-01-17 | Intelligent Security Systems Corporation | Systems and methods for dewarping images |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6122343A (en) * | 1995-04-07 | 2000-09-19 | Technological Resources Pty Limited | Method and an apparatus for analyzing a material |
| RU2218088C1 (ru) * | 2002-05-24 | 2003-12-10 | Щетинин Виктор Васильевич | Цифровой рентгенодиагностический аппарат |
| RU2406975C1 (ru) * | 2009-03-27 | 2010-12-20 | Открытое акционерное общество "Техприбор" | Мультисенсорный анализатор состава и скоростных параметров газожидкостного потока нефтяных скважин |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7433507B2 (en) * | 2003-07-03 | 2008-10-07 | Ge Medical Systems Global Technology Co. | Imaging chain for digital tomosynthesis on a flat panel detector |
| US8116595B2 (en) * | 2007-03-06 | 2012-02-14 | Infimed, Inc. | Universal interface for medical imaging receptors |
| US7923801B2 (en) * | 2007-04-18 | 2011-04-12 | Invisage Technologies, Inc. | Materials, systems and methods for optoelectronic devices |
-
2011
- 2011-07-18 EA EA201100963A patent/EA016381B1/ru active IP Right Revival
-
2012
- 2012-03-29 CN CN2012100890288A patent/CN102891952A/zh active Pending
- 2012-04-13 JP JP2012091963A patent/JP2013024861A/ja active Pending
- 2012-04-20 KR KR1020120041360A patent/KR20130010429A/ko not_active Application Discontinuation
- 2012-05-24 EP EP12169306A patent/EP2549744A2/en not_active Withdrawn
- 2012-07-16 US US13/550,451 patent/US8606052B2/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6122343A (en) * | 1995-04-07 | 2000-09-19 | Technological Resources Pty Limited | Method and an apparatus for analyzing a material |
| RU2218088C1 (ru) * | 2002-05-24 | 2003-12-10 | Щетинин Виктор Васильевич | Цифровой рентгенодиагностический аппарат |
| RU2406975C1 (ru) * | 2009-03-27 | 2010-12-20 | Открытое акционерное общество "Техприбор" | Мультисенсорный анализатор состава и скоростных параметров газожидкостного потока нефтяных скважин |
Non-Patent Citations (1)
| Title |
|---|
| KODAK CCD IMAGE SENSORS. MULTIPLE OUTPUT SENSORS SEAMS CORRECTION. APPLICARION NOTE. Revision 1.0 MTD/PS-1149, 28.04.2009 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EA201100963A1 (ru) | 2012-03-30 |
| US20130020476A1 (en) | 2013-01-24 |
| US8606052B2 (en) | 2013-12-10 |
| EP2549744A2 (en) | 2013-01-23 |
| KR20130010429A (ko) | 2013-01-28 |
| JP2013024861A (ja) | 2013-02-04 |
| CN102891952A (zh) | 2013-01-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EA016381B1 (ru) | Способ сведения и линеаризации амплитудных характеристик мультисенсорных детекторов | |
| EP2867701B1 (en) | Digital positron emission tomography (dpet) energy calibration method | |
| RU2414724C2 (ru) | Способ и устройство для спектральной компьютерной томографии | |
| US8262288B2 (en) | Focal spot position determiner | |
| WO2015190015A1 (ja) | 測距装置 | |
| US8611500B2 (en) | X-ray imaging apparatus and measurement method | |
| KR20180127987A (ko) | X-선 검출기를 위한 시스템 및 방법 | |
| CN105122085A (zh) | 利用调整的能量阈值用于生成能量分辨x射线图像的方法和设备 | |
| JP2017103608A (ja) | 放射線撮像装置及びその制御方法 | |
| JP2020514740A (ja) | 放射線検出器における使用のための画素設計 | |
| CN110072459A (zh) | 用于自校准的自校准ct检测器、系统和方法 | |
| KR20140087755A (ko) | 엑스선 영상 장치 및 그 제어 방법 | |
| JP2019520909A (ja) | 光子カウントコンピュータ断層撮影 | |
| US9453924B2 (en) | Radiation control apparatus and method for acquiring correction data for an imaging device | |
| JP7239612B2 (ja) | 種々の被写体及び放射能分布のための定量的陽電子放射トモグラフィ（ｐｅｔ）再構成における不感時間補正方法 | |
| JP6202855B2 (ja) | 画像処理装置、画像処理装置の制御方法、及びプログラム | |
| EP2702450B1 (en) | System and method for correction of vignetting effect in multi-camera flat panel x-ray detectors | |
| CN111163698A (zh) | 图像处理设备、图像处理方法以及程序 | |
| KR101330117B1 (ko) | 다채널 광전소자를 이용하는 양전자방출 단층촬영장치 | |
| Starck et al. | The use of detective quantum efficiency (DQE) in evaluating the performance of gamma camera systems | |
| JP7246936B2 (ja) | 放射線撮像装置、放射線撮像システムおよび放射線撮像装置の制御方法 | |
| Wang et al. | Novel light-guide-PMT geometries to reduce dead edges of a scintillation camera | |
| US11307312B2 (en) | Image acquisition | |
| HajiRassouliha et al. | The effect of camera settings on image noise and accuracy of subpixel image registration | |
| Friedman et al. | A moving slanted‐edge method to measure the temporal modulation transfer function of fluoroscopic systems |
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): MD |
|
| 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 BY KZ KG 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): RU |
|
| NF4A | Restoration of lapsed right to a eurasian patent |
Designated state(s): RU |