Unsupervised Satellite Image Change Detection Using Quantum-Inspired Mapping and Weighted Hybrid Neighborhood MRFs
تقترح هذه الورقة إطار عمل غير مُشرف للكشف عن التغير في صور الأقمار الاصطناعية يعيد صياغة حقل ماركوف العشوائي المدرك للسياق مع جيران متصلة بثماني اتجاهات وموزونة كمسألة QUBO، محققاً أداءً تنافسياً على معيار OSCD مع الحفاظ على التوافق مع بنيات التلدين الكمي المستقبلية.
تلتقط الأقمار الصناعية التي تدور حول الأرض صوراً للأماكن نفسها باستمرار في أوقات مختلفة، مما يخلق سجلاً مرئياً لكيفية تغير كوكبنا. وتعد هذه الصور حيوية لمراقبة نمو المدن، أو تتبع إزالة الغابات، أو تقييم الأضرار بعد وقوع الكوارث. ومع ذلك، فإن التمييز بين التغيير الحقيقي وبين مجرد خدعة بصرية ناتجة عن الضوء أمر صعب؛ فظلال غيمة ما، أو تغير في زاوية الشمس، أو اختلاف طفيف في مستشعر الكاميرا، يمكن أن يجعل شارعاً هادئاً يبدو وكأنه موقع بناء. ولعقود من الزمن، حاول العلماء حل هذه المشكلة عبر تعليم الحواسيب كيفية رصد الاختلافات، لكن العديد من الطرق تتطلب كميات هائلة من الأمثلة المصنفة لتعلم شكل التغيير، أو أنها تقوم بتنعيم الصورة لدرجة تمحو معها التفاصيل ذاتها التي تحاول البحث عنها.
يقدم نهج جديد، طوره باحثون في جامعة أمير كبير التقنية وجامعة خوارزمي في إيران، مساراً مختلفاً. فبدلاً من الاعتماد على بيانات التدريب أو القواعد الجامدة، ابتكروا نظاماً ينظر إلى النسيج المحلي للصورة ليقرر ما إذا كان التغيير حقيقياً. وقد تعامل الفريق مع المشكلة كأحجية للبحث عن الترتيب الأكثر استقراراً للبكسلات، وهي طريقة متجذرة في مفهوم رياضي يسمى "حقل ماركوف العشوائي" (Markov Random Field). وفي هذا الإطار، لا يتم اتخاذ القرار لكل بكسل بمعزل عن غيره، بل يتأثر بجيرانه، تماماً كما قد يتأثر رأي شخص ما بآراء الأشخاص الواقفين بالقرب منه مباشرة. وقد ترجم الباحثون هذا المنطق إلى صيغة يمكن حلها بواسطة حواسيب مستوحاة من الكم، رغم أنهم اختبروها باستخدام محاكاة كلاسيكية. وكان هدفهم هو بناء نظام يتكيف مع المشهد المحدد الذي يراه، فيقوم بتقوية قواعده في المناطق الصاخبة (التي تحتوي على ضجيج) وتخفيفها في الأماكن الواضحة، وكل ذلك دون الحاجة إلى إنسان ليعلمه ما الذي يجب أن يبحث عنه.
اختبر الباحثون طريقتهم على مجموعة قياسية من صور الأقمار الصناعية تغطي تسع مدن مختلفة، تتراوح من الشبكات الحضرية الكثيفة في هونغ كونغ إلى المناظر الطبيعية المتنوعة في بوردو. وقاموا بتقسيم الصور الكبيرة إلى بلاطات صغيرة سهلة الإدارة وحللوا الاختلافات بين لقطتين لكل موقع. وكان الابتكار الرئيسي في كيفية تعاملهم مع الروابط بين البكسلات؛ فالطرق التقليدية غالباً ما تعامل الجيران من اليسار واليمين والأعلى والأسفل بنفس الطريقة التي تعامل بها الجيران في الزوايا القطرية. أما النظام الجديد فقد أدرك أن الجيران القطريين أبعد جسدياً، وبالتالي هم أقل عرضة لأن يكونوا جزءاً من نفس الميزة المباشرة. ومن خلال تقليل وزن هذه الروابط القطرية، حافظ النظام على حدة الحواف والزوايا التي تميل الطرق الأخرى إلى طمسها. علاوة على على ذلك، قام النظام بحساب قواعده الخاصة لكل بلاطة بناءً على التباين المحلي في الصورة. فإذا كانت بقعة من الصورة فوضوية ومليئة بالضجيج، طبق النظام قاعدة تنعيم أقوى لتصفية التشويش، وإذا كانت البقعة منتظمة، طبق لمسة أخف لتجنب طمس التغييرات الحقيقية.
وعندما قارن الفريق نتائجهم بالطرق الراسخة، أظهر الإطار الجديد ميزة واضحة في تحقيق التوازن في الدقة. وفي الاختبارات، حقق النظام دقة إجمالية بلغت 96.15 بالمائة، مما يعني أنه حدد بدقة الحالة الغالبية العظمى من البكسلات. والأهم من ذلك، أنه وصل إلى درجة (F1-score) بلغت 0.3483، وهي مقياس يوازن بين القدرة على إيجاد التغييرات الحقيقية وبين خطر إطلاق الإنذارات الكاذبة. وكانت هذه الدرجة أعلى من العديد من التقنيات التقليدية التي تعتمد على عتبات إحصائية بسيطة، والتي غالباً ما تعاني للتمييز بين الضجيج والتحول الفعلي. ووجد الباحثون أن طريقتهم كانت جيدة بشكل خاص في كبح ضجيج "الملح والفلفل" (salt-and-pepper noise) الذي يعيب صور الأقمار الصناعية، حيث يتم تصنيف بكسلات معزولة بشكل خاطئ على أنها تغيرت. ومع ذلك، لاحظوا أيضاً أن النظام كان متحفظاً؛ فهو أقل عرضة لتفويت التغيير، لكنه أيضاً فاتته بعض التغييرات الصغيرة والمجزأة التي رصدتها الطرق الأخرى، مفضلاً اليقين على الشمولية.
تؤكد الدراسة أن الجمع بين السياق المحلي وهيكل الجوار الموزون يمكن أن يحسن عملية الكشف عن التغيير غير الخاضعة للإشراف دون الحاجة إلى بيانات تدريب مصنفة. وقد أثبت الباحثون أن نهجهم يعمل بشكل جيد عبر بيئات متنوعة، من شوارع نانت المنظمة إلى تضاريس رين المعقدة، رغم أن الأداء تباين اعتماداً على الخصائص المحددة لكل مدينة. وبينما صُممت الطريقة لتكون متوافقة مع أجهزة الحوسبة الكمومية المستقبلية، فقد تم إنتاج النتيات الحالية باستخدام محاكاة كلاسيكية، مما يظهر أن الصيغة الرياضية نفسها سليمة. ويشير العمل إلى أنه من خلال ترك الصورة تملي قواعدها الخاصة للنعومة واحترام الهندسة الفيزيائية للشبكة، يمكن للحواسيب أن تصبح مراقباً أكثر موثوقية لعالمنا المتغير. ويخطط المؤلفون لتطوير النظام بشكل أكبر من خلال معالجة الحدود بين بلاطات الصور واختبار الطريقة على آلات كمومية فيزيائية، لكن النتيجة الجوهرية تظل قائمة: إن وجود طريقة أذكى وتكيفية للنظر إلى الأرض من الأعلى هو أمر ممكن دون الحاجة إلى معلم.
ملخص تقني: كشف التغير في الصور الساتلية غير الخاضع للإشراف باستخدام خرائط مستوحاة من الكم ومراتب حقول ماركوف العشوائية الهجينة الموزونة
بيان المشكلة يعد كشف التغير (CD) غير الخاضع للإشراف في صور الاستشعار عن بعد ثنائية الزمن أمراً بالغ الأهمية لمراقبة الأرض، ولكنه يظل تحدياً بسبب الضوضاء الإشعاعية، والتغيرات الجوية، والتسجيل غير المثالي للصور. تعتمد النهج التقليدية غالباً على تنظيم ثابت مكانياً، مما قد يؤدي إلى حساسية مفرطة للضوضاء أو تسبب في تنعيم مفرط للحدود في المشاهد غير المتجانسة. علاوة على ذلك، فإن طرق حقل ماركوف العشوائي (MRF) الموجودة غالباً ما تستخدم معاملات تنظيم ثابتة عالمياً وجوارات مكونة من 4 اتصالات أو 8 اتصالات متناحية. تفشل هذه التكوينات الثابتة في التكيف مع إحصائيات المشهد المحلية وتتجاهل التمييز الهندسي بين المسافات العمودية والقطرية، مما يؤدي إلى مقايضات دون المستوى الأمثل بين كبح الإنذارات الكاذبة والحفاظ على الحدود الهيكلية الدقيقة.
المنهجية يقترح المؤلفون إطار عمل غير خاضع للإشراف ولا يتطلب تدريباً، يعيد صياغة حقل ماركوف العشوائي (MRF) المفكك إلى رقع (patches) كمسألة تحسين ثنائي تربيعي غير مقيد (QUBO). تتكون المنهجية من ثلاث مراحل رئيسية:
إعداد البيانات وتفكيكها:
يتم تحويل الصور متعددة الأطياف ثنائية الزمن إلى صور رمادية أحادية النطاق باستخدام أوزان السطوع ITU-R BT.601.
يتم تطبيق مطابقة الهيستوجرام (Histogram matching) على الصورة الثانية بالنسبة للأولى لتحييد الاختلافات الإشعاعية.
تُقسم الصور إلى بلاطات (tiles) غير متداخلة بحجم 32×32 بكسل مع استخدام حشوة انعكاس الحواف المتماثلة لتخفيف آثار الحدود.
يتم إنشاء صورة فرق إشعاعي (ΔI) عبر طرح مطلق لكل بكسل.
صياغة طاقة MRF ورسم خرائط QUBO:
يتم تصنيف مهمة كشف التغير كتقدير الاحتمال الأقصى البعدي (MAP)، من خلال تقليل دالة الطاقة العالمية E(x)=Edata(x)+Espatial(x).
دقة البيانات (Edata): يتم نمذجتها خطياً بناءً على مقدار الفرق بالنسبة لعتبة أساسية ثابتة (λ=0.10).
بيتا التكيفية الواعية بالسياق (βadaptive): بدلاً من معامل تنعيم عالمي، يتم حساب قوة التنعيم المكاني ديناميكياً لكل رقعة. وهي مشتقة من متوسط الرقعة المحلية (μpatch) والانحراف المعياري (σpatch) لصورة الفرق، ومدرجة مقابل عتبة الأساس. تضمن قيمة أرضية قدرها 0.001 الاستقرار العددي.
الجوار الهجين الموزون ذو الـ 8 اتصالات: يتم توسيع طوبولوجيا التفاعل المكاني من 4 اتصالات إلى 8 اتصالات. ومن الأهمية بمكان أن التفاعلات القطرية يتم قياسها بعامل مسافة إقليدية عكسي (1/2≈0.7071) لمراعاة الفصل الفيزيائي الأكبر للجيران القطريين مقارنة بالجيران العموديين.
تحويل QUBO: يتم رسم خريطة دالة طاقة MRF إلى صيغة QUBO (EQUBO(x)=∑Qssxs+∑Qsrxsxr). تدمج الانحيازات الخطية (Qss) دقة البيانات وعقوبات المساحة التكيفية، بينما تمثل الارتباطات التربيعية (Qsr) تفاعلات الجوار الموزونة.
التحسين وإعادة البناء:
على الرغم من أن الصياغة مصممة للتقليب الكمي (Quantum Annealing)، إلا أن الدراسة تستخدم حلاً كلاسيكياً يعتمد على التخمير المحاكي (Simulated Annealing) (مكتبة D-Wave neal) لإيجاد حالة الطاقة الدنيا.
يتم إعادة تجميع ناقلات الحل الثنائية من كل بلاطة لإعادة بناء خريطة التغيير العالمية.
المساهمات الرئيسية ليست المساهمة الأساسية للورقة هي المفاهيم الفردية لـ MRFs أو QUBO أو التخمير، بل هي تكاملها مع آليتين مكانيتين مبتكرتين:
التنظيم الواعي بالسياق: طريقة لتحديد قوة التنعيم المكاني مباشرة من إحصائيات الرقعة المحلية، مما يلغي الحاجة إلى الضبط التجريبي العالمي.
الطوبولوجيا الموزونة بالمسافة غير المتناحية: نموذج جوار ذو 8 اتصالات حيث يتم خفض وزن التفاعلات القطرية صراحةً بناءً على المسافة الإقليدية، بهدف الحفاظ على الدقة الهندسية بشكل أفضل من النماذج المتناحية.
النتائج تم تقييم الإطار باستخدام مجموعة بيانات Onera Satellite Change Detection (OSCD) عبر تسعة مشاهد حضرية متنوعة دون تدريب تحت الإشراف.
الأداء: حققت الطريقة المقترحة متوسط درجة F1 يبلغ 0.3483، ودقة إجمالية (OA) تبلغ 96.15%، ومعامل كابا كوهنز (Cohen's Kappa) يبلغ 0.3323.
التحليل المقارن:
تفوقت الطريقة المقترحة على طرق العتبة التقليدية غير الخاضعة للإشراف (Otsu, Kapur, Kittler-Illingrove, Rosin) من حيث درجة F1.
كشف تحليل حساسية المعلمات الفائقة أن معلمة β العالمية الثابتة تؤدي إلى مقايضة حيث تأتي الدقة العالية على حساب الاستدعاء (تآكل الحدود)، مما يبرر الحاجة إلى النهج التكيفي.
أظهرت المقارنات بين الطوبولوجيات أن تكوين الـ 8 جيران الهجين الموزون وفر توازناً أفضل للمقاييس من متغيرات الـ 4 جيران أو الـ 8 جيران المتناحية، لا سيما في الحفاظ على التفاصيل الهيكلية في المشاهد المعقدة.
القيود الملحوظة: أظهرت الطريقة استدعاءً (Recall) أقل مقارنة ببعض نماذج العتبة البسيطة، مما يشير إلى خريطة تغيير متحفظة تعطي الأولوية لتقليل الإيجابيات الكاذبة على حساب اكتشاف جميع التغييرات المحتملة. تباين الأداء بشكل كبير عبر المدن المختلفة (مثل F1 مرتفع في Rennes، ومنخفض في Bordeaux)، مما يعكس الصعوبة المتأصلة في مجموعة البيانات.
الأهمية والادعاءات يدعي المؤلفون أن هذا الإطار يقدم بديلاً غير خاضع للإشراف ولا يتطلب تدريباً لكشف التغير، ويحتفظ بالتوافق مع أجهزة التقليب الكمي المستقبلية. تكمن الأهمية في إثبات أن الجمع بين التنظيم الواعي بالسياق المحلي والتفاعلات المكانية الموزونة بالمسافة يمكن أن يحسن التوازن بين قمع الضوضاء والحفاظ على الحدود دون الحاجة إلى بيانات مصنفة أو تدريب التعلم العميق.
تصرح الورقة صراحةً بالتواضع بشأن ادعاءاتها:
هي لا تدعي تسارعاً كمياً، حيث أُجريت التجارب باستخدام التخمير المحاكي الكلاسيكي.
تقر بأن التفكيك القائم على الرقع قد يؤدي إلى عدم استمرارية عند حدود البلاطات.
تضع الطريقة كنهج مكمل للسيناريوهات التي يكون فيها التحسين المصاغ صراحةً وغير الخاضع للتدريب مرغوباً، بدلاً من كونها بديلاً للأنظمة القائمة على البيانات والتعلم العميق.
تم تحديد العمل المستقبلي لتقييم الصياغة على أجهزة كمية فيزيائية، ولمعالجة عدم استمرارية الحدود من خلال الرقع المتداخلة.