**بازسازی برش‌نگاری** یا **بازسازی توموگرافی** (به انگلیسی: Tomographic reconstruction) یکی از انواع روش‌های حل معکوس چند‌بعدی است که در آن مسئله برآورد یک سیستم معین با استفاده از تعداد محدودی از پروجکشن‌ها می‌باشد. پایه‌های ریاضی برای تصویربرداری توموگرافی توسط یوهان رادون مطرح شد. یکی از نمونه‌های کاربردی قابل توجه، بازسازی توموگرافی کامپیوتری (CT) است که در آن تصاویر مقطعی بیماران به روش غیر تهاجمی به دست می‌آیند. در پیشرفت‌های اخیر، تبدیل رادون و معکوس آن برای کارهای مربوط به جابجایی جسم واقع گرایانه مورد نیاز برای آزمایش و ارزیابی استفاده از توموگرافی کامپیوتری در امنیت فرودگاه‌ها، استفاده شده‌است.

این مقاله به‌طور کلی به روش‌های بازسازی برای انواع توموگرافی اشاره می‌کند، اما برخی از اصطلاحات و توصیف‌های فیزیکی به‌طور مستقیم به بازسازی توموگرافی کامپیوتری اشعه ایکس، اشاره می‌کنند.

معرفی فرمول
-----------

شکل یک\_ هندسه اشعه موازی استفاده شده در توموگرافی و بازسازی توموگرافی. هر پروجکشن که از توموگرافی تحت یک زاویه خاص به دست آمده، از یک سری از انتگرال‌های خط در طول جسم ساخته شده‌است.
 شکل دوم\_ تصویر توموگرافی حاصل شده از فانتوم جمجمه پلاستیکی. اشعه ایکس‌های فرافکنده شده به وضوح بر روی این مقطع که با استفاده از یک سی تی اسکن گرفته شده مشخص است علاوه بر این اعوجاج‌هایی نیز به دلیل تعداد محدود پروجکشن در زاویه‌ها مشاهده شده‌است.
پروجکشن یک شی که از فرایند اندازه‌گیری توموگرافی در یک زاویه معین \theta  ایجاد شده، از مجموعه‌ای از انتگرال خط ساخته می‌شود (به شکل ۱توجه کنید). مجموعه ای شامل تعداد زیادی از چنین پروجکشن‌هایی تحت زوایای مختلف که در دو بعد (2D) سازماندهی شده‌اند، سینوگرام نامیده می‌شود (به شکل ۳ نگاه کنید). درتوموگرافی کامپیوتری اشعه ایکس، انتگرال خط، نشان دهنده میزان کلی تضعیف پرتو اشعه ایکس می‌باشد که وقتی اشعه در مسیر مستقیم از داخل شیء حرکت می‌کند، گویی این انتگرال در حال تشکیل شدن است. همان‌طور که در بالا ذکر شد، تصویر حاصل یک مدل دوبعدی (یا سه‌بعدی) از ضریب تضعیف است. هدف ما به تصویر کشیدن {\displaystyle \mu (x,y)} است. ساده‌ترین و آسان‌ترین روش برای تجسم روش اسکن، سیستم پروجکشن‌های موازی است، که در اسکنرهای اولیه استفاده شد. برای این بحث، ما فرض می‌کنیم که داده‌ها از یک سری اشعه موازی، در طول یک پروجکشن در زاویه \theta  جمع‌آوری شده‌اند. این روند برای زوایای مختلف تکرار شده‌است. تضعیف به صورن نمایی در بافت اتفاق می‌افتد:

{\displaystyle I=I\_{0}\exp \left({-\int \mu (x,y)\,ds}\right)}
که در آن {\displaystyle \mu (x,y)} ضریب تضعیف به عنوان تابعی از مکان است؛ بنابراین، تضعیف کلی pیک پرتو در یک موقعیت، در پروچکشن ای با زاویه \theta ، توسط انتگرال خط زیر به دست می‌آید:

{\displaystyle p\_{\theta }(r)=\ln \left({\frac {I}{I\_{0}}}\right)=-\int \mu (x,y)\,ds}

با استفاده از سیستم مختصات نشان داده شده در شکل یک، مقدار rکه نقطه {\displaystyle f(x,y)}روی آن تصویر می‌شود، از رابطه زیر به دست می‌آید:

{\displaystyle xcos\theta +ysin\theta =r}

بنابراین معادله بالا را می‌توان به صورت زیر بازنویسی کرد:

{\displaystyle p\_{\theta }(r)=\int \_{-\infty }^{\infty }\int \_{-\infty }^{\infty }f(x,y)\delta (x\cos \theta +y\sin \theta -r)\,dx\,dy}
که در آن {\displaystyle f(x,y)} نشان دهنده {\displaystyle \mu (x,y)}است. این تابع به عنوان تبدیل رادون (یا *سینوگرام*) از شی دوبعدی شناخته می‌شود.

تبدیل فوریه پروجکشن را می‌توان به صورت زیر نوشت:

{\displaystyle P\_{\theta }(\omega )=\int \_{-\infty }^{\infty }\int \_{-\infty }^{\infty }f(x,y)\exp[-j\omega (x\cos \theta +y\sin \theta )]\,dx\,dy=F(\Omega \_{1},\Omega \_{2})}که در آن {\displaystyle \Omega \_{1}=\omega \cos \theta ,\Omega \_{2}=\omega \sin \theta }

{\displaystyle P\_{\theta }(\omega )} نشان دهنده یک مقطع از تبدیل فوریه دوبعدی {\displaystyle f(x,y)} در زاویه \theta است. با استفاده از معکوس تبدیل فوریه، فرمول معکوس تبدیل رادون می‌تواند به راحتی مشتق شود.

{\displaystyle f(x,y)={\frac {1}{2\pi }}\int \limits \_{0}^{\pi }g\_{\theta }(x\cos \theta +y\sin \theta )d\theta }

که در آن {\displaystyle g\_{\theta }(x\cos \theta +y\sin \theta )} مشتق از تبدیل هیلبرت {\displaystyle p\_{\theta }(r)} است.

در تئوری معکوس تبدیل رادون، تصویر اصلی را تولید می‌کند. قضیه پروجکشن مقطعی به ما می‌گوید که اگر ما تعداد بی‌نهایت از پروجکشن‌های یک بعدی جسم داشته باشیم که در بی‌نهایت زاویه گرفته شده‌اند، می‌توانیم شیء اصلی یعنی همان {\displaystyle f(x,y)}را به‌طور کامل بازسازی کنیم. با این حال، در عمل تنها تعداد محدودی پروجکشن در دسترس خواهد بود.

با فرض اینکه {\displaystyle f(x,y)} دارای قطر مؤثر dاست و حد تفکیک مورد نظر {\displaystyle R\_{s}} است، تعداد پروجکشن مورد نیاز برای بازسازی تصویر از رابطه زیر به دست می‌آید: {\displaystyle N>\pi d/R\_{s}}

الگوریتم‌های بازسازی
--------------------

**الگوریتم‌های بازسازی** عملی برای تکمیل روند بازسازی یک جسم سه‌بعدی از پروجکشن‌های آن ساخته شده‌است. این الگوریتم‌ها عمدتاً بر اساس ریاضیات تبدیل رادون، دانش آماری فرایندجمع آوری داده و هندسه سیستم تصویربرداری داده طراحی شده‌اند.

### الگوریتم بازسازی در حوزه فوریه

بازسازی را می‌توان با استفاده از درون یابی انجام داد. فرض کنید N پروجکشن از {\displaystyle f(x,y)}‌در زوایای که به‌طور مساوی فاصله داده می‌شوند و هر کدام سرعت نمونه برداری یکسان دارند، تولید شده‌است. تبدیل فوریه گسسته در هر پروجکشن، نمونه‌برداری را در حوزه فرکانس انجام خواهد داد. ترکیب همه پروجکشن‌های نمونه‌برداری شده در فرکانس، یک خط قطبی در حوزه فرکانس ایجاد می‌کند. خط قطبی نامتراکم خواهد بود بنابرین برای پر کردن نقاط ناشناخته تبدیل فوریه دوبعدی از درون‌یابی استفاده می‌شود و سپس بازسازی می‌تواند از طریق معکوس تبدیل فوریه دیجیتال انجام شود. عملکرد بازسازی ممکن است با طراحی روش‌هایی برای تغییر نامتراکمی خط قطبی و تسهیل اثرگذاری درون‌یابی، بهبود یابند. به عنوان مثال، یک مربع مورب متحدالمرکز در حوزه فرکانس را می‌توان با تغییر زاویه بین هر پروجکشن به صورت زیر بدست آورد:

{\displaystyle \theta '={\frac {R\_{0}}{\max\{|\cos \theta |,|\sin \theta |\}}}}که {\displaystyle R\_{0}}بالاترین فرکانس مورد بررسی است.

مربع مورب متحدالمرکز، کارایی محاسباتی را بهبود می‌بخشد با اینکه تمام موقعیت‌های درون یابی روی شبکه تبدیل فوربه گسسته مستطیلی شکل می‌گیرد. علاوه بر این، خطای درون یابی را نیز کاهش می‌دهد. با این حال، الگوریتم تبدیل فوریه مشکل بزرگ تولید خروجی نویزی را دارد.

### الگوریتم عقبگرد پروجکشن

در عمل بازسازی تصویر توموگرافیک، اغلب از یک نسخه پایدار و گسسته سازی شده تبدیل رادون معکوس استفاده می‌شود، که به عنوان الگوریتم عقبگرد پروجکشن شناخته می‌شود.

در یک سیستم نمونه‌گیری شده گسسته، تبدیل رادون معکوس، مانند زیر است:

{\displaystyle f(x,y)={\frac {1}{2\pi }}\sum \_{i=0}^{N-1}\Delta \theta \_{i}g\_{i}(x\cos \theta \_{i}+y\sin \theta \_{i})}

{\displaystyle g\_{\theta }(t)=p\_{\theta }(t)\cdot k(t)}

که {\displaystyle \Delta \theta } فاصله زاویه ای بین پروجکشن‌ها و {\displaystyle k(t)} کرنل رادون با پاسخ فرکانس {\displaystyle |\omega |}است. نام عقبگرد پروجکشن از این واقعیت سرچشمه گرفته‌است که پروجکشن‌های یک بعدی باید توسط کرنل یک بعدی رادون فیلتر شوند تا بتوان سیگنال دوبعدی را به دست آورد. فیلتر استفاده شده حاوی بهره مقدار ثابت (gain DC) نیست، بنابراین افزودن یک مقدار ثابت (DC) ممکن است مطلوب باشد. بازسازی با استفاده از روش عقبگرد پروجکشن، باعث وضوح بیشتری نسبت به روش درون یابی که در بالا توضیح داده شده‌است، می‌شود. با این حال، باعث ایجاد نویز بیشتر نیز می‌شود، زیرا فیلتر مستعد به افزایش محتوای فرکانس بالا است.

### الگوریتم بازسازی تکراری

الگوریتم تکراری از لحاظ محاسبات پیچیده‌است اما اجازه می‌دهد اطلاعات *پیشین* در مورد سیستم به حساب بیایند. Nتعداد پروجکشن‌ها و {\displaystyle D\_{i}}اپراتور تحریف برای iامین پروجکشن گرفته شده در یک زاویه \theta \_{i}است و {\displaystyle \{\lambda \_{i}\}}مجموعه ای از پارامترها برای بهینه‌سازی تبدیل تکرارها هستند.

{\displaystyle f\_{0}(x,y)=\sum \_{i=1}^{N}\lambda \_{i}p\_{\theta \_{i}}(r)}

{\displaystyle f\_{k}(x,y)=f\_{k-1}(x,y)+\sum \_{i=1}^{N}\lambda \_{i}[p\_{\theta \_{i}}(r)-D\_{i}f\_{k-1}(x,y)]}

یک خانواده جایگزین از الگوریتم‌های بازسازی توموگرافی بازگشتی عبارتند از: تکنیک بازسازی جبری و حداقل تفاضل مجانبی. 

 بازسازی پرتو پنکه ای شپ-لوگان فانتوم با فاصله‌های مختلف سنسور. فاصله کمتر بین سنسورها اجازه بازسازی دقیق تر انجام شود. این شکل با استفاده از MATLAB تولید شد.
### بازسازی پرتو پنکه‌ای

استفاده از یک پرتو پنکه ای هم راستا سازی نشده رایج است زیرا به دست آوردن یک پرتو تابش هم راستا سازی شده دشوار است. پرتوهای پنکه ای، مجموعه ای از انتگرال‌های خط را تولید می‌کنند که مانند پروجکشن‌ها به صورت موازی با یکدیگر نیستند. سیستم پرتو پنکه ای نیاز به محدوده زاویه۳۶۰ درجه ای دارد که موجب محدودیت مکانیکی می‌شود، اما اجازه می‌دهد تا سرعت جمع‌آوری سیگنال سریع تر باشد که این امر در موارد خاص برای مثال در زمینه پزشکی مفید است. روند عقبگرد پروجکشن یک روش دو مرحله ای مشابه را طی می‌کند که با استفاده از محاسبه جمع وزن دار شده عقبگرد پروجکشن‌های گرفته شده از پروجکشن‌های فیلتر شده، بازسازی را انجام می‌دهد.

نرم‌افزار بازسازی توموگرافی
---------------------------

برای بازسازی توموگرافی انعطاف‌پذیر، جعبه ابزارهای منبع بازمانند TomoPy, ODL یا جعبه ابزار ASTRA. در دسترس هستند. TomoPy یک جعبه ابزار منبع باز پایتون برای انجام پردازش داده‌های توموگرافی و امور بازسازی تصویر در منبع پیشرفته فوتون در آزمایشگاه ملی Argonne است. جعبه ابزار TomoPy به‌طور خاصی طراحی شده‌است که برای استفاده آسان است. TomoPy همچنین شامل چندین الگوریتم بازسازی است که می‌تواند در ایستگاه‌های کاری چند هسته‌ای و با امکانات محاسباتی در مقیاس بزرگ اجرا شود. جعبه ابزار ASTRA یک جعبه ابزار MATLAB از پایه‌های پردازنده گرافیکی با عمکرد خوب برای توموگرافی دوبعدی و سه‌بعدی است که ازسال ۲۰۰۹–۲۰۱۴ توسط آزمایشگاه iMinds-Vision Lab University of Antwerp ساخته شده و از سال ۲۰۱۴ به‌طور مشترک توسط iMinds-VisionLab, UAntwerpen و CWI, Amsterdam پیشرفت داده شده‌است. این جعبه ابزار از پرتوهای موازی، پنکه ای و مخروطی با موقعیت منبع/ آشکارساز بسیار قابل انعطاف پشتیبانی می‌کند. تعداد زیادی از الگوریتم‌های بازسازی از طریق مجموعه ابزارTomoPyو ASTRA در دسترس هستند، از جمله این‌ها می‌توان بهFBP, Gridrec, ART , SIRT, SART, BART, CGLS, PML, MLEM و OSEM اشاره کرد. به تازگی، جعبه ابزار ASTRA با TomoPy ادغام شده‌است. با ترکیب جعبه ابزار ASTRA با TomoPy، روش بازسازی مبتنی بر GPU به راحتی در دسترس کاربران قرار گرفت و کاربران جعبه ابزار ASTRA می‌توانند به راحتی داده‌ها را بخوانند و از قابلیت‌های دیگر TomoPy برای فیلتر کردن داده‌ها و تصحیح اعوجاج‌ها استفاده کنند.

نگارخانه
--------

در تصاویر آمده در این قسمت، روند کامل برای یک توموگرافی ساده از جسم و بازسازی توموگرافی بر اساس روش‌های تکراری آمده‌است. 

شکل سوم\_شی فنتوم، دو عدد مربع موربشکل پنجم\_تصویر بازسازی شده از سینوگرام با استفاده از روش‌های تکراری. شی اصلی تقریبا بازسازی شده‌است البته تصویر به دست آمده مقداری اعوجاج هم دارد.شکل چهارم\_سینوگرام یک شی فنتوم-حاصل شده از توموگرافی. پنجاه مقطع پروجکشن طی ۱۸۰ درجه زاویهگرفته شده‌است (نمونه‌گیری با فواصل یکسان)

جستارهای وابسته
---------------

* برش‌نگاری
* سی‌تی اسکن صنعتی

منابع
-----

1. ↑  Megherbi, N. , Breckon, T.P. , Flitton, G.T. , Mouton, A. (October 2013). "Radon Transform based Metal Artefacts Generation in 3D Threat Image Projection". *Proc. SPIE Optics and Photonics for Counterterrorism, Crime Fighting and Defence* (PDF). Vol. 8901. SPIE. pp. 1–7. doi:10.1117/12.2028506. Retrieved 5 November 2013.`{{cite book}}`: نگهداری یادکرد:نام‌های متعدد:فهرست نویسندگان (link)
2. 1 2 3 4 5  Dudgeon and Mersereau (1984). *Multidimensional digital signal processing*. Prentice-Hall.
3. ↑   هرمان GT، اصول توموگرافی کامپیوتری: بازسازی تصویر از طرح، نسخه دوم، Springer، 2009
4. 1 2  R. Mersereau, A. Oppenheim (1974). "Digital reconstruction of multidimensional signals from their projections". *Proceedings of the IEEE*. **62**: 1319–1338. doi:10.1109/proc.1974.9625.
5. ↑  Gursoy D, De Carlo F, Xiao X, and Jacobsen C (2014). "TomoPy: A framework for the analysis of synchrotron tomographic data". *Journal of Synchrotron Radiation*. **22**: 1188–1193. doi:10.1107/S1600577514013939. PMC 4181643.`{{cite journal}}`: نگهداری یادکرد:نام‌های متعدد:فهرست نویسندگان (link)
6. ↑  Van Aarle, W. , Palenstijn, W.J. , De Beenhouwer, J. , Altantzis T. , Bals S. , Batenburg K. J. , and J. Sijbers (October 2015). "The ASTRA Toolbox: a platform for advanced algorithm development in electron tomography". *Ultramicroscopy*. **157**: 35–47. doi:10.1016/j.ultramic.2015.05.002. PMID 26057688.`{{cite journal}}`: نگهداری یادکرد:نام‌های متعدد:فهرست نویسندگان (link)
7. ↑  W. Van Aarle, W J. Palenstijn, J. Cant, E. Janssens, F. Bleichrodt, A. Dabravolski, J. De Beenhouwer, K. J. Batenburg, and J. Sijbers (2016). "Fast and flexible X-ray tomography using the ASTRA toolbox". *Optics Express*. **24**: 35–47. Bibcode:2016OExpr..2425129V. doi:10.1364/OE.24.025129.`{{cite journal}}`: نگهداری یادکرد:نام‌های متعدد:فهرست نویسندگان (link)
8. ↑  Bicer T, Gursoy D, Kettimuthu R, De Carlo F, and Foster I (2016). "Optimization of tomographic reconstruction workflows on geographically distributed resources". *Journal of Synchrotron Radiation*. **23**: 997–1005. doi:10.1107/S1600577516007980. PMC 5315096.`{{cite journal}}`: نگهداری یادکرد:نام‌های متعدد:فهرست نویسندگان (link)
9. ↑  Pelt DM, Gursoy D, Batenburg KJ, De Carlo F, Palenstijna WJ, and Sijbers J (2016). "Integration of TomoPy and the ASTRA toolbox for advanced processing and reconstruction of tomographic synchrotron data". *Journal of Synchrotron Radiation*. **23**: 842–849. doi:10.1107/S1600577516005658. PMC 5315009.`{{cite journal}}`: نگهداری یادکرد:نام‌های متعدد:فهرست نویسندگان (link)

خواندن بیشتر
------------

* Avinash Kak & Malcolm Slaney (1988)، Principles of Computerized Imaging Tomography, IEEE Press،
* Bruyant, PP "الگوریتم‌های بازسازی تحلیلی و تکراری در SPECT" مجله پزشکی هسته ای ۴3 (10): ۱۳۴۳–۱۳۵۸، ۲۰۰۲

پیوند به بیرون
--------------

* Slaney, A. C. Kak and Malcolm. "Principles of Computerized Tomographic Imaging". *Slaney.org*. Retrieved 7 September 2018.
* ابزار Insight ToolKit؛ نرم‌افزار پشتیبانی توموگرافی منبع باز
* "TomoPy — TomoPy 1.1.3 documentation". *Tomopy.readthedocs.org*. Retrieved 7 September 2018.
* ASTRA (همه مقیاس بازسازی تونوگرافی آنتورپ) جعبه ابزار؛ بسیار انعطاف‌پذیر، سریع و نرم‌افزار منبع باز برای بازسازی توموگرافی کامپیوتری
* NiftyRec؛ نرم‌افزار بازسازی توموگرافی باز جامع؛ Matlab و پایتون قابل برنامه‌نویسی
* ابزار توزیع باز و ابزار تجسم باز
* "ITS plc - Electrical Process Tomography For Industrial Visualization". *Itoms.com*. Retrieved 7 September 2018.



| * ن
* ب
* و
تصویربرداری پزشکی |
| --- |
| پرتو ایکس/پرتونگاری | 

|  |  |  |  |  |  |
| --- | --- | --- | --- | --- | --- |
| 2D | 

|  |  |
| --- | --- |
| پزشکی: | 
* پنوموآنسفالوگرافی
* Dental radiography
* سیالوگرافی
* Myelography
* CXR
	+ Bronchography
* AXR
* KUB
* DXA/DXR
* Upper gastrointestinal series/Small-bowel follow-through/Lower gastrointestinal series
* کولانژیوگرافی/کوله‌سیستوگرافی
* ماموگرافی
* Pyelogram
* سیستوگرافی
* Arthrogram
* رحم‌لوله‌نگاری
* Skeletal survey
* آنژیوگرافی
	+ Angiocardiography
	+ Aortography
* Venography
* Lymphogram
* Orbital x-ray

 |
| صنعتی: | 
* آزمون رادیوگرافی

 |

 |
| CT اسکن | 

|  |  |
| --- | --- |
| فنون: | 
* General operation
* Quantitative
* High-resolution
* میکرو سی‌تی
* پرتو الکترونی
* پرتو مخروطی

 |
| اهداف | 
* قلب
	+ calcium scan
	+ angiography
* Abdominal and pelvis
	+ Virtual colonoscopy
* آنژیوگرافی
	+ Coronary
	+ Pulmonary
* Head
* Thyroid
* Whole body imaging
	+ Full-body CT scan

 |

 |
| سایر | 
* فلوئوروسکوپی
* پرتونگاری پانورامیک دندانی
* X-ray motion analysis

 |

 |
| MRI | 
* Brain
	+ عملکردی
* Neurography
* قلبی
	+ perfusion
* Angiography
* Cholangiopancreatography (MRCP)
* Breast
* Sequences
	+ پخش
	+ restriction
	+ Tractography
* Synthetic MRI

 |
| فراصوت | 
* فنون
	+ دوپلر
	+ contrast-enhanced
	+ 3D
	+ endoscopic
	+ دوپلکس
* اکوکاردیوگرافی
	+ Doppler echocardiography
	+ TTE
	+ TEE
* داپلر ترانس کرانیال
* Intravascular
* زنان و زایمان
* Obstetric
* Echoencephalography
* Abdominal ultrasonography
	+ renal
	+ renal tract
* Rectal
* Breast
* کیسه بیضه
* Carotid
* Emergency ultrasound
	+ FAST
	+ pre-hospital

 |
| رادیو-هسته‌ای | 

|  |  |  |  |  |
| --- | --- | --- | --- | --- |
| 2D / جرقه‌نگاری | 

|  |
| --- |
| 
* Cholescintigraphy
* Scintimammography
* اسکن‌های ریوی در پزشکی هسته‌ای
* Radionuclide ventriculography
* اسکن موگا
* Radioisotope renography
* Sestamibi parathyroid scintigraphy
* اسکن تیروئید
* جرقه‌نگاری استخوان
* Immunoscintigraphy
* Dacryoscintigraphy
* DMSA scan
* Gastric emptying scan

 |
| کل بدن: | 
* Octreotide scan
* اسکن گالیوم 67
* Ga-68-DOTATOC
* Indium-111 WBC scan

 |

 |
| 3D / ECT | 

|  |  |
| --- | --- |
| SPECT (پرتو گاما): | 
* Myocardial perfusion imaging

 |
| PET (پوزیترون): | 
* برش نگاری مغز با گسیل پوزیترون(Brain PET)
* Cardiac PET
* PET mammography
* PET-CT
* PET-MRI

 |

 |

 |
| نوری/لیزر | 
* Optical tomography
	+ مقطع‌نگاری همدوسی اپتیکی
* میکروسکوپ هم-کانونی
* Endomicroscopy
* Orthogonal polarization spectral imaging

 |
| دمانگاری | 
* non-contact thermography
* contact thermography
* dynamic angiothermography

 |
| شرایط هدف | 
* سکته حاد
* Pregnancy

 |
| 
* Category

 |