Перетворення Радона

Перетворення Радона. Карти f на (x, y) -область Rf на (α, з) -область.

Перетворення Радона для характеристичної функції двох квадратів, зображених на малюнку нижче. Більш світлі області вказують на більші значення функцій. Чорний колір позначає нуль.

Початкова функція дорівнює одиниці на білій області і нулю в темній області.

В математиці перетворення Радона — це інтегральне перетворення, яке переводить функцію f, визначену на площині, в функцію Rf, визначену на (двовимірному) просторі з прямих на площині, значення якої у певній прямій дорівнює криволінійному інтегралу від цієї функції над цією прямою.

Перетворення було введене в 1917 році Шаблон:НпШаблон:Sfn який також надав формулу зворотного перетворення. Крім того, Радон визначив формули для перетворення в трьох вимірах, в яких інтегрування відбувається на площині (інтегрування по лініях відоме як Шаблон:Нп). Пізніше це було узагальнено до евклідових просторів більш високого розміру, а ширше — в контексті інтегральної геометрії. Комплексний аналог перетворення Радона відомий як Шаблон:Нп. Перетворення Радона широко застосовується для томографії, створення зображення з проєкційних даних, пов'язаних із скануванням поперечного перерізу об'єкта.

Пояснення

Якщо функція $f$ представляє невизначену щільність, тоді перетворення Радона представляє дані проєкції, отримані як вихід томографічного сканування. Отже, зворотне перетворення Радона може бути використане для реконструкції початкової щільності з даних проєкції, і таким чином воно формує математичну основу Шаблон:Нп, також відому як ітеративна реконструкції.

Дані про перетворення Радона часто називають синограмою, оскільки перетворення Радона в центрі точкового джерела є синусоїдою. Отже, перетворення Радона для ряда дрібних предметів графічно постає як кількість розмитих синусоїд з різними амплітудами та фазами.

Перетворення Радона корисна при комп'ютерній томографії (КВТ-сканування), також застосовується у сканерах штрих-кодів, електронній мікроскопії, макромолекулярних комплексах, таких як віруси та білкові комплекси, рефлекторній сейсмології та при вирішенні гіперболічних часткових диференціальних рівнянь.

Визначення

Нехай ƒ(x) = ƒ(x, y) — функція, яка задовольняє трьом умовам регулярності:

1. ƒ(x, y) — неперервна
2. подвійний інтеграл $\int \int \frac{|f(x,y)|}{\sqrt{x^2+y^2}}dxdy$, що визначений на всій площині, збігається
3. для будь-якої довільної точки $(x,y)$ на площині виконується, що $\underset{r\to \mathrm{\infty }}{\lim }\underset{0}{\overset{2\pi }{\int }}f(x+r\cos \varphi ,y+r\sin \varphi )d\varphi =0$

Перетворення Радона, Rƒ, є функцією, визначеною на просторі прямих L в просторі R² криволінійно інтегрованою вздовж кожної такої прямої, як

$Rf(L)=\underset{L}{\int }f(𝐱)|d𝐱|.$

Конкретніше, параметризацію будь-якої прямої L щодо довжини дуги z завжди можна записати як

$(x(z),y(z))=((z\sin \alpha +s\cos \alpha ),(-z\cos \alpha +s\sin \alpha ))$

де s — відстань L від початку і $\alpha$ — кут, який нормальний до L вектор утворює з віссю х. Звідси випливає, що величини (α, s) можна вважати координатами на просторі всіх ліній в R², а перетворення Радона в цих координатах визначають як

$\begin{array}{cc}Rf(\alpha ,s)& =\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(x(z),y(z))dz\\ & =\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f((z\sin \alpha +s\cos \alpha ),(-z\cos \alpha +s\sin \alpha ))dz\end{array}$

Більш загальне визначення

Нехай $f$ - функція на ${ℝ}^n$, інтегровувана по кожній гіперплощині. Нехай також ${\mathbb{P}}^n$ простір усіх гіперплощин у ${ℝ}^n$, наділене відповідною топологією. Перетворення Радона функції $f$ визначається як функція $\hat{f}$ на ${\mathbb{P}}^n$, і задається формулою

$\hat{f}(\xi )=\underset{\xi }{\int }f(x)dm(x),$

де $dm$ - евклідова міра на гіперплощині $\xi .$

Гіперплощину $\xi \in {\mathbb{P}}^n$ можна записати у вигляді $\xi =\{x\in {\mathbb{P}}^n|(x,\omega )=p\},$ де $(\cdot ,\cdot )$ - звичайний скалярний добуток, $\omega =({\omega }_1,...,{\omega }_n)$ - одиничний вектор та $p\in ℝ.$ Пари $(\omega ,p)$ та $(-\omega ,-p)$ приводять до однієї і тої самої гіперплощини $\xi :$ відображення $(\omega ,\pi )\to \xi$ є двократне накриття ${𝕊}^{n-1}\times ℝ$ на ${\mathbb{P}}^n$. Таким чином, простір ${\mathbb{P}}^n$ має канонічну структуру многовиду, відносно якої це накриваюче відображення диференційовуване й регулярне.

Зв'язок з перетворенням Фур'є

Обчислення двовимірного перетворення Радона через два перетворення Фур'є.

Перетворення Радона тісно пов'язане з перетворенням Фур'є. Ми визначаємо тут універсальне перетворення Фур'є як

$\hat{f}(\omega )=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(x)e^{-2\pi ix\omega }dx.$

і для функції для 2-х векторів $𝐱=(x,y)$,

$\hat{f}(𝐰)=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(𝐱)e^{-2\pi i𝐱\cdot 𝐰}dxdy.$

Для зручності позначте ${ℛ}_{\alpha }[f](s)=ℛ[f](\alpha ,s)$. Тоді Шаблон:Нп констатує

$\widehat{{ℛ}_{\alpha }[f]}(\sigma )=\hat{f}(\sigma 𝐧(\alpha ))$

де

$𝐧(\alpha )=(\cos \alpha ,\sin \alpha ).$

Таким чином, двовимірне перетворення Фур'є початкової функції вздовж лінії під кутом нахилу $\alpha$ — це перетворення Фур'є для однієї змінної над перетворенням Радона (отримане під кутом нахила $\alpha$) цієї функції. Цей факт може бути використаний для обчислення як перетворення Радона, так і зворотньої функції.

Результат можна узагальнити для n-мірного випадку

$\hat{f}(r\alpha )=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}ℛf(\alpha ,s)e^{-2\pi isr}ds.$

Подвійне перетворення

Подвійне перетворення Радона є своєрідним примиканням до перетворення Радона. Починаючи з функції g на просторі Σ_n, подвійне перетворення Радона — функція ${ℛ}^{*}g$ на Rⁿ, яка визначається в такий спосіб

${ℛ}^{*}g(x)=\underset{x\in \xi }{\int }g(\xi )d\mu (\xi ).$

Інтеграл тут береться за множину всіх гіперплощин, що проходять через точку х ∈ Rⁿ, а міра d є унікальною мірою ймовірності на множині $\{\xi |x\in \xi \}$ є інваріантом при обертанні навколо точки x.

Конкретно, для двовимірного перетворення Радона подвійне перетворення задається таким чином

${ℛ}^{*}g(x)=\frac{1}{2\pi }\underset{\alpha =0}{\overset{2\pi }{\int }}g(\alpha ,𝐧(\alpha )\cdot 𝐱)d\alpha .$

У контексті обробки зображень подвійне перетворення зазвичай називають зворотньою проєкцієюШаблон:Sfn оскільки воно приймає функцію, визначену в кожній лінії в площині, і «розмазує» або проектує її назад по лінії для створення зображення.

Переплетення власності

Нехай Δ позначає оператор Лапласа на R ⁿ:

$\mathrm{\Delta }=\frac{{\partial }^2}{\partial x_1^2}+\cdots +\frac{{\partial }^2}{\partial x_n^2}.$

Це природний обертальний інваріантний диференціальний оператор другого порядку. На Σ_n, «радіальна» друга похідна

$Lf(\alpha ,s)\equiv \frac{{\partial }^2}{\partial s^2}f(\alpha ,s)$

також інваріантна відносно обертання. Перетворення Радона та його подвійний є переплетеними операторами для цих двох диференціальних операторів у тому сенсі, щоШаблон:Sfn

$ℛ(\mathrm{\Delta }f)=L(ℛf),{ℛ}^{*}(Lg)=\mathrm{\Delta }({ℛ}^{*}g).$

Аналізуючи рішення хвильового рівняння в декількох просторових вимірах, властивість переплетення призводить до поступального представлення Лакса і Філіпса. У візуалізації та чисельному аналізі це використовується для зменшення багатовимірних задач на одновимірні, як метод розмірного розщеплення.

Підходи до реконструкції

Процес реконструкції створює зображення (або функцію $f$ у попередньому розділі) з його проєкційних даних. Реконструкція — обернена задача.

Формула інверсії радону

У двовимірному випадку найбільш часто використовується аналітична формула для відновлення $f$ знаючи його перетворення Радона за допомогою формули відфільтрованої зворотньої проєкції або формули інверсії Радона:

$f(𝐱)=\underset{0}{\overset{\pi }{\int }}(ℛf(\cdot ,\theta )*h)(⟨𝐱,{𝐧}_{\theta }⟩)d\theta$Шаблон:Sfn

де $h$ таке, що $\hat{h}(k)=|k|$Шаблон:Sfn.

Ядро згортки $h$ в деякій літературі згадується як фільтр Рампа.

Нерегулярна похибка

Інтуїтивно, у відфільтрованій формулі зворотного проектування, за аналогією з диференціацією, для якої $\left(\widehat{\frac{d}{dx}f}\right)(k)=ik\hat{f}(k)$, ми бачимо, що фільтр виконує операцію, аналогічну операції взяття похідної. Грубо кажучи, тоді фільтр робить об'єкти більш сингулярними.

Кількісне твердження про нерегулярну похибку радонової інверсії полягає в наступному: Ми маємо $\widehat{{ℛ}^{*}ℛg}(k)=\frac{1}{||𝐤||}\hat{g}(𝐤)$

де ${ℛ}^{*}$ є раніше визначеним Шаблон:Нп до Радонового перетворення.

Таким чином для $g(𝐱)=e^{i⟨{𝐤}_0,𝐱⟩}$,

${ℛ}^{*}ℛg=\frac{1}{||{𝐤}_{\mathrm{𝟎}}||}{e}^{i⟨{𝐤}_0,𝐱⟩}$.

Складний показник $e^{i⟨{𝐤}_0,𝐱⟩}$ таким чином, є власною функцією ${ℛ}^{*}ℛ$ із власним значенням $\frac{1}{||{𝐤}_{\mathrm{𝟎}}||}$. Таким чином, сингулярні значення $ℛ$ є $\sqrt{\frac{1}{||𝐤||}}$. Оскільки ці особливі значення мають прямувати до 0, ${ℛ}^{-1}$ є необмеженим.Шаблон:Sfn

Ітеративні методи реконструкції

Порівняно з методом відфільтрованого зворотного проектування, ітеративна реконструкція коштує великих витрат на обчислення, обмежуючи її практичне використання. Однак через недоброзичливість інверсії радону метод фільтруваної зворотньої проєкції може виявитися нездійсненним при наявності розриву або шуму. Методи ітеративної реконструкції (наприклад, Шаблон:Нп) можуть забезпечити зменшення артефактів металу, зменшення шуму та дози для реконструкції результату, що привертає великий науковий інтерес по всьому світу.

Формули інверсії

Явні та обчислювально ефективні формули інверсії для перетворення Радона та його подвійності. Перетворення Радона в n розмірах може бути інвертованим формулоюШаблон:Sfn

$c_{n}f=(-\mathrm{\Delta }{)}^{(n-1)/2}{R}^{*}Rf$

де

$c_n=(4\pi {)}^{(n-1)/2}\frac{\mathrm{\Gamma }(n/2)}{\mathrm{\Gamma }(1/2)}.$

а потужність Лаплаціана (−Δ)^(n−1)/2 визначається як псевдодиференціальний оператор або, при необхідності, перетворенням Фур'є

$[ℱ(-\mathrm{\Delta }{)}^{(n-1)/2}\varphi ](\xi )=|2\pi \xi {|}^{n-1}(ℱ\varphi )(\xi ).$

Для обчислювальних цілей, потужність Лаплаціани змішується з подвійним перетворенням R^*, щоб отриматиШаблон:Sfn

$c_{n}f=\{\begin{array}{cc}{R}^{*}\frac{d^{n-1}}{d{s}^{n-1}}Rf& n\mathrm{o}\mathrm{d}\mathrm{d}\\ R^{*}{H}_s\frac{d^{n-1}}{d{s}^{n-1}}Rf& n\mathrm{e}\mathrm{v}\mathrm{e}\mathrm{n}\end{array}$

де H_s — перетворення Гільберта відносно змінної s. У двох вимірах оператор H_sd/ds з'являється в обробці зображень як рамповий фільтрШаблон:Sfn. Можна легко довести з теореми Фур'є про зріз і зміни змінних для інтеграції, що для компактно підтримуваної безперервної функції ƒ для двох змінних

через

$f=\frac{1}{2}{R}^{*}{H}_s\frac{d}{ds}Rf.$

Таким чином, в контексті обробки зображень вихідне зображення ƒ може бути відновлено з даних «синограми» Rƒ шляхом застосування рампового фільтру (у $s$ змінній), а потім зворотнього проектування. Оскільки етап фільтрації може бути виконаний ефективно (наприклад, використовуючи методи цифрової обробки сигналу), а один крок зворотнього проектування є просто накопиченням значень у пікселях зображення, що призводить до високоефективного, а отже, широко використовуваного алгоритму.

Явна формула інверсії, отримана останнім методомШаблон:Sfn

$f(x)=\frac{1}{2}(2\pi {)}^{1-n}(-1{)}^{(n-1)/2}\underset{S^{n-1}}{\int }\frac{{\partial }^{n-1}}{\partial s^{n-1}}Rf(\alpha ,\alpha \cdot x)d\alpha$

якщо n непарне, і

$f(x)=(2\pi {)}^{-n}(-1{)}^{n/2}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{1}{q}\underset{S^{n-1}}{\int }\frac{{\partial }^{n-1}}{\partial s^{n-1}}Rf(\alpha ,\alpha \cdot x+q)d\alpha dq$

якщо n парне.

Подвійне перетворення також може бути обернено аналогічною формулою:

$c_{n}g=(-L{)}^{(n-1)/2}R(R^{*}g).$

Перетворення Радона в алгебричній геометрії

В алгебраїчній геометрії перетворення Радона (також відоме як перетворення Брилінського-Радона) будується наступним чином.

Пишуть

${𝐏}^d\stackrel{p_1}{\leftarrow }H\stackrel{p_2}{\to }{𝐏}^{\vee ,d}$

для універсальної гіперплощини, тобто H складається з пар (x, h), де x — точка в d-вимірному проективному просторі ${𝐏}^d$ і h — точка у подвійному проєкційному просторі (іншими словами, x є лінією через початок у (d+1)-вимірному афінному просторі, а h — гіперплощина у цьому просторі), такий, що x міститься в h.

Тоді перетворення Брилінкського — Радона є функтором між відповідними Шаблон:Нп Шаблон:Нп

$Rad:=R{p}_{2,*}{p}_1^{*}:D({𝐏}^d)\to D({𝐏}^{\vee ,d}).$

Основна теорема про це перетворення полягає в тому, що це перетворення проектує еквівалентність категорій Шаблон:Нп на проективному просторі та його подвійному проективному просторі аж до постійних пучків.

Застосування

Шаблон:Нп є найбільш широко використовуваним особливим випадком цієї конструкції, і її отримують шляхом інтеграції з прямими лініями.

Див. також

• Шаблон:Нп
• Шаблон:Нп
• Деконволюція
• Шаблон:Нп
• Шаблон:Нп
• Перетворення Гафа, коли воно пишеться у безперервній формі, дуже схоже, якщо не еквівалентне, перетворенню Радона.Шаблон:Sfn
• Теорема Коші-Крофтона є тісно пов'язаною формулою для обчислення довжини кривих у просторі.
• Швидке перетворення Фур'є

Список літератури

• Шаблон:Citation; Translation: Шаблон:Citation.
• Шаблон:Springer.
• Шаблон:Citation.
• Шаблон:Cite web
• Шаблон:Cite web
• Шаблон:Cite web
• Шаблон:Cite web

Подальше читання

• Шаблон:Cite book
• Шаблон:Citation
• Шаблон:Citation
• Шаблон:Citation
• Шаблон:Springer
• Шаблон:Citation
• Шаблон:Citation
• Шаблон:Citation

Посилання

• Шаблон:MathWorld
• Шаблон:Cite video