МЕХАНИКА И МАШИНОСТРОЕНИЕ
УДК 519.711
RBF-АЛГОРИТМ И ЕГО МОДИФИКАЦИИ ДЛЯ ПОСТРОЕНИЯ ПОВЕРХНОСТНЫХ КОМПЬЮТЕРНЫХ 3-D МОДЕЛЕЙ В МЕДИЦИНСКОЙ ПРАКТИКЕ
© 2011 А.В. Колсанов, А.С. Воронин, Б.И. Яремин, С.С. Чаплыгин
Самарский государственный медицинский университет
Поступила в редакцию 02.10.2011
В работе выполнен сравнительный анализ характеристик модификаций метода RBF применительно к построению ЭБ-моделей органов по результатам 2D и ЭБ-сканирования. Описана схема построения модели на основе метода компактной RBF и приведены результаты моделирования на примере модели позвонка и головного мозга.
Ключевые слова: проекционные данные, поверхность, RBF, реконструкция, ЭБ-модель.
ВВЕДЕНИЕ
Во многих отраслях науки и техники существует практическая необходимость в построении трехмерных компьютерных моделей реальных объектов. Примеры таких задач - исследование органов в медицине, изучение внутренней структуры объектов при неразрушающем контроле, реконструктивная археология, криминалистика, создание систем виртуальной реальности - имитаторов и тренажеров. Исходными данными для реконструкции таких моделей является информация, полученная в результате сканирования реального объекта с использованием одной из известных технологий 2D или 3D-сканирования: рентгеновская компьютерная томография, магниторезонансная томография, ультразвуковое исследование, лазерное трехмерное сканирование и т.д. [1] . Преимущественно результат сканирования представлен в виде точечного описания объекта (поверхности и внутренней структуры). Результат реконструкции на основе вышеуказанного представления - это трехмерная модель в полигональном, поверхностном или ином представлении. Учитывая тот
Колсанов Александр Владимирович, доктор медицинских наук, профессор, заведующий кафедрой оперативной хирургии и клинической анатомии с курсом инновационных технологий E-mail: avkolsanov@mail.ru. Воронин Александр Сергеевич, ассистент кафедры оперативной хирургии и клинической анатомии с курсом инновационных технологий. E-mail: alek.voronin86@yandex.ru.
Яремин Борис Иванович, кандидат медицинских наук, старший преподаватель кафедры оперативной хирургии и клинической анатомии с курсом инновационных технологий. E-mail: boris.yaremin@gmail.com. Чаплыгин Сергей Сергеевич, ассистент кафедры оперативной хирургии и клинической анатомии с курсом инновационных технологий. E-mail: chaplyginss@mail.ru.
факт, что современные средства сканирования объектов предоставляют на выходе количество точек порядка (104 -106), то использование полигонального представления становится нерациональным. Поэтому исследователи, в первую очередь, уделяют внимание получению моделей на основе поверхностного представления. Тем более что поверхностное представление по оценкам экспертов дает более реалистичные модели и позволяет легко устранять такие дефекты входных данных как неполнота, повреждения и др [2].
ЦЕЛЬ ИССЛЕДОВАНИЯ
Для получения подобных моделей существует множество подходов. Задача данного исследования - провести сравнительный анализ существующих методов и оценить возможность их применения в составе программной системы для моделирования объектов в медицинской практике [3].
МАТЕРИАЛЫ И МЕТОДЫ ИССЛЕДОВАНИЯ
Построение модели объекта в трехмерном пространстве на основе проекционных данных
Особенностями моделей объектов непосредственно в медицинской практике являются [3]:
- значительное число проекционных точек (104 -106), что обусловлено особенностями получения проекционных данных, как совокупности точек контуров объектов на 2Б-проекциях;
- повышенные требования к точности модели, что связано с необходимостью
учитывать мелкие множественные артефакты;
- возможность построения модели в реальном времени, что связано с необходимостью учи-
тывать динамику моделируемого объекта (например, в системах визуализации хирургических операций и т.д. [4]).
Для построения вышеуказанных моделей используются сплайны, классическая триангуляция, а -поверхности и др. Авторами уже рассматривались методы моделирования на основе сплайнового представления [5], но как было определено в [6], наиболее эффективным средством моделирования поверхностей в трехмерном пространстве является метод ИВЕ
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
В данной работе ставится задача проанализировать возможность применения различных модификаций метода ИВР применительно к объектам, обладающим ранее перечисленными особенностями (большое число проекционных точек, высокая точность, реальное время).
С точки зрения метода ИВЕ объект (поверхность) представляется в виде:
где Л1 ^ - коэффициенты;
р - функция ИВР, вычисленная для произвольной точки пространства;
С1 N - заданные точки поверхности (проекционные данные).
Коэффициенты £ определяются как вектор-столбец, являющийся решением системы уравнений вида:
Таблица 1. Радиус функции RBF
Тип функции RBF
Piece wise polynomial (RP) <P(r) = \r и
Thin plate spline (TPS) (pir) = \r " In и
Multiquadric (MQ) <?(>•) - л/ l+(£r)2
Inverse miiltiqiiadric (IMQ) <p(r)= ——c
Inverse quadratic (IQ) 1 + (£r)
Gaussian (GS) (p{r) = <? -rV
(р{сг,сх) (р(с2,с2)
<p(cN,cl) <p(cN,c2) где
<P(C2'Cn)
9(cn,cn)
В методах получения 2Б-моделей, например, в методе моделирования рельефа по изолиниям [7] значения вектора-столбца H - это значения, которые принимает искомая модель в заданных точках множества
/'-!«• -i.v,r!;v,|: ¡Г
Для случая трехмерного пространства множество P принимает вид:
В этом случае модель объекта (поверхность) задается неявно в виде функции, принимающей некоторое оговоренное значение в каждой точке поверхности: F( . c) = const В качестве RBF можно использовать различные функции, значение которых определяется расстоянием между двумя точками пространства [8] ( см. табл. 1 ; г = x - c - радиус функции RBF).
Таким образом, процесс построения модели объекта заключается в следующем:
- расчет матрицы значений функции RBF для заданных точек;
- решение системы линейных уравнений для определения коэффициентов Д;
- расчет функции (1) для искомых точек пространства и выделение множества точек, для которых f (x) = const с учетом точности s .
Недостатком такого прямого метода RBF является значительная вычислительная сложность, которая связана с необходимостью решать систему большой размерности, а также необходимостью сохранения больших объемов данных. Обобщенные характеристики временной и пространственной сложности данного алгоритма показаны в табл. 2.
Уменьшение сложности вычислений, описанное, например, в работе [9], достигается за счет уменьшения количества проекционных точек, то есть отбрасывания точек, которые не влияют на результат. Для отбрасывания используется итерационный „жадный" алгоритм [9]. В результате работы алгоритма формируется множество точек, при использовании которых, ошибка в определении поверхности не превышает заданный параметр точности s .
В зависимости от выбранного параметра точности s и набора входных данных, будет наблюдаться разная степень уменьшения сложности вычислений. В соответствии с экспериментальными данными [9], степень уменьшения лежит в пределах 0.14..0.92, при этом любая закономерность отсутствует.
Общее уменьшение времени определяется тем, что при работе „жадного" алгоритма и незначительном результирующем количестве проекционных точек на каждом шаге решается система маленькой размерности. Таким образом данный алгоритм не является универсальным, так как не обеспечивает гарантированного уменьшения вычислительных затрат для произвольных объектов.
Усовершенствованный метод, предложенный в [10] предполагает уменьшение числа операций за счет использования так называемых "компактных" RBF. Такие RBF - это функции, обладающие следующим свойством: функция меньше или равна 1 в окрестности исходной точки и больше 1 при удалении от нее.
Вычисление данной функции до решения системы значительно уменьшает ее размер за счет отбрасывания точек удаленных от заданной и в малой мере влияющих на результат. В зависимости от пространственного расположения заданных точек достигается различная степень снижения вычислительной сложности, а преобладание в матрице системы нулевых элементов позволяет сохранять ее в более компактной форме. На этапе построения матрицы и вычисления поверхности удобной и эффективной структурой для организации хранения данных является k-дерево [11], а результирующая матрица явля-
ется сильно разреженной (см. рис. 1 справа) по сравнению с обычным методом ИБР (см. рис. 2 слева), что позволяет использовать ускоренные методы решения систем, например [12].
Таким образом, средство сокращения вычислительных затрат при построении моделей методом компактной ИБР - это применение эффективной древовидной структуры данных и использование методов ускоренного решения систем сразреженными матрицами.
Еще один известный метод - это метод быстрого вычисления множества потенциалов (РММ), который был адаптирован для решения данной задачи в работе [13]. Метод также позволяет значительно снизить временные и пространственные характеристики алгоритма за счет использования древовидных структур данных и аппроксимации модели с помощью, например, разложения в ряд Тейлора. Сущность метода быстрого вычисления состоит в следующем.
Пространство, в котором расположенные заданные точки, иерархически делится на элементарные подпространства. Такое распределение дает в результате дерево. Для любого подпространства на последнем уровне определяется два множества подпространств: приближенное и отдаленное. Приближенное содержит подпространства, которые расположены на том же уровне, что и искомое и имеют тот же родительский узел. Все другие подпространства составляют отдаленное множество. Процесс формирования модели поверхности с помощью функции вида (1) делится на два этапа: этап подготовки, который выполняется один раз, и этап вычисления, который повторяется для каждой искомой точки. На этапе вычисления часть функции (1), которая определяется точками из приближенных подпространств, вычисляется непосредственно, а часть, которая определяется точками в отдаленных подпространствах - с помощью аппроксимации. В зависимости от вида ИБр которая используется, применяются разные средства аппроксимации [13]. Однако данный метод не является универсальным, и выбор различных функций ИБР требует различных подходов к аппроксимации.
Рис. 1. Вид матрицы системы
Таблица 2. Анализ соответствующих алгоритмов
Этап обработки Алгоритм
Обычный КВЕ КвБс минимизацией Компактный КЕТ Быстрый ЫЗГ
Расчет ГШ Г 0{п2) 0{п2) 0(п 1о§ п) 0(л2)
Решение системы 0{п2) 0{п2) Щп1215) 0(П
Расчет поверхности 0{тп) 0(тп) 0(т 1о§ л) 0(т + п 1о§ и)
Теоретические оценки рассмотренных методов, полученные в результате анализа соответствующих алгоритмов и экспериментальных исследований на наборах тестовых объектов [14], приведены в табл. 2.
Анализируя данные табл. 2, можно сделать следующие выводы:
- на этапе расчета матрицы значений функции RBF наименьшую
вычислительную сложность имеет метод компактной RBF;
- на этапе решения системы наименьшую вычислительную сложность имеет
метод быстрой RBF;
- на этапе расчета поверхности наименьшую вычислительную сложность
имеет метод компактной RBF.
Таким образом, если считать процесс рекон-
струкции последовательным применением рассмотренных ранее этапов, то лучшие результаты с точки зрения минимального времени реконструкции при больших N должен показать метод компактной RBF. Поэтому целесообразно именно его использовать в составе программной системы для реконструкции моделей органов по результатам томографических исследований, структура которой предложена в [15].
ПРАКТИЧЕСКАЯ РЕАЛИЗАЦИЯ МЕТОДА КОМПАКТНОЙ RBF
Для осуществления трехмерного моделирования объектов предлагается следующая схема (см. рис. 2).
Предложенная схема предполагает слияние метода RBF и метода деформационного сопостав-
Рис. 2. Схема получения 3D-модели
Рис. 3. Результаты моделирования
ления трехмерных моделей на основе исследований, изложенных в [16]. Примеры полученных с помощью предложенной схемы и компактного метода RBF моделей показаны на рис. 3 (череп).
В табл. 3 сведены численные результаты экспериментов по построению моделей объектов, представленных на рис. 3. При моделировании использовалась система следующей конфигурации: Intel Core 2 Duo 3.0 ГГц, 2Гб RAM.
В табл. 3 использованы следующие обозначения:
t1 - время расчета матрицы или построения дерева;
t2 - время решения системы;
t3 - время расчета поверхности.
ts - суммарное время построения модели.
Как видно из табл. 3, суммарное время построения модели при количестве точек, достигающем 106...107, достаточно приемлемое, т.е. использование метода компактной RBF в данном случае позволяет приблизиться к реальному времени и получить лучший результат по сравнению, например, с работой.
С помощью предложенной методики ведутся работы по созданию, не имеющего аналогов в мире, аппаратно-программного комплекса "Виртуальный хирург" для 3D моделирования операционного процесса и учебно-методических
модулей для системного обучения врача-хирурга методикам открытой хирургии с небольшим размером операционного поля, методикам эндо-васкулярной хирургии и эндоскопической хирургии на этапах додипломного и последипломного образования. Работы осуществляются при финансовой поддержке Министерства образования и науки Российской Федерации.
ВЫВОДЫ
В результате анализа временных и пространственных характеристик методов построения трехмерных моделей, основанных на использовании ИБР и поверхностном представлении результата, выяснено, что при реконструкции моделей с большим числом точек (4Ы > 10 ) наибольший эффект, с точки зрения минимизации временных затрат, дает алгоритм „компактной ИБР". Данный алгоритм был реализован в составе программной системы для реконструкции трехмерных моделей органов. Результаты экспериментов подтверждают возможность использования метода компактной ИБР при построении 3Б-моделей объектов в реальном времени.
Ближайшей задачей для дальнейших исследований является исследование возможности адаптации метода построения 3Б-моделей на
Таблица 3. Численные результаты экспериментов
Объект Braiii Vertebra
Кол-во точек 10 10* 10J 10° 10' 10 104 10J 10° 10'
ti, с 0.12 0.34 1.05 4.24 24.6 0.16 0.45 1.24 5.48 35.7
t2, с 0.20 0,97 3.65 12.3 67.2 0.31 1.07 5.06 14.6 76.5
t3, с 0.15 0.45 1.34 4.96 26.5 0.1S 0.51 2.04 6.12 30.4
ts. с 0.47 1.76 6.04 21.50 118.3 0.65 2.04 8.34 26.2 142.6
основе RBF к параллельной реализации, что теоретически позволит сократить временные затраты на моделирование сложных динамически изменяющихся во времени объектов.
СПИСОК ЛИТЕРАТУРЫ
1. Farrell E.J. et al. Graphical 3D Medical Image Registration and Quantification // J-Med-Sys, 1997. №21 (3). P. 155-172.
2. Morse B.S. et al. Interpolating implicit surfaces from scattered surface data usingcompactly supported radial basis functions // SMI 2001 International Conference, May 2001. P. 89-98.
3. Бабков В.С. Реконструкция 3D-моделей органов в компьютерной томографии // Научно-практический журнал. "Проблемы моделиванния в автоматизации и проектировании динамических систем". Вып. 52. Донецк: ДонНТУ, 2002. C. 100-105.
4. Xu F, Mueller K. Real-time 3D computed tomographic reconstruction using commodity graphics hardware // Phys. Med. Biol. 52. 2007. P. 3405-3419.
5. Рамаков В. С. Методы обработки данных при томографических исследованиях// Научно-практический журнал Смоленского технического университета. "Информатика, кибернетика и вычислительная техника". Вып. 70. Смоленск: СмолТУ, 2003. р. С. 30-38.
6. DuchonJ. Splines minimizing rotation-invariant semi-norms in Sobolev spaces. In W. Schempp and K. Zeller editors, Constructive Theory of Functions of Several Variables. №571. Berlin: Springer-Verlag, 1977. P. 85-100.
7. PouderoxJ. Adaptive hierarchical RBF interpolation for creating smooth digital elevathion models // Proc. 12th ACM Int. Symp. Advances in Geographical information Systems 2004. ACP Press, 2004. P. 232-240.
8. Larsson E, Fornberg B. A Numerical Study of some Radial Basis Function based Solution Methods for Elliptic PDEs // Computers and Mathematics with Applications. - 2003. - №46. - P. 891-902.
9. Carr J.C. et al. Reconstruction and Representation of 3D Objects with Radial Basis Functions // ACM SIGGRAPH 2001, 12-17 August 2001. P. 67-76.
10. Kojekine N. et al. Software Tools Using CSRBFs for Processing Scattered Data // Computers & Graphics, April 2003, V. 27, №2, 2003. P. 311-319.
11. Bentley J.L. Multidimensional binary search trees used for associative searching // CACM. 1975. №18(9). P. 509-517.
12. Саух С.Е. Метод CR-факторизации неупорядоченных матриц // Сборник трудов международной конференции "Моделирование-2008", 14-16 мая 2008 г. Киев: Институт проблем моделирования в энергетике им. Г.Е. Пухова. 2008. Т. 1. С. 3-10.
13. Beatson R. K. et al. Fast fitting of radialbasis functions: Methods based on preconditioned GMRES iteration // Advances in Computational Math. 1999. №11. P. 253-270.
14. Level of Detail for 3D Graphics [Электронный ресурс], 2008. Режим доступа к ресурсу: http://lodbook.com/ models/ (дата обращения 12.09.2011).
15. Бабков В.С., Ивашкоец Е.В. Проектирование многофункциональной программной системы для реконструкции трехмерных объектов в медицинской практике // Сборник трудов Третьей международной научно-технической конференции молодых ученых и студентов "Информатика и компьютерные технологии" 11-13 декабря 2007 г. Донецк: ДонНТУ, Министерства образования и науки, 2007. С. 285-287.
16. Qiang W, Pan Z, Chun C, Jiajun B. Surface rendering for parallel slice of contours from medical imaging // Computing in science & engineering, January-February 2007, V. 9. №1. 2007. P. 32-37.
RBF-ALGORITHM AND ITS MODIFICATIONS TO THE CONSTRUCTION OF SURFACE 3-D COMPUTER MODELS IN MEDICAL PRACTIK
© 2011 A.V. Kolsanov, A.S. Voronin, B.I. Yaremin, S.S. Chaplygin Samara State Medical University
In work the comparative analyze of RBF-modification characteristic aplied for build 3D-model of organs by 2D and 3D-scan data is proposed. Outline of model build based on compactly supported RBF is described and modelling results for humain brain and vertebra are proposed. Key words: projection data, surface, RBF, reconstruction, 3D-model.
Alexander Kolsanov, Doctor of Medicine, Professor, Head at the Operative Surgery and Clinical Anatomy with a Course of Innovative Technologies Department. E-mail: avkolsanov@mail.ru.
Alexander Voronin, Assistant Lecturer at the Operative Surgery and Clinical Anatomy with a Course of Innovative Technologies Department. E-mail: alek.voronin86@yandex.ru. Boris Yaremin, Candidate of Medicine, Senior Lecturer at the Operative Surgery and Clinical Anatomy with a Course of Innovative Technologies Department. E-mail: boris.yaremin@gmail.com.
Sergey Chaplygin, Assistant Lecturer at the Operative Surgery and Clinical Anatomy with a Course of Innovative Technologies Department. E-mail: chaplyginss@mail.ru.