تمت الموافقة على هذه الدراسة من قبل لجنة أخلاقيات البحوث السريرية في مستشفى شينزين الشعبي (رقم الموافقة LL-KY-2025300-01). وقد أعفت لجنة الأخلاقيات من شرط الحصول على الموافقة المستنيرة لأن هذه الدراسة الاسترجاعية اعتمدت على صور نسيجية مرضية ومواد سريرية موجودة مسبقاً دون أي اتصال مباشر بالمرضى أو تدخل طبي. كما تم تجهيل كافة البيانات السريرية والصور النسيجية المرضية قبل التحليل، ولم يتم استخدام أي معلومات تحدد هوية الأشخاص في هذه الدراسة. وقد تم التعامل مع جميع البيانات المستخدمة في هذا البحث وفقاً للمعايير الأخلاقية للجان البحوث المؤسسية والوطنية.
نظرة عامة على إطار العمل
يظهر الهيكل العام لإطار عمل UAD المقترح، RareCode، في الشكل 1. يعتمد إطار العمل بنية متوازية ذات فرعين؛ فلكل صورة إدخال بحجم 512 × 512، يتم استخراج رقع غير متداخلة بحجم 32 × 32 و 64 × 64 كمدخلات دقيقة وواسعة النطاق على التوالي. يتم تغيير حجم الرقع ذات الـ 64 × 64 إلى 32 × 32 قبل تمريرها إلى الشبكة، بحيث يتشارك كلا الفرعين في نفس حجم إدخال المشفر-المفكك (encoder-decoder) مع الحفاظ على مجالات استقبال مختلفة. تقوم المشفرات في كل فرع برسم خرائط للميزات المستخرجة في فضاء كامن، حيث تفرض وحدة MHC قيود التجزئة. لاحقاً، تقوم المفككات بإعادة بناء الصورة من الميزات المكممة لحساب أخطاء إعادة البناء في المجال المكاني. بعد ذلك، تقيس آلية CAR درجة الشذوذ على المستوى الدلالي من خلال تحليل ترددات تنشيط دليل الرموز (codebook). وأخيراً، يدمج النموذج درجات إعادة البناء ودرجات CAR من خلال دمج موزون للحصول على درجات الشذوذ على مستوى الصورة. يتم تدريب إطار عمل RareCode بشكل شامل (end-to-end) باستخدام العينات الطبيعية فقط، دون الحاجة إلى تعليقات توضيحية للشذوذ. وقد استُخدم تصميم الفرعين لالتقاط كل من الميزات الخلوية المحلية والهياكل الغدية الأوسع، بينما استُخدمت وحدة MHC لنمذجة هذه الميزات باستخدام أدلة رموز بسعات مختلفة.
بنية المشفر-المفكك (Encoder-decoder architecture)
يقوم RareCode ببناء مشفرات ومفككات مستقلة هيكلياً للفروع ذات المقياس الدقيق والمقياس الخشن. يستخدم كل فرع نفس تصميم المشفر-المفكك ولكن دون مشاركة المعلمات. يتكون المشفر من أربع كتل تلافيفية، تشتمل كل منها على تلافيف 3 x 3، وتطبيع الدفعة (batch normalization)، وتنشيط ReLU، وإسقاط (dropout). يزداد عرض القناة من 64 إلى 128، ثم إلى 256، ويتم تسطيح خريطة الميزات النهائية وإسقاطها في تضمين كامن ذي 64 بُعداً. أما المفكك فيعكس تصميم المشفر بطبقة إسقاط متصلة بالكامل تليها أربع طبقات تلافيف منقولة (transposed convolution layers)، مما يعيد بناء كل رقعة إلى حجم الإدخال الأصلي للشبكة. يتم حساب خطأ إعادة البناء كمتوسط مربع الخطأ بين الرقع المدخلة والرقع المعاد بناؤها. ومن خلال الجمع بين بنية المشفر-المفكك هذه وقيود دفتر الرموز المنفصلة، يتعلم RareCode تمثيلات مدمجة للأنسجة الطبيعية ويقيس الانحرافات أثناء عملية الاستنتاج.
دفتر رموز هرمي متعدد المقاييس
لالتقاط السمات المرضية عند مستويات متعددة من التجريد - بدءًا من أنماط الأنسجة الواسعة وصولاً إلى التباينات الخلوية الموضعية - تستخدم وحدة MHC عدد K من دفاتر الرموز المنفصلة
بسعات متفاوتة n1، n2، ...، nK. في تكوين FourScales النهائي، يحتوي كل فرع على أربعة دفاتر رموز تتضمن 64 و 128 و 256 و 512 إدخالاً على التوالي، ولكل إدخال في دفتر الرموز تضمين (embedding) ذو 64 بُعداً. وبالنسبة لكل سمة رقعة مشفرة، يتم اختيار أقرب إدخال في دفتر الرموز وفقاً للمسافة الإقليدية. ثم يتم دمج المخرجات المكممة من دفاتر الرموز المختلفة باستخدام أوزان قابلة للتعلم يتم توحيدها عبر دفاتر الرموز. ومن المتوقع أن تقوم دفاتر الرموز الأصغر، بسبب قيود الضغط القوية، بتشفير أنماط نموذجية خشنة الحبيبات تجرد التباينات المحلية، بينما تحافظ دفاتر الرموز الأكبر على سمات أدق حبيبات تلتقط خصائص هيكلية أكثر تحديداً. لكل دفتر رموز
، حيث يمثل ek(i)
Rd متجه التضمين رقم i في دفتر الرموز رقم k، يتم تعيين السمات المستمرة إلى فضاء دفتر الرموز المنفصل عبر البحث عن أقرب جار (المعادلتين 1 و 2):


لتحقيق دمج تكيفي للميزات متعددة المقاييس، نقوم بإدخال أوزان قابلة للتعلم. وبعد إجراء تسوية softmax، يتم إجراء جمع موزون للمخرجات المكممة من كل كتاب رموز (المعادلة 3):

حيث تشير σ(·) إلى دالة softmax، مما يضمن استيفاء الأوزان للمعادلة Σk σ(wk) = 1. يتيح هذا التصميم للنموذج ضبط نسب المساهمة لدفاتر الرموز (codebooks) تلقائيًا عند مستويات دقة مختلفة بناءً على المحتوى الدلالي لميزات الإدخال.
لتحسين تعلم دفتر الرموز (codebook)، نعتمد دالة خسارة VQ القياسية (المعادلة 4):

حيث تشير sg[·] إلى عملية إيقاف التدرج (stop-gradient)، ويمثل β = 0.25 معامل وزن خسارة الالتزام. يحفز الحد الأول نواقل دفتر الرموز (codebook vectors) للتحرك نحو مخرجات المشفر، بينما يحفز الحد الثاني مخرجات المشفر على البقاء متسقة مع نواقل دفتر الرموز المقابلة لها. وتظهر الآلية الكاملة للكمية المتجهة متعددة المقاييس (multi-scale VQ) موصحة في الشكل 2.
تسجيل ندرة تنشيط دفتر الرموز
يقيس CAR الشذوذ على المستوى الدلالي بناءً على إحصائيات تنشيط دفتر الرموز المقدرة من عينات التدريب الطبيعية. بعد تدريب النموذج، تم تمرير جميع صور التدريب الطبيعية عبر المشفر ووحدة MHC دون استخدام تعزيز البيانات. وبالنسبة لكل فرع وكل دفتر رموز، تم تسجيل مؤشر التعيين لكل ميزة رقعة، وحساب عدد التعيينات لكل إدخال في دفتر الرموز. ثم تم تطبيع هذه الحسابات للحصول على توزيع تكرار التنشيط لدفتر الرموز هذا (المعادلة 5):

خلال مرحلة الاستدلال، تمت معالجة كل صورة اختبار عبر نفس إجراء استخراج الرقع، والترميز، وتعيين دفتر رموز الجار الأقرب. وبالنسبة لكل إدخال معين في دفتر الرموز، يتم حساب درجة الندرة على أنها اللوغاريتم السالب لتردد تنشيطه في مجموعة التدريب الطبيعية (المعادلة 6):

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

تُشتق هذه الدرجة من المسافة الإقليدية بين كل متجه سمات مُشفّر وأقرب إدخال له في دفتر الرموز. ويتم تقدير توزيع المسافة من عينات تدريب طبيعية، ثم تُحول مسافات عينات الاختبار إلى درجات مئينية. وتشير الدرجات المئينية الأكبر إلى أن السمة المُشفرة أكثر صعوبة في التمثيل باستخدام نماذج دفتر الرموز الطبيعية المُتعلمة.
تجمع درجة CAR النهائية بين ندرة التنشيط ومسافة التكميم عبر جميع كتب الشفرات باستخدام أوزان كتاب الشفرات القابلة للتعلم (المعادلة 8):

تُطبق نفس أوزان دليل الرموز الموحدة المستخدمة لدمج الميزات متعددة المقاييس على دمج مستوى الدرجات، وذلك للحفاظ على الاتساق بين تعلم التمثيلات وتسجيل الشذوذ. يتم حساب درجة CAR الناتجة على مستوى الصورة، ثم تُدمج مع درجة إعادة البناء للحصول على درجة الشذوذ النهائية. وتتضح آلية تسجيل CAR الشاملة في الشكل 3.
حساب درجة الشذوذ النهائية
نقوم بدمج أخطاء إعادة البناء في المجال المكاني مع درجات CAR على المستوى الدلالي. وقبل عملية الدمج، يتم تطبيع كلا الدرجتين بشكل منفصل باستخدام تطبيع الحد الأدنى والأقصى القائم على المئوية (percentile-based min-max normalization) لتقليل تأثير القيم المتطرفة الشديدة. ومن ثم تُعرَّف درجة الشذوذ النهائية على مستوى الصورة وفقاً لـ (المعادلة 9):

حيث يشير Srecon إلى متوسط خطأ إعادة البناء المربع المعير، وهو ما يعكس جودة إعادة بناء العينة في حيز البكسل؛ وSCAR هو درجة CAR الموصوفة سابقاً، والتي تقيس درجات الشذوذ على المستوى الدلالي. ويتحكم المعامل الفائق α
[0,1] في الوزن النسبي بين المكونين، مع تحديد القيم المثلى بناءً على مجموعات التحقق.
هدف التدريب
يتم تدريب إطار عمل RareCode من البداية إلى النهاية باستخدام العينات العادية فقط. وتتكون خسارة التدريب الإجمالية من المكونات التالية (المعادلة 10):

يتم تعريف كل مصطلح على النحو التالي: فقد إعادة البناء:
، وهو يقيس متوسط مربع الخطأ بين الرقع المدخلة p والرقع المعاد بناؤها D(zq) بعد الكمية. يتم تطبيق هذا الفقد بشكل منفصل على فرع الرقع بحجم 32 × 32 (
) وفرع الرقع بحجم 64 × 64 (
) لضمان التعلم الفعال للميزات على كلا المقياسين المكانيين.
وصف مجموعة البيانات
لتقييم RareCode، أُجريت تجارب على مجموعة بيانات صور نسيجية مرضية لسرطان القول والمستقيم (CRC) مُصنفة إكلينيكياً، وتضم صوراً لمقاطع أنسجة القول والمستقيم مصبوغة بالهيماتوكسيلين والإيوزين (H&؛هـ)، والتي تم الحصول عليها من مستشفى شينزين الشعبي (Shenzhen People's Hospital). وجميع الصور الواردة في مجموعة البيانات هذه مستمدة من حالات سريرية حقيقية. تم تحويل الصور إلى صيغة رقمية بتكبير 40x وبدقة أصلية تبلغ 1024 × ١٠٢٤ بكسل. ولضمان موثوقية وجودة التوصيف، تمت مراجعة جميع تسميات فئات العينات وتأكيدها بشكل مشترك من قبل اثنين أو أكثر من استشاريي علم الأمراض وفقاً لبروتوكول التعمية المزدوجة.
بناءً على الخصائص النسيجية المرضية، تم تصنيف مجموعة البيانات إلى ثلاث فئات: الأنسجة الطبيعية (1,226 صورة)، والسرطان (1,215 صورة)، والالتهاب (1,214 صورة). وفي سياق الكشف الثنائي عن الشذوذ (UAD)، عُدّت عينات السرطان والالتهاب كحالات شاذة، بينما استُخدمت العينات الطبيعية كتوزيع مرجعي. ويعكس هذا الصياغة مهمة الفرز المستهدفة المتمثلة في فصل الحالات التي تتطلب مراجعة مرضية إضافية عن تلك التي تعتبر طبيعية من الناحية المورفولوجية.
جميع H&تم تغيير حجم الصور المصبوغة بـ H&E إلى 512 × 512 بكسل قبل استخدامها كمدخلات للشبكة. واستُخدم التحقق المتقاطع الطبقي خماسي الطيات (Stratified 5-fold cross-validation) لأغراض التقييم. وفي كل طية، قُسِّمت العينات الطبيعية أولاً إلى مجموعات فرعية للتدريب، والتحقق، والاختبار المستبعد. تدرَّب RareCode على عينات التدريب الطبيعية فقط. أما المجموعة الفرعية للتحقق، والتي احتوت أيضاً على عينات طبيعية فقط، فقد استُخدمت لاختيار النموذج، وضبط المعلمات الفائقة (hyperparameter tuning)، واختيار وزن الدمج. αاستُخدمت مجموعة اختبار مستبعدة، تتكون من عينات طبيعية غير مرئية وعينات شاذة تشمل السرطان والالتهابات، لغرض التقييم النهائي للأداء فقط، ولم تُستخدم في تدريب النموذج، أو اختيار المعلمات الفائقة، أو α الاختيار. تظهر الفئات التمثيلية، وبنية مجموعة البيانات، وبروتوكول التحقق في الشكل 4.
تفاصيل التنفيذ
تم تنفيذ إطار عمل RareCode المقترح ونماذج المقارنة المرجعية باستخدام PyTorch. وأُجريت جميع التجارب على محطة عمل مجهزة بوحدة معالجة رسومات واحدة من نوع NVIDIA GeForce RTX 4060 Ti (16 GB). أما بالنسبة لتكوين بنية الشبكة، ومن أجل التقاط السمات بمدى مكاني مختلف، تقوم البنية ثنائية الفرع بمعالجة رقع غير متداخلة بمقياسين: 32 × 32 بكسل (مجال استقبال أصغر، لالتقاط أنماط النسيج المحلية) و 64 × 64 بكسل (مجال استقبال أكبر، لالتقاط سياق مكاني أوسع)، على التوالي. وقد تم تعيين بُعد تضمين السمات المستمرة الناتج عن مشفرات كلا الفرعين بشكل موحد عند 64. وفي وحدة MHC، قمنا بتكوين أربعة دفاتر رموز بسعات مختلفة لكل فرع، بكميات نماذج أولية منفصلة تبلغ 64 و 128 و 256 و 512، على التوالي.
خلال عملية التحسين، تم تدريب الشبكة من الطرف إلى الطرف لمدة 50 حقبة (epochs) وبحجم دفعة قدره 8. وقد استخدمنا مُحسن AdamW لتحديثات الأوزان، بمعدل تعلم أولي قدره 5 × واضمحلال للوزن من أجل منع الإفراط في التجهيز (overfitting). بالإضافة إلى ذلك، ولضمان تقارب سلس للتدريب، استُخدم جدول زمني لمعدل التعلم يعتمد على التلدين الجيبي (cosine annealing)، مما أتاح تحسيناً أدق في المراحل المتأخرة من التدريب. كما تم اختيار أحجام الرقع (patch sizes) وأحجام دفتر الرموز (codebook sizes) لتحقيق التوازن بين التمثيل متعدد المقاييس والتكلفة الحسابية. وقد تم دعم التكوين النهائي لـ FourScales من خلال تحليل الاستئصال (ablation analysis)، كما تم اختيار α بناءً على مجموعة التحقق قبل التقييم باستخدام مجموعة الاختبار.
بالنسبة لاستئصال تعقيد وحدة فك التشفير (decoder-complexity ablation)، تم تنفيذ متغير وحدة فك تشفير معقدة استناداً إلى تكوين FourScales. استخدم هذا المتغير نفس تقسيم البيانات، وأحجام دفتر الرموز (codebook sizes)، وعصور التدريب، والمحسن، ومعدل التعلم، واختيار α بناءً على التحقق، ومقاييس التقييم المستخدمة في نموذج FourScales. تم استبدال وحدة فك التشفير الأساسية بوحدة فك تشفير من نمط EMCAD تحتوي على كتل تلافيفية بعمق متعدد المقاييس (multi-scale depthwise convolution blocks) ووحدات انتباه الكتل التلافيفية (CBAM). بعد كل مرحلة من مراحل زيادة العينات بالتلافيف المنقولة (transposed-convolution upsampling) الثلاث الأولى، تم تطبيق تلافيف بعمق بمقاسات 3 × 3 و 5 × 5 و 7 × 7 على التوازي، تلاها دمج نقطي (pointwise fusion) بمقاس 1 × 1، وتطبيع الدفعة (batch normalization)، وتنشيط ReLU، واتصال متبقٍ (residual connection)، وانتباه CBAM. تضمن CBAM انتباهاً للقنوات بناءً على واصفات التجميع المتوسطي والأقصى (average- and max-pooling descriptors)، وانتباهاً مكانياً باستخدام تلافيف بمقاس 7 × 7. واستخدمت المرحلة النهائية التلافيف المنقولة وتنشيط sigmoid لإعادة بناء رقع (patches) بمقاس 32 × 32.