С.В. МАЛЬЦЕВ1, В.И. ВОРОНОВ2, Т.В. КОЛЕСНИЧЕНКО1
1Казанская государственная медицинская академия, 420012, г. Казань, ул. Бутлерова, д. 36
2Казанский национальный исследовательский технический университет им. А.Н. Туполева (КНИТУ-КАИ), 420111, г. Казань, ул. К. Маркса, д. 10
Мальцев Станислав Викторович — доктор медицинских наук, профессор, заведующий кафедрой педиатрии с курсом поликлинической педиатрии, тел. (843) 236-20-84, e-mail: [email protected]1
Воронов Виктор Иванович — доктор технических наук, профессор кафедры радиоэлектронных и квантовых устройств, тел. (843) 238-40-67, e-mail: [email protected]2
Колесниченко Татьяна Викторовна — соискатель кафедры педиатрии с курсом поликлинической педиатрии, тел. (843) 236-20-84, e-mail: [email protected]1
Оценка минеральной плотности кости (МПК) и выявление остеопении и остеопороза у детей младших возрастных групп имеет особое значение, поскольку патология костной массы в этом возрасте играет определяющую роль в течение всего последующего развития ребенка. В статье представлено описание программы (X-Ray Analysis), разработанной авторами для измерения МПК по цифровым рентгенограммам отсканированных пленочных снимков, получаемых на рентгеновских аппаратах общего назначения. Результаты измерений рентгенограмм позволили выявить отличия МПК между двумя группами обследованных новорожденных.
Ключевые слова: минеральная плотность кости у детей, остеопения, остеопороз, компьютерный анализ цифровых рентгенограмм.
S.V. MALTSEV1, V.I. VORONOV2, T.V. KOLESNICHENKO1
1Kazan State Medical Academy, 49 Butlerova St., Kazan, Russian Federation 420012
2TupolevKazan National Research Technical University (KNRTU-KAI), 10 K. Marx St., Kazan, Russian Federation 420111
Using computer analysis of film X-ray patterns for determining of bone mineral density in new-born children
Maltsev S.V. — D. Med. Sc., Professor, Head of the Department of Pediatrics with a course of outpatient pediatrics, tel. (843) 236-20-84, e-mail: [email protected]1
Voronov V.I. — Doctor of Engineering, Professor of the Department of radio-electronic and quantum devices, tel. (843) 238-40-67, e-mail: [email protected]2
Kolesnichenko T.V. — Candidate of the Department of pediatrics with a course outpatient pediatrics, tel. (843) 236-20-84, e-mail: [email protected]1
Evaluation of bone mineral density (BMD) and detection of osteopenia and osteoporosis in children of younger age groups are of particular importance, because pathology of bone mass at that age plays a crucial role throughout the evolution of the child. The article presents the description of the program (X-Ray Analysis), developed by the authors to measure BMD by digitized radiographs of scanned film images from X-ray equipment of general purpose. The measurement results of radiographs revealed differences of BMD between two groups of surveyed newborns.
Key words: bone mineral density in children, osteopenia, osteoporosis, computer analysis of digital radiographs.
Одной из проблем педиатрии продолжает оставаться проблема остеопении и остеопороза, которые встречаются у детей всех возрастов. Особое значение эта проблема имеет у детей младших возрастных групп, поскольку недостаток или избыток костной массы в этом возрасте играет определяющую роль в течение всего последующего развития ребенка. Примерно 80% пиковой костной массы генетически детерминировано и только 20% определяются факторами образа жизни. Учитывая, что все периоды детства (особенно младенчество и подростковый период) являются важными для формирования костной ткани, возникает необходимость иметь оптимальные методы измерения минеральной плотности кости (МПК) для клинических и исследовательских целей, расширению возможностей костной денситометрии у детей [1].
Главным показателем остеопении и остеопороза считается величина минеральной плотности кости, которая является основным параметром известных критериев костной патологии у детей (Z-критерий) и взрослых (T-критерий). Для определения МПК в основном используются следующие методы [2]: двухэнергетическая рентгеновская абсорбциометрия (dual-energy X-ray absoptiometry — DEXA), метод компьютерной количественной томографии (ККТ), метод рентгеновской денситометрии и метод ультрасонометрии (количественное ультразвуковой исследование). Анализ публикаций свидетельствует о том, что все эти методы имеют ряд недостатков. В частности, определение МПК каждым из методов требует наличия современного и достаточно дорогостоящего оборудования. Кроме того, данные по плотности костной ткани, полученные различными методами, не всегда сравнимы между собой, поскольку имеют различный физический смысл, различные численные значения, а также несовпадающую размерность. Так, например, размерность МПК, получаемая методами DEXA и рентгеновской денситометрии выражается в г/см2 (количество эквивалента минеральной плотности кости — гидроксиапатита кальция, приходящегося на единицу площади поперечного сечения просвечиваемого материала). В методе ККТ определяется по сути дела физическая объемная плотность костной ткани, характеризуемая размерностью г/см3 (масса костной ткани, приходящая на единицу объема). Ультрасонометрия оперирует скоростью распространения ультразвуковой волны в кортикальном слое кости и измеряемые данные имеют размерность м/c (метр в секунду). Именно по этой причине в 1994 году рабочей группой ВОЗ было решено диагностировать снижение МПК на основании универсального критерия (Z-критерия у детей), выраженного в единицах стандартного отклонения (SD), единого для всех методов костной денситометрии.
Следует также отметить, что при диагностике костной патологии у детей имеются известные ограничения по дозе рентгеновского излучения так, что сравнение в этом отношении методов DEXA и метода рентгеновской денситометрии оказывается в пользу последнего, поскольку пациент при использовании метода DEXA подвергается облучению в течение достаточно длительного периода времени (3-4 минуты [3]), а при использовании метода рентгеновской денситометрии эта операция выполняется одномоментно.
В клинической практике в нашей стране, особенно в сельских лечебных заведениях, неоснащенных современными аппаратами для DEXA, ККТ или ультрасонометрии, но имеющих стандартное рентгеновское оборудование, часто требуется оценка минеральной плотности костной ткани по рентгеновским снимкам, что можно реализовать на основе денситометрических измерений. В этом случае МПК оценивается по степени почернения участков рентгеновской фотопленки в характерных областях интереса. Однако степень почернения зависит не только от плотности материала кости, но и от значительного числа дополнительных факторов. Главными из них являются интенсивность рентгеновского излучения источника, характерные свойства самой рентгеновской фотопленки, в частности, ее контрастность и чувствительность к рентгеновскому излучению, условия последующей обработки в проявляюще-закрепляющих растворах, химический состав и чистота растворов, температура обработки и ряд других факторов, меняющихся от снимка к снимку.
На практике проконтролировать влияние всех вышеупомянутых факторов, действующих на пленку в конкретном сеансе получения рентгеновского изображения, не представляется возможным. Поэтому для их, по сути дела компенсационного учета, используют различные калибровочные элементы — тест-фантомы, располагаемые при получении рентгенограмм рядом с диагностируемым объектом [4]. В рентгеновских аппаратах общего назначения в качестве калибровочных элементов используют тест-фантомы в виде ступенчатых клиньев, выполненных из материала, поглощающего рентгеновское излучение. Высота каждой ступени в направлении прохождения рентгеновского излучения определяет степень почернения рентгеновской пленки для области, соответствующей этой ступени. Эквивалентное значение МПК для каждой ступени при известном материале клина можно найти либо расчетным методом [5], либо экспериментально на аттестованных денситометрах. Таким образом, сравнивая степень почернения ступеней с почернением областей на рентгеновском изображении диагностируемого медицинского объекта можно определить МПК для этих областей.
Современная компьютерная техника позволяет проводить сравнение степени почернения участков тест-фантома с областями на рентгеновском снимке с помощью компьютерных программ. С этой целью снимок сканируется, оцифровывается и вводится в программу. Однако стандартные медицинские программы, с помощью которых можно выполнить оценку МПК по снимкам, получаемым с помощью рентгеновской аппаратуры общего назначения, производителями такого оборудования не выпускаются. Тем не менее, информация по использованию подобных программных средств имеется в ряде публикаций. В частности, в работах С.В. Мальцева и Н.В. Зиатдиновой [6, 7] с помощью компьютерной программы проводилась оценка МПК у детей с костной патологией, обусловленной различными формами рахита. Информация о разработке подобной программы представлена также в работе [8]. Однако в указанных работах сопоставление степени почернения различных участков изображения с соответствующим почернением ступеней тест-фантома осуществлялось, по сути дела, визуально, что снижает точность получаемых результатов.
Нами разработана программа — X-Ray Analysis, которая позволяет в автоматическом режиме определить значения МПК не только в конкретной точке, но и вдоль трассы, указанной оператором на снимке, а также в области произвольной формы, ограниченной кусочно-ломанной замкнутой линией. Первый вариант программы для использования в OS Windows 95/98 описан в [9]. В настоящее время более функциональная версия программы, разработанной в системе LabView для OSWindows 7, используется на кафедре педиатрии с курсом поликлинической педиатрии КГМА при проведении научных исследований.
Рисунок.
Результат работы программы при измерении МПК вдоль линейной трассы
Методика оценки МПК с помощью данной программы сводится к следующему. При получении рентгеновского изображения объекта тест-фантом располагают рядом с диагностируемым объектом. На рентгеновской пленке одновременно фиксируется изображение объекта и изображение тест-фантома. Поэтому все факторы, влияющие на денситометрические свойства пленки (в том числе и на степень почернения), одинаковым образом будут влиять и на изображение объекта, и на изображение тест-фантома. После стандартной процедуры химической обработки и сушки рентгеновская пленка сканируется, что позволяет получить компьютерное изображение диагностируемого объекта и тест-фантома. На первом этапе работы X-Ray Analysis осуществляется ввод изображения в программу и выполняется ее калибровка. Для этого оператор с помощью манипулятора мышь проводит на изображении тест-фантома калибровочную трассу. Яркости пикселов изображения вдоль трассы считываются программой, сравниваются с соответствующими эквивалентными значениями МПК тест-фантома для этой трассы и запоминаются. На следующем этапе работы — этапе измерений — оператор проводит линейную трассу в интересующей области, либо вычерчивает замкнутую ломаную, охватывающую эту область. Программа автоматически пересчитывает яркости пикселов в значения МПК и выдает на экран график этих значений вдоль трассы, либо среднюю величину МПК в области. Пользовательский интерфейс программы в режиме измерений проиллюстрирован на рисунке. В левой части рисунка показано цифровое рентгеновское изображение диагностируемого объекта и трасса (в данном случае прямая линия), проведенная оператором. В правой части представлен график значений МПК вдоль этой трассы. График вычерчивается автоматически по окончании проведения трассы. При использовании курсоров, включаемых оператором на поле графика, программой рассчитывается и величина кортико-диафизарного индекса, являющегося, как известно, относительным показателем толщины кортикального слоя.
С помощью данной программы нами было выполнено измерение МПК у 120 новорожденных детей гестационного возраста 32-40 недель (50 — доношенных, 70 — недоношенных новорожденных). Целью исследования было изучение МПК у новорожденных детей на основе компьютерного анализа пленочных рентгенограмм. Оценка МПК проводилась на уровне эпифиза и диафиза плечевой кости, диафиза ключицы. Данные области изучения МПК были выбраны в связи с наименьшим их «перекрытием» тенями других костных образований, тенями внутренних органов, что повышает точность измерения МПК.
Были выявлены отличия МПК между двумя группами обследованных новорожденных. В частности, уже предварительный анализ показал, что у недоношенных детей средние значения МПК в области эпифиза плечевой кости находятся в диапазоне от 0,5 до 0,9 г/см2 в зависимости от гестационного возраста, у доношенных детей диапазон составил от 1,0 до 1,4 г/см2. В области диафиза плечевой кости диапазон МПК составил от 0,8 до 1,1 г/см2 у недоношенных детей, от 1,1 до 1,5 г/см2 — у доношенных детей. В области диафиза ключицы МПК находилась в диапазоне от 0,8 до 1,0 г/см2 у недоношенных детей и от 1,0 до 1,3 г/см2 у доношенных.
В заключении следует отметить, что оценка МПК у детей продолжает оставаться актуальной проблемой и использование компьютерного анализа делает метод рентгеновской денситометрии более доступным для широкого использования даже при наличии рентгеновской аппаратуры общего назначения.
ЛИТЕРАТУРА
1. Практическое руководство по костной денситометрии у детей. Сокращенныйпереводматериалов National Osteoporosis Society. — Novembre 2009.
2. Казначеева Т.В., Осипова А.А. Современные методы определения минеральной плотности костной ткани // Проблемы репродукции. — 2007. — № 6. — С. 57-61
3. Баранова О.В., Малевич Э.Е. Современные возможности лучевых методов диагностики остеопороза // Медицинские новости. — 2009. — № 10. — С. 12-16.
4. Чиж Г.В., Полойко Ю.Ф. Оценка качественных характеристик медицинских рентгеновских пленок. Инструкция по применению // Министерство здравоохранения республики Беларусь. — 2003. — 26 с.
5. Гайдарова Т.А. и др. Метод использования рентгеноморфометрии для оценки минеральной плотности костной ткани альвеолярной кости // Бюллетень ВСНЦ СО РАМН. — 2006. — № 5 (51). — С. 46-49
6. Зиатдинова Н.В. Клинико-диагностическое значение определения цинка, хрома и стронция при различных формах рахита у детей: автореф. дис. … к.м.н. — Казань, 2001.
7. Мальцев С.В., Зиатдинова Н.В., Фофанов В.Б. Рентгеноденситометрия костной ткани у детей с различными формами рахита // Казанский медицинский журнал. — 2003. — Т. 84, № 1. — С. 41-42
8. Ломовцева О.В. Проектирование специализированной компьютерной системы определения минеральной плотности костной ткани по рентгеновскому снимку [Электронный ресурс] –– Режим доступа: http://masters.donntu.edu.ua/2006/kita/lomovceva/diss/index.htm
9. Воронов В.И., Воронова Т.В. Компьютерная программа обработки цифровых рентгеновских снимков/ Зарегистрирована ОФАП, ID 6635, N 50200601393 от 8.09.2006,
см. также «Инновации в науке и образовании, издание» ФГНУ «Гос. коор. центр информационных технологий». — 2006. — № 7 (18). — с. 29.
REFERENCES
1. Prakticheskoe rukovodstvo po kostnoy densitometrii u detey. Sokrashchennyy perevod materialov National Osteoporosis Society [A practical guide to bone densitometry in children. Abridged translation of materials National Osteoporosis Society]. 2009, November.
2. Kaznacheeva T.V., Osipova A.A. Modern methods for the determination of bone mineral density. Problemy reproduktsii, 2007, no. 6, pp. 57-61 (in Russ.).
3. Baranova O.V., Malevich E.E. Modern possibilities of radiological methods of diagnosing osteoporosis. Meditsinskie novosti, 2009, no. 10, pp. 12-16 (in Russ.).
4. Chizh G.V., Poloyko Yu.F. Otsenka kachestvennykh kharakteristik meditsinskikh rentgenovskikh plenok. Instruktsiya po primeneniyu [Evaluation of qualitative characteristics of medical X-ray films. Instructions for use]. Minsk: Ministerstvo zdravookhraneniya respubliki Belarus, 2003. 26 p.
5. Gaydarova T.A. et al. The method rentgenomorfometrii use to assess bone mineral density of alveolar bone. Byulleten’ VSNTs SO RAMN, 2006, no. 5 (51), pp. 46-49 (in Russ.).
6. Ziatdinova N.V. Kliniko-diagnosticheskoe znachenie opredeleniya tsinka, khroma i strontsiya pri razlichnykh formakh rakhita u detey: avtoref. dis. … k.m.n. [Clinical and diagnostic value of the determination of zinc, chromium and strontium in various forms of rickets in children. Synopsis of the dis. Of PhD med. sci.]. Kazan, 2001.
7. Mal’tsev S.V., Ziatdinova N.V., Fofanov V.B. Roentgen bone tissue in children with various forms of rickets. Kazanskiy meditsinskiy zhurnal, 2003, vol. 84, no. 1, pp. 41-42 (in Russ.).
8. Lomovtseva O.V. Proektirovanie spetsializirovannoy komp’yuternoy sistemy opredeleniya mineral’noy plotnosti kostnoy tkani po rentgenovskomu snimku [Designing specialized computer system determine bone mineral density by an X-ray], available at:
http://masters.donntu.edu.ua/2006/kita/lomovceva/diss/index.htm
9. Voronov V.I., Voronova T.V. Komp’yuternaya programma obrabotki tsifrovykh rentgenovskikh snimkov [Computer software processing digital radiographs]. Zaregistrirovana OFAP, ID 6635, № 50200601393 ot 8.09.2006, sm. takzhe “Innovatsii v nauke i obrazovanii”, izdanie FGNU “Gos. koor. tsentr informatsionnykh tekhnologiy”, 2006, no. 7 (18), p. 29.