Revisiting PSF models: unifying framework and high-performance implementation
Liu, Stergiopoulou, Chuah et al.
Localization microscopy often relies on detailed models of point spread functions. For applications such as deconvolution or PSF engineering, accurate models for light propagation in imaging systems with high numerical aperture are required. Different models have been proposed based on 2D Fourier transforms or 1D Bessel integrals. The most precise ones combine a vectorial description of the electric field and precise aberration models. However, it may be unclear which model to choose, as there is no comprehensive comparison between the Fourier and Bessel approaches yet. Moreover, many existing libraries are written in Java (e.g. our previous PSF generator software) or MATLAB, which hinders the integration into deep learning algorithms. In this work, we start from the original Richards-Wolf integral and revisit both approaches in a systematic way. We present a unifying framework in which we prove the equivalence between the Fourier and Bessel strategies and detail a variety of correction factors applicable to both of them. Then, we provide a high-performance implementation of our theoretical framework in the form of an open-source library that is built on top of PyTorch, a popular library for deep learning. It enables us to benchmark the accuracy and computational speed of different models, thus allowing for an in-depth comparison of the existing models for the first time. We show that the Bessel strategy is optimal for axisymmetric beams while the Fourier approach can be applied to more general scenarios. Our work enables efficient PSF computation on CPU or GPU, which can then be included in simulation and optimization pipelines.
academic
PSF मॉडल्स का पुनर्विचार: एकीकृत ढांचा और उच्च-प्रदर्शन कार्यान्वयन
बिंदु प्रसार फलन (PSF) मॉडल स्थानीयकरण सूक्ष्मदर्शी के मूल उपकरण हैं। यह पेपर उच्च संख्यात्मक एपर्चर (NA) इमेजिंग प्रणालियों में प्रकाश प्रसार मॉडलिंग की समस्या के लिए एक एकीकृत सैद्धांतिक ढांचा प्रस्तुत करता है, जो 2D फूरियर रूपांतरण और 1D बेसेल इंटीग्रल आधारित दोनों विधियों की समानता को प्रमाणित करता है। अनुसंधान ने PyTorch आधारित एक उच्च-प्रदर्शन ओपन-सोर्स लाइब्रेरी विकसित की है, जो विभिन्न मॉडल्स का पहली बार व्यवस्थित बेंचमार्किंग करती है। प्रयोग दर्शाते हैं कि बेसेल रणनीति अक्षीय-सममित बीम के लिए इष्टतम है, जबकि फूरियर विधि अधिक सामान्य परिदृश्यों के लिए उपयुक्त है। यह कार्य CPU/GPU दक्ष संगणना का समर्थन करता है और गहन शिक्षा और अनुकूलन प्रवाह में निर्बाध रूप से एकीकृत हो सकता है।
बिंदु प्रसार फलन (PSF) प्रकाशीय सूक्ष्मदर्शी की एक मौलिक अवधारणा है, जो इमेजिंग प्रणाली की आवेग प्रतिक्रिया का वर्णन करती है। उच्च संख्यात्मक एपर्चर (NA) इमेजिंग प्रणालियों में, सटीक PSF मॉडल निम्नलिखित अनुप्रयोगों के लिए महत्वपूर्ण है:
एकल-अणु स्थानीयकरण सूक्ष्मदर्शी (SMLM): अणु स्थानीयकरण के लिए सटीक PSF की आवश्यकता
विकनवोल्यूशन सूक्ष्मदर्शी: छवि पुनर्प्राप्ति के लिए सटीक PSF मॉडल
PSF इंजीनियरिंग: विशेष PSF डिजाइन के माध्यम से सुपर-रेजोल्यूशन इमेजिंग
सैद्धांतिक स्तर: दो मुख्य विधियां (फूरियर रूपांतरण विधि और बेसेल फलन विधि) मौजूद हैं, लेकिन उनके बीच संबंध स्पष्ट नहीं है, व्यवस्थित तुलना की कमी है
कार्यान्वयन स्तर: मौजूदा लाइब्रेरीज मुख्यतः Java या MATLAB में लिखी गई हैं (जैसे लेखकों का पूर्व PSF Generator), जो आधुनिक गहन शिक्षा ढांचे में एकीकृत करना कठिन है
अनुप्रयोग स्तर: सटीकता और कम्प्यूटेशनल गति की व्यवस्थित बेंचमार्किंग की कमी है, उपयोगकर्ताओं को उपयुक्त मॉडल चुनना कठिन है
सैद्धांतिक एकीकृत ढांचा: Richards-Wolf इंटीग्रल से शुरू करते हुए, यह प्रमाणित किया गया है कि फूरियर (कार्टेशियन) और बेसेल (गोलीय) विधियां अनिवार्य रूप से एक ही प्रसार इंटीग्रल के विभिन्न पैरामीटरकरण हैं
सार्वभौमिक सुधार कारक: कई भौतिक सुधार कारकों (Gibson-Lanni गोलीय विपथन, विलोपन कारक, Fresnel संचरण गुणांक, मनमाना चरण विकृति आदि) को व्यवस्थित रूप से प्राप्त किया गया है, और दोनों विधियों पर एकीकृत रूप से लागू किया गया है
उच्च-प्रदर्शन PyTorch कार्यान्वयन: ओपन-सोर्स लाइब्रेरी psf-generator विकसित की गई है, जो चार प्रसारकों (अदिश/सदिश × कार्टेशियन/गोलीय) को लागू करती है, CPU/GPU संगणना और स्वचालित अवकलन का समर्थन करती है
व्यवस्थित बेंचमार्किंग: विभिन्न PSF मॉडल्स की सटीकता और कम्प्यूटेशनल गति की पहली बार व्यापक तुलना, व्यावहारिक अनुप्रयोग के लिए चयन मार्गदर्शन प्रदान करती है
पारिस्थितिकी तंत्र एकीकरण: napari ग्राफिकल इंटरफेस प्लगइन और chromatix ऑप्टिकल सिमुलेशन ढांचा एकीकरण प्रदान करता है, ओपन-सोर्स समुदाय अपनाने को बढ़ावा देता है
समानता प्रमाण: पहली बार कार्टेशियन और गोलीय पैरामीटरकरण को एक ही इंटीग्रल के विभिन्न प्रतिनिधित्व के रूप में कठोरता से प्रमाणित किया गया है, क्षेत्र में लंबे समय से मौजूद सैद्धांतिक अस्पष्टता को समाप्त किया गया है
सुधार कारक सामान्यीकरण: पूर्व में केवल बेसेल विधि पर लागू सुधार कारकों (जैसे Gibson-Lanni, विलोपन कारक) को फूरियर विधि तक विस्तारित किया गया है
कस्टम FFT कार्यान्वयन: chirp Z रूपांतरण के आधार पर मनमानी पिक्सेल आकार के 2D FFT को लागू किया गया है, स्थानीयकरण सूक्ष्मदर्शी में अत्यंत छोटे पिक्सेल आकार के नमूनाकरण समस्या को हल किया गया है
कुशल संख्यात्मक एकीकरण: गोलीय विधि Simpson नियम का उपयोग करके 4-वें क्रम की सटीकता को लागू करती है, और torch.vmap के माध्यम से बैच प्रसंस्करण को वेक्टरीकृत करती है
अवकलनीय बेसेल फलन: PyTorch मूल रूप से समर्थित नहीं बेसेल फलन स्वचालित अवकलन कार्यक्षमता को विस्तारित किया गया है
सैद्धांतिक एकीकरण: Richards-Wolf इंटीग्रल से शुरू करते हुए, यह प्रमाणित किया गया है कि फूरियर (कार्टेशियन) और बेसेल (गोलीय) विधियां एक ही प्रसार इंटीग्रल के विभिन्न पैरामीटरकरण हैं, लंबे समय से मौजूद सैद्धांतिक अस्पष्टता को समाप्त किया गया है
Leutenegger et al. (2006): तेज़ फोकस क्षेत्र संगणना, कार्टेशियन विधि का आधुनिक कार्यान्वयन
Aguet (2009): सदिश PSF मॉडल की व्यवस्थित पेपर (डॉक्टरेट थीसिस)
Kirshner et al. (2013): लेखकों की पूर्व PSF Generator Java संस्करण
Miora et al. (2024): नवीनतम PSF संगणना विधि सारांश और MATLAB कार्यान्वयन
समग्र मूल्यांकन: यह सिद्धांत और अभ्यास को कसकर एकीकृत करने वाला एक उत्कृष्ट पेपर है। सैद्धांतिक रूप से, यह पहली बार दोनों मुख्य PSF मॉडलिंग विधियों की समानता को कठोरता से प्रमाणित करता है, एक एकीकृत ढांचा स्थापित करता है; व्यावहारिक रूप से, यह उच्च गुणवत्ता का PyTorch कार्यान्वयन और व्यवस्थित बेंचमार्किंग प्रदान करता है। पेपर की ओपन-सोर्स भावना और पारिस्थितिकी तंत्र एकीकरण (napari, chromatix) वैज्ञानिक समुदाय के लिए महत्वपूर्ण योगदान प्रदर्शित करते हैं। मुख्य कमी प्रायोगिक डेटा और मौजूदा सॉफ्टवेयर के साथ सीधी तुलना की कमी है। यह कार्य कम्प्यूटेशनल सूक्ष्मदर्शी क्षेत्र का एक महत्वपूर्ण उपकरण बन जाएगा, विशेष रूप से गहन शिक्षा और भौतिक मॉडल के संयोजन के अनुप्रयोग में व्यापक संभावनाएं हैं।