$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
استخدمت هذه الدراسة حصريا مجموعات بيانات التصوير المقطعي المحوسب والتصوير بالرنين المغناطيسي المتاحة للجمهور والتي تم تحديد هويتها. لم يكن هناك أشخاص أحياء من البشر أو. لذلك ، لم تكن هناك حاجة إلى موافقة مجلس المراجعة المؤسسية (IRB) أو لجنة الأخلاقيات.
نظرة عامة على الطريقة
يقدم هذا البروتوكول خط أنابيب قابلا للتكرار لتقليل تقليل الصور الطبية الموفرة للطاقة. فهو يجمع بين تقنيات المعالجة المسبقة ، بما في ذلك مرشحات الشحذ وتجميع K-means ، مع التشفير الذاتي المستند إلى الشبكة العصبية التلافيفية (CNN) لإزالة الضوضاء من الصور. تعمل هذه الطريقة المتكاملة على تحسين جودة الصورة مع تقليل وقت التدريب واستهلاك الطاقة للأجهزة ، ودعم التشخيص الطبي المستدام19،20،21،22،23. يلخص الشكل 5 الإطار الشامل.

الشكل 5: بنية إطار عمل تقليل الضوضاء المقترحة. يوضح هذا الشكل خط الأنابيب الكامل: الحصول على البيانات ، والمعالجة المسبقة (الشحذ + تجزئة الوسائل K) ، متبوعا بتقليل التشفير الذاتي المستند إلى CNN. ويؤكد على النهج الهجين الذي يهدف إلى تقليل استخدام الطاقة مع الحفاظ على الدقة الهيكلية. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.
يتم تلخيص سير العمل الكامل في الشكل 5 ، والذي يوضح التسلسل من إعداد مجموعة البيانات إلى المعالجة المسبقة (الشحذ و K-means) إلى تقليل الضوضاء من المشفر الذاتي ل CNN إلى الاختبار.
إعداد البرامج والبيئة
الخيار أ - Google Colab (موصى به لقابلية التكرار): انتقل إلى colab.research.google.com، انقر فوق دفتر ملاحظات جديد. حدد وقت التشغيل > تغيير نوع وقت التشغيل > مسرع الأجهزة: وحدة معالجة الرسومات > حفظ. انقر على ملف > تحميل لتحميل دفتر الملاحظات المقدم ومجموعة البيانات .zip (أو اربط Google Drive بالنقر على الملفات > تحميل Drive). في خلية التعليمات البرمجية الأولى ، قم بتثبيت التبعيات (بالفعل في دفتر الملاحظات): pip install opencv-python scikit-image scikit-learn tensorflow matplotlib numpy pandas. انقر فوق وقت التشغيل > تشغيل الكل. تأكد من رؤية: (i) ملخص البيئة مع إصدارات الحزمة، (ii) اسم GPU (في إخراج الخلية)، و (iii) إنشاء مجلد تجربة يعمل/YYYY-MM-DD/.
الخيار ب - محلي (كوندا): افتح موجه Anaconda / Terminal وقم بتشغيل:

ضع مجموعة البيانات الخاصة بك ضمن البيانات / raw / والبرامج النصية ضمن التعليمات البرمجية /. تشغيل: python code / train_autoencoder.py --data_root data --img_size 256 --k_values 3 5 --epochs 100 --batch 32 --lr 0.001. قم بتأكيد لرؤية: نسخة مطبوعة للبيئة، ووحدة معالجة الرسومات المكتشفة، وتشغيل دليل السجل/YYYY-MM-DD/.
إعداد مجموعة البيانات
تم الحصول على مجموعات البيانات المتاحة للجمهور من فحوصات التصوير المقطعي المحوسب والتصوير بالرنين المغناطيسي من مستودعات التصوير الطبي. احتوت مجموعات البيانات هذه على ضوضاء نموذجية لعمليات المسح السريرية في العالم الحقيقي وتم إلغاء تحديدها وفقا للمعايير المؤسسيةوالأخلاقية 8،9. تم تغيير حجم كل صورة إلى 256 × 256 بكسل وتخزينها بتنسيق PNG أو DICOM للتوافق. تم تقسيم مجموعات البيانات بشكل عشوائي على النحو التالي: 70٪ للتدريب ، و 15٪ للتحقق من الصحة ، و 15٪ للاختبار ، مما يضمن التوزيع الطبقي لطرق التصوير ومستويات الضوضاء⁶. تم قياس إحصائيات الضوضاء باستخدام متوسط تباين البكسل قبل تقليل الضوضاء.
تنظيم المجلدات: قم بإنشاء مجلد كما هو موضح أدناه:

لإضافة صور ، انسخ صور التصوير المقطعي المحوسب / التصوير بالرنين المغناطيسي غير المحددة إلى البيانات / raw / CT والبيانات / raw / التصوير بالرنين المغناطيسي (PNG أو JPG أو DICOM). توحيد حجم الصورة باستخدام Python (موصى به): قم بتشغيل خلية دفتر الملاحظات تغيير الحجم والتحويل. يقوم بتحميل كل صورة ، ويحول التدرج الرمادي إذا لزم الأمر ، ويغير حجمه إلى 256 × 256 .

البديل هو استخدام واجهة المستخدم الرسومية (بديل ImageJ / Fiji) كما هو موضح. انقر فوق ملف > استيراد تسلسل الصور > (أو صور مفردة) ، ثم صورة > اكتب > 8 بت ، صورة > ضبط حجم > ... > 256 × 256 ، ثم قم بتقديم > حفظ باسم > PNG في البيانات / preproc /.
إنشاء تقسيمات: قم بتشغيل خلية دفتر الملاحظات Train/Val/Test split (70/15/15)؛ يقوم بخلط أسماء الملفات ونسخها إلى بيانات / تقسيمات / قطار | val | اختبار /. تطبع وحدة التحكم الأعداد (على سبيل المثال ، القطار: 700 ، القابلية: 150 ، الاختبار: 150). نقطة التفتيش (مراقبة): افتح البيانات / الانقسامات / التدريب / وتحقق من أن الصور 256 × 256 وتدرج الرمادي (قناة واحدة). احتفظ ببيان CSV صغير (splits.csv) يسرد مسار الملف والطريقة والتسمية المقسمة.
المعالجة المسبقة للصور
استخدم نواة الشحذ 3 × 3. قم بإجراء تحسين الصورة باستخدام نواة التوضيح. تم استخدام نواة الالتفاف الشحذ لتعزيز الحدود التشريحية قبل التجزئة:

يبرز هذا المرشح الهياكل المهمة من خلال تعزيز المكونات عالية التردد4،5. تم تجميع كل صورة باستخدام مكتبة OpenCV الخاصة ب Python (cv2.filter2D ()) لإنشاء الصورة المحسنة. يوضح الشكل 6 المقارنات قبل وبعد.

الشكل 6: مخطط دقة التحقق من الصحة لقيم مختلفة من k. مخطط شريطي يوضح دقة التحقق من الصحة التي تم تحقيقها لقيم مختلفة من K المستخدمة في خطوة تجميع K-means (K = 2 ، 3 ، 4 ، 5). تبلغ الدقة ذروتها عند K = 3 ، مما يشير إلى الفصل الأمثل بين الهياكل التشريحية ومناطق الضوضاء. المقياس: قيم الدقة الطبيعية (0-1). الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.
قم بتطبيق (Python / OpenCV) باستخدام الكود أدناه.

أو استخدم واجهة المستخدم الرسومية البديلة (ImageJ/Fiji) بالنقر فوق معالجة > عوامل التصفية > Convolve. الصق مصفوفة 3 × 3 أعلاه وانقر فوق موافق > ملف > حفظ باسم > PNG في البيانات / preproc / enhanced /. راقب نقطة التفتيش حيث تبدو الحواف وحدود الأعضاء أكثر وضوحا. إذا شوهدت هالات شديدة الشحذ ، فقم بتقليل الوزن المركزي للنواة من 5 إلى 4.5 وأعد التشغيل. يوضح الشكل 7 النتائج قبل وبعد تقليل الضوضاء. تظهر الصور السابقة بوضوح ضوضاء وضبابية ، بينما تعرض الصور اللاحقة وضوحا محسنا ، وحدودا تشريحية أكثر وضوحا ، وتباينا محسنا ، مما يدل على فعالية خط الأنابيب المقترح لتقليل الضوضاء.

الشكل 7: بنية شبكة Convolution NAutoencoder. يوضح بنية التشفير الذاتي: طبقة الإدخال ، التشفير (طبقات التحويل العشوائي + البركة) ، عنق الزجاجة ، وحدة فك التشفير (أخذ العينات + طبقات التحويل المحمولة). يتم تصنيف كل طبقة بالحجم والوظيفة. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.
تجميع K-Means للتجزئة
تمت إعادة تشكيل الصور المحسنة باستخدام NumPy (image.reshape(-1, 1)) وتجميعها باستخدام sklearn.cluster.KMeans(n_clusters=3 أو 5). تمت إعادة تشكيل الإخراج المجزأ مرة أخرى إلى 2D (np.reshape (clustered_array ، image.shape)) لتصور المناطق التشريحية مقابل مناطق الضوضاء14. يسرد الجدول 1 إعدادات التجزئة.
| رقم الريال | قيمة ك | دقة |
| 1 | 3 | 0.761 |
| 2 | 5 | 0.869 |
| 3 | 7 | 0.75 |
الجدول 1: قيم k مع دقة التحقق من الصحة الخاصة بها. يعرض هذا الجدول دقة التحقق من الصحة التي تم الحصول عليها لقيم مختلفة لمعلمة التجميع k المستخدمة في خطوة تجزئة K-means لخط أنابيب إزالة الضوضاء. توضح النتائج أن k = 5 ينتج عنه أعلى دقة ، مما يوجه الاختيار الأمثل لمعلمات التجميع.
تجزئة K-Means: إعادة التشكيل والكتلة (Python / scikit-learn) باستخدام الكود أدناه.

قم بإجراء معاينة الألوان (اختياري) وتعيين التسميات إلى الألوان لمراقبة الجودة المرئية وحواف التراكب من الصورة الحادة. اختر K عن طريق الجري باستخدام K = 3 و K = 5 ؛ حساب دقة التحقق من الصحة في اتجاه مجرى النهر (الجدول 1 يسرد الإعدادات; يوضح الشكل 6 الدقة مقابل K). هذه نقطة تفتيش. مراقبة وحدات البكسل المهيمنة على الضوضاء التي تشكل مجموعات مميزة ؛ يجب أن تظل المناطق التشريحية متجاورة. إذا ظهرت بقع صغيرة ، فقم بتطبيق الحدود التشريحية.
تقليل الضوضاء المستندة إلى الشبكة العصبية
وصف الهندسة المعمارية: انظر الشكل 8 للاطلاع على مخطط وحدة فك تشفير عنق الزجاجة. تم تطوير جهاز تشفير ذاتي قائم على CNN باستخدام TensorFlow / Keras. تضمنت الهندسة المعمارية: طبقة الإدخال: صورة ذات تدرج رمادي 256 × 256 ، جهاز تشفير بثلاث طبقات تلافيفية (النواة: 3 × 3 ، خطوة: 1 ، تنشيط ReLU) ، كل منها متبوعا بأقصى تجميع ؛ عنق الزجاجة كتمثيل كامن كثيف ، وحدة فك ترميز بثلاث طبقات لأخذ العينات مع الالتفافات المنقولة ، الإخراج كطبقة منشطة بالسيني تنتج صورة خالية من الضوضاء. يوفر الجدول 2 معلمات معمارية كاملة للطبقة.

الشكل 8: النتائج المرئية لعملية إزالة الضوضاء. يقدم صورة صاخبة أصلية وصورة معالجة مسبقا ومخرجات نهائية غير ضوئية جنبا إلى جنب للتقييم النوعي. يوضح الحفاظ على الحواف وتقليل القطع الأثرية. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.
| رقم الريال | سمات التحويل البرمجي | قيم سمات التحويل البرمجي |
| 1 | محسن | آدم |
| 2 | ضياع | إنتروبيا متقاطعة فئوية |
| 3 | المقاييس | دقة |
الجدول 2: سمات التجميع للشبكة العصبية. يوضح هذا الجدول بالتفصيل معلمات التجميع الرئيسية المستخدمة لتدريب نموذج التشفير الذاتي التلافيفي. تم تنفيذ هذه الإعدادات في TensorFlow وتتضمن المحسن ووظيفة الخسارة ومقياس التقييم وعدد الفترات وحجم الدفعة.
تكوين التدريب: تضمنت معلمات التدريب: وظيفة الخسارة كمتوسط الخطأ التربيعي (MSE) ، والمحسن مثل آدم (معدل التعلم = 0.001) ، وحجم الدفعة 32 ، والعصور 100 مع التوقف المبكر (الصبر = 10).
فيما يلي لمحة موجزة عن التدريب. ابدأ التدريب بالانتقال إلى Colab والنقر فوق وقت التشغيل > تشغيل الكل ؛ تأكد من إدراج وحدة معالجة الرسومات (على سبيل المثال ، Tesla T4). استخدم رمز python / train_autoencoder.py --data_root البيانات --epochs 100 --batch 32 --lr 0.001 --early_stop 10. المراقبة باستخدام وحدة التحكم تطبع العصر train_loss و val_loss والوقت/الحقبة. التوقف المبكر وراء الصبر = 10 حقب دون تحسن. يتم حفظ أفضل نموذج في runs/.../checkpoints/best.h5. هذه نقطة تفتيش ، لاحظ باستخدام model.summary () الذي يعرض مكدس الطبقة. يجب أن يتطابق عدد المعلمات مع الجدول 2. إذا تلقيت أخطاء OOM، فقم بتقليل الدفعة إلى 16 أو قم بتعيين حجم الإدخال إلى 224 × 224.
نقاط التفتيش المتوسطة واستكشاف الأخطاء وإصلاحها
بعد التوضيح والتجزئة ، قم بمعاينة ثلاث صور لكل تقسيم وتأكد من (i) تحسين الحواف ، (ii) محاذاة أقنعة الكتلة مع علم التشريح. أثناء التدريب ، تأكد من انخفاض فقدان التحقق من الصحة وعدم تباعده. إذا بدا الإخراج غير الضوضي مفرطا في سلاسة ، فقم بتقليل مركز النواة (4.5) ، أو قم بزيادة فترات التدريب بمقدار 10 بمعدل تعلم أقل (على سبيل المثال ، 5e-4).
تقييم الأداء: تم استخدام المقاييس الكمية التالية: ذروة نسبة الإشارة إلى الضوضاء (PSNR) ، ومقياس مؤشر التشابه الهيكلي (SSIM) ، ودقة التحقق من تصنيف الصور غير الضوضية.
قامت الطريقة المقترحة بتحسين PSNR من 21.52 إلى 28.14 ديسيبل و SSIM من 0.76 إلى 0.86 مقارنة بنماذج خط الأساس المعروضة في 7,16. تم تسجيل كفاءة الطاقة من خلال مراقبة استخدام وحدة معالجة الرسومات (سجلات NVIDIA-SMI) ووقت التدريب. يلخص الجدول 3 النتائج المخفضة للضوضاء.
| رقم الريال | متري | نموذج خط الأساس |
| 1 | متوسط وقت العصر (ثانية) | 25.8 |
| 2 | إجمالي طاقة التدريب (كيلوواط ساعة) | 0.52 |
| 3 | استخدام وحدة معالجة الرسومات (٪) | 85% |
| 4 | وقت الاستدلال لكل صورة (مللي ثانية) | 18.7 |
| 5 | دقة التحقق من الصحة (٪) | 76.19% |
| 6 | PSNR (ديسيبل) | 21.52 |
| 7 | SSIM | 0.7619 |
الجدول 3: مقاييس تقييم الأداء لخط الأساس مقابل الطريقة المقترحة. يقارن هذا الجدول خط أنابيب تقليل الضوضاء المقترح بنموذج خط الأساس عبر العديد من مقاييس الأداء ، بما في ذلك وقت التدريب واستخدام وحدة معالجة الرسومات ومقاييس الجودة (PSNR و SSIM ودقة التحقق من الصحة). تشير النتائج إلى تحسين كفاءة الطاقة وجودة الصورة بالطريقة المقترحة.
حساب PSNR / SSIM واستخدام الطاقة. بالنسبة إلى PSNR / SSIM (صورة scikit) ، استخدم الكود أدناه.

لاستخدام طاقة الحوسبة/وحدة معالجة الرسومات، قم بتشغيل الخلية التي تسجل nvidia-smi --query-gpu=power.draw,utilization.gpu --format=csv -l 1 لتشغيل /.../gpu_log.csv أثناء التدريب (متوفرة في دفتر الملاحظات). في محطة طرفية ثانية، قم بتشغيل التعليمات البرمجية التالية:

تحليل CSV لحساب متوسط الطاقة (W) والتكامل بمرور وقت التدريب للطاقة المقدرة (Wh = kWh). يجمع الكمبيوتر المحمول PSNR و SSIM ودقة التحقق من الصحة ووقت العصر والطاقة في results_table3.csv للتضمين المباشر.
تقييم الاستدامة ومحاكاة التطبيب عن بعد
لتقييم استدامة أجهزة الجهاز ، قمنا بتطبيق ضوضاء Gaussian و Poisson لمحاكاة تدهور الصورة من الأجهزة القديمة. أعاد النموذج هذه المدخلات المتدهورة إلى جودة شبه تشخيصية ، مما يتحقق من قوتها27،28،29. في محاكاة التطبيب عن بعد ، تم إرسال الصور غير الضوضوية عن طريق محاكاة عرض النطاق الترددي 256 كيلوبت في الثانية باستخدام مآخذ بايثون. تم الاحتفاظ بالوضوح البصري ، ودعم التشخيص عن بعد30،31،32.
لمحاكاة شيخوخة الجهاز ، قم بتطبيق الضوضاء الغاوسية (σ = 10-30) وضوضاء بواسون لتنظيف الصور (محاكاة تقادم الخلية) وحفظها في البيانات / محاكاة / قديمة /. تشغيل النموذج المدرب على المدخلات القديمة / ؛ حفظ المخرجات في النتائج/aged_denoised/. لاحظ أن الجودة المرئية يجب أن تصل إلى دقة شبه تشخيصية ؛ قارن PSNR / SSIM مع خط الأساس.
لاختبار عرض النطاق الترددي للتطبيب عن بعد ، قم بضغط الصور غير الضوئية إلى PNG / JPEG بنسبة 85٪ -95٪ وأرسل ارتباطا محاكاة 256 كيلوبت في الثانية باستخدام اختبار مقبس Python (telemed_sim.py) المقدم قياس وقت الرحلة ذهابا وإيابا وتجزئة سلامة الملفات. تحقق من عدم تقديم أي عناصر تشخيصية. لاحظ أنه يتم الحفاظ على الوضوح المرئي وأحجام الملفات المناسبة لمهام سير العمل ذات النطاق الترددي المنخفض.
نقطة نهاية البروتوكول
عند الانتهاء ، يجب أن يكون لديك: (1) صور حادة في preproc / enhanced / ، (ii) صور مجزأة في preproc / مجزأة / ، (3) جهاز تشفير ذاتي مدرب مع best.h5 في عمليات التشغيل /.../نقاط التفتيش / ، (4) مخرجات غير ضوضاء للاختبار والصور القديمة في النتائج / (v) ملف مقاييس مترجم (results_table3.csv) يلخص PSNR و SSIM ودقة التحقق من الصحة ووقت العصر والطاقة المقدرة.