Optimal mitigation of random telegraph noise for improved photometry at high frame rates

ملخص تقني: التخفيف الأمثل لضجيج التردد العشوائي (RTN) لتحسين القياس الضوئي عند معدلات الإطارات العالية

بيان المشكلة
يُعد ضجيج التردد العشوائي (RTN)، المعروف أيضاً باسم ضجيج "الملح والفلفل" أو ضجيج "الفشار"، مصدراً كبيراً لضجيج القراءة في مستشعرات الصور ذات أشباه الموصلات المكملة (CMOS) المستخدمة في علم الفلك. وبينما توفر مستشعرات CMOS العلمية (sCMOS) بشكل عام ضجيج قراءة أقل من أجهزة (CCD) التقليدية، فإن وجود الـ RTN يقلل من هذه الميزة، لا سيما في التطبيقات ذات معدل الإطارات العالي والتدفق المنخفض. ينشأ الـ RTN من التقاط وتحرير الإلكترونات بواسطة مصائد الحوامل بالقرب من الواجهة بين (Si/SiO2) في ترانزستورات تتبع المصدر. في البكسلات التي تعاني من الـ RTN، يتذبذب مستوى الانحياز بين حالات منفصلة، مما يخلق ذيلاً طويلاً في توزيع ضجيج القراءة. يؤدي هذا إلى زيادة جذر متوسط المربع (RMS) لضجيج القراءة بشكل كبير، حتى لو ظل وسيط ضجيج القراءة منخفضاً، مما يقلل من نسبة الإشارة إلى الضجيج (SNR) للمصادر الخافتة. تفشل استراتيجيات التخفيف الموجودة، مثل حجب البكسلات المعيبة أو تطبيق مرشحات الوسيط البسيطة، في كثير من الأحيان في مراعاة توزيعات مستويات الانحياز المحددة لكل بكسل، وقد تؤدي إلى تدهور دقة القياس الضوئي عندما تكون دالة انتشار النقطة (PSF) ناقصة العينات، أو قد تسبب انحيازات منهجية.

المنهجية
استقصى المؤلفون ثلاثة مستشعرات CMOS متوفرة تجارياً تُستخدم في الأدوات الفلكية: Sony IMX455 (في كاميرات QHY600M)، وGpixel GSENSE400 (في كاميرات QHY42)، وFairchild Imaging HWK4123 (في كاميرات Hamamatsu ORCA-Quest 2).

1. التوصيف: جمع الفريق آلاف إطارات الانحياز (bias frames) لقياس ضجيج القراءة وتحديد بكسلات الـ RTN. استخدموا اختبار "أندرسون-دارلينج" للنمطية للكشف عن توزيعات الانحياز غير الغاوسية. بالنسبة للبكسلات التي فشلت في هذا الاختبار، قاموا بملاءمة نموذج غاوسي ثلاثي لتوصيف توزيع مستوى الانحياز، واستخراج المعلمات للغاوسي المركزي (المتوسط $\mu_b$، والانحراف المعياري $\sigma_r$) والقمم الجانبية التي تمثل قفزات الـ RTN (السعات $B_1, B_2$، ومسافة القفزة $d$).
2. خوارزمية التصحيح: تم تطوير خوارزمية جديدة تعتمد على لغة بايثون لتحديد وتصحيح قفزات الـ RTN لكل بكسل وكل إطار على حدة.
   • توليد المرجع: تقدر الخوارزمية مستوى الإشارة المتوقع للبكسل باستخدام الوسيط الزمني (للبيانات السلاسل الزمنية) أو الوسيط المكاني (للإطارات الفردية)، مع تعديله وفقاً لمعلمات الانحياز الخاصة بالبكسل.
   • تحديد القفزة: باستخدام دالة التوزيع التراكمي (CDF) لـ "بواسون-غاوس"، تحسب الخوارزمية عتبات الأهمية ($\alpha = 0.003$) لتحديد ما إذا كانت إشارة البكسل تنحرف بشكل كبير عن التوزيع الغاوسي المركزي المتوقع، مع مراعاة الانحراف نحو اليمين المتأصل في إحصاءات بواسون.
   • التصحيح: إذا تم تحديد قفزة، يتم إزاحة الإشارة بمقدار القفزة $d$ (محولة إلى إلكترونات). تتضمن الخوارزمية خطوة تصحيح الانحياز لمنع التضخم الاصطناعي لقيم التدفق المتوسط الناتج عن التصحيح غير المتماثل للقفزات الموجبة مقابل السالبة.
3. التحقق من الصحة: قارن المؤلفون خوارزميتهم باستراتيجية "الحجب" البسيطة (معاملة بكسلات الـ RTN كبكسلات معيبة) باستخدام محاكاة مونت كارلو والملاحظات الفلكية الحقيقية. شمل الاختبار الواقعي القياس الضوئي للمجال النجمي وحدث احتجاب كويكب تم رصده عبر شبكة تلسكوبات "لاس كومبريس" (LCO) الصغيرة.

المساهمات الرئيسية

• تكميم تأثير الـ RTN: تحدد الدراسة مساهمة الـ RTN في أرضية ضجيج القراءة الإجمالية. بالنسبة لـ Sony IMX455، يزيد الـ RTN من جذر متوسط مربع (RMS) ضجيج القراءة بنسبة $\gtrsim 20\%$. بالنسبة لـ GSENSE400 وHWK4123، تكون الزيادات أصغر ($\gtrsim 4\%$ و$\gtrsim 2.5\%$ على التوالي)، وإن كانت لا تزال ذات صلة بالتطبيقات ذات الضجيج الفائق الانخفاض.
• تطوير الخوارزمية: يقدم المؤلفون خوارزمية تصحيح قوية لكل بكسل، تتكيف مع توزيع مستوى الانحياز المحدد لكل بكسل RTN وتدفق الفوتونات المحلي، متجنبة عيوب الطرق ذات العتبة الثابتة.
• التحليل المقارن: تقدم الورقة مقارنة مفصلة بين خوارزمية التصحيح المقترحة وحجب البكسلات، وتوضح أن الخوارزمية تتفوق على الحجب، خاصة عندما تكون الـ PSF ناقصة العينات، أو عندما تكون البكسلات المحجوبة قريبة من مركز المصدر، أو عندما يكون ضجيج القراءة مقارباً لضجيج الفوتونات (shot noise).
• البرمجيات مفتوحة المصدر: أطلق المؤلفون مجموعة برمجيات RTN-Corrector التي تتضمن سكربتات لتحديد بكسلات الـ RTN (rtn_fitter) وتطبيق التصحيحات (rtn_fixer)، مما يجعل هذه التقنيات متاحة للمجتمع الفلكي.

النتائج

• المحاكاة: أظهرت محاكاة مونت كارلو أن خوارزمية التصحيح حسنت نسبة الإشارة إلى الضجيب (SNR) باستمرار مقارنة بالحجب. وُجد أن الحجب يؤدي إلى تدهور الـ SNR في الحالات التي تقع فيها بكسلات الـ RTN بالقرب من قلب المصدر، خاصة في حالات الـ PSF ناقصة العينات.
• القياس الضوئي النجمي: في ملاحظات مجال نجمي باستخدام IMX455، أدى كل من خوارزمية التصحيح والحجب إلى تحسين SNR للمخططات الضوئية بمتوسط $>5\%$ للمصادر الخافتة ($< 3 e^-$/pix/frame). أظهرت خوارزمية التصحيح تبايناً أقل في الأداء وحافظت على دقة قياس ضوئي مطلقة، بينما مال الحجب إلى التقليل من تقدير السطوع.
• احتجاب الكويكب: خلال احتجاب الكويكب الطروجان "أنشيس" (Anchises)، قللت خوارزمية التصحيح من بقايا جذر متوسط المربع (RMS residuals) للمخطط الضوئي بنسبة $\approx 7\%$ وصححت المبالغة المنهجية في تقدير عمق الاحتجاب الناتجة عن عدم توازن قفزات الـ RTN السالبة أثناء الحدث.
• خصوصيات المستشعر: أشارت الدراسة إلى أن IMX455 أظهر توزيعات RTN متماثلة (قفزات موجبة وسالبة)، بينما أظهر GSENSE400 وHWK4123 توزيعات غير متماثلة للغاية (قفزات موجبة في الغالب) ومقادير قفزات أصغر. أدى تصميم HWK4123 ذو السعة الفائقة الانخفاض إلى جعل قفزات الـ RTN أصغر من فوتون واحد في كثير من الحالات، مما حد من فعالية تصحيح ما بعد المعالجة لهذا المستشعر تحديداً.

الأهمية والادعاءات
تزعم الورقة أن التخفيف من الـ RTN ضروري لتحقيق الإمكانات الكاملة لمستشعرات sCMOS، لا سيما للقياس الضوئي عالي السرعة والملاحظات للأجسام الخافتة حيث يهيمن ضجيج القراءة. ويؤكد المؤلفون أن خوارزميتهم توفر طريقة موثوقة لاستعادة دقة القياس الضوئي التي قد تضيع بسبب الـ RTN، مما يقدم بديلاً متفوقاً للحجب البسيط للبكسلات. ويشددون على أنه بينما يكون الحجب فعالاً في أنظمة محددة (PSF جيدة العينات، وضجيج القراءة مهيمن)، فإن خوارزمية التصحيح توفر تحسينات أكثر اتساقاً عبر مجموعة أوسع من الظروف، بما في ذلك الـ PSF ناقصة العينات وتغير سطوع المصدر. ويُقدم العمل كأداة عملية لعلماء الفلك الذين يستخدمون كواشف تعاني من الـ RTN، مما يوفر دفعة فورية لجودة البيانات دون الحاجة إلى تغييرات في الأجهزة. كما يشير المؤلفون إلى أن تخفيف الـ RTN هو أمر بالغ الأهمية بشكل خاص للتلسكوبات الفضائية، حيث يُعرف أن الأضرار الإشعاعية تزيد من انتشار الـ RTN.