Рассматривается актуальная задача определения максимального подмножества множества точек в R
по ограниченному набору их двумерных проекций. Данная задача естественно возникает в приложениях физической гидроаэродинамики для оптической диагностики реальных течений жидкости и газа путем измерения мгновенных полей скорости в объеме потока. Предложены методы восстановления указанного множества и определения достаточности измерений для однозначного решения обратной задачи.
Discrete scattering centers reconstruction problem on a limited set of projections.pdf В настоящее время томографические методы широко применяются для исследования внутренней структуры объектов. В основном они используются в медицинской диагностике, в промышленности и строительстве для контроля качества всевозможных изделий и инженерных сооружений. Математический аппарат вычислительной томографии изложен в монографиях [1, 2]. В последнее время все чаще стали появляться примеры использования томографии и в других областях науки и техники. Например, оптическая томография рассеивающих сред сравнительно недавно стала применяться в хорошо зарекомендовавшем себя, методе цифровой трассерной визуализации - Particle Image Velocimetry (PIV) [3]. Благодаря возможности регистрировать мгновенные пространственные распределения скорости метод цифровой трассерной визуализации используется в фундаментальных научных исследованиях для изучения потоков, содержащих вихревые структуры, информация о которых частично теряется при применении одноточечных методов диагностики. К подобным течениям можно отнести турбулентные течения, включая струи, следы, слои смешения и многое другое. Применение полевых методов дает возможность получения информации о динамике структур, их масштабах, расчета дифференциальных характеристик, пространственных и пространственно-временных корреляций, а также статистических характеристик потока. Кроме того, измерительные системы на основе данного метода применяются в прикладных исследованиях для авиастроительной индустрии, автомобилестроения, энергетики, химической промышленности, биологии и медицины. Основным отличием томографической цифровой трассерной визуализации (Tomo PIV) для экспериментальной гидроаэродинамики от медицинских томографических исследований является малое количество проекций, используемое для восстановления объемного изображения измерительной области потока. Вследствие записи небольшого числа проекций задача получения приближения исходного распределения интенсивности рассеянного света на трассерах усложняется, однако форма и разреженность полезного сигнала (локализованные в пространстве узкие пики интенсивности), в конечном итоге, позволяют реконструировать объемное распределение интенсивности требуемого качества, подходящее для дальнейшего корреляционного анализа и оценки смещения трассеров в потоке. В данном случае для томографической реконструкции измерительной области потока применяются итерационные алгебраические методы [4, 5]. На рис. 1 представлена схема расчета поля скорости в измерительном объеме. Альтернативными методами получения трехкомпонентной трехмерной картины поля течения являются голографический, сканирующий и 3D Particle Tracking Velocimetry (3D PTV) [6]. Последний метод, 3D PTV, по принципу действия наиболее близок к томографической цифровой трассерной визуализации. Однако он требует поддержания сравнительно низкой концентрации трассеров для успешной идентификации отдельных частиц на всех проекциях путем триангуляции и отслеживания их перемещения для оценки скорости. Перечислим некоторые известные сложности Tomo PIV-измерений. Во-первых, концентрация частиц в области измерения должна быть достаточной для получения хорошего пространственного разрешения оценки скорости. C другой стороны, она не должна превышать критического значения ~0,2 частиц/пикс, выше которого томографическая реконструкция существенно затрудняется вследствие затенения частицами друг друга. При этом стандартные алгоритмы реконструкции ведут к появлению ложных образов трассеров в восстанавливаемом объеме вследствие возрастания неопределенности при недостатке проекционных данных [7, 8]. Таким образом, в эксперименте необходимо поддерживать определенную концентрацию трассеров на проекции, которая бы позволила увеличить вероятность получения истинной конфигурации рассеивающих центров и, в итоге, надежной оценки поля скорости потока. В данной работе предложен и математически обоснован алгоритм получения максимального множества точек, проекции которого в плоскости, заданной светочувствительными матрицами камер, совпадают с зарегистрированными проекциями. Обсуждаются вопросы необходимости и достаточности проекционных данных для надежного определения внутренней структуры среды, состоящей из набора дискретных рассеивающих центров. 3D корреляционная функция Реконструированное объемное изображение Вектор скорости для расчетной области Рис. 1. Иллюстрация принципов работы алгоритма оценки поля скорости в методе Tomo PIV Томографические трассерные проекции 1. Метод реконструкции Сформулируем задачу в ее математической постановке. Пусть X = {, x2,..., xn} - множество попарно различных точек в R3. Обозначим через Рп = {p1, p2,..., pm} - проекцию точек из множества X на плоскость п. Очевидным образом выполняется неравенство m < n и, в общем случае, m = m(X,п). Здесь и далее, будем считать, что любая точка из исходного множества не совпадает с любой ее проекцией. Пусть далее задан набор плоскостей п1,п2,...,пки набор проекций точек из множества Хна каждую из плоскостей Р Рп2,--,РПк. Сформулируем следующую обратную задачу: для заданных множеств точек проекций Р Р,2,...,РПк на плоскости п1,п2,...,пк восстановить исходное множество X = {x1, x2,..., xn}. Очевидным образом, исходное множество точек X для конкретного набора плоскостей и проекций не всегда может быть восстановимо. В частности, задача не имеет однозначного решения, если количество плоскостей равно 2. Без потери общности можно считать, что плоскости п J попарно непараллельны, так как в случае параллельных плоскостей множества проекций совпадают с точностью до параллельного переноса. Рассмотрим плоскость пj и проекции исходного множества точек на нее. Тогда ' - tJ 1nj = £1 J,1 Pj,1 + C1J,2РJ,2 + - + c1 J,mjPj,mj , x2 - tJ,2nj = c2j,1 PJ,1 + c2 j,2PJ,2 + - + c2 j,mj Pj,mj , xn - tJ,nnj = CnJ,1 pj ,1 + cnj,2 pj,2 + .•• + cnj,mjpj ,m j ■ Здесь nj - вектор нормали к плоскости пj , pj к - проекция некоторой точки x{ исходного множества X на плоскость пj . Неизвестные коэффициенты tj n е R , а cJj к таковы, что для каждой строки I все, кроме одного, равны 0, а один равен 1. Система (1) представляет собой однородную систему линейных уравнений из 3n уравнений и 3(2n + mj) неизвестных x, t, c (каждый со своим набором индексов). Записав систему (1) к раз (по числу проекций), получим систему уравнений (2), решения которой будут определять решение обратной задачи по восстановлению множества X через набор проекций на плоскости п1,п2,...,пк . Система (2) совместна по построению, т.е. решение существует, но, в общем случае, не единственно. Как показано на рис. 2, даже в двумерном случае (аналогичные примеры можно построить и для R3) без обладания какой-либо дополнительной информацией об исходном множестве X достаточно сложно, да и не всегда возможно, однозначно восстановить множество по его проекциям. Добавление третьей плоскости, на которую выполняется проекция исходного множества, расположенной под углом 45° к каждой плоскости (как показано на рис. 3), не уменьшает неопределенности - допустимы 2 конфигурации точек исходного множества. Изменив угол (сделав его отличным от 45 и 90°) между плоскостями, получим возможность однозначного восстановления множества по его проекциям. Х2 Х3 Рис. 3. Сохранение неопределенности x3 x1 Р1 Р4 x3 Р3 Р2 Рис. 2. Набор проекций p1,p2,p3,p4 на 2 координатные плоскости, определяющий 2 набора исходных точек x1, x2 , x3 и x1, x2 , x3 Таким образом, мы доказали следующую теорему: Теорема 1. Пусть X = {x1, x2,..., xn} - множество точек в R3, P , Pn,..., P проекции X на набор плоскостей п1, п2,...,nk (Pn = {p1, p2,..., pm}). Пусть также система (2) имеет единственное решение, такое, что tj n е R , все, кроме одного, c{j k для каждой строки l равны 0, а один равен 1. Тогда исходное множество X единственным образом определяется по множеству проекций. Замечание 1. Отметим, что теорема 1 явным образом не дает конструктивного алгоритма для вычисления исходного множества. Кроме того, априори требуется знание количества точек в исходном множестве. Практически более интересной является задача определения характеристик исходного множества без знания какой-либо дополнительной информации о нем. Как было уже замечено, обратная задача по восстановлению множества по набору проекций в общем случае не имеет единственного решения, поэтому переформулируем ее в несколько ином контексте. Определение. В условиях предыдущих обозначений назовем множество Y = {, y2,..., ys} подходящим, если проекция этого множества на плоскость пj совпадает с множеством Pj для любого J и любое расширение множества Y нарушает это условие. Тем самым множество Y является максимальным множеством, проекции которого на заданные плоскости пj совпадают с множествами P-. Очевидно, что X с Y . Дадим конструктивный алгоритм восстановления подходящего множества Y. Рассмотрим xq - некоторую точку из исходного множества X. Очевидным образом найдутся минимум 2 точки проекции pi и pj такие, что xq = Pi + n1t1 = Pj + n2^ (2) где n1 и n2 - нормали к плоскостям п1 и п2; t1 и t2 - некоторые неизвестные. Запишем теперь эти уравнения на t1 и t2 вида (2) для всех возможных pi и Pj . Решения этих уравнений определят множество Z = {pi + n1t} (t - всевозможные решения вышеописанных уравнений). Верно следующее включение X с Y с Z . Для получения множества Y из множества Z необходимо исключить те точки, которые не являются проекциями на все оставшиеся плоскости п3, —, пк . Чтобы упростить задачу исключения точек, попадающих в Z \ Y, достаточно применить аффинное преобразование, переводящее плоскости п1, п2 и п (3 < q < к) во взаимно ортогональные (координатные) плоскости. Тогда точка (x, y, z) будет принадлежать множеству Y, если точка (x, y, 0) принадлежит множеству проекций на п1, точка (x, 0, z) - проекций на п2, а точка (0, y, z) - проекций на п. В противном случае, ее нужно будет исключить из множества Z и продолжить ту же схему до перебора всех возможных вариантов. Замечание 1. Для построения исходного множества Z также можно использовать аффинное преобразование, приводящее плоскости п1 и п2 во взаимно ортогональные. Замечание 2. Последовательное исключение лишних точек из Z обеспечит большую скорость сходимости алгоритма, чем полный перебор. Компьютерная реализация алгоритма поиска также может быть хорошо оптимизирована. Тем самым мы предложили алгоритм восстановления подходящего множества для множества проекций. Более точное приближение множества Y к множеству X можно выполнять за счет увеличения количества проекций, что очень трудоемко в условиях реальных практических измерений. Заключение Рассмотрена обратная задача восстановления максимального множества точек в R3 по набору его двумерных проекций. Доказана теорема, позволяющая дать ответ на вопрос о достаточности изменений для однозначного определения множества по набору его двумерных проекций. Предлагается итеративный алгоритм для восстановления подходящего (максимального) множества по его набору проекций.
Хермен Г. Т. Восстановление изображений по проекциям. Основы реконструктивной томографии. М.: Мир, 1983.
Лаврентьев М.М., Зеркаль С.М., Трофимов О.Е. Численное моделирование в томографии и условно-корректные задачи. Новосибирск: Изд-во ИДМИ НГУ, 1999.
Raff el M., Willet C.E., Were ley S.T., Kompenhans J. Particle Image Velocimetry. Practical Guide. Berlin: Springer, 2007.
Elsinga G.E., Scarano F., Wieneke B., van Oudheusden B.W. Tomagraphic particle image ve-locimetry // Exp. Fluids. 2006. V. 41. P. 933-947.
Бильский А.В., Ложкин В.А., Маркович Д.М. и др. Оптимизация и тестирование томографического метода измерения скорости в объеме потока // Теплофизика и аэромеханика. 2011. Т. 18. № 4. C. 555-566.
Maas H.G., Gruen A.,Papantoniou,D. Particle tracking velocimetry in tree-dimentional flows // Exp. Fluids. 1993. V. 15. P. 133-146.
Novara M., Scarano F. Performances of motion tracking enhanced Tomo-PIV on turbulent shear flows // Exp. Fluids. 2012. V. 52. P. 1027-1041.
Wieneke B. Iterative reconstruction of volumetric particle distribution // Proc. 9th International Symposium on Particle Image Velocimetry, July 21-23, 2011, Kobe, Japan.