प्रस्तावना: छायाओं के माध्यम से अदृश्य गहराइयों को देखना
हमारे दैनिक जीवन में कार्य-कारण का संबंध अत्यंत सीधा और सहज होता है: यदि हमें किसी वस्तु का आकार और प्रकाश स्रोत की स्थिति ज्ञात हो, तो दीवार पर पड़ने वाली छाया की रूपरेखा की गणना ज्यामितीय प्रकाशिकी के माध्यम से सरलता से की जा सकती है। इसी प्रकार, यदि पृथ्वी के आंतरिक भूगर्भीय घनत्व और लोचदार गुणों का ज्ञान हो, तो भूकंप के केंद्र से सतह तक पहुंचने वाली भूकंपीय तरंगों का सटीक अनुकरण तरंग समीकरण की सहायता से किया जा सकता है। इस प्रक्रिया को, जहां कारण (आंतरिक भौतिक पैरामीटर) से प्रभाव (बाहरी प्रेक्षण डेटा) प्राप्त किया जाता है, अग्रगामी समस्या (Forward Problem) कहा जाता है। अग्रगामी समस्याएं स्वभाव से स्थिर और अद्वितीय हल वाली होती हैं।
परंतु आधुनिक चिकित्सा विज्ञान, भूभौतिकी और खगोल विज्ञान में हमारे समक्ष इसके ठीक विपरीत चुनौती होती है:
- चिकित्सीय सीटी स्कैन (CT Scan): रोगी के शरीर को चीरे बिना, केवल विभिन्न कोणों से शरीर से गुजरने वाली एक्स-रे किरणों के क्षीणन डेटा (1D छाया) से आंतरिक अंगों, कोमल ऊतकों और सूक्ष्म ट्यूमर का 3D चित्र कैसे तैयार किया जाए?
- भूभौतिकीय अन्वेषण: सतह पर दर्ज भूकंपीय परावर्तन तरंगों के आधार पर जमीन से कई किलोमीटर नीचे तेल और गैस के भंडारों या विवर्तनिक दरारों का मानचित्रण कैसे किया जाए?
- खगोलीय रेडियो इंटरफेरोमेट्री: इवेंट होराइजन टेलीस्कोप (EHT) ने अत्यंत विरल डेटा से सुपरमैसिव ब्लैक होल के अभिवृद्धि चक्र की पहली छवि कैसे बनाई?
बाहरी प्रेक्षित डेटा से आंतरिक भौतिक कारण को ज्ञात करने के इस गणितीय कार्य को प्रतिलोम समस्या (Inverse Problem) कहा जाता है।
flowchart LR
subgraph Forward ["अग्रगामी समस्या (Forward Problem)"]
direction TB
A["भौतिक कारण f(x, y)<br/>(आंतरिक ऊतक घनत्व)"] -->|भौतिक ऑपरेटर A| B["प्रेक्षित डेटा g(s, θ)<br/>(एक्स-रे प्रक्षेपण छाया)"]
end
subgraph Inverse ["प्रतिलोम समस्या (Inverse Problem)"]
direction TB
C["बाहरी प्रेक्षण g(s, θ)<br/>(शोरयुक्त माप)"] -->|प्रतिलोम ऑपरेटर A⁻¹ / नियमितीकरण| D["पुनर्निर्मित संरचना f*(x, y)<br/>(सटीक आंतरिक क्रॉस-सेक्शनल छवि)"]
end
1. प्रतिलोम समस्या क्या है?：हाडामार्ड सुस्थितता और अस्वस्थता की त्रासदी
प्रतिलोम समस्याओं को हल करना अग्रगामी समस्याओं की तुलना में अत्यधिक कठिन क्यों है?
1.1 ऑपरेटर समीकरण का स्वरूप
$$\mathcal{A} f = g$$जहाँ $f \in \mathcal{X}$ आंतरिक अवशोषण गुणांक है, $g \in \mathcal{Y}$ मापा गया प्रक्षेपण डेटा है, और $\mathcal{A}$ एक रैखिक समाकल ऑपरेटर है।
1.2 हाडामार्ड के सुस्थितता मानदंड (1902)
फ्रांसीसी गणितज्ञ जाक हाडामार्ड ने कहा कि कोई समस्या सुस्थित (well-posed) तभी कहलाती है जब वह निम्नलिखित तीन शर्तों को पूरा करे:
- अस्तित्व (Existence): प्रत्येक डेटा $g$ के लिए कम से कम एक हल $f$ अवश्य हो।
- अद्वितीयता (Uniqueness): हल अद्वितीय हो ($\ker(\mathcal{A}) = \{0\}$)।
- स्थिरता (Stability): इनपुट डेटा में सूक्ष्म परिवर्तन होने पर हल में कोई अनियंत्रित विस्फोट न हो: $\|\delta f\| \le C \|\delta g\|$।
यदि इनमें से कोई भी शर्त पूरी नहीं होती, तो समस्या अस्वस्थ (ill-posed) कहलाती है।
1.3 कॉम्पैक्ट ऑपरेटर और उच्च-आवृत्ति शोर का संकट
$$g(s) = \int_a^b K(s, t) f(t) \, dt$$$$f = \sum_{n=1}^\infty \frac{1}{\sigma_n} \langle g, u_n \rangle v_n$$माप में उपस्थित मामूली सा शोर भी $\frac{1}{\sigma_n} \to \infty$ के कारण हजारों गुना बढ़ जाता है, जिससे बिना नियमितीकरण के पूरी छवि नष्ट हो जाती है।
2. समाकल ज्यामिति: जोहान रेडॉन (1917)
ऑस्ट्रियाई गणितज्ञ जोहान रेडॉन ने 1917 में एक ज्यामितीय प्रश्न का समाधान प्रस्तुत किया:
“यदि हमें किसी समतल में सभी संभव सरल रेखाओं पर किसी फलन के रेखा-समाकल ज्ञात हों, तो क्या मूल फलन को अद्वितीय रूप से पुनर्निर्मित किया जा सकता है?”
रेडॉन रूपांतरण की परिभाषा
$$x \cos\theta + y \sin\theta = s$$$$(\mathcal{R}f)(s, \theta) = \iint_{\mathbb{R}^2} f(x, y) \, \delta(x \cos\theta + y \sin\theta - s) \, dx dy$$चिकित्सा इमेजिंग में $(\mathcal{R}f)(s, \theta)$ के 2D निरूपण को साइनोग्राम (Sinogram) कहा जाता है।
3. सीटी स्कैन का भौतिक आधार: बीयर-लैम्बर्ट नियम
$$\frac{dI}{dt} = - f(\mathbf{x}) I \implies I_{\text{det}} = I_0 \exp\left( - \int_L f(\mathbf{x}) \, dt \right)$$$$g(s, \theta) = \ln\left( \frac{I_0}{I_{\text{det}}} \right) = \int_L f(\mathbf{x}) \, dt = (\mathcal{R}f)(s, \theta)$$अतः सीटी स्कैन की समस्या गणितीय रूप से रेडॉन रूपांतरण के प्रतिलोमन के समतुल्य है।
4. फूरियर स्लाइस प्रमेय (Fourier Slice Theorem)
$$\left( \mathcal{F}_1 [\mathcal{R}f(\cdot, \theta)] \right)(\omega) = \left( \mathcal{F}_2 f \right)(\omega \cos\theta, \; \omega \sin\theta)$$यह दर्शाता है कि कोण $\theta$ पर 1D प्रक्षेपण का फूरियर रूपांतरण, मूल छवि के 2D फूरियर रूपांतरण के उसी कोण पर मूल बिंदु से गुजरने वाले स्लाइस के बिल्कुल बराबर होता है।
5. फ़िल्टर्ड बैकप्रोजेक्शन (Filtered Backprojection, FBP)
$$f(x, y) = \int_0^\pi \left[ \int_{-\infty}^\infty G(\omega, \theta) \cdot |\omega| \cdot e^{i 2\pi \omega (x\cos\theta + y\sin\theta)} \, d\omega \right] d\theta$$- रैंप फ़िल्टरिंग (Ramp Filter): $|\omega|$ से गुणा करके साधारण बैकप्रोजेक्शन के धुंधलेपन ($1/r$) को समाप्त किया जाता है।
- बैकप्रोजेक्शन: फ़िल्टर किए गए डेटा को सभी कोणों पर छवि ग्रिड पर वापस फैलाकर जोड़ा जाता है।
6. अस्वस्थता का समाधान: तिखोनोव नियमितीकरण
$$\min_{\mathbf{f}} \left\{ \|\mathbf{A}\mathbf{f} - \mathbf{g}\|_2^2 + \alpha \|\mathbf{L}\mathbf{f}\|_2^2 \right\}$$$$\mathbf{f}_\alpha = (\mathbf{A}^T \mathbf{A} + \alpha \mathbf{L}^T \mathbf{L})^{-1} \mathbf{A}^T \mathbf{g}$$फ़िल्टर कारक $\frac{\sigma_i^2}{\sigma_i^2 + \alpha}$ शून्य के निकट वाले सिंगुलर मानों को नियंत्रित करके शोर के विस्फोट को रोकता है।
7. आधुनिक क्षितिज: टोटल वेरिएशन और कम्प्रेस्ड सेंसिंग
- टोटल वेरिएशन (TV): ग्रेडिएंट के $L_1$ मानक $\|\nabla \mathbf{f}\|_1$ का उपयोग करके किनारों को धुंधला किए बिना शोर को समाप्त करता है।
- कम्प्रेस्ड सेंसिंग (Compressed Sensing): विरलता (sparsity) के सिद्धांत का उपयोग करके बहुत कम प्रक्षेपण कोणों से भी पूर्ण गुणवत्ता वाली छवि का पुनर्निर्माण करता है, जिससे एक्स-रे विकिरण का खतरा 80% से अधिक कम हो जाता है।
- डीप अनफोल्डिंग (Deep Unfolding): रेडॉन ऑपरेटर को न्यूरल नेटवर्क में समाहित करके अत्यंत तेज़ और भौतिक रूप से विश्वसनीय पुनर्निर्माण प्रदान करता है।
निष्कर्ष
प्लेटो की गुफा के रूपक की भांति, जहां केवल छायाएं दिखाई देती हैं, प्रतिलोम समस्याओं का गणित मानव बुद्धि की असीम शक्ति का प्रमाण है: केवल बाहरी छायाओं के विश्लेषण से हम प्रकृति और मानव शरीर के भीतर छिपी जटिल वास्तविकताओं को पूर्ण स्पष्टता के साथ उजागर कर सकते हैं।