आंतरिक पुनर्निर्माण

From Vigyanwiki

डिजिटल इमेजिंग में पुनरावृत्त पुनर्निर्माण में, आंतरिक पुनर्निर्माण (जिसे दृश्य के सीमित क्षेत्र (एलएफवी) पुनर्निर्माण के रूप में भी जाना जाता है) छवि डेटा को दृश्य के एक छोटे क्षेत्र तक सीमित करने के कारण होने वाली ट्रंकेशन कलाकृतियों को सही करने की एक तकनीक है। पुनर्निर्माण एक ऐसे क्षेत्र पर केंद्रित है जिसे रुचि का क्षेत्र (आरओआई) कहा जाता है। यद्यपि आंतरिक पुनर्निर्माण को दंत या कार्डियक एक्स-रे कंप्यूटेड टोमोग्राफी छवियों पर प्रयुक्त किया जा सकता है, यह अवधारणा सीटी तक सीमित नहीं है। इसे कई विधियों में से एक के साथ प्रयुक्त किया जाता है।

विधियों

प्रत्येक विधि का उद्देश्य सदिश का समाधान करना है निम्नलिखित समस्या में:

Two diagrams
किसी वस्तु को दिखाने वाली छवि का रुचि का क्षेत्र (आरओआई)।

मान लीजिए कि रुचि का क्षेत्र (आरओआई) है और , के बाहर का क्षेत्र है। मान लीजिए कि , , , ज्ञात आव्यूह हैं; और मूल छवि के अज्ञात सदिश हैं, जबकि और प्रतिक्रियाओं के सदिश माप हैं ( ज्ञात है और अज्ञात है)। क्षेत्र के अंदर है और , क्षेत्र में,, क्षेत्र के बाहर है। , के अनुरूप माप में एक क्षेत्र के अंदर है। इस क्षेत्र को , ( के रूप में दर्शाया गया है, जबकि g है क्षेत्र F के बाहर। यह से मेल खाता है और इसे , (). के रूप में दर्शाया गया है

सीटी छवि-पुनर्निर्माण उद्देश्यों के लिए, .

आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, आव्यूह , , , सम्मिश्र ऑपरेटर (गणित) के अतिरिक्त छवि पुनर्निर्माण के लिए प्रयुक्त किया जाता है।

नीचे सूचीबद्ध पहली आंतरिक-पुनर्निर्माण विधि एक्सट्रपलेशन है। यह एक स्थानीय टोमोग्राफी विधि है जो ट्रंकेशन कलाकृतियों को समाप्त करती है किंतु एक अन्य प्रकार की कलाकृतियों का परिचय देती है: एक कटोरा प्रभाव एक सुधार को अनुकूली एक्सट्रपलेशन विधि के रूप में जाना जाता है, चूँकि नीचे दी गई पुनरावृत्त एक्सट्रपलेशन विधि भी पुनर्निर्माण परिणामों में सुधार करती है। कुछ स्थितियों में, आंतरिक पुनर्निर्माण के लिए स्पष्ट पुनर्निर्माण पाया जा सकता है। नीचे दी गई स्थानीय विपरीत विधि स्थानीय टोमोग्राफी विधि को संशोधित करती है, और स्थानीय टोमोग्राफी के पुनर्निर्माण परिणाम में सुधार कर सकती है; पुनरावृत्तीय पुनर्निर्माण पद्धति को आंतरिक पुनर्निर्माण पर प्रयुक्त किया जा सकता है। उपरोक्त विधियों में, एक्सट्रपलेशन को अधिकांशतः प्रयुक्त किया जाता है।

एक्सट्रपलेशन विधि

Six views of an image
1) शीप-लोगन फैंटम के प्रक्षेपण 2) काटे गए प्रक्षेपण (शून्य एक्सट्रपलेशन) 3) स्थिर, 4) घातीय और 5) द्विघात एक्सट्रपलेशन 6) 4 और 5 का मिश्रित एक्सट्रपलेशन

:

, , , ज्ञात आव्यूह हैं; और अज्ञात सदिश हैं; एक ज्ञात सदिश है, और एक अज्ञात सदिश है। हमें सदिश को जानना होगा। और मूल छवि हैं, जबकि और प्रतिक्रियाओं के माप हैं। सदिश रुचि के क्षेत्र , (). के अंदर है। सदिश , क्षेत्र के बाहर है। बाहरी क्षेत्र को , () कहा जाता है और , के अनुरूप माप में एक क्षेत्र के अंदर है। इस क्षेत्र को , () दर्शाया गया है। सदिश का क्षेत्र (क्षेत्र के बाहर) भी Y से मेल खाता है और इसे , (). के रूप में दर्शाया गया है। सीटी छवि पुनर्निर्माण में यह है

आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, आव्यूह , , , एक सम्मिश्र ऑपरेटर के अतिरिक्त छवि पुनर्निर्माण के लिए प्रयुक्त किया जाता है।

बाहरी क्षेत्र में प्रतिक्रिया का अनुमान हो सकता है; उदाहरण के लिए, मान लें कि यह है।

Eight views of an image
ए) शेप-लोगान हेड फैंटम बी) फैंटम की फसल सी) एक्सट्रपलेशन के बिना पुनर्निर्माण डी) स्थिरांक के साथ पुनर्निर्माण, (ई) द्विघात और (एफ) मिश्रित एक्सट्रपलेशन

का एक समाधान के रूप में लिखा जाता है, और इसे एक्सट्रपलेशन विधि के रूप में जाना जाता है। परिणाम इस बात पर निर्भर करता है कि एक्सट्रपलेशन फलन कितना अच्छा है। एक अधिकांशतः पसंद है

दो क्षेत्रों की सीमा पर.[1][2][3][4] एक्सट्रपलेशन विधि को अधिकांशतः प्राथमिक ज्ञान के साथ जोड़ दिया जाता है,[5][6] और एक एक्सट्रपलेशन विधि जो गणना समय को कम करती है, नीचे दिखाई गई है।

अनुकूली एक्सट्रपलेशन विधि

मान लें कि एक समान्य समाधान, और ऊपर वर्णित एक्सट्रपलेशन विधि से प्राप्त किया गया है। बाहरी क्षेत्र में प्रतिक्रिया की गणना निम्नानुसार की जा सकती है:

पुनर्निर्मित छवि की गणना निम्नानुसार की जा सकती है:

यह मान लिया है कि

आंतरिक क्षेत्र की सीमा पर; समस्या का समाधान करता है, और इसे अनुकूली एक्सट्रपलेशन विधि के रूप में जाना जाता है। जिसमे अनुकूली एक्सट्रपलेशन फलन है।[7][8][9][10][5]

पुनरावृत्त एक्सट्रपलेशन विधि

यह माना जाता है कि एक समान्य समाधान, और , नीचे वर्णित एक्सट्रपलेशन विधि से प्राप्त किया गया है:

या

पुनर्निर्माण के रूप में प्राप्त किया जा सकता है

यहाँ एक एक्सट्रपलेशन फलन है, और यह माना जाता है

इस समस्या का एक समाधान है.[11]

स्थानीय टोमोग्राफी

बहुत छोटे फिल्टर वाली स्थानीय टोमोग्राफी को लैम्ब्डा टोमोग्राफी के रूप में भी जाना जाता है।[12][13]

स्थानीय व्युत्क्रम विधि

स्थानीय व्युत्क्रम विधि स्थानीय टोमोग्राफी की अवधारणा का विस्तार करती है। बाहरी क्षेत्र में प्रतिक्रिया की गणना इस प्रकार की जा सकती है:

सामान्यीकृत व्युत्क्रम को संतोषजनक मानें

परिभाषित करना

जिससे

इस तरह,

उपरोक्त समीकरण को इस प्रकार हल किया जा सकता है

,

ध्यान में रख कर

, का सामान्यीकृत व्युत्क्रम है , अर्थात।

समाधान को इस प्रकार सरल बनाया जा सकता है

.

गणित का प्रश्न

आव्यूह को आव्यूह के स्थानीय व्युत्क्रम के रूप में जाना जाता है। इस प्रकार के अनुरूप इसे स्थानीय व्युत्क्रम विधि के रूप में जाना जाता है। [11]


पुनरावृत्तीय पुनर्निर्माण विधि

यहां एक लक्ष्य फलन परिभाषित किया गया है, और यह विधि लक्ष्य को पुनरावृत्त रूप से प्राप्त करती है। यदि लक्ष्य फलन किसी प्रकार का सामान्य हो सकता है, तो इसे न्यूनतम मानक विधि के रूप में जाना जाता है।

,

का विषय है

और ज्ञात है,

जहां , और न्यूनतमकरण और के भार स्थिरांक हैं एक प्रकार का आदर्श है. अधिकांशतः उपयोग किए जाने वाले मानदंड , , , कुल भिन्नता (टीवी) मानदंड या उपरोक्त मानदंडों का संयोजन हैं। इस विधि का एक उदाहरण उत्तल सेट (पीओसीएस) विधि पर प्रक्षेपण है।[14][15]

विश्लेषणात्मक समाधान

विशेष परिस्थितियों में, आंतरिक पुनर्निर्माण को एक विश्लेषणात्मक समाधान के रूप में प्राप्त किया जा सकता है; जिसमे का समाधान है जो ऐसे स्थितियों में स्पष्ट है.[16][17][18]


तेज़ एक्सट्रपलेशन

एक्सट्रपोलेटेड डेटा अधिकांशतः धनात्मक-निश्चित कर्नेल में कनवल्शन होता है। डेटा को एक्सट्रपलेशन के बाद इसका आकार N गुना बढ़ जाता है, जहां N = 2 ~ 3. यदि डेटा को ज्ञात कर्नेल फलन में संयोजित करने की आवश्यकता है, तो संख्यात्मक गणना log(NN गुना बढ़ जाएगी, यहां तक ​​कि तेज़ फास्ट फूरियर ट्रांसफॉर्म (FFT) के साथ भी एक एल्गोरिदम उपयोग है, जो एक्सट्रपलेटेड डेटा के भाग से योगदान की विश्लेषणात्मक गणना करता है। मूल कनवल्शन गणना की तुलना में गणना समय छोड़ा जा सकता है; इस एल्गोरिथ्म के साथ, एक्सट्रपोलेटेड डेटा का उपयोग करके कनवल्शन की गणना में उल्लेखनीय वृद्धि नहीं होती है। इसे तीव्र एक्सट्रपलेशन के रूप में जाना जाता है।[19]

विधियों की तुलना

एक्सट्रपलेशन विधि ऐसी स्थिति में उपयुक्त होती है

और
अथार्त एक छोटी ट्रिमिंग कलाकृतियों की स्थिति है ।

अनुकूली एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है

और
अथार्त एक सामान्य ट्रंकेशन कलाकृतियों की स्थिति यह विधि बाहरी क्षेत्र के लिए एक समान्य समाधान भी प्रस्तुत करती है।

पुनरावृत्त एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है

और
अथार्त एक सामान्य ट्रंकेशन कलाकृतियों की स्थिति यद्यपि यह विधि अनुकूली पुनर्निर्माण की तुलना में उत्तम आंतरिक पुनर्निर्माण प्राप्त करती है, किंतु बाहरी क्षेत्र में इसका परिणाम नहीं मिलता है।

ऐसी स्थिति के लिए स्थानीय टोमोग्राफी उपयुक्त है

और
अर्थात सबसे बड़ी ट्रंकेशन कलाकृतियों की स्थिति। चूँकि इस पद्धति में कोई ट्रिमिंग कलाकृतियाँ नहीं हैं, पुनर्निर्माण में एक निश्चित त्रुटि ( के मान से स्वतंत्र) है।

स्थानीय व्युत्क्रम विधि स्थानीय टोमोग्राफी के समान ऐसी स्थिति में उपयुक्त है

और
अथार्त सबसे बड़ी ट्रंकेशन कलाकृतियों की स्थिति। चूँकि इस विधि के लिए कोई खंडन कलाकृतियाँ नहीं हैं, पुनर्निर्माण में एक निश्चित त्रुटि ( के मान से स्वतंत्र) है जो स्थानीय टोमोग्राफी से छोटी हो सकती है।

पुनरावृत्तीय पुनर्निर्माण विधि बड़ी गणनाओं के साथ अच्छा परिणाम प्राप्त करती है। यद्यपि विश्लेषणात्मक विधि स्पष्ट परिणाम प्राप्त करती है, यह केवल कुछ स्थितियों में ही कार्यात्मक होती है। तेज़ एक्सट्रपलेशन विधि अन्य एक्सट्रपलेशन विधियों के समान परिणाम प्राप्त कर सकती है, और गणना को कम करने के लिए उपरोक्त आंतरिक पुनर्निर्माण विधियों पर प्रयुक्त किया जा सकता है।

यह भी देखें

टिप्पणियाँ

  1. M.M. Seger, Rampfilter implementation on truncated projection data. Application to 3D linear tomography for logs, Proceedings SSAB02, Symposium on Image Analysis, Lund, Sweden, 7–8 March 2002. Editor Astrom.
  2. F .Rashid-Farrokhi, K.J.R. Liu, C.A. Berenstein and D.Walnut, Wavelet-based Multiresolution Local Tomography, IEEE Transactions on Image Processing 6 (1997), 1412–1430.
  3. M. Nilsson, Local Tomography at a Glance, Licentiate Theses in Mathematical Sciences 2003:3 ISSN 1404-028X, ISBN 91-628-5741-X, LUTFMA-2007-2003. Printed in Sweden by KFS AB Lund, 2003.
  4. P.S. Cho, A.D. Rudd and R.H. Johnson, Cone-beam CT from width truncated projections, Computerized Medical Imaging and Graphics 20(1) (1996), 49–57, 49–57.
  5. 5.0 5.1 J. Hsieh, E. Chao, J. Thibault, B. Grekowicz, A. Horst, S. McOlash and T.J. Myers, A novel reconstruction algorithm to extend the CT scan fieldofview, Medical Phys 31 (2004), 2385–2391.
  6. K.J. Ruchala, G.H. Olivera, J.M. Kapatoes, P.J. Reckwerdt and T.R. Mack, Methods for improving limited fieldofview radiotherapy reconstructions using imperfect a priori images, Med Phys 29 (2002), 2590–2605.
  7. M. Nassi, W.R. Brody, B.P.Medoff and A.Macovski, Iterative reconstruction reprojection: an algorithm for limited data cardiac computed tomography, IEEE trans Biomed Engineering 295 (1982), 333–340.
  8. J.H. Kim, K.Y. KWAK, S.B. Park and Z.H. Cho, Projection space iteration reconstruction reprojection, IEEE transaction on Medical Imaging 4 (1983), 139–143
  9. P.S.Cho, A.D. Rudd and R.H. Johnson, Conebeam CT from width truncated projections Computerized, Medical Imaging and Graphics 20 (1996), 49–57.
  10. B. Ohnesorge, T. Flohr, K. Schwarz, J.P. Heiken and K.T. Bae, 2000 Efficient correction for CT image artifacts caused by objects extending outside the scan field of view, Med Phys 27, 39–46.
  11. 11.0 11.1 Shuangren Zhao, Kang Yang, Dazong Jiang, Xintie Yang, Interior reconstruction using local inverse, J Xray Sci Technol. 2011; 19(1): 69-90
  12. A. Faridani, E.L. Ritman and K.T. Smith, Local tomography, SIAM J APPL MATH 52 (1992), 459–484.
  13. A. Katsevich, 1999 Cone beam local tomography, SIAM J APPL MATH 59, 2224–2246.
  14. Ye. Yangbo, Yu. 1 Hengyong 2 and GeWang, Exact Interior Reconstruction from Truncated Limited Angle Projection Data, International Journal of Biomedical Imaging (2008), 1–6.
  15. L. Zeng, B. Liu, L. Liu and C. Xiang, A new iterative reconstruction algorithm for 2D exterior fanbeam CT, Journal of XRay Science and Technology 18 (2010), 267–277.
  16. Y. Zou and X. Pan, 2004, Exact image reconstruction on PIlines from minimum data in helical conebeam CT, Physics in Medicine and Biology 49(6), 941–959.
  17. M. Defrise, F. Noo, R. Clackdoyle and H. Kudo, Truncated Hilbert transform and image reconstruction from limited tomographic data. IOPscience.iop.org, 2006
  18. F. Noo, R. Clackdoyle and J.D. Pack, A twostep Hilbert transform method for 2D image reconstruction, Phys Med Biol 49 (2004), 3903–3923.
  19. S Zhao, K Yang, X Yang, Reconstruction from truncated projections using mixed extrapolations of exponential and quadratic functions, Journal of X-ray Science and Technology, 2011, 19(2) pp 155–72