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

बाएं से: मूल छवि, धुंधली छवि, वीनर विसंवलन का उपयोग करके छवि को धुंधला करना।

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

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

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

परिभाषा

एक प्रणाली दी गई है:

$\ 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{aligned}G(f)&={\frac {1}{H(f)}}\left[{\frac {1}{1+1/(|H(f)|^{2}\mathrm {SNR} (f))}}\right]\end{aligned}}$

यहाँ, $\ 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{aligned}\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{aligned}}$

यदि हम द्विघात का विस्तार करें, तो हमें निम्नलिखित मिलता है:

${\begin{aligned}\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{aligned}}$

हालाँकि, हम यह मान रहे हैं कि नॉइज़ सिग्नल से स्वतंत्र है, इसलिए:

$\ \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.