Томография плазменно-пылевых структур тема автореферата и диссертации по физике, 01.04.04 ВАК РФ

УДК 533 9 082 5

На правах рукописи

Бульба Артем Владимирович

ТОМОГРАФИЯ ПЛАЗМЕННО-ПЫЛЕВЫХ СТРУКТУР

Специальность 01 04 04 - физическая электроника

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Петрозаводск 2007

003158756

Работа выполнена на кафедре информационно-измерительных систем и физической электроники физико-технического факультета Петрозаводского государственного университета

Научный руководитель

Официальные оппоненты

Ведущая организация

доктор физ -мат наук, профессор Хахаев Анатолий Диамидович

доктор физ -мат наук, профессор Герасимов Геннадий Николаевич

кандидат физ -мат наук, доцент Карасев Виктор Юрьевич

Объединенный институт высоких температур РАН

Защита состоится 007 г на заседании диссертационного со-

вета К 212 190 01 при Петрозаводском государственном университете по адресу г Петрозаводск, пр Ленина, д 33

Отзывы направлять по адресу

185910, г Петрозаводск, пр Ленина 33, ПетрГУ, ученому секретарю диссертационного совета К 212 190 01

С диссертацией можно ознакомиться в библиотеке Петрозаводского государственного университета

Автореферат разослан 2007 г

Ученый секретарь диссертационного совета К 212 190 01

доктор физ -мат наук, профессор г ,

Фофанов Анатолий Дмитриевич Чг^/Г

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы

При определенных условиях в плазме, в частности, в тлеющем разряде постоянного тока в газах из введенных в плазму макрочастиц размером от нескольких до сотен микрометров возникают самоупорядоченные плазменно-пылевые структуры В случае высокой степени упорядоченности структуры называют "плазменным кристаллом"

При диагностике пылевых образований особую роль играет выяснение геометрических параметров формируемой в плазме структуры пылевых частиц Традиционно пространственная диагностика плазменно-пылевых структур сводится к анализу их изображений, полученных в лазерном "ноже", т е в каждый момент времени регистрируется одно горизонтальное или вертикальное изображение структуры, не отражающее распределение частиц в пространстве, если речь не касается исследования монослойных плазменных образований Достоверность получаемых таким способом результатов зависит от толщины "ножа" и характера структуры, в некоторых случаях возможны значительные ошибки как в оценке среднего расстояния между частицами, так и в виде парной корреляционной функции, характеризующей степень упорядоченности структуры В некоторых экспериментах плазменно-пылевые объекты наблюдались путем горизонтального сканирования вертикальным лазерным "ножом" При этом осуществлялась запись положения каждой частицы с помощью телекамеры Полная длительность такого сканирования составляла более десятка секунд, что совершенно неприемлемо в случае осцилляции частиц Использование же классической стереоскопической системы машинного зрения позволяет осуществлять реконструкцию только приблизительно, поскольку остается не решенной задача сопряжения характерных точек изображений

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

Цель работы, задачи

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

позволяющий получать полную и достоверную информацию о распределении частиц в плазменно-пылевых образованиях

Для достижения этой цели необходимо решить ряд задач

1 Проанализировать существующие оптические методы пространственной диагностики плазменно-пылевых структур

2 Адаптировать метод оптической томографии применительно к задаче пространственной диагностики пылевых структур

3 Сформулировать ограничения и условия применимости метода

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

5 Разработать необходимое ПО Проверить работоспособность и эффективность разработанного алгоритма с использованием модельных и экспериментальных проекционных данных

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

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

Научная новизна

1 Развит метод оптической томографии применительно к задаче пространственной диагностики плазменно-пылевых образований

2 Сформулированы ограничения и условия применимости развитого метода

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

4 Впервые предложенным методом в тлеющем разряде постоянного тока проведено исследование зависимости межчастичных расстояний трехмерной пылевой структуры от условий разряда

Научно-практическая значимость работы

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

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

Изучение плазменных кластеров (как двумерных, так и трехмерных) представляет собой отдельную, очень интересную, развивающуюся область исследований Здесь поставлено множество вопросов, в решение которых представленная методика может сыграть существенную роль, как инструмент контроля геометрических параметров

- одной из наиболее интересных проблем является образование центров кристаллизации, переход от "пылевой" молекулы к молекулярному кластеру и, затем, к твердому состоянию

- проблема непрерывного перехода малых пылевых кластеров с небольшим количеством пылинок к трехмерным или двумерным "кристаллам" с многими пылинками

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

- исследование двойных смесей (кластеров с пылинками различного размера) и неоднородных кластеров

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

На защиту выносятся:

1 Метод оптической томографии, адаптированный применительно к задаче пространственной диагностики плазменно-пылевых структур, для которого сформулированы ограничения и условия применимости

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

3 Набор программных модулей "Tomography 2", основанный на предложенном алгоритме, в реализации которых учитывалось влияние

используемых оптических систем регистрации и зашумленность проекционных данных

4 Зависимости межчастичных расстояний трехмерной пылевой структуры из полидисперсного порошка с диаметром частиц от нескольких единиц до нескольких десятков микрон от токов и давлений в тлеющем разряде постоянного тока в Аг

Апробация работы

Содержание работы докладывалось на Международной конференции Оптических Обществ Дании, Финляндии, Норвегии и Швеции «Северная оптика 2003» (г Эспоо, Финляндия, 16-18 июня 2003 г), Международной конференции «Образование и обучение в оптике и фотонике (ЕТОР)» (г Туссон, USA, штат Аризона, 5-9 октября 2003 г), IV международном симпозиуме по теоретической и прикладной плаз-мохимии ISTAPC-2005 (г Иваново, Россия, 13-18 мая 2005 г), V международной конференции «Физика плазмы и плазменные технологии» РРРТ-5 (г Минск, Белоруссия, 18-22 сентября 2006 г), IV Курчатовской молодежной научной школе (г Москва, 20-22 ноября 2006 г), Всероссийской (с международным участием) конференции "Физика низкотемпературной плазмы - 2007" (г Петрозаводск, 24-28 июня 2007 г)

Успешно пройдена процедура государственной регистрации в «Национальном информационном фонде неопубликованных документов» разработки, предъявленной в отраслевой фонд алгоритмов и программ программа для определения пространственного расположения набора точечных объектов методом оптической реконструктивной томографии "Tomography 2" (номер государственной регистрации 50200701083 от 28 05 2007) Отраслевым фондом алгоритмов и программ выдано свидетельство об отраслевой регистрации разработки № 8356

Федеральной службой по интеллектуальной собственности, патентам и товарным знакам (РОСПАТЕНТ) выдано свидетельство об официальной регистрации программы для ЭВМ "Tomography 2" (объекту регистрации присвоен регистрационный № 2007612629 от 20 06 2007)

Пройдена процедура государственной регистрации в «Национальном информационном фонде неопубликованных документов» разработки, предъявленной в отраслевой фонд алгоритмов и программ программа для обучения алгоритмам реконструктивной томографии "Tomography" (номер государственной регистрации 50200700997 от 15 05 2007) Отраслевым фондом алгоритмов и программ выдано свидетельство об отраслевой регистрации разработки №8286

Федеральной службой по интеллектуальной собственности, патентам и товарным знакам (РОСПАТЕНТ) выдано свидетельство об официальной регистрации программы для ЭВМ "Tomography" (объекту регистрации присвоен регистрационный № 2007611674 от 20 04 2007)

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

Список работ, опубликованных по результатам диссертации, приведен в конце автореферата

Структура и объем диссертации

Диссертация состоит из введения, пяти глав, заключения и приложения Работа содержит 120 страниц, 74 рисунка, 6 таблиц Список использованной литературы включает 54 наименования

Благодарности

Работа выполнена, благодаря финансовой поддержке по гранту RUX0-000013-PZ-06 (НОЦ «Плазма») Американского фонда гражданских исследований и развития, а также Министерства образования и науки РФ и правительства Республики Карелии

СОДЕРЖАНИЕ РАБОТЫ

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

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

Впервые трехмерная плазменно-пылевая структура была исследована в [18] с использованием горизонтально и вертикально ориентированных лазерных ножей толщиной 90 мкм Эксперимент проводился в ВЧ - разряде в криптоне при давлении 1 4 Тор В разряд инжектировались почти монодисперсные (9 4±0 3 мкм) сферические полимерные пылинки Объемная структура наблюдается путем сканирования лазерным "ножом" При этом осуществлялась запись положения каждой частицы по "х", "у" и "г" координатам с помощью цифровой телекамеры, а данные представляли собой набор плоских горизонтальных изображений, по которым вычислялось объемное представление Эта методика применялась при исследовании пылевых структур в [19-24], однако ее временное разрешение пока все еще недостаточно Так, в [23] скорость развертки лазерного пучка и телекамеры была равна 0 3 мм/с, а полная длительность развертки - 18 3с, что совершенно неприемлемо в случае осцилляции частиц и, например, для трехмерных наблюдений фазовых переходов

В работе [25] предлагается использовать стереоскопический подход в решении задачи трехмерного восстановления плазменного "кристалла"

Классическое стереозрение включает два процесса совмещение деталей, наблюдаемых двумя камерами и восстановление их трехмерного прообраза [26, 27] Последний процесс относительно прост прообраз соответствующих точек можно найти как точку пересечения лучей, проходящих через эти точки и центры соответствующих диафрагм камер [28] Следовательно, если характерную точку изображения можно наблюдать и однозначно идентифицировать на всех камерах в любой момент времени, стереозрение реализуется просто В то же время, каждая сцена состоит из сотен тысяч пикселей (или десятков и сотен характерных объектов), так что для определения точных соответствий между точками различных изображений сцены требуются некоторые методы, исключающие ошибочные измерения глубины Таким образом, центральной проблемой компьютерного стереозрения является поиск сопряженных точек

Традиционно [28] задача поиска сопряженных точек решается двумя известными способами

1 учетом эпиполярного условия,

2 корреляционным методом нахождения пиксельных соответствий путем сравнения профилей яркости [27]

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

приблизительной, т к метод позволяет только несколько сократить размер зоны поиска

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

В [29] частицы освещены двумя параллельными световыми лучами лазеров (686 и 656 нм), чьи интенсивности предварительно были про-модулированы по глубине кристалла Рассеянный свет регистрируется двумя камерами, расположенными под углом в 68° относительно друг друга В этом случае сопряженные точки находятся корреляционным методом

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

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

Во второй главе рассмотрены принципы и некоторые алгоритмы томографии

Суть всех типов томографии едина по интегральной информации (проекциям) полученной от некоторого слоя вещества, нужно определить локальную информацию А(х,у) в том же сечении [30], а затем по сечениям восстановить объем ^х,у,г), где "г" - координата, перпендикулярная слоям

Современную томографию можно без преувеличения считать всеволновой В зависимости от характера конкретной задачи восстановление структуры объекта может основываться на регистрации пучков электронов, ионов, нейтронов, фотонов во всем диапазоне электромагнитного спектра, звуковых волн [31, 32] Исследуемый объект при этом облучается извне или сам является излучателем В различных типах томографии суммарная информация качественно различна, однако есть

тенденция описывать все виды томографии с точки зрения единого подхода, т к конечная цель одна [33]

Рассмотрим простейший случай. Пусть двухмерное сечение исследуемого объекта "1" (рис 1) описывается функцией плотности f(x,y) Вдоль подвижной оси "v" сканируем слой объекта просвечивающим излучением "4", используя излучатель "3" (рис 1), и фиксируем ослабление луча детектором "2", получая проекцию при заданном угле в Следует отметить, что объект может излучать и сам (в этом случае излучатель не потребуется)

Проекция представляет собой отображение N-мерной функции в (N-1 )-мерную функцию, получаемое путем ее интегрирования в заданном направлении В рассматриваемом двухмерном случае проекция рв(у) функции "f" будет определена как

pe(v)= \f(x,y)du, (1)

Итак, задача восстановления заключается в решении интегрального уравнения (1) Самый естественный и мощный способ дает преобразование Фурье Воспользовавшись основными свойствами преобразования, а также теоремой о проекциях и сечениях [32] получаем

Л*. У) = ±? )р-'{Р{рв(у)}\У\№, (2)

о

71

Если в (2) принять = 1, тогда /(х, у) = У)(10, (3)

о

(3) представляет собой метод "обратного проецирования", позволяющий оценивать плотность в любой точке сечения посредством сло-

Y

/ \ 1 ч V

\4 %) X

\ ^-Дг v>

л

Рис 1 Схема сканирования слоя

жения лучевых сумм рв(у') для всех лучей, проходящих через искомую точку

Предположим, в объеме восстановления есть набор светящихся точечных объектов "А", "В" и пусть известны их двухмерные проекции "Р1", "Р2", "РЗ" при трех различных углах (рис 2)

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

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

В третьей главе описан и протестирован предлагаемый алгоритм Чтобы применить томографические методы к пылевой плазме, мы должны осветить анализируемую область плазмы с рассеивающими частицами (Б1, БЗ) расширенным лазерным пучком, установить фоторегистрирукмцие системы в горизонтальной плоскости, и зарегистрировать несколько двумерных изображений (проекций Р1, Р2, РЗ) этой области под различными углами наблюдения , (Рис 3)

X

Рис 3 Схема томографирования 01, 02, ОЗ - оптические оси

При использовании оптической системы с большим увеличением можно применить видоизмененный алгоритм обратной проекции Трехмерный массив "V", соответствующий всем-элементам области восстановления структуры, первоначально обнуляется Рассчитываются координаты точек, через которые пройдет прямая, восстановленная из каждой точки плоскости проекции ("Р1", "Р2", "РЗ", рис 3) и проходящая через центр оптической системы, и ко всем элементам массива области восстановления, через которые проходит эта прямая, добавляется "1" В результате за координаты реальных частиц принимаются те элементы массива, в которых появилось значение "3" (при трех проекциях) Однако, в реальном эксперименте фоторегистрирующие системы не откалиброваны, главные оптические оси не находятся в одной плоскости (оптические системы сдвинуты в пространстве на неизвестные "8Ь1Й_У", "БЫй г" - рис 3), межугловые расстояния известны приблизительно (в диапазоне 0±А0), а следы частиц на проекциях являются не точками, а пятнами из-за ограниченной глубины резкости оптической системы

Остается возможность варьировать неизвестные значения в некотором диапазоне и осуществлять проверку точности восстановления, например, вычисляя проекции с восстановленного объекта и накладывая вычисленные проекции на экспериментально полученные Остановить такой итерационный процесс можно при достижении некоего ми-

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

1 воспользоваться априорной информацией об исследуемом объекте,

2 отказаться от принципа заполнения и анализа трехмерного массива "V"

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

Отказаться от заполнения и анализа трехмерного массива "V" возможно, решив задачу аналитически Лучи, идущие от каждого следа на проекции через объем восстановления, будем задавать пересечением двух плоскостей При поиске точки пересечения лучей сделаем допуск на возможность некоторого расхождения лучей (будем рассматривать объем, в котором могут пересекаться искомые прямые), поскольку в реальном эксперименте следы частиц на проекциях являются не точками, а пятнами и треками Из рисунка 4 видно, что прямая "А1К1", задаваемая плоскостями "А1А5К1" и "А1А4К1" (где "Ь" - малое произвольное смещение относительно следа "А1"), проходит через след на проекции "Р1" и оптический центр линзы "К1"

Для более компактного описания алгоритма введем следующие операторы ¥{} - итерационный оператор, С{} - оператор априорных ограничений,

М{} - оператор вычисления координат объема пересечения, Н{} -оператор вычисления проекций с восстановленного объекта,

Тогда итерационный алгоритм можно представить системой (4,5) Задаем начальные приближения (выбираем начальную точку у стенки газоразрядной трубки на прямой "В1К2", рис 3) и осуществляем восстановление на приближении "w" с вычислением проекций (4) Далее проводим проверку точности восстановления (6) и, если отклонение меньше предыдущего (итерация "w-1"), запоминаем параметры ("shift", "AQ" для каждой проекции) С учетом априорных ограничений продолжаем варьировать параметры с целью минимизации отклонения (6) и вычисления координат частиц плазменно-пылевой структуры (5), т е будем программно изменять координаты фоторегистрирующих систем "PI", "РЗ" (при неизменном расстоянии между трубкой и фотоматрицей "L") посредством одновременного смещения прямых "А1К1" и "С1КЗ" вдоль прямой, идущей из следа "В1" частицы "D1"

Н

м

mm

]Гдг, , ]Гдг, , ,

. '=i Л V i=i )г v '=i

рл

(wf

де,

shif т

Н М'

(4)

М< mm

F

лв

shift

II f

Щрч-*

= х„

(5)

где ^Аг, — сумма расстояний между точками множества, образую-

щих объем пересечения, п=6, "Р1" - следов на проекции "Р1", "р3" -следов на проекции "РЗ", "т" — порядковый номер следа на проекции "Р2" и восстановленной частицы, к е \ ..рхръ, ^^ = (х,у)]

к т

- координаты следа г на проекции у, ' - указывает на приближение, = (и,у?У™* - вычисленная проекция с координат

x[y'J при заданных AQ и смещении "shift", "shift" - смещение фоторе-

гистрирующих систем, (w) - w-я итерация

В результате выполнения итерационного процесса будет найдено оптимальное пространственное расположение частиц восстанавливаемого объекта с минимальным отклонением "R" (6) Варьирование углов в диапазоне " ± AQ " для каждой проекции можно реализовать как

дополнительный процесс восстановления с проверкой по уравнению (6) после каждого малого смещения по прямой "В1К2" (рис 3)

* = (6)

j=i i=i

где "р"- число проекций, "N" - количество частиц, vt =(u,v)t} - координаты следа "/" на проекции "j ", V(} = (u,v] - вычисленная проекция Фактически "R" показывает среднюю ошибку позиционирования каждой восстановленной частицы

Далее в этой главе на математических моделях проиллюстрирована работоспособность и эффективность предложенного алгоритма Рассмотрены ограничения, накладываемые на метод, вследствие применения оптических систем с большим увеличением

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

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

и 1 [34], с целью выделения множеств смежных объектов и координат их центров, необходимых в предложенном алгоритме томографии

Однако в экспериментах реальный сигнал на выходе прибора содержит кроме полезного сигнала еще и шумы

Так, в нашем случае шум возникает из-за недостатка освещения и использования лазера (шум Эрланга), что приводит к появлению несуществующих деталей или разорванности объектов

Хотя человеческое зрение без труда ликвидирует образовавшиеся пробелы в разрывах и несуществующие детали, автоматические системы сталкиваются с серьезными трудностями при распознавании

Поскольку исследуемый объект освещается лазером (532 нм), из полученных изображений можно сразу отсечь два канала, отвечающих за красный и синий цвета

Известно, что резкие перепады яркости на изображении (например, связанные с шумом) вносят значительный вклад в высокочастотную часть его Фурье-преобразования Следовательно, сглаживание достигается в частотной области ослаблением высокочастотных компонент определенного диапазона Фурье-образа данного изображения [34], что позволяет исключить или существенно уменьшить яркость деталей, имеющих размеры меньше интересующих, а также ликвидировать образовавшиеся пробелы в разрывах В качестве же фильтра в работе применяется гауссов фильтр низких частот, отражающий наши априорные предположения о характере искомого сигнала

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

Описанный подход предобработки экспериментальных данных реализован в виде модуля 'ТййаЬоп" в составе программы "Тото-§гар1у 2"

В пятой главе описаны условия эксперимента и результаты применения разработанного алгоритма и программных средств к задаче пространственной диагностики плазменно-пылевых образований методом оптической томографии

Экспериментальная установка представляет собой вакуумную установку с вертикально установленной газоразрядной трубкой с внутренним диаметром 2 6 см (внешний 3 см) Промежуток между электродами 45 см Макрочастицы - полидисперсный порошок А12Оз (1-60 мкм), инжектируются из контейнера, расположенного в верхней части труб-

кн. Исследования проводились в постоянном токе тлеющего разряда в атмосфере Аг.

8 установке использовались три идентичные фотокамеры "Canon EOS 10D" с объективами "Canon ЕГ 100 F2.8 Macro" и удлинительными кольцами. Камеры устанавливались по уровню в штативы па одинаковом расстоянии от трубки. Синхронизация фотокамер осуществлялась, благодаря возможности одновременного электронного управления несколькими затворами. Сгруктура была освещена расширенным лучом полупроводникового лазера с длиной волны 532 нм и мощностью 51 мВт. При этом луч направлялся не перпендикулярно трубке, а под уг лом - 45", чтобы избежать бликов.

Для обработки результатов эксперимента предложенным алгоритмом использовалась разработанная программа 'Tomography 2", В качестве входных параметров используются: L расстояние от газоразрядной трубки до светочувствительной матрицы (на корпусе фотоаппаратов Canon EOS 10D есть метка, указывающая на расположение матрицы); R tr радиус газоразрядной трубки; физические размеры пикселя светочувствительной матрицы; масштаб проекций; углы поворота фотоаппаратов; диапазон угловой погрешности AQ\ начальные смещения проекций (задаются нулевыми).

• - ---W— •

. » * » »

* »

г ' i- ■ f * ♦ • •

•> ' > }: ' - * * • • •

% 4 * • # • •

• » i ь - 4 ♦

» « 4 » t i — - ( 1 • ' i # • * •

f • * —__ , » • y-

4 * 1 » *

* » * t * f / a *

и ш - - ♦

Рис.5. Проекции (а, Ь, с) и восстановленная структура (d).

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

рис 5 Газ Аг, Р=29 Па, 1=0 31 мА, масштаб 105 pix на 250 мкм, время экспозиции 1/500 с , диафрагма 4, L=388 мм, размер пикселя светочувствительной матрицы 7 25 мкм

Серые размытые пятна (рис 5) — частицы плазменно-пылевой структуры, черные одинаковые кружки - проекции, вычисленные с восстановленной структуры и наложенные поверх изображений частиц плазменно-пылевой структуры для проверки точности восстановления

Результат обработки можно сразу оценить визуально, т к восстановленная структура выводится в 3-х мерном виде (рис 5d), благодаря использованию разработанного модуля "Result of reconstruction" в составе программы "Tomography 2"

В качестве начальных приближенных значений в программу были введены следующие углы проекции 0°, -36 03°, -160 02° (относительно первой проекции). В результате обработки оптимальными признаны углы -1 05°, -31 63° , -158 53° Для смещений оптимальные результаты shift Х=6 pix (вторая проекция), shift Y=29 pix (вторая проекция), shift Х=2 pix (третья проекция), shift Y=67 pix (третья проекция)

Рис 6 Парная корреляционная функция

Восстановленная структура содержит 40 частиц Минимальное отклонение "Л" составляет 2 83 р1х Среднее межчастичное расстояние 298±7 (мкм) На рис 6 представлена парная корреляционная функция восстановленной трехмерной структуры с гЕтах=364 мкм

Предложенным методом исследованы зависимости параметров трехмерной пылевой структуры из полидисперсного порошка А1203 с диаметром частиц от нескольких единиц до нескольких десятков

микрон от токов (в диапазоне 0 1-12 мА, с шагом 0 1 мА) и давлений (27-80 Па, с шагом 13 Па) в тлеющем разряде постоянного тока в Аг

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

Применение томографии позволило обнаружить разницу между горизонтальными (гху) и вертикальными (г2) составляющими межчастичного расстояния (гср) Выявлено, что "г2" в пределах погрешности не зависит от токов и давлений "гху" показывало зависимость от токов при всех давлениях линейную при давлениях 27 Па (рис 7), 40 Па и 53 Па, не линейную, которую можно аппроксимировать квадратичной зависимостью, при, например, 67 Па (рис. 8), 80 Па Зависимость "гср" от тока проявляется как линейно, при давлении 67 Па, так и не

линейно при 80 Па (рис 9) Зависимость погрешности не обнаружена

от давления в пределах

Рис 7 Зависимость горизонтальной составляющей от тока при Р=27 Па

Зависимость гху(1), Аг, 67 Па

0,4 0,5

0,7 0,8

♦ Эксперимент —"Аппроксимация

1,1

Рис 8 Зависимость горизонтальной составляющей от тока при Р=67 Па

В заключении сформулированы основные результаты, полученные автором

1 Адаптирован метод оптической томографии применительно к задаче пространственной диагностики плазменно-пылевых образований

2 Сформулированы ограничения и условия применимости адаптированного метода

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

4 Предложенным методом проведено исследование зависимости средних межчастичных расстояний при различных токах (0,1-1,2 мА, с шагом 0,1 мА) и давлениях (27-80 Па, с шагом 13 Па) в тлеющем разряде постоянного тока в Аг для полидисперсного порошка А1203

В приложении даны описания интерфейсов разработанного программного обеспечения и базовые понятия теории множеств

Публикации по материалам диссертации

1 Bulba А V , Luizova L А , Khakhaev A D The computer training

program on tomography // Education and Training in Optics and

Photonics (ETOP) Conference (October 5-9, 2003, Tucson, Arizona) Tucson, 2003 Pp 89-97

2 Bulba A V , Luizova L A, Khakhaev A D Dusty plasma tomography // Northern Optics 2003 The joint conference of the Optical Societies of Denmark, Finland, Norway and Sweden (16-18 June 2003, Espoo, Finland) Espoo, 2003 P 49

3 Luizova L , Khakhaev A , Podryadchikov S., Scherbina A., Bulba A Investigation of ordered dusty structures in the glow discharge // Proceedings of the 30th EPS Conference on Controlled Fusion and Plasma Physics (July 7-11,2003, St Petersburg, Russia) St Petersburg, 2003

4 Бульба А В Томография плазменного кристалла // Восьмая Санкт-Петербургская ассамблея молодых ученых и специалистов СПб, изд СПбГУ, 2003 С 27

5 Бульба А В , Луизова JI А Компьютерная обучающая программа по томографии // Технологии Microsoft в теории и практике программирования (Санкт-Петербургский государственный политехнический университет, 4-5 марта 2004 г) СПб, изд СПбГПУ, 2004 С 80

6 Bulba А V , Luizova L А, Khakhaev А , Podryadchikov S F , Shty-kov A S , Scherbina A I The Influence of Physical-Chemical Characteristics of Plasma-Forming Gas and Macroparticle Matter in Complex Plasma on Ordered Structure Self-Organization // New Vistas in Dusty Plasmas AIP Conference proceeding V 799, Melville, New York, 2005,pp 359-362

7 Bulba A V, Luizova L A, Khakhaev A D Research Possibilities of 3D Dasty Plasma Structures // New Vistas m Dusty Plasmas AIP Conference proceeding V 799, Melville, New York, 2005, pp 355-358

8 Bulba A V , Ekimov D A , Luizova L A , Khakhaev A D Investigation of 3D dusty plasma structures // V International Conference «Plasma Physics and Plasma Technology» (18-22 September 2006, Minsk, Belarus) Contributed papers Minsk, 2006 V 1 Pp 396-399

9 Бульба А. В, Луизова Л A, Подрядчиков С Ф, Хахаев А Д, Щербина А И Самоорганизация и рост пылевых структур в тлеющем разряде//Химия высоких энергий 2006 том 40 №2 С 155-160

10 Анисимов А Л , Бульба А В , Луизова Л А , Хахаев А Д , Штыков А С Невозмущающие методы оптической диагностики гетерогенной плазмы // Химия высоких энергий 2006 том 40 № 3 С 233-237

11. Бульба А В, Луизова Л А, Хахаев А Д Обучение принципам компьютерной томографии для решения задач диагностики плазмы

11 Электронный журнал «Исследовано в России» 2006 С 677-686 http //¿hurnal аре relarn ru/articles/2006/ 068 pdf

12 Бульба А В , Екимов Д А , Луизова Л А , Хахаев А Д Пространственная диагностика плазменно-пылевых структур // IV Курчатовская молодежная научная школа (20-22 ноября 2006 г, Москва, РНЦ Курчатовский Институт) Сборник трудов Москва, 2006

13 Бульба А В Программа для определения пространственного расположения набора точечных объектов методом оптической реконструктивной томографии «Tomography 2» // Компьютерные учебные программы и инновации (Телеграф отраслевого фонда алгоритмов и программ) - 2007, № 5 С 26

14 Бульба А В Программа для обучения алгоритмам реконструктивной томографии «Tomography» // Компьютерные учебные программы и инновации (Телеграф отраслевого фонда алгоритмов и программ) -2007, №5 С 10

15 Бульба А В Исследование пространственного распределения точечных масс методом оптической томографии // Вестник ИРГТУ 2007 №3 С 116-125

16 Бульба А В, Луизова Л А, Хахаев А Д Томография плазменно-пылевых структур // Материалы семинара-школы молодых ученых, студентов и аспирантов "Физика низкотемпературной плазмы 2004" Изд Скандинавия, 2005 С 116-126

17 Бульба А В, Луизова Л А, Пискунов А А, Соловьев А В Исследование структуры плазменно-пылевых образований // Материалы Всероссийской (с международным участием) конференции "Физика низкотемпературной плазмы — 2007" (24—28 июня 2007 г), изд ПетрГУ, 2007, Т 2, С 214-218

Цитированная литература

18 Chu J H , Lin I Coulomb Lattice m a Weakly Ionized Colloidal Plasma // Physica A 1994 V 205 P 183

19 Thomas H, Morfill G E , Demmel V, Goree J Plasma crystal coulomb crystallization in a dusty plasma//Phys Rev Lett 1994 V 73 P 652

20 Нефедов А П, Петров О Ф , Фортов В Е Кристаллические структуры в плазме с сильным взаимодействием макрочастиц // УФН 1997 Т 167 № 11 С 1215-1226

21 Фортов В Е , Петров О Ф , Молотков В И Пылевая плазма в условиях микрогравитации эксперименты на международной космической станции // Электронный научный журнал «ИССЛЕДО-

ВАНО В РОССИИ» http //zhurnal аре relarn m/articles/2002/204 pdf С 2244-2247

22 Ваулина О С , Нефедов А П , Петров О Ф , Фортов В Е Формирование упорядоченных структур заряженных макрочастиц в фотоэмиссионной ловушке//ЖЭТФ 2000 Т 118 С 351

23 Zuzic М , Ivlev А V , Goree J , et al Three-Dimensional Strongly Coupled Plasma Crystal under Gravity Conditions // Phys Rev Lett 2000 V 85, № 19 P 4064

24 Arp О , Block D , Piel A Structure and trapping of three-dimensional dust clouds in a capacitively coupled rf-discharge // New Vistas in Dusty Plasmas AIP Conference proceeding V 799, Melville, New York, 2005, pp 20-21

25 Kading S , Melzer A, Arp О , Block D , Piel A Stereoscopic Investigations of 3D Coulomb Balls // New Vistas in Dusty Plasmas AIP Conference proceeding V 799, Melville, New York, 2005, pp 335-338

26 ХорнБ К П Зрение роботов - М Мир, 1989

27 Форсайт Д , Понс Ж Компьютерное зрение Современный подход Пер с англ - М Издательский дом "Вильяме", 2004

28 Абрамов В В , Киричук В С и др Реконструкция трехмерных поверхностей по двум проекциям при отслеживании камерой заданной точки сцены // Автометрия, 1998, № 5, с 3

29 Antipov S N, Khrustalyov Yu V , Petrov О F , Statsenko KB, Va-siliev M M, Levchenko V D Structure and melting of 3D anisotropic dust crystals in dc glow discharges // 33rd EPS Conference on Plasma Phys Rome, 19-23 June 2006 ECA Vol 301, D-5 028 (2006)

30 Сизиков В С Математические методы обработки результатов измерений Учебник для вузов - СПб Политехника, 2001

31 Пикалов В В , Преображенский Н Г Вычислительная томография и физический эксперимент //УФН 1983 Т 143 вып 3 С 469-498

32 Терновой К С , Синков М В Введение в современную томографию -Киев Науковадумка 1983

33 Тихонов А Н, Арсенин В Я , Тимонов А А Математические задачи компьютерной томографии М Наука 1987

34 Гонсалес Р , Вудс Р Цифровая обработка изображений - М Техносфера 2006

Подписано в печать 19 09 07 Формат 60x84 '/16 1 уч -изд л Тираж 100 экз Изд № 220

Государственное образовательное учреждение высшего профессионального образования ПЕТРОЗАВОДСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

Отпечатано в типографии Издательства ПетрГУ 185910, г Петрозаводск, пр Ленина, 33

ВВЕДЕНИЕ.

1 ПРОСТРАНСТВЕННАЯ ДИАГНОСТИКА ПЛАЗМЕННО-ПЫЛЕВЫХ СТРУКТУР.

1.1 Суть явления плазменно-пылевых образований.

1.2 Методы пространственной диагностики пылевых образований.

1.2.1 Сканирование лазерным ножом.

1.2.2 Стереоскопический подход.

1.2.2.1 Поиск сопряженных точек по эпиполярному условию.

1.2.2.2 Корреляционный критерий сходства.

2 ПРИНЦИПЫ ТОМОГРАФИИ.

2.1 Определение томографии.

2.2 Классификация методов вычислительной томографии.

2.3 Математическая постановка задачи томографии.

2.4 Восстановление с использованием Фурье преобразований.

2.5 Метод обратного проецирования.

3 АЛГОРИТМ ТОМОГРАФИРОВАНИЯ ПЫЛЕВЫХ СТРУКТУР.

3.1 Постановка задачи.

3.2 Математическое описание алгоритма.

3.3 Оценка факторов, ограничивающих применение томографии.

3.4 Тестирование алгоритма на математических моделях.

4 ПРЕДОБРАБОТКА ИЗОБРАЖЕНИЙ ПЫЛЕВЫХ СТРУКТУР.

Введение.

4.1 Цветовое представление изображений.

4.1.1 Цветовые модели.

4.1.2 Цветовая модель RGB.

4.2 Фильтрация изображений.

4.2.1 Фильтрация в частотной области.

4.2.2 Гауссов фильтр низких частот.

4.2.3 Пороговая обработка.

4.3 Морфологическая обработка изображений.

4.3.1 Введение.

4.3.2 Выделение связных компонент.

5 ВОССТАНОВЛЕНИЕ ПРОСТРАНСТВЕННОГО РАСПРЕДЕЛЕНИЯ ЧАСТИЦ

ПЛАЗМЕННО-ПЫЛЕВОЙ СТРУКТУРЫ.

5.1 Экспериментальная установка.

5.1.1 Вакуумная часть экспериментальной установки.

5.1.2 Электрическая часть экспериментальной установки.

5.1.3 Оптическая часть экспериментальной установки.

5.2 Исследование пространственной структуры пылевых образований.

Актуальность работы

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

При диагностике пылевых образований особую роль играет выяснение геометрических параметров формируемой в плазме структуры пылевых частиц. Традиционно пространственная диагностика плазменно-пылевых структур сводится к анализу их изображений, полученных в лазерном "ноже", т. е. в каждый момент времени регистрируется одно горизонтальное или вертикальное изображение структуры, не отражающее распределение частиц в пространстве, если речь не касается исследования монослойных плазменных образований. Достоверность получаемых таким способом результатов зависит от толщины "ножа" и характера структуры; в некоторых случаях возможны значительные ошибки как в оценке среднего расстояния между частицами, так и в виде парной корреляционной функции, характеризующей степень упорядоченности структуры. В некоторых экспериментах плазменно-пылевые объекты наблюдались путём горизонтального сканирования вертикальным лазерным "ножом". При этом осуществлялась запись положения каждой частицы с помощью телекамеры. Полная длительность такого сканирования составляла более десятка секунд, что совершенно неприемлемо в случае осцилляции частиц. Использование же классической стереоскопической системы машинного зрения позволяет осуществлять реконструкцию только приблизительно, поскольку остаётся не решённой задача сопряжения характерных точек изображений.

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

Цель работы, задачи

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

1. Проанализировать существующие оптические методы пространственной диагностики плазменно-пылевых структур.

2. Адаптировать метод оптической томографии применительно к задаче пространственной диагностики пылевых структур.

3. Сформулировать ограничения и условия применимости метода.

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

5. Разработать необходимое ПО. Проверить работоспособность и эффективность разработанного алгоритма с использованием модельных и экспериментальных проекционных данных.

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

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

Научная новизна

1. Развит метод оптической томографии применительно к задаче пространственной диагностики плазменно-пылевых образований.

2. Сформулированы ограничения и условия применимости развитого метода.

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

4. Впервые предложенным методом в тлеющем разряде постоянного тока проведено исследование зависимости межчастичных расстояний трехмерной пылевой структуры от условий разряда.

Научно-практическая значимость работы

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

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

Изучение плазменных кластеров (как двумерных, так и трёхмерных) представляет собой отдельную, очень интересную, развивающуюся область исследований. Здесь поставлено множество вопросов, в решение которых представленная методика может сыграть существенную роль, как инструмент контроля геометрических параметров:

- одной из наиболее интересных проблем является образование центров кристаллизации, переход от "пылевой" молекулы к молекулярному кластеру и, затем, к твердому состоянию.

- проблема непрерывного перехода малых пылевых кластеров с небольшим количеством пылинок к трёхмерным или двумерным "кристаллам" с многими пылинками.

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

- исследование двойных смесей (кластеров с пылинками различного размера) и неоднородных кластеров.

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

Основные положения, выносимые на защиту:

1. Метод оптической томографии, адаптированный применительно к задаче пространственной диагностики плазменно-пылевых структур, для которого сформулированы ограничения и условия применимости.

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

3. Набор программных модулей "Tomography 2", основанный на предложенном алгоритме, в реализации которых учитывалось влияние используемых оптических систем регистрации и зашумленность проекционных данных.

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

Апробация работы

Содержание работы докладывалось на Международной конференции Оптических Обществ Дании, Финляндии, Норвегии и Швеции «Северная оптика 2003» (г. Эспоо, Финляндия, 16-18 июня 2003 г.), Международной конференции «Образование и обучение в оптике и фотонике (ЕТОР)» (г. Туссон, USA, штат Аризона, 5-9 октября 2003 г.), IV международном симпозиуме по теоретической и прикладной плазмохимии ISTAPC-2005 (г. Иваново, Россия, 13-18 мая 2005 г.), V международной конференции «Физика плазмы и плазменные технологии» РРРТ-5 (г. Минск, Белоруссия, 18-22 сентября 2006 г.), IV Курчатовской молодежной научной школе (г. Москва, 20 - 22 ноября 2006 г.), Всероссийской (с международным участием) конференции "Физика низкотемпературной плазмы - 2007" (г. Петрозаводск, 24—28 июня 2007 г.).

Успешно пройдена процедура государственной регистрации в «Национальном информационном фонде неопубликованных документов» разработки, предъявленной в отраслевой фонд алгоритмов и программ: программа для определения пространственного расположения набора точечных объектов методом оптической реконструктивной томографии «Tomography 2» (номер государственной регистрации: 50200701083 от 28.05.2007). Отраслевым фондом алгоритмов и программ выдано свидетельство об отраслевой регистрации разработки №8356.

Федеральной службой по интеллектуальной собственности, патентам и товарным знакам (РОСПАТЕНТ) выдано свидетельство об официальной регистрации программы для ЭВМ «Tomography 2» (объекту регистрации присвоен регистрационный № 2007612629 от 20.06.2007).

Пройдена процедура государственной регистрации в «Национальном информационном фонде неопубликованных документов» разработки, предъявленной в отраслевой фонд алгоритмов и программ: программа для обучения алгоритмам реконструктивной томографии «Tomography» (номер государственной регистрации: 50200700997 от 15.05.2007). Отраслевым фондом алгоритмов и программ выдано свидетельство об отраслевой регистрации разработки №8286.

Федеральной службой по интеллектуальной собственности, патентам и товарным знакам (РОСПАТЕНТ) выдано свидетельство об официальной регистрации программы для ЭВМ «Tomography» (объекту регистрации присвоен регистрационный № 2007611674 от 20.04.2007).

По результатам работы опубликованы следующие статьи:

1. Bulba А.V., Luizova L.A., Khakhaev A.D. The computer training program on tomography // Education and Training in Optics and Photonics (ETOP) Conference (October 5-9, 2003, Tucson, Arizona). Tucson, 2003. Pp. 89-97.

2. Bulba A.V., Luizova L.A., Khakhaev A.D. Dusty plasma tomography // Northern Optics 2003. The joint conference of the Optical Societies of Denmark, Finland, Norway and Sweden (16-18 June 2003, Espoo, Finland). Espoo, 2003. P. 49.

3. Luizova L., Khakhaev A., Podryadchikov S., Scherbina A., Bulba A. Investigation of ordered dusty structures in the glow discharge // Proceedings of the 30th EPS Conference on Controlled Fusion and Plasma Physics ( July 7-11, 2003, St. Petersburg, Russia). St. Petersburg, 2003.

4. Бульба A.B. Томография плазменного кристалла // Восьмая Санкт-Петербургская ассамблея молодых ученых и специалистов. СПб, изд. СПбГУ, 2003. С. 27.

5. Бульба А.В., Луизова Л.А. Компьютерная обучающая программа по томографии // Технологии Microsoft в теории и практике программирования (Санкт-Петербургский государственный политехнический университет, 4-5 марта 2004 г.). СПб, изд. СПбГПУ, 2004. С. 80.

6. Bulba А. V., Luizova L. A., Khakhaev A., Podryadchikov S. F., Shty-kov A. S., Scherbina A. I. The Influence of Physical-Chemical Characteristics of Plasma-Forming Gas and Macroparticle Matter in Complex Plasma on Ordered Structure Self-Organization // New Vistas in Dusty Plasmas. AIP Conference proceeding .V. 799, Melville, New York, 2005, pp. 359-362.

7. Bulba A. V., Luizova L. A., Khakhaev A. D. Research Possibilities of 3D Dasty Plasma Structures // New Vistas in Dusty Plasmas. AIP Conference proceeding. V. 799, Melville, New York, 2005, pp. 355-358.

8. Bulba A. V., Ekimov D.A., Luizova L.A., Khakhaev A.D. Investigation of 3D dusty plasma structures // V International Conference «Plasma Physics and Plasma

Technology» (18-22 September 2006, Minsk, Belarus). Contributed papers. Minsk, 2006. V. l.pp. 396-399.

9. Бульба A.B., Луизова Л.А., Подрядчиков С.Ф., Хахаев А.Д., Щербина А.И. Самоорганизация и рост пылевых структур в тлеющем разряде // Химия высоких энергий. 2006. том 40. №2. С. 155-160.

10. Анисимов А.Л., Бульба А.В., Луизова Л.А., Хахаев А.Д., Штыков А.С. Невозмущающие методы оптической диагностики гетерогенной плазмы // Химия высоких энергий. 2006. том 40. №3. С. 233-237.

11. Бульба А.В., Луизова Л.А., Хахаев А.Д. Обучение принципам компьютерной томографии для решения задач диагностики плазмы // Электронный журнал «Исследовано в России». 2006. С. 677 - 686. http://zhurnal.ape.relarn.ru/articles/2006/ 068.pdf

12. Бульба А.В., Екимов Д.А., Луизова Л.А., Хахаев А.Д. Пространственная диагностика плазменно-пылевых структур // IV Курчатовская молодежная научная школа (20 - 22 ноября 2006 г., Москва, РНЦ Курчатовский Институт). Сборник трудов. Москва, 2006.

13. Бульба А.В. Программа для определения пространственного расположения набора точечных объектов методом оптической реконструктивной томографии «Tomography 2». // Компьютерные учебные программы и инновации (Телеграф отраслевого фонда алгоритмов и программ). - 2007, №5. С. 26.

14. Бульба А.В. Программа для обучения алгоритмам реконструктивной томографии «Tomography». // Компьютерные учебные программы и инновации (Телеграф отраслевого фонда алгоритмов и программ). - 2007, №5. С. 10.

15. Бульба А.В. Исследование пространственного распределения точечных масс методом оптической томографии // Вестник ИРГТУ. 2007. №3. С.116-125.

16. Бульба А.В., Луизова Л.А., Хахаев А.Д. Томография плазменно-пылевых структур // Материалы семинара-школы молодых ученых, студентов и аспирантов "Физика низкотемпературной плазмы 2004". 2005. С. 116-126.

17. Бульба А. В., Луизова Л. А., Пискунов А. А., Соловьев А. В. Исследование структуры плазменно-пылевых образований // Материалы Всероссийской (с международным участием) конференции "Физика низкотемпературной плазмы - 2007" (24-28 июня 2007 г), изд. ПетрГУ, 2007, Т2, С.214-218.

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

Вклад автора

Адаптирован метод оптической томографии применительно к задаче пространственной диагностики плазменно-пылевых образований.

Сформулированы ограничения и условия применимости адаптированного метода.

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

Предложенным методом проведено исследование зависимости средних межчастичных расстояния при различных токах и давлениях в тлеющем разряде постоянного тока в Аг для полидисперсного порошка AI2O3.

Объём и структура диссертации

Диссертация состоит из введения, пяти глав, заключения и приложения. Работа содержит 120 страниц, 74 рисунка, 6 таблиц. Список использованной литературы включает 54 наименования.

ЗАКЛЮЧЕНИЕ

1. Адаптирован метод оптической томографии применительно к задаче пространственной диагностики плазменно-пылевых образований.

2. Сформулированы ограничения и условия применимости адаптированного метода.

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

4. Предложенным методом проведено исследование зависимости средних межчастичных расстояния при различных токах (0,1-1,2 мА, с шагом 0,1 мА) и давлениях (27 - 80 Па, с шагом 13 Па) в тлеющем разряде постоянного тока в Аг для полидисперсного порошка AI2O3.

Работа выполнена в Научно-образовательном центре по фундаментальным проблемам приложений физики низкотемпературной плазмы (НОЦ «Плазма») - грант RUX0-000013-PZ-06 Американского фонда гражданских исследований и развития, при под держке Министерства образования и науки РФ и Правительства Республики Карелии.

1. Нефёдов А.П. Плазменно-пылевые структуры в низкотемпературной плазме // Материалы пленарных докладов конференции по физике низкотемпературной плазмы ФНТП-2001. Петрозаводский государственный университет, 2001, С.5.

2. Pollick Е. L., Hansen J.P. Statistical Mechanics of Dense Ionized Matter // Phys Rev.A.-1973.-V. 8, №6. P. 3110-3122.

3. Ichimaru S. Strongly Coupled Plasmas: High-Density Classical Plasmas and Degenerate Electron Liquids//Rev. Mod. Phys.- 1982.-V. 54. P. 1017-1059.

4. Луизова JI.A., Хахаев А.Д. Проблемы и перспективы исследования упорядоченных структур в плазме // НОЦ "Плазма", Петрозаводск, 2002.

5. Ikezi Н. Coulomb Solid of Small Particles in Plasmas // Phys. Fluids.-1986.- V.29.-P.1764-1766.

6. Hayashi Y., Tachibana K. // Jpn. Appl. Phys. 1994, V. 33., P.301.

7. Fortov V.E., Nefedov A.P., Petrov O.F. et al. // Phys. Let. A. 1996, V. 219., P.89.

8. Fortov V.E., Molotkov V.I., Nefedov A.P., Petrov O.F. // Phys. Plasmas. 1999, V. 6, P.1759.

9. Fortov V.E., Nefedov A.P., Sinel"shchikov V. et al. // Phys. Lett. A. 2000, V. 267, P. 179.

10. Chu J.H., Lin I. Coulomb Lattice in a Weakly Ionized Colloidal Plasma // Physica A. 1994. V. 205. P. 183.

11. Chu J.H. Direct Observation of Coulomb Crystals and Liquids in Strongly Coupled Rf Dusty Plasmas // Phys. Rev. Lett. 1994, V. 72, P. 4009.

12. Thomas H., Morfill G.E., Demmel V., Goree J. Plasma crystal: coulomb crystallization in a dusty plasma // Phys. Rev. Lett. 1994, V.73, P. 652.

13. Нефедов А.П., Петров О.Ф., Фортов B.E. Кристаллические структуры в плазме с сильным взаимодействием макрочастиц // УФН. 1997, Т. 167, № 11, С. 1215 1226.

14. Ваулина О.С., Нефедов А.П., Петров О.Ф., Фортов В.Е. Формирование упорядоченных структур заряженных макрочастиц в фотоэмиссионной ловушке // ЖЭТФ. 2000. Т. 118. С. 351.

15. Zuzic M., Ivlev A.V., Goree J., et al Three-Dimensional Strongly Coupled Plasma Crystal under Gravity Conditions // Phys. Rev. Lett. 2000, V. 85, №19, P.4064.

16. ArpO., Block D., PielA. Structure and trapping of three-dimensional dust clouds in a capacitively coupled rf-discharge // New Vistas in Dusty Plasmas. AIP Conference proceeding. V. 799, Melville, New York, 2005, pp. 20-21.

17. Грузман И.С., Киричук B.C., Косых В.П., Перетяган Г.И., Спектор А.А. Цифровая обработка изображений в информационных системах: Учебное пособие.-Новосибисрк: Изд-во НГТУ, 2000. 168.

18. KadingS., MelzerA., Агр О., Block D., PielA. Stereoscopic Investigations of 3D Coulomb Balls // New Vistas in Dusty Plasmas. AIP Conference proceeding. V. 799, Melville, New York, 2005, pp. 335-338.

19. Kading S., Melzer A. Stereoscopic investigations of 2D and 3D dust crystals // Fourth International Conference on the Physics of Dasty Plasmas (13-17 June 2005, Orleans): book of abstracts. Orleans, 2005. II. 3.

20. Edward Thomas, Jeremiah Williams Experimental Measurements of Velocity Dissipation and Neutral-Drag Effects during the Formation of a Dusty Plasma // Physical Review Letters, V. 95, Issue 5, id. 055001, jul., 2005.

21. О.Ф.Петров, В.И. Молотков, O.C. Ваулина, В.Е.Фортов. Диагностика пылевой плазмы. ЭНТП, серия Б, том V-1, часть 1. М.: Изд. Янус-К, 2006. С. 136-147.

22. Хорн Б.К.П. Зрение роботов. -М.: Мир, 1989.

23. Форсайт Д., Понс Ж. Компьютерное зрение. Современный подход. : Пер. с англ. -М.: Издательский дом "Вильяме", 2004.

24. Абрамов В.В., Киричук B.C. и др. Реконструкция трехмерных поверхностей по двум проекциям при отслеживании камерой заданной точки сцены. // Автометрия, 1998, №5, с.З.

25. Annaratone B.M., Antonova Т., Goldbeck D.D., Tomas H.M., Morfill G.E. // Plasma Physics and Controlled Fusion 46. 2004, B495.

26. Antonova Т., Annaratone B.M., Goldbeck D.D., Yaroshenko V. Interaction among particles in 3D plasma clusters // New Vistas in Dusty Plasmas. AIP Conferenceproceeding. V. 799, Melville, New York, 2005, pp. 299-302.

27. Antonova Т., Annaratone B.M., Thomas H.M., Morfill G.E. The structure of a plasma cluster as seen by an injected particle // 32nd EPS Conference on Plasma Phys.Tarragona, 27 June 1 Jule 2005 ECA Vol.29C, P- 2.132.

28. Битюцкий О.И., Перетяган Г.И. Поиск и локализация реперных фрагментов при совмещении повторных снимков // Автометрия. 1988, №3.

29. Васильева И. А. Основы спектральной диагностики газа с конденсированной дисперсной фазой УФН, 1993, Т.163, №8, С.47-88.

30. Пикалов В.В., Преображенский Н.Г. Вычислительная томография и физический эксперимент. // УФН. 1983. Т.143. вып. 3. С. 469^198.

31. Мельникова Т. С, Пикалов В. В., Преображенский Н. Г. О локальной диагностике оптически плотной асимметричной плазмы. // Опт. и спектр. 1980. Т. 48, С. 474.

32. Вишняков Г.Н., Левин Г.Г. Оптическая томография фазовых объектов. // Оптика и спектроскопия. 1982. Т.53, Вып. 4, С.731-735.

33. Пикалов В.В., Преображенский Н. Г. Реконструктивная томография в газодинамике и физике плазмы. Новосибирск: Наука. 1987.

34. Пикалов В.В., Мельникова Т.С. Томография плазмы. Новосибирск: Наука. Сибирская издательская фирма РАН 1995.

35. Резник. А.Н., Юрасова Н.В. Ближнепольная СВЧ томография биологических сред. // ЖТФ. 2004. Т.74. Вып. 4, С. 108-116.

36. Филонин О.В. Малоракурсная томография крупногабаритных объектов. // Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ», 2006, С. 1455, http://zhurnal.ape.relarn.ru/articles/2006/157.pdf

37. Наттерер Ф. Математические аспекты компьютерной томографии. М.: Мир. 1990.

38. Терновой К.С., Синков М.В. Введение в современную томографию. Киев: Наукова думка. 1983.

39. Тихонов А.Н., Арсенин В.Я., Тимонов А.А. Математические задачи компьютерной томографии. М.: Наука. 1987.

40. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука. 1986.

41. Тихонов А.Н., Гончаровский А.В., Степанов В.В., Ягола А.Г. Численные методы решения некорректных задач. М.: Наука. 1990.

42. Календер В. Компьютерная томография. М: Техносфера, 2006. 344с.

43. Сизиков B.C. Математические методы обработки результатов измерений: Учебник для вузов. СПб: Политехника, 2001.

44. Хермен Г. Восстановление изображений по проекциям: основы реконструктивной томографии. Пер. с англ. М.: Мир. 1983.

45. Денисова Н.В., Пикалов В.В., Баландин А.Л. Модифицированный метод максимума энтропии в томографии плазмы. // Оптика и спектроскопия. 1996, Т.81, №1, С.43-48.

46. Bracewell R.N. Strip integration in radioastronomy // Aust. J. Phys., 1956, 9, P. 198-217.

47. De Roster D., Klug A. Reconstruction of three-dimensional structures from electron micrographs//Nature, 1968,217, P. 130-134.

48. Rowley P.D. Quantitative interpretation of three-dimensional weakly refractive phase objects using holographic interferometry // J. Opt. Soc. Amer., 1969, 59, P. 1496-1498.

49. Бульба A.B., Луизова Л.А., Хахаев А.Д. Обучение принципам компьютерной томографии для решения задач диагностики плазмы // Электронный журнал «Исследовано в России». 2006, С. 677 686. http://zhurnal.ape.relarn.ru/articles/2006/ 068.pdf

50. Выгодский М. Я. Справочник по высшей математике. М: государственное издательство физико-математической литературы, 1958.

51. Гонсалес Р., Вудс Р. Цифровая обработка изображений. М.: Техносфера. 2006.

52. Лайонс Р. Цифровая обработка сигналов. М.: Бином. 2006.