आंतरिक पुनर्निर्माण: Difference between revisions
No edit summary |
|||
| (8 intermediate revisions by 4 users not shown) | |||
| Line 1: | Line 1: | ||
[[डिजिटल इमेजिंग]] में पुनरावृत्त पुनर्निर्माण में, आंतरिक पुनर्निर्माण (जिसे दृश्य के सीमित क्षेत्र (एलएफवी) पुनर्निर्माण के रूप में भी जाना जाता है) छवि डेटा को दृश्य के एक छोटे क्षेत्र तक सीमित करने के कारण होने वाली ट्रंकेशन कलाकृतियों को सही करने की एक तकनीक है। पुनर्निर्माण एक ऐसे क्षेत्र पर केंद्रित है जिसे रुचि का क्षेत्र (आरओआई) कहा जाता है। यद्यपि आंतरिक पुनर्निर्माण को दंत या कार्डियक [[एक्स-रे कंप्यूटेड टोमोग्राफी]] छवियों पर | [[डिजिटल इमेजिंग]] में पुनरावृत्त पुनर्निर्माण में, '''आंतरिक पुनर्निर्माण''' (जिसे दृश्य के सीमित क्षेत्र (एलएफवी) पुनर्निर्माण के रूप में भी जाना जाता है) छवि डेटा को दृश्य के एक छोटे क्षेत्र तक सीमित करने के कारण होने वाली ट्रंकेशन कलाकृतियों को सही करने की एक तकनीक है। पुनर्निर्माण एक ऐसे क्षेत्र पर केंद्रित है जिसे रुचि का क्षेत्र (आरओआई) कहा जाता है। यद्यपि आंतरिक पुनर्निर्माण को दंत या कार्डियक [[एक्स-रे कंप्यूटेड टोमोग्राफी]] छवियों पर प्रयुक्त किया जा सकता है, यह अवधारणा सीटी तक सीमित नहीं है। इसे कई विधियों में से एक के साथ प्रयुक्त किया जाता है। | ||
== | ==विधियों == | ||
प्रत्येक विधि का उद्देश्य | प्रत्येक विधि का उद्देश्य सदिश <math>x | ||
</math> का समाधान करना है निम्नलिखित समस्या में: | |||
:<math> | :<math> | ||
| Line 18: | Line 20: | ||
</math> | </math> | ||
[[File:ROI and object.png|thumb|alt=Two diagrams|किसी वस्तु को दिखाने वाली छवि का रुचि का क्षेत्र (आरओआई)।]] | [[File:ROI and object.png|thumb|alt=Two diagrams|किसी वस्तु को दिखाने वाली छवि का रुचि का क्षेत्र (आरओआई)।]]मान लीजिए कि <math>X</math> रुचि का क्षेत्र (आरओआई) है और <math>Y</math>, <math>X</math> के बाहर का क्षेत्र है। मान लीजिए कि <math>A</math>, <math>B</math>, <math>C</math>, <math>D</math> ज्ञात आव्यूह हैं; <math>x</math> और <math>y</math> मूल छवि के अज्ञात सदिश हैं, जबकि <math>f</math> और <math>g</math> प्रतिक्रियाओं के सदिश माप हैं (<math>f</math> ज्ञात है और <math>g</math> अज्ञात है)। <math>X</math> क्षेत्र <math>x \in X</math> के अंदर है और <math>Y</math>, क्षेत्र <math>Y</math> में,<math>y \in Y</math>, क्षेत्र <math>X</math> के बाहर है। <math>F</math>, <math>X</math> के अनुरूप माप में एक क्षेत्र के अंदर है। इस क्षेत्र को <math>F</math>, (<math>f \in F</math> के रूप में दर्शाया गया है, जबकि g है क्षेत्र F के बाहर। यह <math>Y</math> से मेल खाता है और इसे <math>G</math>, (<math>g \in G</math>). के रूप में दर्शाया गया है | ||
मान लीजिए <math>A</math>, <math>B</math>, <math>C</math>, <math>D</math> ज्ञात आव्यूह हैं; <math>x</math> और <math>y</math> | |||
सीटी छवि-पुनर्निर्माण उद्देश्यों के लिए, <math> C = 0 </math>. | सीटी छवि-पुनर्निर्माण उद्देश्यों के लिए, <math> C = 0 </math>. | ||
आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, | आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, आव्यूह <math>A</math>, <math>B</math>, <math>C</math>, <math>D</math> सम्मिश्र [[ऑपरेटर (गणित)]] के अतिरिक्त छवि पुनर्निर्माण के लिए प्रयुक्त किया जाता है। | ||
नीचे सूचीबद्ध पहली आंतरिक-पुनर्निर्माण विधि [[एक्सट्रपलेशन]] है। यह एक स्थानीय टोमोग्राफी विधि है जो ट्रंकेशन कलाकृतियों को समाप्त करती है किंतु एक अन्य प्रकार की कलाकृतियों का परिचय देती है: एक कटोरा प्रभाव एक सुधार को अनुकूली एक्सट्रपलेशन विधि के रूप में जाना जाता है, चूँकि नीचे दी गई पुनरावृत्त एक्सट्रपलेशन विधि भी पुनर्निर्माण परिणामों में सुधार करती है। कुछ स्थितियों में, आंतरिक पुनर्निर्माण के लिए स्पष्ट पुनर्निर्माण पाया जा सकता है। नीचे दी गई स्थानीय विपरीत विधि स्थानीय टोमोग्राफी विधि को संशोधित करती है, और स्थानीय टोमोग्राफी के पुनर्निर्माण परिणाम में सुधार कर सकती है; पुनरावृत्तीय पुनर्निर्माण पद्धति को आंतरिक पुनर्निर्माण पर प्रयुक्त किया जा सकता है। उपरोक्त विधियों में, एक्सट्रपलेशन को अधिकांशतः प्रयुक्त किया जाता है। | |||
=== एक्सट्रपलेशन विधि === | === एक्सट्रपलेशन विधि === | ||
[[File:Different extrapolations of the Sheep-Logan phantoms.jpg|thumb|alt=Six views of an image|1) शीप-लोगन फैंटम के प्रक्षेपण 2) काटे गए प्रक्षेपण (शून्य एक्सट्रपलेशन) 3) स्थिर, 4) घातीय और 5) द्विघात एक्सट्रपलेशन 6) 4 और 5 का मिश्रित एक्सट्रपलेशन]]:<math> | [[File:Different extrapolations of the Sheep-Logan phantoms.jpg|thumb|alt=Six views of an image|1) शीप-लोगन फैंटम के प्रक्षेपण 2) काटे गए प्रक्षेपण (शून्य एक्सट्रपलेशन) 3) स्थिर, 4) घातीय और 5) द्विघात एक्सट्रपलेशन 6) 4 और 5 का मिश्रित एक्सट्रपलेशन]]:<math> | ||
| Line 44: | Line 42: | ||
</math> | </math> | ||
<math>A</math>, <math>B</math>, <math>C</math>, <math>D</math> ज्ञात आव्यूह हैं; <math>x</math> और <math>y</math> अज्ञात सदिश हैं; <math>f</math> एक ज्ञात | <math>A</math>, <math>B</math>, <math>C</math>, <math>D</math> ज्ञात आव्यूह हैं; <math>x</math> और <math>y</math> अज्ञात सदिश हैं; <math>f</math> एक ज्ञात सदिश है, और <math>g</math> एक अज्ञात सदिश है। हमें सदिश <math>x</math> को जानना होगा। <math>x</math> और <math>y</math> मूल छवि हैं, जबकि <math>f</math> और <math>g</math> प्रतिक्रियाओं के माप हैं। सदिश <math>x</math> रुचि के क्षेत्र <math>X</math>, (<math>x \in X</math>). के अंदर है। सदिश <math>y</math>, क्षेत्र <math>X</math> के बाहर है। बाहरी क्षेत्र को <math>Y</math>, (<math>y \in Y</math>) कहा जाता है और <math>f</math>, <math>X</math> के अनुरूप माप में एक क्षेत्र के अंदर है। इस क्षेत्र को <math>F</math>, (<math>f \in F</math>) दर्शाया गया है। सदिश <math>g</math> का क्षेत्र (क्षेत्र <math>F</math> के बाहर) भी Y से मेल खाता है और इसे <math>G</math>, (<math>g \in G</math>). के रूप में दर्शाया गया है। सीटी छवि पुनर्निर्माण में यह है | ||
सीटी छवि पुनर्निर्माण में | |||
:<math> C = 0 </math> | :<math> C = 0 </math> | ||
आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, | आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, आव्यूह <math>A</math>, <math>B</math>, <math>C</math>, <math>D</math> एक सम्मिश्र ऑपरेटर के अतिरिक्त छवि पुनर्निर्माण के लिए प्रयुक्त किया जाता है। | ||
बाहरी क्षेत्र में प्रतिक्रिया का | बाहरी क्षेत्र में प्रतिक्रिया का अनुमान <math>G</math> हो सकता है; उदाहरण के लिए, मान लें कि यह <math>g_{ex}</math> है। | ||
: <math> | : <math> | ||
\begin{bmatrix} | \begin{bmatrix} | ||
| Line 64: | Line 62: | ||
</math> | </math> | ||
[[File:ExtrapolatedReconstruction.png|thumb|alt=Eight views of an image|ए) शेप-लोगान हेड फैंटम बी) फैंटम की फसल सी) एक्सट्रपलेशन के बिना पुनर्निर्माण डी) स्थिरांक के साथ पुनर्निर्माण, (ई) द्विघात और (एफ) मिश्रित एक्सट्रपलेशन]] | [[File:ExtrapolatedReconstruction.png|thumb|alt=Eight views of an image|ए) शेप-लोगान हेड फैंटम बी) फैंटम की फसल सी) एक्सट्रपलेशन के बिना पुनर्निर्माण डी) स्थिरांक के साथ पुनर्निर्माण, (ई) द्विघात और (एफ) मिश्रित एक्सट्रपलेशन]] | ||
<math>x</math> का एक समाधान <math>x_0</math> के रूप में लिखा जाता है, और इसे एक्सट्रपलेशन विधि के रूप में जाना जाता है। परिणाम इस बात पर निर्भर करता है कि एक्सट्रपलेशन फलन <math>g_{ex}</math> कितना अच्छा है। एक अधिकांशतः पसंद है | |||
: <math>g_{ex}|_{\partial G}=f |_{\partial F}</math> | : <math>g_{ex}|_{\partial G}=f |_{\partial F}</math> | ||
दो क्षेत्रों की सीमा पर.<ref>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.</ref><ref>F .Rashid-Farrokhi, K.J.R. Liu, C.A. Berenstein and D.Walnut, Wavelet-based Multiresolution Local Tomography, | दो क्षेत्रों की सीमा पर.<ref>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.</ref><ref>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.</ref><ref>M. Nilsson, Local Tomography at a Glance, Licentiate Theses in Mathematical Sciences 2003:3 {{ISSN|1404-028X}}, | IEEE Transactions on Image Processing 6 (1997), 1412–1430.</ref><ref>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.</ref><ref>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.</ref> | {{ISBN|91-628-5741-X}}, LUTFMA-2007-2003. Printed in Sweden by KFS AB Lund, 2003.</ref><ref>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.</ref> एक्सट्रपलेशन विधि को अधिकांशतः प्राथमिक ज्ञान के साथ जोड़ दिया जाता है,<ref name="ReferenceA">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, | extend the CT scan fieldofview, | ||
Medical Phys 31 (2004), 2385–2391.</ref><ref>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.</ref> और एक एक्सट्रपलेशन विधि जो गणना समय को कम करती है, नीचे दिखाई गई है। | Medical Phys 31 (2004), 2385–2391.</ref><ref>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.</ref> और एक एक्सट्रपलेशन विधि जो गणना समय को कम करती है, नीचे दिखाई गई है। | ||
| Line 76: | Line 75: | ||
=== अनुकूली एक्सट्रपलेशन विधि === | === अनुकूली एक्सट्रपलेशन विधि === | ||
एक | मान लें कि एक समान्य समाधान, <math>x_0</math> और <math>y_0</math> ऊपर वर्णित एक्सट्रपलेशन विधि से प्राप्त किया गया है। बाहरी क्षेत्र <math>g_1</math> में प्रतिक्रिया की गणना निम्नानुसार की जा सकती है: | ||
:<math> | :<math> | ||
| Line 97: | Line 96: | ||
</math> | </math> | ||
यह मान लिया है कि | यह मान लिया है कि | ||
:<math>f |_{\partial F}=(g_1+g_{1ex})|_{\partial G}</math> आंतरिक क्षेत्र की सीमा पर; <math>x_1</math> समस्या का समाधान करता है, और इसे अनुकूली एक्सट्रपलेशन विधि के रूप में जाना जाता है। | :<math>f |_{\partial F}=(g_1+g_{1ex})|_{\partial G}</math> | ||
<math>g_{1ex}</math> अनुकूली एक्सट्रपलेशन | |||
आंतरिक क्षेत्र की सीमा पर; <math>x_1</math> समस्या का समाधान करता है, और इसे अनुकूली एक्सट्रपलेशन विधि के रूप में जाना जाता है। जिसमे <math>g_{1ex}</math>अनुकूली एक्सट्रपलेशन फलन है।<ref>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.</ref><ref>J.H. Kim, K.Y. KWAK, S.B. Park and Z.H. Cho, Projection space iteration reconstruction reprojection, | tomography, IEEE trans Biomed Engineering 295 (1982), 333–340.</ref><ref>J.H. Kim, K.Y. KWAK, S.B. Park and Z.H. Cho, Projection space iteration reconstruction reprojection, | ||
IEEE transaction | IEEE transaction | ||
| Line 104: | Line 104: | ||
CT from width truncated projections Computerized, Medical Imaging | CT from width truncated projections Computerized, Medical Imaging | ||
and Graphics 20 (1996), 49–57.</ref><ref>B. Ohnesorge, T. Flohr, K. Schwarz, J.P. Heiken and K.T. Bae, 2000 Efficient correction for CT image artifacts caused | and Graphics 20 (1996), 49–57.</ref><ref>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.</ref><ref name="ReferenceA"/> | by objects extending outside the scan field of view, Med Phys 27, 39–46.</ref><ref name="ReferenceA" /> | ||
=== पुनरावृत्त एक्सट्रपलेशन विधि === | === पुनरावृत्त एक्सट्रपलेशन विधि === | ||
यह माना जाता है कि एक | यह माना जाता है कि एक समान्य समाधान, <math>x_0</math> और <math>y_0</math>, नीचे वर्णित एक्सट्रपलेशन विधि से प्राप्त किया गया है: | ||
:<math> | :<math> | ||
| Line 143: | Line 141: | ||
\end{bmatrix} | \end{bmatrix} | ||
</math> | </math> | ||
यहाँ <math>g_{1ex}</math> एक एक्सट्रपलेशन | यहाँ <math>g_{1ex}</math> एक एक्सट्रपलेशन फलन है, और यह माना जाता है | ||
:<math>(f-f_1)|_{\partial F}=g_{1ex}|_{\partial G}</math> | :<math>(f-f_1)|_{\partial F}=g_{1ex}|_{\partial G}</math> | ||
<math>x_1</math> | इस समस्या का <math>x_1</math> एक समाधान है.<ref name="shuangrenZhao 2011">Shuangren Zhao, Kang Yang, Dazong Jiang, Xintie Yang, ''Interior reconstruction using local inverse'', ''J Xray Sci Technol''. 2011; 19(1): 69-90 | ||
</ref> | </ref> | ||
=== स्थानीय टोमोग्राफी === | === स्थानीय टोमोग्राफी === | ||
बहुत छोटे फिल्टर वाली स्थानीय टोमोग्राफी को लैम्ब्डा टोमोग्राफी के रूप में भी जाना जाता है।<ref>A. Faridani, E.L. Ritman and K.T. Smith, Local tomography, SIAM J APPL MATH 52 (1992), 459–484.</ref><ref>A. Katsevich, 1999 Cone beam local tomography, SIAM J APPL MATH 59, 2224–2246.</ref> | बहुत छोटे फिल्टर वाली स्थानीय टोमोग्राफी को लैम्ब्डा टोमोग्राफी के रूप में भी जाना जाता है।<ref>A. Faridani, E.L. Ritman and K.T. Smith, Local tomography, SIAM J APPL MATH 52 (1992), 459–484.</ref><ref>A. Katsevich, 1999 Cone beam local tomography, SIAM J APPL MATH 59, 2224–2246.</ref> | ||
=== स्थानीय व्युत्क्रम विधि === | === स्थानीय व्युत्क्रम विधि === | ||
स्थानीय व्युत्क्रम विधि स्थानीय टोमोग्राफी की अवधारणा का विस्तार करती है। बाहरी क्षेत्र में प्रतिक्रिया की गणना इस प्रकार की जा सकती है: | स्थानीय व्युत्क्रम विधि स्थानीय टोमोग्राफी की अवधारणा का विस्तार करती है। बाहरी क्षेत्र में प्रतिक्रिया की गणना इस प्रकार की जा सकती है: | ||
| Line 159: | Line 153: | ||
f = A x + B y | f = A x + B y | ||
</math> | </math> | ||
सामान्यीकृत व्युत्क्रम | सामान्यीकृत व्युत्क्रम <math>B^+</math> को संतोषजनक मानें | ||
:<math>B B^+ B = B </math> | :<math>B B^+ B = B </math> | ||
परिभाषित करना | परिभाषित करना | ||
:<math>Q=[I-BB^+]</math> | :<math>Q=[I-BB^+]</math> | ||
जिससे | |||
:<math>QB = 0</math> | :<math>QB = 0</math> | ||
इस तरह, | इस तरह, | ||
| Line 175: | Line 169: | ||
:<math>QQ = Q </math> | :<math>QQ = Q </math> | ||
:<math>QQQ = Q </math> | :<math>QQQ = Q </math> | ||
<math>Q</math> | <math>Q</math>, <math>Q</math> का सामान्यीकृत व्युत्क्रम है , अर्थात। | ||
:<math>Q^+ = Q </math> | :<math>Q^+ = Q </math> | ||
समाधान को इस प्रकार सरल बनाया जा सकता है | समाधान को इस प्रकार सरल बनाया जा सकता है | ||
:<math>x_1 = A^+ Q f</math>. | :<math>x_1 = A^+ Q f</math>. | ||
गणित का | गणित का प्रश्न | ||
आव्यूह <math>A^+Q = A^+ [I-BB^+]</math> को आव्यूह <math>{\begin{bmatrix} | |||
A & B \\ | A & B \\ | ||
C & D \\ | C & D \\ | ||
\end{bmatrix}} </math> | \end{bmatrix}} </math> के स्थानीय व्युत्क्रम के रूप में जाना जाता है। इस प्रकार <math>A</math> के अनुरूप इसे स्थानीय व्युत्क्रम विधि के रूप में जाना जाता है। <ref name="shuangrenZhao 2011" /> | ||
=== पुनरावृत्तीय पुनर्निर्माण विधि === | === पुनरावृत्तीय पुनर्निर्माण विधि === | ||
यहां एक लक्ष्य | यहां एक लक्ष्य फलन परिभाषित किया गया है, और यह विधि लक्ष्य को पुनरावृत्त रूप से प्राप्त करती है। यदि लक्ष्य फलन किसी प्रकार का सामान्य हो सकता है, तो इसे न्यूनतम मानक विधि के रूप में जाना जाता है। | ||
:<math> \min( R\|x\| + S\|y\|+T\|g\|)</math>, | :<math> \min( R\|x\| + S\|y\|+T\|g\|)</math>, | ||
| Line 207: | Line 201: | ||
\end{bmatrix} | \end{bmatrix} | ||
</math> | </math> | ||
और | और <math> f</math> ज्ञात है, | ||
जहां <math>R</math>, <math>S</math> और <math> T </math> न्यूनतमकरण और <math>\|\cdot\|</math> के भार स्थिरांक हैं एक प्रकार का आदर्श है. अधिकांशतः उपयोग किए जाने वाले मानदंड <math>L_0</math>, <math>L_1</math>, <math>L_2</math>, <math>L_{+\infty}</math> कुल भिन्नता (टीवी) मानदंड या उपरोक्त मानदंडों का संयोजन हैं। इस विधि का एक उदाहरण उत्तल सेट (पीओसीएस) विधि पर प्रक्षेपण है।<ref>Ye. Yangbo, Yu. 1 Hengyong 2 and GeWang, Exact Interior Reconstruction from Truncated Limited Angle | |||
Projection | Projection | ||
Data, International Journal of Biomedical Imaging (2008), 1–6. | Data, International Journal of Biomedical Imaging (2008), 1–6. | ||
</ref> | </ref><ref>L. Zeng, B. Liu, L. Liu and C. Xiang, A new iterative reconstruction algorithm for 2D exterior fanbeam | ||
<ref>L. Zeng, B. Liu, L. Liu and C. Xiang, A new iterative reconstruction algorithm for 2D exterior fanbeam | |||
CT, Journal of | CT, Journal of | ||
XRay | XRay | ||
Science and Technology 18 (2010), 267–277.</ref> | Science and Technology 18 (2010), 267–277.</ref> | ||
===विश्लेषणात्मक समाधान=== | ===विश्लेषणात्मक समाधान=== | ||
विशेष परिस्थितियों में, आंतरिक पुनर्निर्माण को एक विश्लेषणात्मक समाधान के रूप में प्राप्त किया जा सकता है; | विशेष परिस्थितियों में, आंतरिक पुनर्निर्माण को एक विश्लेषणात्मक समाधान के रूप में प्राप्त किया जा सकता है; जिसमे <math>x</math> का समाधान है जो ऐसे स्थितियों में स्पष्ट है.<ref>Y. Zou and X. Pan, 2004, Exact image reconstruction on PIlines | ||
from minimum data in helical conebeam | from minimum data in helical conebeam | ||
CT, Physics | CT, Physics | ||
| Line 233: | Line 223: | ||
===तेज़ एक्सट्रपलेशन=== | ===तेज़ एक्सट्रपलेशन=== | ||
एक्सट्रपोलेटेड डेटा | एक्सट्रपोलेटेड डेटा अधिकांशतः [[सकारात्मक-निश्चित कर्नेल|धनात्मक-निश्चित कर्नेल]] में [[कनवल्शन]] होता है। डेटा को एक्सट्रपलेशन के बाद इसका आकार N गुना बढ़ जाता है, जहां N = 2 ~ 3. यदि डेटा को ज्ञात कर्नेल फलन में संयोजित करने की आवश्यकता है, तो संख्यात्मक गणना log(''N'')·''N'' गुना बढ़ जाएगी, यहां तक कि तेज़ [[फास्ट फूरियर ट्रांसफॉर्म]] (FFT) के साथ भी एक एल्गोरिदम उपयोग है, जो एक्सट्रपलेटेड डेटा के भाग से योगदान की विश्लेषणात्मक गणना करता है। मूल कनवल्शन गणना की तुलना में गणना समय छोड़ा जा सकता है; इस [[ कलन विधि |एल्गोरिथ्म]] के साथ, एक्सट्रपोलेटेड डेटा का उपयोग करके कनवल्शन की गणना में उल्लेखनीय वृद्धि नहीं होती है। इसे तीव्र एक्सट्रपलेशन के रूप में जाना जाता है।<ref> | ||
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</ref> | 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</ref> | ||
==विधियों की तुलना== | |||
== | |||
एक्सट्रपलेशन विधि ऐसी स्थिति में उपयुक्त होती है | एक्सट्रपलेशन विधि ऐसी स्थिति में उपयुक्त होती है | ||
:: <math> |x| > |y| </math> और <math> |X| > |Y|</math> | :: <math> |x| > |y| </math> और <math> |X| > |Y|</math> | ||
: | : अथार्त एक छोटी ट्रिमिंग कलाकृतियों की स्थिति है । | ||
अनुकूली एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है | अनुकूली एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है | ||
:: <math> |x| \sim |y| </math> और <math> |X| \sim |Y|</math> | :: <math> |x| \sim |y| </math> और <math> |X| \sim |Y|</math> | ||
: | : अथार्त एक सामान्य ट्रंकेशन कलाकृतियों की स्थिति यह विधि बाहरी क्षेत्र के लिए एक समान्य समाधान भी प्रस्तुत करती है। | ||
पुनरावृत्त एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है | पुनरावृत्त एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है | ||
:: <math> |x| \sim |y| </math> और <math> |X| \sim |Y|</math> | :: <math> |x| \sim |y| </math> और <math> |X| \sim |Y|</math> | ||
: | : अथार्त एक सामान्य ट्रंकेशन कलाकृतियों की स्थिति यद्यपि यह विधि अनुकूली पुनर्निर्माण की तुलना में उत्तम आंतरिक पुनर्निर्माण प्राप्त करती है, किंतु बाहरी क्षेत्र में इसका परिणाम नहीं मिलता है। | ||
ऐसी स्थिति के लिए स्थानीय टोमोग्राफी उपयुक्त है | ऐसी स्थिति के लिए स्थानीय टोमोग्राफी उपयुक्त है | ||
:: <math> |x| \ll |y| </math> और <math> |X| \ll |Y|</math> | :: <math> |x| \ll |y| </math> और <math> |X| \ll |Y|</math> | ||
: | :अर्थात सबसे बड़ी ट्रंकेशन कलाकृतियों की स्थिति। चूँकि इस पद्धति में कोई ट्रिमिंग कलाकृतियाँ नहीं हैं, पुनर्निर्माण में एक निश्चित त्रुटि (<math>|y|</math> के मान से स्वतंत्र) है। | ||
स्थानीय व्युत्क्रम विधि | स्थानीय व्युत्क्रम विधि स्थानीय टोमोग्राफी के समान ऐसी स्थिति में उपयुक्त है | ||
:: <math> |x| \ll |y| </math> और <math> |X| \ll |Y|</math> | :: <math> |x| \ll |y| </math> और <math> |X| \ll |Y|</math> | ||
: | :अथार्त सबसे बड़ी ट्रंकेशन कलाकृतियों की स्थिति। चूँकि इस विधि के लिए कोई खंडन कलाकृतियाँ नहीं हैं, पुनर्निर्माण में एक निश्चित त्रुटि (<math>|y|</math> के मान से स्वतंत्र) है जो स्थानीय टोमोग्राफी से छोटी हो सकती है। | ||
पुनरावृत्तीय पुनर्निर्माण विधि बड़ी गणनाओं के साथ अच्छा परिणाम प्राप्त करती है। यद्यपि विश्लेषणात्मक विधि | पुनरावृत्तीय पुनर्निर्माण विधि बड़ी गणनाओं के साथ अच्छा परिणाम प्राप्त करती है। यद्यपि विश्लेषणात्मक विधि स्पष्ट परिणाम प्राप्त करती है, यह केवल कुछ स्थितियों में ही कार्यात्मक होती है। तेज़ एक्सट्रपलेशन विधि अन्य एक्सट्रपलेशन विधियों के समान परिणाम प्राप्त कर सकती है, और गणना को कम करने के लिए उपरोक्त आंतरिक पुनर्निर्माण विधियों पर प्रयुक्त किया जा सकता है। | ||
==यह भी देखें== | ==यह भी देखें == | ||
{{Wiktionary|extrapolation}} | {{Wiktionary|extrapolation}} | ||
* | *पूर्वानुमान | ||
*[[न्यूनतम बहुपद एक्सट्रपलेशन]] | *[[न्यूनतम बहुपद एक्सट्रपलेशन]] | ||
*[[मल्टीग्रिड विधि]] | *[[मल्टीग्रिड विधि]] | ||
*[[भविष्यवाणी अंतराल]] | *[[भविष्यवाणी अंतराल|पूर्वानुमान अंतराल]] | ||
*[[प्रतिगमन विश्लेषण]] | *[[प्रतिगमन विश्लेषण]] | ||
*[[रिचर्डसन एक्सट्रपलेशन]] | *[[रिचर्डसन एक्सट्रपलेशन]] | ||
| Line 274: | Line 262: | ||
*[[मृत गणना]] | *[[मृत गणना]] | ||
*[[छवि पुनर्निर्माण]] | *[[छवि पुनर्निर्माण]] | ||
*[[स्थानीय उलटा]] | *[[स्थानीय उलटा|स्थानीय व्युत्क्रम]] | ||
*[[सामान्यीकृत उलटा]] | *[[सामान्यीकृत उलटा|सामान्यीकृत व्युत्क्रम]] | ||
*एक्सट्रपोलेशन | *एक्सट्रपोलेशन | ||
==टिप्पणियाँ== | ==टिप्पणियाँ == | ||
{{reflist}} | {{reflist}} | ||
[[Category:Created On 25/07/2023]] | [[Category:Created On 25/07/2023]] | ||
[[Category:Machine Translated Page]] | |||
[[Category:Pages with script errors]] | |||
[[Category:चिकित्सा छवि]] | |||
Latest revision as of 15:05, 6 September 2023
डिजिटल इमेजिंग में पुनरावृत्त पुनर्निर्माण में, आंतरिक पुनर्निर्माण (जिसे दृश्य के सीमित क्षेत्र (एलएफवी) पुनर्निर्माण के रूप में भी जाना जाता है) छवि डेटा को दृश्य के एक छोटे क्षेत्र तक सीमित करने के कारण होने वाली ट्रंकेशन कलाकृतियों को सही करने की एक तकनीक है। पुनर्निर्माण एक ऐसे क्षेत्र पर केंद्रित है जिसे रुचि का क्षेत्र (आरओआई) कहा जाता है। यद्यपि आंतरिक पुनर्निर्माण को दंत या कार्डियक एक्स-रे कंप्यूटेड टोमोग्राफी छवियों पर प्रयुक्त किया जा सकता है, यह अवधारणा सीटी तक सीमित नहीं है। इसे कई विधियों में से एक के साथ प्रयुक्त किया जाता है।
विधियों
प्रत्येक विधि का उद्देश्य सदिश का समाधान करना है निम्नलिखित समस्या में:
मान लीजिए कि रुचि का क्षेत्र (आरओआई) है और , के बाहर का क्षेत्र है। मान लीजिए कि , , , ज्ञात आव्यूह हैं; और मूल छवि के अज्ञात सदिश हैं, जबकि और प्रतिक्रियाओं के सदिश माप हैं ( ज्ञात है और अज्ञात है)। क्षेत्र के अंदर है और , क्षेत्र में,, क्षेत्र के बाहर है। , के अनुरूप माप में एक क्षेत्र के अंदर है। इस क्षेत्र को , ( के रूप में दर्शाया गया है, जबकि g है क्षेत्र F के बाहर। यह से मेल खाता है और इसे , (). के रूप में दर्शाया गया है
सीटी छवि-पुनर्निर्माण उद्देश्यों के लिए, .
आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, आव्यूह , , , सम्मिश्र ऑपरेटर (गणित) के अतिरिक्त छवि पुनर्निर्माण के लिए प्रयुक्त किया जाता है।
नीचे सूचीबद्ध पहली आंतरिक-पुनर्निर्माण विधि एक्सट्रपलेशन है। यह एक स्थानीय टोमोग्राफी विधि है जो ट्रंकेशन कलाकृतियों को समाप्त करती है किंतु एक अन्य प्रकार की कलाकृतियों का परिचय देती है: एक कटोरा प्रभाव एक सुधार को अनुकूली एक्सट्रपलेशन विधि के रूप में जाना जाता है, चूँकि नीचे दी गई पुनरावृत्त एक्सट्रपलेशन विधि भी पुनर्निर्माण परिणामों में सुधार करती है। कुछ स्थितियों में, आंतरिक पुनर्निर्माण के लिए स्पष्ट पुनर्निर्माण पाया जा सकता है। नीचे दी गई स्थानीय विपरीत विधि स्थानीय टोमोग्राफी विधि को संशोधित करती है, और स्थानीय टोमोग्राफी के पुनर्निर्माण परिणाम में सुधार कर सकती है; पुनरावृत्तीय पुनर्निर्माण पद्धति को आंतरिक पुनर्निर्माण पर प्रयुक्त किया जा सकता है। उपरोक्त विधियों में, एक्सट्रपलेशन को अधिकांशतः प्रयुक्त किया जाता है।
एक्सट्रपलेशन विधि
:
, , , ज्ञात आव्यूह हैं; और अज्ञात सदिश हैं; एक ज्ञात सदिश है, और एक अज्ञात सदिश है। हमें सदिश को जानना होगा। और मूल छवि हैं, जबकि और प्रतिक्रियाओं के माप हैं। सदिश रुचि के क्षेत्र , (). के अंदर है। सदिश , क्षेत्र के बाहर है। बाहरी क्षेत्र को , () कहा जाता है और , के अनुरूप माप में एक क्षेत्र के अंदर है। इस क्षेत्र को , () दर्शाया गया है। सदिश का क्षेत्र (क्षेत्र के बाहर) भी Y से मेल खाता है और इसे , (). के रूप में दर्शाया गया है। सीटी छवि पुनर्निर्माण में यह है
आंतरिक पुनर्निर्माण की अवधारणा को सरल बनाने के लिए, आव्यूह , , , एक सम्मिश्र ऑपरेटर के अतिरिक्त छवि पुनर्निर्माण के लिए प्रयुक्त किया जाता है।
बाहरी क्षेत्र में प्रतिक्रिया का अनुमान हो सकता है; उदाहरण के लिए, मान लें कि यह है।
का एक समाधान के रूप में लिखा जाता है, और इसे एक्सट्रपलेशन विधि के रूप में जाना जाता है। परिणाम इस बात पर निर्भर करता है कि एक्सट्रपलेशन फलन कितना अच्छा है। एक अधिकांशतः पसंद है
दो क्षेत्रों की सीमा पर.[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(N)·N गुना बढ़ जाएगी, यहां तक कि तेज़ फास्ट फूरियर ट्रांसफॉर्म (FFT) के साथ भी एक एल्गोरिदम उपयोग है, जो एक्सट्रपलेटेड डेटा के भाग से योगदान की विश्लेषणात्मक गणना करता है। मूल कनवल्शन गणना की तुलना में गणना समय छोड़ा जा सकता है; इस एल्गोरिथ्म के साथ, एक्सट्रपोलेटेड डेटा का उपयोग करके कनवल्शन की गणना में उल्लेखनीय वृद्धि नहीं होती है। इसे तीव्र एक्सट्रपलेशन के रूप में जाना जाता है।[19]
विधियों की तुलना
एक्सट्रपलेशन विधि ऐसी स्थिति में उपयुक्त होती है
- और
- अथार्त एक छोटी ट्रिमिंग कलाकृतियों की स्थिति है ।
अनुकूली एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है
- और
- अथार्त एक सामान्य ट्रंकेशन कलाकृतियों की स्थिति यह विधि बाहरी क्षेत्र के लिए एक समान्य समाधान भी प्रस्तुत करती है।
पुनरावृत्त एक्सट्रपलेशन विधि ऐसी स्थिति के लिए उपयुक्त है
- और
- अथार्त एक सामान्य ट्रंकेशन कलाकृतियों की स्थिति यद्यपि यह विधि अनुकूली पुनर्निर्माण की तुलना में उत्तम आंतरिक पुनर्निर्माण प्राप्त करती है, किंतु बाहरी क्षेत्र में इसका परिणाम नहीं मिलता है।
ऐसी स्थिति के लिए स्थानीय टोमोग्राफी उपयुक्त है
- और
- अर्थात सबसे बड़ी ट्रंकेशन कलाकृतियों की स्थिति। चूँकि इस पद्धति में कोई ट्रिमिंग कलाकृतियाँ नहीं हैं, पुनर्निर्माण में एक निश्चित त्रुटि ( के मान से स्वतंत्र) है।
स्थानीय व्युत्क्रम विधि स्थानीय टोमोग्राफी के समान ऐसी स्थिति में उपयुक्त है
- और
- अथार्त सबसे बड़ी ट्रंकेशन कलाकृतियों की स्थिति। चूँकि इस विधि के लिए कोई खंडन कलाकृतियाँ नहीं हैं, पुनर्निर्माण में एक निश्चित त्रुटि ( के मान से स्वतंत्र) है जो स्थानीय टोमोग्राफी से छोटी हो सकती है।
पुनरावृत्तीय पुनर्निर्माण विधि बड़ी गणनाओं के साथ अच्छा परिणाम प्राप्त करती है। यद्यपि विश्लेषणात्मक विधि स्पष्ट परिणाम प्राप्त करती है, यह केवल कुछ स्थितियों में ही कार्यात्मक होती है। तेज़ एक्सट्रपलेशन विधि अन्य एक्सट्रपलेशन विधियों के समान परिणाम प्राप्त कर सकती है, और गणना को कम करने के लिए उपरोक्त आंतरिक पुनर्निर्माण विधियों पर प्रयुक्त किया जा सकता है।
यह भी देखें
- पूर्वानुमान
- न्यूनतम बहुपद एक्सट्रपलेशन
- मल्टीग्रिड विधि
- पूर्वानुमान अंतराल
- प्रतिगमन विश्लेषण
- रिचर्डसन एक्सट्रपलेशन
- स्थिर विश्लेषण
- रुझान अनुमान
- प्रक्षेप
- एक्सट्रपलेशन डोमेन विश्लेषण
- मृत गणना
- छवि पुनर्निर्माण
- स्थानीय व्युत्क्रम
- सामान्यीकृत व्युत्क्रम
- एक्सट्रपोलेशन
टिप्पणियाँ
- ↑ 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.
- ↑ 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.
- ↑ 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.
- ↑ 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.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.
- ↑ 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.
- ↑ 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.
- ↑ 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
- ↑ 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.
- ↑ 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.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
- ↑ A. Faridani, E.L. Ritman and K.T. Smith, Local tomography, SIAM J APPL MATH 52 (1992), 459–484.
- ↑ A. Katsevich, 1999 Cone beam local tomography, SIAM J APPL MATH 59, 2224–2246.
- ↑ 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.
- ↑ 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.
- ↑ 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.
- ↑ M. Defrise, F. Noo, R. Clackdoyle and H. Kudo, Truncated Hilbert transform and image reconstruction from limited tomographic data. IOPscience.iop.org, 2006
- ↑ F. Noo, R. Clackdoyle and J.D. Pack, A twostep Hilbert transform method for 2D image reconstruction, Phys Med Biol 49 (2004), 3903–3923.
- ↑ 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