ГОСТ Р МЭК 62220-1-2-2010. Национальный стандарт Российской Федерации. Изделия медицинские электрические. Характеристики устройств для получения цифровых рентгеновских изображений. Часть 1-2. Определение квантовой эффективности регистрации. Детекторы, используемые при маммографии

(утв. и введен в действие Приказом Росстандарта от 12.11.2010 № 452-ст)

Утвержден и введен в действие Приказом Федерального агентства

ИЗДЕЛИЯ МЕДИЦИНСКИЕ ЭЛЕКТРИЧЕСКИЕ

ХАРАКТЕРИСТИКИ УСТРОЙСТВ ДЛЯ ПОЛУЧЕНИЯ ЦИФРОВЫХ РЕНТГЕНОВСКИХ ИЗОБРАЖЕНИЙ

ЧАСТЬ 1-2

ОПРЕДЕЛЕНИЕ КВАНТОВОЙ ЭФФЕКТИВНОСТИ РЕГИСТРАЦИИ ДЕТЕКТОРЫ, ИСПОЛЬЗУЕМЫЕ ПРИ МАММОГРАФИИ

Medical electrical equipment. Characteristics of digital X-ray imaging devices. Part 1-2. Determination of the detective quantum efficiency. Detectors used in mammography

IEC 62220-1-2:2007 Medical electrical equipment — Characteristics of digital X-ray imaging devices — Part 1-2: Determination of the detective quantum efficiency — Detectors used in mammography (IDT)

ГОСТ Р МЭК 62220-1-2-2010

Группа Ф32

ОКС 11.040.50

ОКП 94 4220

Предисловие

Сведения о стандарте

3. Утвержден и введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 12 ноября 2010 г. № 452-ст.

4. Настоящий стандарт идентичен международному стандарту МЭК 62220-1-2:2007 «Изделия медицинские электрические. Характеристики цифровых приемников рентгеновского изображения. Часть 1-2. Определение квантовой эффективности регистрации. Детекторы, используемые при маммографии» (МЭК 62220-1-2:2007 «Medical electrical equipment — Characteristics of digital X-ray imaging devices — Part 1-2: Determination of the detective quantum efficiency — Detectors used in mammography).»

При применении настоящего стандарта рекомендуется использовать вместо ссылочных международных стандартов соответствующие им национальные стандарты Российской Федерации, сведения о которых приведены в дополнительном Приложении ДА.

Информация об изменениях к настоящему стандарту публикуется в ежегодно издаваемом информационном указателе «Национальные стандарты», а текст изменений и поправок — в ежемесячно издаваемых информационных указателях «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячно издаваемом информационном указателе «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет.

Введение

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

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

Примечание. Несмотря на то, что DQE широко используется для описания устройств для получения изображения, связь между этим физическим параметром и восприятием изображения наблюдателем еще не вполне понятна [1], [3].

DQE уже широко используется изготовителями для описания устройств для получения цифровых рентгеновских изображений. DQE также требуется регулирующим организациям [таким как Администрация продуктов питания и напитков (FDA)] для процедуры допуска продукции к употреблению. Однако до настоящего времени нет стандарта ни для условий измерения, ни для процедуры измерения, что ведет к тому, что значения, полученные из разных источников, нельзя сравнивать.

Поэтому настоящий стандарт создан для регламентации процедуры измерений и формата декларации соответствия квантовой эффективности регистрации для устройств для получения цифровых рентгеновских изображений.

Настоящий стандарт будет полезен изготовителям, потребителям, дистрибьюторам и регулирующим организациям. Это второй документ из серии из трех стандартов:

• Часть 1, которая будет предназначена для рентгенографии, за исключением маммографии и рентгеноскопии;
• Часть 1-2 (настоящий стандарт), предназначенная для маммографии;
• Часть 1-3, которая будет предназначена для детекторов динамических изображений.

Эти стандарты могут рассматриваться как первая часть группы стандартов МЭК 62220, описывающих параметры устройств для получения цифрового рентгеновского изображения.

1. ОБЛАСТЬ ПРИМЕНЕНИЯ

Настоящий стандарт устанавливает метод определения квантовой эффективности детектора для устройств для получения цифрового рентгеновского изображения в функции воздушной кермы и пространственной частоты в диапазоне рабочих условий, установленных изготовителем. Настоящий стандарт МЭК 62220 рекомендован к применению изготовителями, а также в хорошо оборудованных испытательных лабораториях.

Настоящий стандарт распространяется на устройства для получения цифровых рентгеновских изображений, применяемые для получения маммографических изображений, таких как системы цифровой рентгенографии (CR), системы на основе прямых и непрямых плоских панелей, сканирующие системы (на основе CCD или счета фотонов). Настоящий стандарт не распространяется:

• на устройства для получения цифровых рентгеновских изображений, применяемые в общей рентгенографии и в детальной рентгенографии;
• компьютерную томографию;
• устройства для получения динамических изображений (где получают серию изображений, например, при рентгеноскопии или при исследовании сердца).

Примечание. Вышеперечисленные устройства исключены, поскольку многие их параметры (например, качество пучка излучения, геометрия, зависимость во времени и т.д.) отличаются от соответствующих параметров при маммографии. Некоторые из этих устройств будут рассмотрены в отдельных стандартах (МЭК 62220-1 и МЭК 62220-1-3) подобно тому, как это сделано по отношению к другим параметрам, например чувствительности и контрасту, в стандартах МЭК и ИСО.

2. НОРМАТИВНЫЕ ССЫЛКИ

В настоящем стандарте использованы нормативные ссылки на следующие стандарты. Для датированных ссылок применяют только указанные изделия. Для недатированных ссылок применяют самые последние изделия (включая любые дополнения и изменения).

МЭК 60336 Изделия медицинские электрические. Излучатели медицинские рентгенодиагностические. Характеристики фокусных пятен

МЭК 60601-2-45:2005 Изделия медицинские электрические. Часть 2-45. Частные требования безопасности к маммографическим рентгеновским аппаратам и маммографическим устройствам для стереотаксиса

МЭК 61267:2001 Медицинская диагностическая рентгеновская аппаратура. Радиационные условия при определении характеристик

МЭК 62220-1:2003 Изделия медицинские электрические. Характеристики устройств для получения цифрового рентгеновского изображения

ИСО 12232:1998 Фотография. Электронные камеры для регистрации статических изображений. Определение чувствительности по ИСО.

3. ТЕРМИНЫ И ОПРЕДЕЛЕНИЯ

В настоящем стандарте применены термины по МЭК 60788, а также следующие термины с соответствующими определениями:

3.1. Функция преобразования: участок уровня выходного сигнала с большой площади (исходные данные) устройства для получения цифрового рентгеновского изображения по отношению к числу экспонирующих квантов на единицу площади (Q) в плоскости поверхности детектора.

Примечания

1. При расчете Q следует умножать измеренное значение воздушной кермы за вычетом обратного рассеяния на значение, указанное в таблице 2, столбец 4.

2. Многие калибровочные лаборатории, такие как национальные институты метрологии, калибруют измерители излучения в единицах воздушной кермы.

3.2. Квантовая эффективность регистрации DQE(u, v): отношение двух функций NPS: в числителе NPS входного сигнала на поверхности детектора цифрового рентгеновского детектора после прохождения определяющего фильтра, указанного для функции передачи системы, а в знаменателе — измеренный NPS выходного сигнала (исходные данные).

Примечание. Вместо двумерной квантовой эффективности регистрации часто публикуют сечение двумерной квантовой эффективности регистрации вдоль указанной линии в плоскости частот.

3.3. Поверхность детектора: доступная поверхность, ближайшая к плоскости приемника изображения.

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

[МЭК 62220-1:2003, определение 3.4, модифицировано]

3.4. Устройство для получения цифрового рентгеновского изображения: устройство, состоящее из цифрового рентгеновского детектора, включая защитные слои, установленные для практического использования электроники усиления и оцифровки и компьютера, и обеспечивающего исходные данные (DN) изображения.

[МЭК 62220-1:2003, определение 3.5]

3.5. Матрица изображения: сборка матричных элементов, преимущественно в декартовой системе координат.

[МЭК 62220-1:2003, определение 3.6, модифицировано]

3.6 Лаг-эффект: влияние предыдущего изображения на текущее.

[МЭК 62220-1:2003, определение 3.7]

3.7. Линеаризованные данные: исходные данные, к которым приложена обратная функция преобразования.

Примечание. Линеаризованные данные прямо пропорциональны воздушной керме.

[МЭК 62220-1:2003, определение 3.8]

3.8. Функция передачи модуляции MTF(u, v): модуль обычно комплексной функции оптической передачи, выраженной как функция пространственных частот u и v.

[МЭК 62220-1:2003, определение 3.9]

3.9 Шум: флюктуации от ожидаемого значения стохастического процесса.

[МЭК 62220-1:2003, определение 3.10]

3.10. Спектр мощности шума NPS-W(u, v): модуль преобразования Фурье от автоковариантной функции шума. Мощность шума, содержащегося в двумерном интервале пространственной частоты, в функции двумерной частоты.

Примечание. В литературе спектр мощности шума часто называют «Спектр Винера» в честь математика Норберта Винера.

[МЭК 62220-1:2003, определение 3.11]

3.11. Исходные данные DN: необработанные данные, подвергшиеся коррекции, указанной в настоящем стандарте.

[МЭК 62220-1:2003, определение 3.12]

3.12 Поток фотонов Q: среднее число фотонов на единицу площади.

[МЭК 62220-1:2003, определение 3.13]

3.13. Необработанные данные: сигнал пикселя, считанный непосредственно после аналого-цифрового преобразования с устройства для получения цифрового рентгеновского изображения или полученный с системы счета фотонов без какой-либо программной коррекции.

[МЭК 62220-1:2003, определение 3.14, модифицировано]

3.14. Пространственная частота u или v: обратная величина периода повторения пространственного явления. Размерность пространственной частоты обратно длине.

[МЭК 62220-1:2003, определение 3.15]

4. ТРЕБОВАНИЯ

4.1 Условия работы

Устройство для получения цифрового рентгеновского изображения должно храниться и работать в соответствии с рекомендациями изготовителя. Условия работы должны быть такими же, как предусмотренные для работы в клинике, и поддерживаться в течение всего времени испытаний.

Условия окружающей среды в помещении, где работает устройство для получения цифрового рентгеновского изображения, должны быть указаны вместе с результатами испытаний.

4.2 Рентгеновский аппарат

Для всех испытаний по настоящему стандарту должно использоваться РПУ постоянного напряжения (МЭК 60601-2-45). Процентная пульсация должна быть равна или меньше 4%.

Номинальное значение фокусного пятна (МЭК 60336) должно быть не больше 0,4.

Для измерения воздушной кермы должен использоваться калиброванный измеритель излучения. Погрешность измерений (коэффициент покрытия 2) [2] должна быть меньше 5%.

Примечания

1. Термины «погрешность» и «коэффициент покрытия» определены в справочнике ИСО по выражению погрешности измерений [2].

4.3 Качество излучения

Качество излучения должно быть RQA-M 2 по МЭК 61267, если это подходит для клинического использования этого детектора. Может быть использовано другое качество излучения, которое используется при клиническом использовании устройства для получения цифрового рентгеновского изображения, например RQA-M 1, RQA-M 3, RQA-M 4, или качество излучения при материале анода, отличном от молибдена (см. таблицу 1).

Таблица 1

Качество излучения для определения квантовой эффективности регистрации и соответствующие параметры

Характеристика стандартного качества излучения (МЭК 61267)	Толщина фильтра, мм	Номинальное анодное напряжение, кВ	Номинальный первый слой половинного ослабления (СПО), мм AI	Дополнительный фильтр, мм AI
Mo/Mo (RQA-M 1)	0,032	25	0,56	2
Mo/Mo (RQA-M 2)	0,032	28	0,60	2
Mo/Mo (RQA-M 3)	0,032	30	0,62	2
Mo/Mo (RQA-M 3)	0,032	35	0,68	2
Mo/Rh	0,025	28	0,65	2
Rh/Rh	0,025	28	0,74	2
W/Rh	0,050	28	0,75	2
W/AI	0,500	28	0,83	2

4.4 Тест-объект

Минимальные размеры пластины: a = 120 мм, f = 60 мм.

Область интереса (ROI) при определении MTF b x c, 25 x 50 мм (указана внутренней штрихпунктирной линией).

После облучения на детекторе (наружная штрихпунктирная линия) не менее 100 x 100 мм.

Рисунок 1. Тест-объект

Пластина из нержавеющей стали используется как краевой тест-объект. Поэтому край пластины должен быть тщательно отполирован под углом 90° к плоскости пластины. При облучении края рентгеновскими лучами в контакте с безэкранной пленкой изображение края на пленке не должно быть размыто более чем на 5 мкм.

4.5 Геометрия

Геометрия установки для измерений должна соответствовать рисунку 2. В этой геометрической конфигурации рентгеновский аппарат используется таким же образом, как при нормальном диагностическом применении. Расстояние между фокусным пятном рентгеновской трубки и поверхностью детектора должно быть от 600 до 700 мм. Если по техническим причинам расстояние нельзя установить в этом диапазоне, может быть выбрано другое расстояние, но это должно быть указано при декларации результата.

Рисунок 2. Геометрия при экспонировании устройства для получения цифрового рентгеновского изображения с целью

или функции передачи модуляции за тест-объектом

Тест-объект должен быть размещен непосредственно перед поверхностью детектора. Центр края тест-объекта размещают на расстоянии 60 мм от центра поверхности корпуса детектора. Облучаемая площадь поверхности детектора должна быть 100 x 100 мм, центр этой площади должен находиться на расстоянии 60 мм от центра поверхности корпуса детектора.

В установке по рисунку 2 диафрагма B1 и дополнительный фильтр должны быть расположены близко от фокусного пятна рентгеновской трубки. Диафрагму B2 следует использовать, но она может не применяться в том случае, если доказано, что ее наличие не изменяет результатов измерений.

Детектор-монитор следует применять для того, чтобы подтвердить точность рентгеновского генератора. Детектор-монитор R1 должен быть размещен за той частью пучка, которая проходит через диафрагму B2. Точность детектора-монитора (стандартное отклонение 16) должна быть лучше 2%. Соотношение между показанием, считываемым с монитора, и воздушной кермой на поверхности детектора должно быть откалибровано для каждого используемого качества излучения. При калибровке этого соотношения должно быть обращено внимание на то, чтобы на показания измерителя излучения не влияло обратное рассеяние излучения от оборудования, находящегося сзади измерителя излучения. В любом случае необходимо проверить, что детектор-монитор не влияет на измерения функции преобразования, MTF и спектра мощности шума. Для минимизации обратного рассеяния от слоев, расположенных за детектором, следует обеспечить расстояние до других объектов не менее 250 мм.

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

Эта геометрия используется и для единообразного облучения поверхности детектора при определении функции преобразования и спектра мощности шума, и для облучения поверхности детектора за тест-объектом (см. 4.6.6). При всех измерениях должна облучаться одна и та же область поверхности детектора.

4.6 Условия облучения

4.6.1 Общие условия

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

Уровень экспозиции должен быть выбран таким же, как используемый при работе цифрового рентгеновского детектора в клинической практике, для которой он предназначен. Этот уровень называется «опорным», и он должен быть нормирован изготовителем. Должны быть выбраны, по крайней мере, два дополнительных уровня облучения, один равен двум «опорным» уровням и один равен 1/2 «опорного» уровня. При изменении уровня экспозиции не должно быть допущено изменение системы уставок (увеличение мощности и т.д.).

Для того чтобы перекрыть диапазон различных клинических исследований, могут быть выбраны дополнительные уровни. Для этих дополнительных уровней могут быть выбраны другие уставки системы, и их следует поддерживать постоянными в течение процедуры испытаний.

Изменение воздушной кермы должно достигаться путем изменения анодного тока, или время облучения должно быть подобно тому, которое применяется при клиническом использовании цифрового рентгеновского детектора. Лаг-эффекты должны быть устранены (см. 4.6.3).

Радиационные условия должны быть указаны вместе с результатами испытаний (см. раздел 7).

4.6.2 Измерение воздушной кермы

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

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

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

Примечания

1. Необходимо принять во внимание ослабление в воздухе.

2. Если коллиматор перед пациентом является многощелевым, экспозиция должна быть интегрирована в течение сканирования. Многощелевая диафрагма создает негомогенное радиационное поле на поверхности детектора, поэтому для получения правильного результата необходимо провести более длительное сканирование по детектору излучения.

Если используется детектор-монитор, следует начертить график кривой в функции расстояния d между фокусным пятном и детектором излучения по следующему уравнению:

f(d) = \sqrt{\frac{\text{показание детектора} - \text{монитора}}{\text{показание детектора излучения}}} .

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

Если детектор-монитор не используют, то вычерчивают график зависимости квадратного корня от обратной величины показания измерителя излучения в функции от расстояния между измерителем излучения и фокусным пятном. Экстраполяцию и т.д. проводят так же, как изложено в предыдущем абзаце.

4.6.3 Исключение лаг-эффектов

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

Это влияние можно разделить на две компоненты — аддитивную (аддитивная ветвь) и мультипликативную (изменение прироста). Величины обеих компонент должны быть определены. Для общей информации см. [10], [11] и [12].

При определении возможных лаг-эффектов цифровой детектор должен работать в соответствии со спецификациями изготовителя. Между двумя последовательными экспозициями необходимо выдерживать минимальный интервал времени (определяемый с помощью испытаний по Приложению A) во избежание искажающего влияния лаг-эффектов на измерение квантовой эффективности регистрации.

Примечание. На лаг-эффекты могут влиять следующие факторы: время облучения по отношению к считыванию, метод очистки от остатков предыдущего облучения, время от стирания до повторного облучения, время от считывания до повторного облучения, включение «фиктивного» считывания для стирания эффектов от предыдущего облучения.

Для определения значения лаг-эффектов должны быть использованы процедуры испытаний, приведенные в Приложении A.

4.6.4 Облучение для определения функции преобразования

Уставки устройства для получения цифрового рентгеновского изображения должны быть такими же, как при экспонировании тест-объекта. Облучение должно быть выполнено в геометрии рисунка 2, но без тест-объекта в пучке. Воздушную керму измеряют в соответствии с 4.6.2. Функция преобразования должна быть определена при уровне воздушной кермы от нулевой до значения, на 20% большего, чем максимальный испытанный уровень воздушной кермы.

Коэффициент преобразования при нулевом уровне воздушной кермы должен быть определен по темному изображению, которое реализуется при тех же условиях, как рентгеновское изображение. Минимальный уровень рентгеновской воздушной кермы должен быть не больше, чем 1/5 опорного уровня воздушной кермы.

В зависимости от процедуры расчета (см. 6.3.1) меняется число различных экспозиций; если нужно проверить только линейность функции преобразования, достаточно пяти экспозиций, распределенных равномерно в желаемом диапазоне. Если нужно определить полную функцию преобразования, воздушная керма должна изменяться таким образом, чтобы максимальное увеличение логарифмической (по основанию 10) воздушной кермы было не больше 0,1.

4.6.5 Облучение для определения спектра мощности шума

Уставки устройства для получения цифрового рентгеновского изображения должны быть такими же, как при экспонировании тест-объекта. Облучение должно быть выполнено в геометрии рисунка 2, но без тест-объекта в пучке. Воздушную керму измеряют в соответствии с 4.6.2.

Для приблизительного расчета спектра мощности шума, который используется в дальнейшем для расчета DQE, используют прямоугольную площадку размерами приблизительно 50 x 50 мм, расположенную в центре облучаемой области размерами 100 x 100 мм.

Для этой цели набор входных данных должен состоять, по меньшей мере, из 4 млн независимых пикселей, образующих одно или несколько независимых изображений, каждое из которых содержит не менее 256 пикселей в любом пространственном направлении. Если необходимо более чем одно изображение, все изображения должны быть получены при одних и тех же качестве излучения и уровне воздушной кермы. Стандартное отклонение излучений, использованных для получения разных изображений, должно быть не больше 10% среднего значения.

Примечание. Минимальное число требуемых пикселей изображения определяется требуемой точностью, что определяет минимальное число областей интереса. Для получения точности двумерного спектра мощности шума 5% необходимо минимум 960 областей интереса (частично перекрывающихся), содержащих 16 млн независимых пикселей изображения в данных размерах области интереса. Процесс усреднения, проводимый в дальнейшем для получения одномерного среза, снижает минимальное число необходимых пикселей изображения до 4 млн, что все еще достаточно для достижения необходимой точности.

Должно быть обращено внимание на то, чтобы не было корреляции между последовательными изображениями (лаг-эффект; см. 4.6.2). При проведении облучений не допускается изменение уставок системы.

Изображения для определения спектра мощности шума должны быть сделаны при уровнях воздушной кермы по 4.6.1.

4.6.6 Облучение с тест-объектом в пучке излучения

Облучение следует проводить в геометрии рисунка 2. Тест-объект размещают непосредственно на поверхности детектора. Тест-объект должен быть размещен таким образом, чтобы край был наклонен на угол \alpha относительно оси столбцов пикселей или рядов пикселей, должен быть в интервале от 1,5° до 3°.

Примечание. Метод наклона тест-объекта относительно рядов или столбцов матрицы изображения является общим с другими стандартами (ИСО 15529 и ИСО 12223) и изложен в ряде публикаций, связанных с получением предварительных образцов функции передачи модуляции.

Должно быть сделано, по меньшей мере, два облучения с тест-объектом в пучке излучения, один с тест-объектом, ориентированным приблизительно вдоль столбцов, и один с тест-объектом, ориентированным приблизительно вдоль рядов матрицы изображения. Как известно, в системах CR резкость зависит от ориентации края относительно направления перемещения лазерного луча при сканировании. Поэтому в системах CR должно быть сделано четыре облучения с тест-объектом в пучке излучения, с поворотом тест-объекта на 90° между последовательными экспозициями. Положение других компонентов не должно изменяться. Для нового положения следует провести новую настройку тест-объекта.

Изображения для определения MTF следует производить при одном из трех уровней воздушной кермы (см. 4.6.1).

5. КОРРЕКТИРОВКИ НЕОБРАБОТАННЫХ ДАННЫХ

Допускаются следующие линейные и не зависящие от изображения корректировки данных для определения функции преобразования, спектра мощности шума и функции передачи модуляции.

Все нижеследующие коррекции следует проводить, как при нормальном клиническом использовании:

• замена необработанных данных от плохих или дефективных пикселей подходящими данными;
• коррекция плоского поля, включающая в себя:

коррекцию неоднородности радиационного поля;

коррекцию контраста индивидуальных пикселей;

коррекцию сигнала индивидуальных пикселей;

коррекцию изменения скорости в течение сканирования;

• коррекция геометрической дисторсии.

1. Некоторые детекторы осуществляют линейную обработку изображений в соответствии с их физической концепцией. Так как эта обработка является линейной и независимой от изображения, эти операции допускаются как исключение.

2. Коррекция изображения считается независимой от изображения, если такой же коррекции подвергаются все изображения независимо от их содержания.

6. ОПРЕДЕЛЕНИЕ КВАНТОВОЙ ЭФФЕКТИВНОСТИ РЕГИСТРАЦИИ

6.1 Определение и формула DQE(u, v)

Квантовую эффективность регистрации DQE(u, v) в зависимости от частоты определяют по следующей формуле:

DQE(u,v) = G^{2}MTF^{2}(u,v) \frac{W_{in}(u,v)}{W_{out}(u,v)} .(1)

Источником этой формулы является [4], формула 2.153.

В настоящем стандарте спектр мощности шума на выходе W_{out}(u,v) и функцию передачи модуляции MTF(u, v) устройства для получения цифрового рентгеновского изображения следует рассчитывать по линеаризованным данным. Линеаризованные данные рассчитывают из исходных данных с помощью обратной функции преобразования и выражают числом квантов экспозиции на единицу площади. Сигнал детектора при нулевой пространственной частоте G в формуле (1) является частью функции преобразования и не нуждается в отдельном определении.

Поэтому в настоящем стандарте рабочей формулой для определения квантовой эффективности регистрации DQE(u, v) в зависимости от частоты является

где MTF(u, v) — предварительная функция передачи модуляции устройства для получения цифрового рентгеновского изображения;

W_{in}(u,v) — спектр мощности шума радиационного поля на поверхности детектора, определяемый в соответствии с 6.2;

W_{out}(u,v) — спектр мощности шума на выходе устройства для получения цифрового рентгеновского изображения, определяемый в соответствии с 6.3.2.

6.2 Параметры, используемые при расчетах

W_{in(u,v)} = K_{a} \cdot SNR_{in}^{2},(3)

где K_{a} — измеренная воздушная керма, мкГр;

SNR_{in}^{2} — квадрат отношения сигнал — шум к воздушной керме, 1/мм2мкГр, значение, приведенное в таблице 2, столбец 4.

В настоящем стандарте следует применять значения SNR_{in}^{2} по таблице 2.

Радиационный параметр SNR_{in}^{2} для применения в настоящем стандарте (дополнительный фильтр 2 мм AI)

Качество излучения	Толщина фильтра, мм	Номинальное анодное напряжение, кВ	Расчетное значение SNR_{in}^{2} , 1/мм2мкГр
Mo/Mo (RQA-M 1)	0,032	25	4639
Mo/Mo (RQA-M 2)	0,032	28	4981
Mo/Mo (RQA-M 3)	0,032	30	5303
Mo/Mo (RQA-M 3)	0,032	35	6325
Mo/Rh	0,025	28	5439
Rh/Rh	0,025	28	5944
W/Rh	0,050	28	5975
W/AI	0,500	28	6575

6.3 Определение различных параметров по изображениям

6.3.1 Линеаризация данных

6.3.2 Спектр мощности шума

Однородно экспонированная площадь цифрового рентгеновского детектора должна быть поделена на квадратные площадки, называемые областями интереса (ОИ). Каждая область интереса для расчета индивидуального образца спектра мощности шума должна быть размерами 256 x 256 пикселей. Эти площадки должны перекрываться на 128 пикселей в горизонтальном и вертикальном направлениях (см. рисунок 3). Первая площадка должна находиться в верхнем левом углу изображения, тогда следующая площадка получается передвижением на 128 пикселей вправо в горизонтальном направлении, создавая вторую площадку, которая на половину перекрывает первую. Следующая площадка получается перемещением второй вновь на 128 пикселей. Это повторяется до конца первой горизонтальной «полосы». Начиная вновь с левой стороны изображения и одновременно смещаясь на 128 пикселей в вертикальном направлении, получают вторую горизонтальную «полосу». Смещение в вертикальном направлении создает следующие полосы, пока вся площадь размерами примерно 50 x 50 мм не будет покрыта областями интереса.

Размер областей интереса должен быть n = 256

Рисунок 3. Расположение областей интереса

Может быть проведено направленное удаление с помощью приложения двумерного полинома второго порядка к линеаризованным данным каждого полного изображения, используемого для расчета спектра, и вычитания этой функции [ S(x_{i},y_{i}) , см. формулу (4)] из линеаризованных данных. Для каждой области интереса вычисляют двумерное преобразование Фурье без применения какого-либо окна.

Двумерное преобразование Фурье применяют в соответствии с формулой (4). Исходя из формулы 3.44 в (см. [4]), получают рабочую формулу для определения спектра мощности шума

W_{out}(u_{n},v_{k}) = \frac{\Delta x\Delta y}{M \cdot 256 \cdot 256} \sum_{m = 1}^{M} | \sum_{i = 1}^{256} \sum_{j = 1}^{256} (I(x_{i}y_{j}) - S(x_{i},y_{j}))exp(- 2\pi i(u_{n}x_{i} + v_{k}y_{j}))|^{2},(4)

где \Delta x , \Delta y — размер пикселя соответственно в горизонтальном или вертикальном направлении;

M — число областей интереса;

I(x_{i}y_{j}) — линеаризованные данные;

S(x_{i},y_{j}) — двумерный полином, применяемый опционно.

Усредненный двумерный спектр мощности шума получается при усреднении по всем измеренным образцам спектра.

Для получения одномерных срезов двумерного спектра мощности шума вдоль осей в плоскости пространственных частот используют 15 рядов или столбцов двумерного спектра по каждой оси. Однако усреднение проводят только по семи рядам или столбцам спектра мощности шума по обеим сторонам соответствующих осей (общее число 14), не учитывая измерения по самим осям. Для всех точек должны быть рассчитаны точные пространственные частоты с учетом расстояния от начала по радиусу. Должно быть проведено сглаживание путем усреднения данных по 14 рядам и столбцам, которые находятся в интервале частот 2f_{int}(f - f_{int} \leq f \leq f + f_{int}) вокруг пространственных частот, которые должны быть указаны в отчете (см. раздел 7). f_{int} вычисляют по формуле

f_{int} = \frac{0,001}{\text{шаг пикселя }(\text{мм})} .

Примечание. Тот факт, что диапазон частот устанавливается в зависимости от шага пикселей, всегда обеспечивает постоянство числа используемых точек данных независимо от шага пикселей, а следовательно, постоянство точности.

Чтобы определить, влияют ли эффекты квантования на спектр мощности шума, должны быть рассчитаны вариации исходных данных (DN), которые используются для расчета спектра мощности шума, для одного изображения. Если вариации больше 0,25 (см. ИСО 12232), можно считать, что шумы квантования пренебрежимо малы. Если вариации меньше 0,25, данные признают не подходящими для определения спектра мощности шума.

Примечание. В общем случае вариации исходных данных больше, чем четверть интервала квантования. Только если число бит при квантовании очень мало, вариации могут быть меньше. Для расчета вариаций квантования, т.е. 1/12, предполагается, что аналоговые величины, которые оцифровываются, имеют равномерное или прямоугольное распределение по отношению к каждому интервалу квантования [2].

Если спектр мощности шума определяется вдоль диагонали (45° по отношению к горизонтальной или вертикальной оси), усреднение с отдельных образцов должно проводиться аналогично тому, как описано в предыдущих абзацах, но включая значения вдоль диагонали. Эти измерения для 45° могут также потребовать усреднения по сечениям, смежным с 45°, для улучшения точности определения NPS.

6.3.3 Определение функции передачи модуляции

Для тех систем, в которых коррекция неоднородности радиационного поля не предусмотрена при нормальном клиническом использовании, должна быть проведена двумерная коррекция изображения тест-объекта. В радиационном поле однородно экспонированного изображения область интереса выбирают, по крайней мере, в 1,5 больше, чем область интереса при определении MTF, показанная на рисунке 1. По отношению к линеаризованным данным однородно экспонированного изображения применяют двумерный полином второго порядка [S(x_{i},y_{j})] . Проводят направленное удаление в соответствующих областях интереса в изображении тест-объекта по отношению к линеаризованным данным в соответствии с формулой 5

I_{\text{кор}}(x_{i},y_{j})/S(x_{i},y_{j}) \cdot S_{\text{ср}},(5)

где S_{\text{ср}} — усредненное значение пикселя линеаризованных данных в области интереса однородно экспонированного изображения.

Предварительная функция передачи модуляции должна быть определена вдоль двух взаимно перпендикулярных осей, параллельных рядам или столбцами матрицы изображения. Предварительную функцию передачи модуляции определяют как усредненную из тех, которые получены из изображений объекта, повернутого приблизительно на 180°, если такой поворот возможен.

Для определения MTF должна быть использована полная длина функции распределения края (ESF), определенной для области интереса, показанной на рисунке 1.

Определяют целое число № линий (т.е. рядов или столбцов), приводящих к боковому сдвигу края в направлении линии, наиболее близко соответствующей расстоянию пикселей. Для этого могут быть использованы разные методы. Один из них состоит в том, что определяют угол \alpha между краем и столбцами или рядами матрицы изображения и вычисляют N = округленно (1/tg\alpha) , где «округление» обозначает выбор ближайшего целого числа. № должно быть целым числом.

Для создания образцового краевого профиля или ESF используют значения пикселей из линеаризованных данных (см. 6.3.1) № последовательных линий (т.е. рядов или столбцов), пересекающих край. Значение первого пикселя первой линии дает первую точку данных в образцовой ESF, первый пиксель во второй линии — вторую точку данных, первый пиксель в N-й линии — N-ю точку данных. Эту процедуру повторяют для других пикселей в № последовательных линиях, например, значение второго пикселя в первой линии дает (N + 1)-ю точку данных, второй пиксель во второй линии — (N + 2)-ю точку данных и т.д.

Для расчета усредненной ESF эту процедуру повторяют для других групп из № последовательных линий вдоль края. Проводят усреднение всех функций распределения края и на основе этой усредненной образцовой ESF рассчитывают MTF.

Расстояние в образцовой ESF считается постоянным и равным длине участка пикселей \Delta x , деленной на N, т.е. получается ESF(x_{n}) при x_{n} = n(\Delta x/N) . Образцовая ESF дифференцируется на отрезке [-1, 0, 1] или [-0,5, 0, 0,5], при этом получается образцовая функция распределения линии (LSF). Эффект спектрального сглаживания от дифференцирования с конечным элементом может быть скорректирован [6]. Вычисляют преобразование Фурье функции распределения линии, и модуль этого преобразования дает MTF. MTF нормализован по значению при нулевой частоте. Поскольку расстояние от индивидуальных пикселей до края рассчитывают вдоль направления линии, а не в направлении, перпендикулярном к краю, может быть проведена коррекция шкалы по оси частот (коэффициент шкалы 1/cos\alpha ).

Примечание. Если коррекцию с коэффициентом 1/cos\alpha не проводят, ошибка коэффициента шкалы составляет 0,1%.

Для получения MTF при пространственных частотах, которые должны быть указаны в декларации соответствия (см. раздел 7), должны быть получены данные в интервале частот 2f_{int}^{\text{мм}^{- 1}}(f - f_{int} \leq f \leq f + f_{int}) , (определение f_{int} см. 6.3.2) вокруг этих пространственных частот.

7. ФОРМАТ ДЕКЛАРАЦИИ СООТВЕТСТВИЯ

При указании квантовой эффективности регистрации должны быть также указаны следующие параметры:

• качество излучения в соответствии с таблицей 1, в случае использования другого качества излучения — материал мишени, материал и толщина фильтра, анодное напряжение, слой половинного ослабления (СПО), в мм AI, и использованное значение SNR_{in}^{2} ;

• расстояние между фокусным пятном и поверхностью детектора, если оно меньше 600 или больше 700 мм;
• отклонения от рекомендуемой геометрии (см. 4.5);
• метод, использованный для определения MTF, и его обоснование в случае использования краевого метода, отличного от стандартизованного;
• указания всех расстояний при получении изображения, например, расстояние между тест-объектом и плоскостью детектора;
• условия окружающей среды.

Результаты измерения DQE должны быть приведены в виде чисел, сведенных в таблицу. Значения DQE должны быть указаны для пространственных частот 0,5\text{мм}^{- 1} , 1\text{мм}^{- 1} , 1,5\text{мм}^{- 1} вплоть до наибольшей пространственной частоты, которая ниже частоты Найквиста. В таблицу могут быть добавлены другие параметры. Дополнительно результаты измерения могут быть представлены в виде графика как значения DQE(u, v) в функции пространственных частот при воздушной керме в качестве параметра, по обеим осям используется линейная шкала.

В общем случае значения DQE(u, v) должны быть приведены для обеих осей, горизонтальной и вертикальной. Если отношение DQE(u, o)/DQE(o, v)|u = v находится в диапазоне от 0,9 до 1,1, значения DQE(u, v) для обеих осей могут быть усреднены и указаны как действительные для обеих осей.

Дополнительно могут быть приведены значения DQE вдоль диагональной оси, при этом должно быть ясно указано, что DQE относится к диагональной оси.

8. ТОЧНОСТЬ

Неопределенность DQE следует определять в соответствии с инструкциями GUM [2], используя формулу (2).

Неопределенность представленных значений DQE (коэффициент покрытия 2 в соответствии с [2]) должна быть меньше, чем

\Delta(DQE(u))/(DQE(u)) = \pm 0,1,

причем выбирается большее значение.

Неопределенность следует указывать в листах данных.

Приложение A (обязательное)

ОПРЕДЕЛЕНИЕ ЛАГ-ЭФФЕКТОВ

A.1. Испытания для определения аддитивных лаг-эффектов

1. По методу, изложенному в 4.6.6, проводят облучение краевого тест-объекта. Следует убедиться, что объект правильно установлен в пучке в соответствии с 4.6.6. Облучение должно быть проведено при «опорном» уровне воздушной кермы по 4.6.1.

2. Создают изображение, полученное в результате облучения по перечислению 1), следуя методу, предложенному изготовителем.

3. Следуют всем этапам, которые являются частями предлагаемого изготовителем метода при работе цифрового рентгеновского детектора в промежутке между облучениями.

4. Не проводя облучения поверхности детектора, создают второе изображение по методу перечисления 2.

5. Записывают промежуток времени между первым (после облучения) и вторым (без облучения) считываниями с цифрового рентгеновского детектора. Большее из этого промежутка и времени, определенного в A.2, должно быть минимальным временем между снимками, используемым при определении функции преобразования, спектра мощности шума и MTF.

6. На изображении по перечислению 2) (после облучения) измеряют усредненное значение линеаризованных данных в прямоугольном участке, включающем, по крайней мере, 1000 пикселей, смежном, но не перекрывающем площадь изображения объекта высокого контраста по перечислению 2) (область интереса 2, рисунок A.1).

Примечание. 1000 пикселей — минимальное значение, необходимое для обеспечения того, что относительная разница 0,005 детектируется с доверительной вероятностью 95% при вероятности детектирования 80%. Предпочтительнее использовать 10000 пикселей.

Рисунок A.1. Определение областей интереса (ОИ)

7. На изображении по перечислению 4) (без облучения) измеряют усредненное значение линеаризованных данных в прямоугольном участке, включающем, по крайней мере, 1000 пикселей, смежном, но не перекрывающем площадь изображения объекта высокого контраста по перечислению 2) (область интереса 2, рисунок A.1).

8. На изображении по перечислению 4) (без облучения) измеряют усредненное значение линеаризованных данных в прямоугольном участке, включающем, по крайней мере, 100 пикселей, в пределах области, покрытой изображением объекта высокого контраста (область интереса 1, рисунок A.1).

9. Испытание признают удовлетворительным, если разность измерений по перечислениям 7) и 8), разделенная на измерение по перечислению 6), меньше 0,005.

Следует отметить, что наличие лаг-эффектов за тест-объектом, даже если их значение меньше 0,5%, может негативно влиять на определение MTF.

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

1. По методу, изложенному в 4.6.1, проводят облучение без объекта в пучке, используя «опорный» уровень воздушной кермы.

2. Создают изображение, полученное в результате облучения по перечислению 1) (изображение 1, после облучения без тест-объекта),следуя методу, предложенному изготовителем.

3. Следуют всем этапам, которые являются частями предлагаемого изготовителем метода при работе цифрового рентгеновского детектора в промежутке между облучениями.

4. По методу, изложенному в 4.6.6, проводят облучение краевого тест-объекта. Следует убедиться, что объект правильно установлен в пучке в соответствии с 4.6.6. облучение должно быть проведено при наиболее высоком уровне воздушной кермы, использованном при измерениях.

5. Создают изображение, полученное в результате облучения по перечислению 4) (изображение 2), следуя методу, предложенному изготовителем.

6. Следуют всем этапам, которые являются частями предлагаемого изготовителем метода при работе цифрового рентгеновского детектора в промежутке между облучениями.

7. По методу, изложенному в 4.6.1, проводят второе облучение без объекта в пучке, используя нормальный уровень воздушной кермы.

8. Создают изображение, полученное в результате облучения по перечислению 7) (изображение 3, после облучения без тест-объекта), следуя методу, предложенному изготовителем. Записывают промежуток времени между вторым (после облучения с тест-объектом) и третьим (после облучения без тест-объекта) считываниями с цифрового рентгеновского детектора. Больший из этого промежутка и времени, определенного в A.1, должен быть минимальным временем между снимками, используемым при определении функции преобразования, спектра мощности шума и MTF.

9. На изображениях 1 и 3 соответственно измеряют усредненное значение линеаризованных данных в прямоугольном участке, включающем, по крайней мере, 1000 пикселей в пределах области, покрытой изображением объекта высокого контраста (область интереса 1, рисунок A.1).

10. На изображениях 1 и 3 соответственно измеряют усредненное значение линеаризованных данных в прямоугольном участке, включающем, по крайней мере, 1000 пикселей, смежном, но не перекрывающем площадь изображения объекта высокого контраста (область интереса 2, рисунок A.1).

Это обеспечивает, что вклад лаг-эффекта меньше, чем 0,5% действующей воздушной кермы.

Приложение B (рекомендуемое)

РАСЧЕТ ВХОДНОГО СПЕКТРА МОЩНОСТИ ШУМА

Входной спектр мощности шума равен входящему потоку фотонов (см. [4], формула 2.134)

W_{in}(u,v) = Q,(B.1)

где Q — поток фотонов, т.е. число квантов экспозиции на единицу площади, 1/мм2.

Q зависит от спектра рентгеновского излучения и уровня воздушной кермы:

Q = K_{a} \cdot \int (\Phi(E)/K_{a})dE = K_{a} \cdot SNR_{in}^{2},(B.2)

где K_{a} — воздушная керма, мкГр;

E — энергия рентгеновского излучения, кэВ;

\Phi(E)/K_{a} — спектральный поток рентгеновского излучения на воздушную керму, 1/(мм2кэВмкГр);

Значения, приведенные в таблице 2, рассчитаны с использованием компьютерной программы BOONE [7], основанной на методе Дж. М. Буне [8]. Использованы коэффициенты взаимодействия из NIST [9]. При использовании других программ полученные значения могут незначительно отличаться.

Приложение ДА (справочное)

Таблица ДА.1

Обозначение ссылочного международного стандарта	Степень соответствия	Обозначение и наименование соответствующего национального стандарта
МЭК 60336:2005	IDT	ГОСТ Р МЭК 60336-2010 Излучатели медицинские рентгенодиагностические. Характеристики фокусных пятен
МЭК 61267:2001	-	<*>
МЭК 62220-1:2003	IDT	ГОСТ Р МЭК 62220-1-2006 Изделия медицинские электрические. Характеристики цифровых приемников рентгеновского изображения. Часть 1. Определение квантовой эффективности регистрации
ИСО 12232:1998	-	<*>
<*> Соответствующий национальный стандарт отсутствует. До его утверждения рекомендуется использовать перевод на русский язык данного международного стандарта. Перевод данного международного стандарта находится в Федеральном информационном фонде технических регламентов и стандартов.
Примечание. В настоящей таблице использовано следующее условное обозначение степени соответствия стандартов:
- IDT — идентичные стандарты.

ТЕРМИНЫ

МЭК 60788	МР-...-...-...
Сокращенный термин	МР-...-...c
Термин, определяемый в настоящем стандарте	3.xx
Автоматическое управление экспозиционной дозой	МР-36-46
Анодное напряжение	МР-36-02
Анодный ток	МР-36-07
Блок источника излучения	МР-20-05
Воздушная керма	МР-13-11
Время облучения	МР-36-11
Детектор облучения	МР-51-01
Диафрагма	МР-37-29
Дополнительный фильтр	МР-35-02
Измеритель излучения	МР-50-01
Исходные данные DN	3.11
Качество излучения	МР-13-28
Квантовая эффективность регистрации, DQE(u, v)	3.2
Лаг-эффект	3.6
Линеаризованные данные	3.7
Матрица изображения	3.5
Необработанные данные	3.13
Номинальное значение	МР-20-14
Облучение	МР-12-09
Отсеивающий растр	МР-32-06
Пиксель	МР-32-60
Плоскость приемника изображения	МР-37-15
Поверхность детектора	3.3
Поток фотонов	3.12
Пространственная частота u или v	3.14
Процентная пульсация	МР-36-17
Пучок излучения	МР-37-05
Рассеянное излучение	МР-11-13
Рентгеновская трубка	МР-22-03
Рентгеновский аппарат	МР-20-20
Рентгеновский генератор	МР-20-17
РПУ постоянного напряжения	МР-21-06
Слой половинного ослабления	МР-13-42
Спектр мощности шума (NPS), W(u, v)	3.10
Тест-объект	МР-71-04
Усилитель рентгеновского изображения	МР-32-39
Устройство для получения цифрового рентгеновского изображения	3.4
Фокусное пятно	МР-20-13С
Функция передачи модуляции MTF(u, v)	3.8
Функция преобразования	3.1
Шум	3.9

БИБЛИОГРАФИЯ

Литература, на которую даются ссылки

1. Доклад ИКРУ 54:1996 Медицинские изображения. Оценка качества изображения

2. БИПМ, МЭК, ИФСС, ИСО, ИУПАК, ИУПАП и ОИМЛ: 1993, Руководство по выражению неопределенности измерений

3. Метц У.Ц., Вагнер Р.Ф., Дои К., Броун Д.Г., Нишикава Р.М., Майерс К.Дж. О консенсусе при квантовой оценке медицинских систем изображения. Med. Phys., 1995, 22, с. 1057 — 1061

4. Справочник по медицинским изображениям. Т.1: Физика и психофизика. Издатели: Дж. Бойтель, Х.Л. Кундель, Р.Л. Ван Меттер, SPIE 2000

5. Тапиовара М.Дж., Вагнер Р.Ф. Анализ SNR и DQE рентгеновских изображений широкого спектра. Phys. Med. Biol. 1985, 30, с. 519 — 529 и перепечатка, Phys. Med. Biol., 1986, 31,с. 195

6. Каннингхе Н.А., Фенстер А. Метод определения функции передачи модуляции с помощью краевых профилей с коррекцией на различия определенных элементов. Med. Phys., 1987, 14, с. 533 — 537

7. Пакет программ BOONE, версия 1.2, 2005 (Л. Батц, Siemens Medical Solutions, Эрланген, Германия)

8. Буне Дж.М., Фьюелл Т.Р., Дженнингс Р.Дж. Спектральные модели молибденового, родиевого и вольфрамового анода с использованием интерполирующих полиномов и приложение к маммографии. Med. Phys., 1997, 24(12), с. 1863 — 1874

9. Веб-сайт НИСТ: http//physics.nist.gov/PhysRetData/XrayMassCoef/cover.html

10. Гранфорс П.Р., Ауфрихтер Р. DQE(f) рентгеновского детектора с плоской панелью на основе аморфного кремния: факторы, влияющие на параметры детектора, и методология измерения Proc. SPIE, 2000, 3977, с. 2 — 13

11. Мензер Б., Бастьяенс Р.Дж.М., Насцетти А., Овердик М., Симон М. Линейные системные модели лаг-эффекта в плоских динамических рентгеновских детекторах. Proc. SPIE, 2005, 5745, с. 430 — 441

12. Овердик М., Сольф Т., Вишманн Х.А. Временные артефакты в плоских динамических рентгеновских детекторах. Proc. SPIE, 2001, 4320, с. 47 — 58

Другая литература, представляющая интерес

Дейнти Дж.К., Шоу Р. Наука об изображениях. Academic Press, Лондон, 1974, глава 5, с. 153

---