वीनर डिकोनवोल्यूशन

गणित में, वीनर विसंवलन, विसंवलन में निहित शोर समस्याओं के लिए विनीज़ निस्यन्दक का एक अनुप्रयोग है। यह आवृत्ति कार्यक्षेत्र में काम करता है, उन आवृति पर डिकंवॉल्व्ड शोर के प्रभाव को कम करने का प्रयास करता है जिनमें संकेत-से-शोर अनुपात खराब होता है।

वीनर विसंवलन विधि का छवि विसंवलन अनुप्रयोगों में व्यापक उपयोग होता है, क्योंकि अधिकांश दृश्य छवियों का आवृत्ति वर्णक्रम काफी अच्छी तरह से व्यवहार किया जाता है और आसानी से अनुमान लगाया जा सकता है।

वीनर विसंवलन का नाम नॉर्बर्ट वीनर के नाम पर रखा गया है।

परिभाषा
एक प्रणाली दी गई:


 * $$\ y(t) = (h*x)(t) + n(t)$$

जहाँ $$*$$ संवलन को दर्शाता है और:


 * $$\ x(t)$$ समय$$\ t $$ पर कुछ मूल संकेत (अज्ञात) है।
 * $$\ h(t)$$ रैखिक समय-अपरिवर्तनीय प्रणाली की ज्ञात आवेग प्रतिक्रिया है।
 * $$\ n(t)$$ कुछ अज्ञात योगात्मक शोर है, जो$$\ x(t)$$ से स्वतंत्र है।
 * $$\ y(t)$$ हमारा देखा हुआ संकेत है।

हमारा लक्ष्य कुछ $$\ g(t)$$ खोजना है ताकि हम निम्नलिखित नुसार $$\ x(t)$$ का अनुमान लगा सकें:


 * $$\ \hat{x}(t) = (g*y)(t)$$

जहाँ $$\ \hat{x}(t)$$ का एक अनुमान $$\ x(t)$$ है जो माध्य वर्ग त्रुटि को न्यूनतम करता है


 * $$\ \epsilon(t) = \mathbb{E} \left| x(t) - \hat{x}(t) \right|^2$$,

साथ ही$$\ \mathbb{E}$$ अपेक्षित मूल्य को दर्शाता है। वीनर विसंवलन निस्यन्दक $$\ g(t)$$ प्रदान करता है। निस्यन्दक को आवृति कार्यक्षेत्र में सबसे आसानी से वर्णित किया गया है:


 * $$\ G(f) = \frac{H^*(f)S(f)}{ |H(f)|^2 S(f) + N(f) }$$

जहाँ:


 * $$\ G(f)$$ और $$\ H(f)$$ के फूरियर रूपांतरण$$\ g(t)$$ और$$\ h(t)$$ हैं,
 * $$\ S(f) = \mathbb{E}|X(f)|^2 $$ मूल संकेत का औसत पावर वर्णक्रमीय घनत्व $$\ x(t)$$ है,
 * $$\ N(f) = \mathbb{E}|V(f)|^2 $$ शोर का औसत शक्ति वर्णक्रमीय घनत्व $$\ n(t)$$ है,
 * $$X(f)$$, $$Y(f)$$, और $$V(f)$$ के फूरियर रूपांतरण $$x(t)$$, और $$y(t)$$हैं, क्रमश,
 * अधिलेख $${}^*$$ जटिल संयुग्म को दर्शाता है।

निस्यंदन संचालन या तो समय-कार्यक्षेत्र में, जैसा कि ऊपर किया गया है, या आवृति कार्यक्षेत्र में किया जा सकता है:


 * $$\ \hat{X}(f) = G(f)Y(f)$$

और फिर उलटा फूरियर रूपांतरण $$\ \hat{X}(f)$$ प्राप्त करने के लिए $$\ \hat{x}(t)$$ निष्पादित करना

ध्यान दें कि छवियों की स्तिथि में, तर्क $$\ t $$ और $$\ f $$ द्वि-आयामी हो जाते हैं; हालाँकि परिणाम वही है।

व्याख्या
वीनर निस्यन्दक का संचालन तब स्पष्ट हो जाता है जब उपरोक्त निस्यन्दक समीकरण को फिर से लिखा जाता है:



\begin{align} G(f) & = \frac{1}{H(f)} \left[ \frac{ 1 }{ 1 + 1/(|H(f)|^2 \mathrm{SNR}(f))} \right] \end{align} $$ यहाँ, $$\ 1/H(f)$$ मूल प्रणाली का उलटा है, $$\ \mathrm{SNR}(f) = S(f)/N(f)$$ संकेत-से-शोर अनुपात है, और $$\ |H(f)|^2 \mathrm{SNR}(f)$$ शोर वर्णक्रमीय घनत्व के लिए शुद्ध निस्यन्दक किए गए संकेत का अनुपात है। जब शून्य शोर (यानी अनंत संकेत-से-शोर) होता है, तो वर्गाकार कोष्ठक के अंदर का शब्द 1 के बराबर होता है, जिसका अर्थ है कि वीनर निस्यन्दक प्रणाली का उलटा है, जैसा कि हम उम्मीद कर सकते हैं। हालाँकि, जैसे-जैसे कुछ आवृत्तियों पर शोर बढ़ता है, संकेत-से-शोर अनुपात कम हो जाता है, इसलिए वर्ग कोष्ठक के अंदर का शब्द भी कम हो जाता है। इसका मतलब यह है कि वीनर निस्यन्दक उनके निस्यन्दक किए गए संकेत-से-शोर अनुपात के अनुसार आवृत्तियों को क्षीण करता है।

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

व्युत्पत्ति
जैसा कि ऊपर उल्लेख किया गया है, हम मूल संकेत का एक अनुमान तैयार करना चाहते हैं जो माध्य वर्ग त्रुटि को कम करता है, जिसे निम्नलिखित रूप से व्यक्त किया जा सकता है:


 * $$\ \epsilon(f) = \mathbb{E} \left| X(f) - \hat{X}(f) \right|^2$$

$$\epsilon$$ की पिछली परिभाषा के समतुल्य, फूरियर रूपांतरण के लिए प्लांचरेल प्रमेय या पार्सेवल के प्रमेय का उपयोग करके प्राप्त किया जा सकता है।

यदि हम अभिव्यक्ति में इसके लिए $$\ \hat{X}(f)$$ स्थानापन्न करते हैं, उपरोक्त को पुनर्व्यवस्थित किया जा सकता है



\begin{align} \epsilon(f) & = \mathbb{E} \left| X(f) - G(f)Y(f) \right|^2 \\ & = \mathbb{E} \left| X(f) - G(f) \left[ H(f)X(f) + V(f) \right] \right|^2 \\ & = \mathbb{E} \big| \left[ 1 - G(f)H(f) \right] X(f) - G(f)V(f) \big|^2 \end{align} $$ यदि हम द्विघात का विस्तार करें, तो हमें निम्नलिखित मिलता है:



\begin{align} \epsilon(f) & = \Big[ 1-G(f)H(f) \Big] \Big[ 1-G(f)H(f) \Big]^*\, \mathbb{E}|X(f)|^2 \\ & {} - \Big[ 1-G(f)H(f) \Big] G^*(f)\, \mathbb{E}\Big\{X(f)V^*(f)\Big\} \\ & {} - G(f) \Big[ 1-G(f)H(f) \Big]^*\, \mathbb{E}\Big\{V(f)X^*(f)\Big\} \\ & {} + G(f) G^*(f)\, \mathbb{E}|V(f)|^2 \end{align} $$ हालाँकि, हम यह मान रहे हैं कि शोर संकेत से स्वतंत्र है, इसलिए:


 * $$\ \mathbb{E}\Big\{X(f)V^*(f)\Big\} = \mathbb{E}\Big\{V(f)X^*(f)\Big\} = 0$$

शक्ति वर्णक्रमीय घनत्व$$\ S(f) $$ और $$\ N(f) $$ को प्रतिस्थापित करने पर, हमारे पास निम्न है:



\epsilon(f) = \Big[ 1-G(f)H(f) \Big]\Big[ 1-G(f)H(f) \Big]^ * S(f)  + G(f)G^*(f)N(f) $$ न्यूनतम त्रुटि मान ज्ञात करने के लिए, हम इसके संबंध में विर्टिंगर व्युत्पन्न $$\ G(f)$$ की गणना करते हैं और इसे शून्य के बराबर निर्धारित करें।



\frac{d\epsilon(f)}{dG(f)} = 0 \Rightarrow G^*(f)N(f) - H(f)\Big[1 - G(f)H(f)\Big]^* S(f) = 0 $$ इस अंतिम समानता को वीनर निस्यन्दक देने के लिए पुन: व्यवस्थित किया जा सकता है।

यह भी देखें

 * सूचना क्षेत्र सिद्धांत
 * विखंडन
 * वीनर फिल्टर
 * प्वाइंट स्प्रेड फ़ंक्शन
 * अंधा विखंडन
 * फूरियर रूपांतरण

संदर्भ

 * Rafael Gonzalez, Richard Woods, and Steven Eddins. Digital Image Processing Using Matlab. Prentice Hall, 2003.

बाहरी संबंध

 * Comparison of different deconvolution methods.