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

Перетворення Радона
Названо на честь	Johann Radond
Формула	${\displaystyle \mathrm {R} f(L)=\int _{L}f(\mathbf {x} )\,|\mathrm {d} \mathbf {x} |}$
Підтримується Вікіпроєктом	Вікіпедія:Проєкт:Математика
Перетворення Радона у Вікісховищі

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

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

Якщо функція ${\displaystyle f}$ представляє невизначену щільність, тоді перетворення Радона представляє дані проєкції, отримані як вихід томографічного сканування. Отже, зворотне перетворення Радона може бути використане для реконструкції початкової щільності з даних проєкції, і таким чином воно формує математичну основу томографічної реконструкції^[en], також відому як ітеративна реконструкції.

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

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

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

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

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

${\displaystyle Rf(L)=\int _{L}f(\mathbf {x} )\,|d\mathbf {x} |.}$

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

${\displaystyle (x(z),y(z))={\Big (}(z\sin \alpha +s\cos \alpha ),(-z\cos \alpha +s\sin \alpha ){\Big )}\,}$

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

${\displaystyle {\begin{aligned}Rf(\alpha ,s)&=\int _{-\infty }^{\infty }f(x(z),y(z))\,dz\\&=\int _{-\infty }^{\infty }f{\big (}(z\sin \alpha +s\cos \alpha ),(-z\cos \alpha +s\sin \alpha ){\big )}\,dz\end{aligned}}}$

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

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

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

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

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

${\displaystyle {\hat {f}}(\omega )=\int _{-\infty }^{\infty }f(x)e^{-2\pi ix\omega }\,dx.}$

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

${\displaystyle {\hat {f}}(\mathbf {w} )=\int \limits _{-\infty }^{\infty }\int \limits _{-\infty }^{\infty }f(\mathbf {x} )e^{-2\pi i\mathbf {x} \cdot \mathbf {w} }\,dx\,dy.}$

Для зручності позначте ${\displaystyle {\mathcal {R}}_{\alpha }[f](s)={\mathcal {R}}[f](\alpha ,s)}$. Тоді теорема Фур'є про зріз функції^[en] констатує

${\displaystyle {\widehat {{\mathcal {R}}_{\alpha }[f]}}(\sigma )={\hat {f}}(\sigma \mathbf {n} (\alpha ))}$

де

${\displaystyle \mathbf {n} (\alpha )=(\cos \alpha ,\sin \alpha ).}$

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

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

${\displaystyle {\hat {f}}(r\alpha )=\int _{-\infty }^{\infty }{\mathcal {R}}f(\alpha ,s)e^{-2\pi isr}\,ds.}$

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

${\displaystyle {\mathcal {R}}^{*}g(x)=\int _{x\in \xi }g(\xi )\,d\mu (\xi ).}$

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

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

${\displaystyle {\mathcal {R}}^{*}g(x)={\frac {1}{2\pi }}\int _{\alpha =0}^{2\pi }g(\alpha ,\mathbf {n} (\alpha )\cdot \mathbf {x} )\,d\alpha .}$

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

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

${\displaystyle \Delta ={\frac {\partial ^{2}}{\partial x_{1}^{2}}}+\cdots +{\frac {\partial ^{2}}{\partial x_{n}^{2}}}.}$

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

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

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

${\displaystyle {\mathcal {R}}(\Delta f)=L({\mathcal {R}}f),\quad {\mathcal {R}}^{*}(Lg)=\Delta ({\mathcal {R}}^{*}g).}$

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

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

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

${\displaystyle f(\mathbf {x} )=\int _{0}^{\pi }({\mathcal {R}}f(\cdot ,\theta )*h)(\left\langle \mathbf {x} ,\mathbf {n} _{\theta }\right\rangle )d\theta }$^[9]

де ${\displaystyle h}$ таке, що ${\displaystyle {\hat {h}}(k)=|k|}$^[10].

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

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

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

де ${\displaystyle {\mathcal {R}}^{*}}$ є раніше визначеним примиканням^[en] до Радонового перетворення.

Таким чином для ${\displaystyle g(\mathbf {x} )=e^{i\left\langle \mathbf {k} _{0},\mathbf {x} \right\rangle }}$,

${\displaystyle {\mathcal {R}}^{*}{\mathcal {R}}g={\frac {1}{||\mathbf {k_{0}} ||}}e^{i\left\langle \mathbf {k} _{0},\mathbf {x} \right\rangle }}$.

Складний показник ${\displaystyle e^{i\left\langle \mathbf {k} _{0},\mathbf {x} \right\rangle }}$ таким чином, є власною функцією ${\displaystyle {\mathcal {R}}^{*}{\mathcal {R}}}$ із власним значенням ${\displaystyle {\frac {1}{||\mathbf {k_{0}} ||}}}$. Таким чином, сингулярні значення ${\displaystyle {\mathcal {R}}}$ є ${\displaystyle {\sqrt {\frac {1}{||\mathbf {k} ||}}}}$. Оскільки ці особливі значення мають прямувати до 0, ${\displaystyle {\mathcal {R}}^{-1}}$ є необмеженим.^[10]

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

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

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

де

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

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

${\displaystyle \left[{\mathcal {F}}(-\Delta )^{(n-1)/2}\phi \right](\xi )=|2\pi \xi |^{n-1}({\mathcal {F}}\phi )(\xi ).}$

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

${\displaystyle c_{n}f={\begin{cases}R^{*}{\frac {d^{n-1}}{ds^{n-1}}}Rf&n{\rm {\ odd}}\\R^{*}H_{s}{\frac {d^{n-1}}{ds^{n-1}}}Rf&n{\rm {\ even}}\end{cases}}}$

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

через

${\displaystyle f={\frac {1}{2}}R^{*}H_{s}{\frac {d}{ds}}Rf.}$

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

Явна формула інверсії, отримана останнім методом^[4]

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

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

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

якщо n парне.

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

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

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

Пишуть

${\displaystyle \mathbf {P} ^{d}{\stackrel {p_{1}}{\gets }}H{\stackrel {p_{2}}{\to }}\mathbf {P} ^{\vee ,d}}$

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

Тоді перетворення Брилінкського — Радона є функтором між відповідними похідними категоріями^[en] еталевих пучків^[en]

${\displaystyle Rad:=Rp_{2,*}p_{1}^{*}:D(\mathbf {P} ^{d})\to D(\mathbf {P} ^{\vee ,d}).}$

Основна теорема про це перетворення полягає в тому, що це перетворення проектує еквівалентність категорій перекручених пучків^[en] на проективному просторі та його подвійному проективному просторі аж до постійних пучків^[15].

Рентгенівська трансформація^[en] є найбільш широко використовуваним особливим випадком цієї конструкції, і її отримують шляхом інтеграції з прямими лініями.

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

• Lokenath Debnath; Dambaru Bhatta (19 квітня 2016). Integral Transforms and Their Applications. CRC Press. ISBN 978-1-4200-1091-6.
• Deans, Stanley R. (1983), The Radon Transform and Some of Its Applications, New York: John Wiley & Sons
• Helgason, Sigurdur (2008), Geometric analysis on symmetric spaces, Mathematical Surveys and Monographs, т. 39 (вид. 2nd), Providence, R.I.: American Mathematical Society, doi:10.1090/surv/039, ISBN 978-0-8218-4530-1, MR 2463854
• Herman, Gabor T. (2009), Fundamentals of Computerized Tomography: Image Reconstruction from Projections (вид. 2nd), Springer, ISBN 978-1-85233-617-2
• Hazewinkel, Michiel, ред. (2001), transform Radon transform, Математична енциклопедія, Springer, ISBN 978-1-55608-010-4
• Natterer, Frank (June 2001), The Mathematics of Computerized Tomography, Classics in Applied Mathematics, т. 32, Society for Industrial and Applied Mathematics, ISBN 0-89871-493-1
• Natterer, Frank; Wübbeling, Frank (2001), Mathematical Methods in Image Reconstruction, Society for Industrial and Applied Mathematics, ISBN 0-89871-472-9
• Kiehl, Reinhardt; Weissauer, Rainer (2001), Weil conjectures, perverse sheaves and l'adic Fourier transform, Springer, doi:10.1007/978-3-662-04576-3, ISBN 3-540-41457-6, MR 1855066

• Weisstein, Eric W. Перетворення Радона(англ.) на сайті Wolfram MathWorld.
• Analytical projection (the Radon transform) (video). 10 вересня 2015.