CC BY
39
алгоритм значительно превосходит оба этих алгоритма (рис. 2в).
Изолинии были построены для 20 уровней интенсивности из диапазона [min I (x, y) +10,
L x y
min I ( x, y) + 701. Эти значения подобраны в ходе
х y J
экспериментов.
Для реализации предложенного алгоритма было разработано интерактивное приложение на языке C#. Основными функциями приложения являются:
- загрузка и отображение изображений МРТ в формате DICOM;
- предложенный полуавтоматический алгоритм для выделения очагов поражения (отдельные окна для выделения прямоугольных регионов интереса, обучения классификатора, тестирования);
- дополнительные функции для предварительной обработки изображений (эквализация гистограммы, гамма-коррекция и др.).
На основании изложенного можно сделать следующие выводы. Предложенный подход к выделению областей интереса на изображениях основан на классификации изолиний, он может быть особенно полезен при обработке изображений, полученных с помощью медицинского оборудования. Разработанный алгоритм применен к задаче выделения очагов поражения на изображениях МРТ мозга. В работе проведен анализ существующих
алгоритмов для решения подобных задач, перечислены их недостатки. Показано значительное превосходство предложенного алгоритма над базовыми алгоритмами сегментации, основанными на выделении краев, разрастании регионов, кластеризации. Кроме того, проведено сравнение алгоритма с одним из новых полуавтоматических алгоритмов выделения очагов поражения [5], которое показало общее превосходство предложенного алгоритма почти на 10 %.
Авторы выражают благодарность А.В. Петряйки-ну, А.С. Крылову и Т.А. Ахадову за ценные консультации.
Литература
1. Zijdenbos A.P., Dawant B.M., Margolin R.A., Palmer A.C. Morphometric Analysis of White Matter Lesions in MR Images: Method and Validation, IEEE Transactions on Medical Imaging. 1994. Vol. 13, № 4, pp. 716-724.
2. Udupa J.K., Samarasekera L.W.S., Miki Y. [et al.]. Multiple Sclerosis Lesion Quantification Using Fuzzy-Connectedness Principles. IEEE Transactions on Medical Imaging. 1997. Vol. 16, № 5, pp. 598-609.
3. Scully M., Anderson B., Lane T. [et al.]. An Automated Method for Segmenting White Matter Lesions through Multi-level Morphometric Feature Classification with Application to Lupus, Frontiers in Human Neuroscience. 2010. Vol. 4, № 27.
4. Hillary F.G., Biswal B.B. Automated Detection and Quantification of Brain Lesions in Acute Traumatic Brain Injury Using MRI. Brain Imaging and Behavior. 2009. Vol. 3, pp. 111-122.
5. Xie Y., Tao X. White Matter Lesion Segmentation Using Machine Learning and Weakly Labeled MR Images. Proceedings of SPIE. 2011. Vol. 7962, pp. 79622G. URL: http://spiedigitallibra-ry.org/proceedings/resource/2/psisdg/7962/1/79622G_1 (дата обращения: 14.09.2011).
УДК 004.925.86
ПОСТРОЕНИЕ КРИВОЙ СКОЛЬЖЕНИЯ КОНИЧЕСКОГО ИНСТРУМЕНТА ПРИ МНОГОКООРДИНАТНОЙ ОБРАБОТКЕ
А.И. Будник; Е.И. Кац, к.т.н.
(Уральский федеральный университет им. Б.Н. Ельцина, г. Екатеринбург,
Ь^тка1еха^г@таИ ты)
Рассматривается задача построения кривой скольжения на поверхности конического инструмента в определенный момент времени. Разработан алгоритм построения кривой скольжения на поверхности конического инструмента. Моделирование такой кривой необходимо, в первую очередь, для моделирования поверхности (объема), заметаемой инструментом в процессе многокоординатной фрезерной обработки.
Ключевые слова: многокоординатная обработка, геометрическое моделирование, кривая скольжения, граничные точки, конический инструмент.
Решение задачи построения кривой скольжения на поверхности конического инструмента в определенный момент времени необходимо, в первую очередь, для моделирования поверхности (объема), заметаемой инструментом в процессе многокоординатной фрезерной обработки. Эта задача, в свою очередь, актуальна при разработке САМ-систем и систем моделирования обработки.
Для решения рассматриваемой задачи рядом исследователей предложены различные методы
(см. например [1-3]). Настоящая статья посвящена разработке подхода, в наибольшей степени соответствующего требованиям имеющегося пакета геометрического моделирования фрезерной обработки на станках с ЧПУ Получено решение, обеспечивающее приемлемое быстродействие при требуемой точности построения.
Особенностью предлагаемого подхода является использование линейности изменения станочных координат положения инструмента при обработке.
В любой момент времени вектор мгновенной скорости в пространстве изменяется линейно от одной точки к другой, что позволяет существенно упростить поиск точек на кривой скольжения. Кроме того, удалось построить алгоритм поиска таким образом, чтобы отсекаемая часть аппроксимированной поверхности инструмента оставалась выпуклой.
Постановка задачи: пусть на входе имеются контур конического инструмента, полная производная матрицы скорости по времени, матрица инструмента. Изменение всех координат считается линейным по времени [4]. Необходимо построить кривую скольжения на поверхности инструмента в данный момент времени.
При рассмотрении задачи контур инструмента задан в плоскости xOz в положительной области по оси Ox. Инструмент может состоять из нескольких простых инструментов. Сам контур инструмента изначально не имеет экстремумов вдоль оси Ox. Контур непрерывен. Он начинается и оканчивается на оси Oz.
Для построения кривой скольжения необходимо дать ее определение. Кривая скольжения - это замкнутый контур на поверхности инструмента, в каждой точке которого скалярное произведение вектора внешней нормали к поверхности и вектора мгновенной скорости равно нулю [1]. Такое скалярное произведение в точках на поверхности режущей части инструмента является положительным (назовем эту область положительной). В результате задача построения кривой скольжения переходит в задачу построения кривой, ограничивающей положительную область поверхности инструмента. Для решения задачи необходимо научиться находить точки кривой скольжения.
Как известно, при движении сферической поверхности точки, лежащие на окружности большого радиуса в плоскости, перпендикулярной поступательной составляющей вектора скорости в центре сферы, будут образовывать кривую скольжения [5]. Назовем их граничными точками (grazing points). Поскольку инструмент, как и его контур, имеет вертикальную ориентацию относительно своей оси и существует понятие режущей стороны инструмента, есть смысл говорить о левой и правой сторонах кривой скольжения в направлении движения инструмента. Для построения кривой скольжения возьмем окружность с центром на оси инструмента, перпендикулярную этой оси, то есть набор точек на поверхности инструмента, расположенных на определенной высоте. Так как вектор скорости d в граничной точке должен быть перпендикулярным вектору внешней нормали n,
пара граничных точек слева и справа должна располагаться в плоскости, перпендикулярной вектору d, а вектор п лежать в этой плоскости. Из рисунка 1 видно, что точки М0 и Ы\ лежат на кривой скольжения, так как удовлетворяют ее определению. При одной точке пересечения плоскости, перпендикулярной вектору d и окружности, взятой на некоторой высоте, можно говорить о петле кривой скольжения. Тогда кривая скольжения либо заканчивается в этой точке по причине замкнутости, либо будет состоять из нескольких компонент связности, каждая из которых удовлетворяет определению, данному выше.
Рассмотрим случай с конической поверхностью (рис. 1). Для нахождения граничных точек на некоторой высоте по оси Oz возьмем окружность с центром Е на оси конуса, лежащую на его боковой поверхности и одновременно в плоскости, перпендикулярной оси Oz. Точки М0 и М\ имеют ту же скорость, что и точка Q. В этих точках внешние нормали перпендикулярны вектору скорости d, а значит, являются граничными точками.
O
Рис. 1
O
Рис. 2
При варианте с конусом наличие лишь одной граничной точки на определенной высоте по Oz возможно в двух случаях (рис. 2). Для этого необходимо, чтобы вектор скорости d в точке Q центра сферы имел тот же угол наклона к оси Oz, что и контур конуса. Тогда получаются две петли на конусе, причем будет справедливо следующее утверждение:
d ■ k
-q- = cos q , (1)
k
где d - вектор мгновенной скорости в точке Q;
k - орт оси Oz; q - угол полураствора контура. Так как скорость вычисляется по формуле d (t) = P(t) = P0M (t), (2)
где d (t) - вектор мгновенной скорости в момент
времени t; P (t) - производная закона движения точки по t; P0 - точка в системе координат инструмента; M (T) - полная производная матрицы ско-
z
рости по времени, то можно утверждать, что вектор скорости линейно зависит от координат точки.
Представим в параметрическом виде вектор скорости d через векторы скорости и0 и и1:
(( и 0 + ( и1 -и 0 ) •{ ) •к ) ( и 0 + ( и1 -и 0 ) •{ )2
= cos q
(3)
где и0, и1 - векторы скоростей в точках Q0 и Q1; t - параметр вдоль оси инструмента (рис. 3). Корнями уравнения
т =
\б - бо\ б - бо
(4)
где Q, Q0, Ql - центры сфер, касательных к поверхности инструмента, служат значения параметров t, используя которые, можно найти точки Q для вершин петель, образованных кривой скольжения.
Используя рисунок 3, получим следующее выражение для вычисления координаты г экстремумов петель:
д) •( г, -2р) +(г„ +Ц, •tgq) -(Ц, •tgq-z0 •tg2 д) " (^ д)
(5)
где г0, х1 - координаты г вершин А и В контура конуса; Я0 - радиус нижнего основания по оси Oz.
На элементе инструмента может располагаться одна петля в случае выхода одного из найденных г за пределы этого элемента или при вырождении уравнения (3) в линейное. У цилиндра, очевидно, С0Б2д=1 и, следовательно, совпадут точки Q и Е. В итоге возможны два варианта: либо плоскость, перпендикулярная d, пересекает выбранную на определенной высоте окружность в двух точках, либо содержит ее целиком. Для построения кривой скольжения необходимо найти граничные точки, лежащие на ней.
Рассмотрим несколько способов нахождения граничных точек на кривой скольжения с последующим ее построением.
Первый способ - нахождение граничных точек для конкретных значений г, взятых с некоторым фиксированным шагом по высоте инструмента. Отклонение аппроксимирующей кривой от истинной кривой скольжения при этом не учитывается. Количество точек определено лишь количеством шагов вдоль оси инструмента и видом кривой.
Второй способ - поиск граничных точек с учетом отклонения аппроксимирующей ломаной от истинной кривой скольжения. Для каждой найденной точки аппроксимирующей кривой в направлении касательного вектора к кривой подбирается такой шаг, чтобы отклонение касательного вектора от истинной кривой скольжения не превышало допустимого значения. Для полученного значения г конца вектора ищется следующая граничная точка. В первом и втором способах используется метод касательной сферы [5].
Третий способ основан на поиске граничных точек с учетом линейности изменения вектора мгновенной скорости при движении по отрезку. Таким образом можно находить граничные точки на прямых, образующих боковую поверхность конического инструмента (рис. 4), используя определение кривой скольжения. На концах образующих вычисляются векторы мгновенных скоростей. Из линейного уравнения находится точка на образующей, в которой скалярное произведение нормали к поверхности и вектора мгновенной скорости равно нулю. При этом учитывается отклонение при построении образующих конуса. Количество точек ограничено габаритами инструмента, значением максимально допустимого отклонения и геометрией самой кривой скольжения. Для большей точности образующие подбираются так, чтобы точно найти экстремум петли граничной кривой.
Такой способ позволяет просто решить еще одну важную задачу. Нахождение кривой скольжения нужно, в первую очередь, для построения заметаемой поверхности, частью которой будет часть поверхности инструмента, ограниченная кривой скольжения. При аппроксимации этой поверхности многогранником (сеткой из треугольников) важно сохранить выпуклость. Если находить граничные точки одним из первых двух способов, а затем просто соединять их с точками на направляющей, то выпуклость, вообще говоря, будет нарушена. Третий же способ лишен этого недостатка, так как, если соединять найденные точки с концами соответствующих образующих, выпуклость гарантирована.
На основании полученных результатов разработаны программы моделирования кривой скольжения
на поверхности конического инструмента при многокоординатной обработке. Программный модуль добавлен в состав симулятора обработки на станках с ЧПУ Полученные результаты служат основой для построения поверхности, заметаемой инструментом при многокоординатной обработке.
Литература
1. Wang W.P., Wang K.K. Geometric modeling for swept volume of moving solids. IEEE Computer Graphics and Applications. 1986. № 6 (12), pp. 8-17.
2. Abdel-Malek K., Seaman W., Yeh H.-J. NC-verification of up to 5 axis machining processes using manifold stratification. ASME Journal of Manufacturing Science and Engineering. 2000. № 122, pp. 1-11.
3. Weinert K., Du S.-J., Damm P., Stautner M. Swept volume generation for the simulation of machining processes. International Journal of Machine Tools and Manufacture. 2004. № 44 (6), pp. 617-28.
4. Будник А.И., Кац Е.И. Построение траектории движения инструмента при многокоординатной обработке // Программные продукты и системы. 2009. № 3.
5. Eyyup Aras. Generating cutter swept envelopes in five-axis milling by two-parameter families of spheres. Computer-Aided Design. 2009. № 41, pp. 95-105.
УДК 004.925.86
КРИВАЯ СКОЛЬЖЕНИЯ НА ИНСТРУМЕНТЕ ПРОИЗВОЛЬНОЙ ФОРМЫ ПРИ МНОГОКООРДИНАТНОЙ ОБРАБОТКЕ
А.И. Будник (Уральский федеральный университет им. Б.Н. Ельцина, г. Екатеринбург,
Ь^тка1еха^г@таИ ги)
Построение кривой скольжения при многокоординатной фрезерной обработке аналитическим методом позволяет рассчитывать на более точное представление обработанной поверхности. В статье рассматриваются все возможные варианты форм кривой в зависимости от формы инструмента. На основе алгоритма построения создан программный модуль в составе симулятора NCManager.
Ключевые слова: многокоординатная обработка, геометрическое моделирование, кривая скольжения, граничные точки, заметаемая поверхность.
Для построения кривой скольжения существуют различные методы, предложенные, например, в работах [1, 2]. Решение задачи в дальнейшем необходимо для построения поверхности, заметенной инструментом при геометрическом моделировании фрезерной обработки на станках с ЧПУ В данной статье рассматриваются особенности построения кривой скольжения на поверхности произвольного по форме инструмента.
Основная постановка задачи: на входе имеются контур инструмента, полная производная матрицы скорости по времени, матрица инструмента. Необходимо построить кривую скольжения на инструменте при движении в данный момент.
При рассмотрении задачи контур инструмента задан в плоскости xOz в положительной области по оси Ох. Инструмент может состоять из нескольких простых инструментов, контур каждого из которых представляет набор определенных элементов. В качестве элементов контура могут быть дуга окружности и прямая. В результате инструмент как тело вращения может представлять собой совокупность нескольких поверхностей вращения (рис. 1), таких как конус (цилиндр - как частный случай), сфера, тор, плоскость.
Кривой скольжения у сферы всегда будет являться окружность. Из постановки задачи понятно, что сфера может быть представлена максимум четвертью своей поверхности. Значит, кривая
скольжения будет выглядеть как дуга, не превышающая четверти окружности.
Рассмотрим применение способа нахождения точек кривой скольжения методом касательной сферы (рис. 2) для таких видов элементов инструмента, как тор, сфера и плоскость.
Рис. 1
Рис. 2
На инструменте фрагмент поверхности тора может быть выпуклым как внутрь, так и наружу. Рассмотрим оба варианта расположения контура элемента инструмента. На рисунке 3 а изображен контур выпуклого наружу фрагмента тора.
Для поиска центра ( сферы, касающейся выпуклого тора в точках на высоте точки Е, используется следующее выражение, полученное по рисунку 3 а:
