إعادة أخذ العينات باستخدام لانكزوس

استيفاء لانكزوس بنصف قطر 1
استيفاء لانكزوس بنصف قطر 2
استيفاء لانكزوس بنصف قطر 3
رسم بياني جزئي لإشارة منفصلة (نقاط سوداء) وتقريب لانكزوس الخاص بها (منحنى أزرق متصل)، مع معامل الحجم a يساوي 1 (أعلى)، 2 (وسط)، و3 (أسفل). كما يظهر نسختان من نواة لانكزوس، مُزاحتان ومُقاستان، تُقابلان العينتين 4 و11 (منحنيات متقطعة).

يُعدّ ترشيح لانكزوس وإعادة تشكيل لانكزوس تطبيقين لصيغة رياضية محددة. يمكن استخدامها كمرشح تمرير منخفض أو لاستكمال قيمة الإشارة الرقمية بسلاسة بين عيناتِها . في الحالة الأخيرة، تُحوّل كل عينة من الإشارة المعطاة إلى نسخة مُترجمة ومُقاسة من نواة لانكزوس ، وهي دالة sinc مُقسّمة بواسطة الفص المركزي لدالة sinc ثانية أطول. ثم يُحسب مجموع هاتين النواتين المُترجمة والمُقاسة عند النقاط المطلوبة.

تُستخدم إعادة التشكيل لانكزوس عادةً لزيادة معدل أخذ العينات للإشارة الرقمية، أو لإزاحتها بمقدار جزء من فاصل أخذ العينات. كما تُستخدم غالبًا في الاستيفاء متعدد المتغيرات ، على سبيل المثال لتغيير حجم الصورة الرقمية أو تدويرها . وقد اعتُبرت "الحل الأمثل" من بين العديد من المرشحات البسيطة لهذا الغرض. [ 1 ]

اخترع المرشح كلود دوشون ، الذي أطلق عليه اسم كورنيليوس لانكزوس بسبب استخدام دوشون لتقريب سيجما في بناء المرشح، وهي تقنية ابتكرها لانكزوس. [ 2 ]

تعريف

نواة لانكزوس

نوافذ لانكزوس لـ a = 1، 2، 3
نواة لانكزوس للحالات a = 1 و 2 و 3، مع أطياف تردداتها. سيكون لمرشح sinc قطع عند التردد 0.5.

يُحدد تأثير كل عينة إدخال على القيم المُستكملة بواسطة نواة إعادة بناء المرشح L ( x ) ، والتي تُسمى نواة لانكزوس. وهي دالة sinc المُعَيَّرة sinc( x ) ، مُطبَّقة (مُضروبة) في نافذة لانكزوس .أو نافذة sinc ، وهي الفص المركزي لدالة sinc الممتدة أفقيًا sinc( x / a ) لـ axa .

ل(x)={منذ(x)منذ(x/أ)لو -أ<x<أ،0خلاف ذلك.{\displaystyle L(x)={\begin{cases}\operatorname {sinc} (x)\operatorname {sinc} (x/a)&{\text{إذا كان}}\ -a<x<a,\\0&{\text{فيما عدا ذلك}}.\end{cases}}}

وبعبارة أخرى،

ل(x)={1لو x=0،أالخطيئة(πx)الخطيئة(πx/أ)π2x2لو -أx<أ و x0،0خلاف ذلك.{\displaystyle L(x)={\begin{cases}1&{\text{إذا كان}}\ x=0,\\{\dfrac {a\sin(\pi x)\sin(\pi x/a)}{\pi ^{2}x^{2}}}&{\text{إذا كان}}\ -a\leq x<a\ {\text{و}}\ x\neq 0,\\0&{\text{فيما عدا ذلك}}.\end{cases}}}

المعامل a هو عدد صحيح موجب، عادةً 2 أو 3، يحدد حجم النواة. تحتوي نواة لانكزوس على 2a - 1 فصًا: فص موجب في المركز، و a - 1 فصًا متناوبًا بين السالب والموجب على كل جانب.

صيغة الاستيفاء

بالنظر إلى إشارة أحادية البعد مع عينات s i ، بالنسبة للقيم الصحيحة لـ i ، يتم الحصول على القيمة S ( x ) المستوفاة عند وسيط حقيقي عشوائي x عن طريق الالتفاف المنفصل لتلك العينات مع نواة لانكزوس: [ 3 ]

S(x)=أنا=x-أ+1x+أsأنال(x-أنا)،{\displaystyle S(x)=\sum _{i=\lfloor x\rfloor -a+1}^{\lfloor x\rfloor +a}s_{i}L(xi),}

حيث يمثل a معامل حجم المرشح، وx{\displaystyle \lfloor x\rfloor }هي دالة الجزء الصحيح . حدود هذا المجموع هي بحيث تكون النواة صفرًا خارجها.

ملكيات

طالما أن المعامل a عدد صحيح موجب، فإن نواة لانكزوس متصلة في كل مكان، ومشتقتها معرفة ومتصلة في كل مكان (حتى عند x = ± a ، حيث تؤول كلتا دالتي sinc إلى الصفر). لذلك، ستكون الإشارة المُعاد بناؤها S ( x ) متصلة أيضًا، ومشتقتها متصلة.

تكون نواة لانكزوس صفرًا عند كل وسيط صحيح x ، باستثناء x = 0 ، حيث تكون قيمتها 1. لذلك، فإن الإشارة المعاد بناؤها تقوم باستيفاء العينات المعطاة بدقة: سيكون لدينا S ( x ) = s i لكل وسيط صحيح x = i .

إعادة أخذ العينات لانكزوس هي أحد أشكال الطريقة العامة التي طورها لانكزوس لمواجهة ظاهرة جيبس ​​عن طريق ضرب معاملات متسلسلة فورييه المقتطعة بـsأنانج(πك/م){\displaystyle \mathrm {sinc} (\pi k/m)}، أينك{\displaystyle k}هو معامل الدليل وم{\displaystyle m}يمثل هذا عدد المعاملات التي نحتفظ بها. [ 4 ] [ 5 ] وينطبق المنطق نفسه في حالة الدوال المقتطعة إذا أردنا إزالة تذبذبات جيبس ​​في طيفها.

الاستيفاء متعدد الأبعاد

بداية صورة بالأبيض والأسود. نسخة أصلية منخفضة الجودة مع تشوهات JPEG .
أُعيدت معاينة الصورة نفسها بخمسة أضعاف عدد العينات في كل اتجاه، باستخدام إعادة معاينة لانكزوس. وتمت إزالة تشوهات البكسلة عن طريق تغيير دالة نقل الصورة.

نواة مرشح لانكزوس في بعدين هي

ل(x،y)=ل(x)ل(y).{\displaystyle L(x,y)=L(x)L(y).}

تقييم

المزايا

نافذة لانكزوس المنفصلة واستجابتها الترددية ؛ انظر دالة النافذة للمقارنة مع النوافذ الأخرى.

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

تمت مقارنة مرشح لانكزوس مع طرق استيفاء أخرى للإشارات المنفصلة، ​​وخاصةً الإصدارات الأخرى ذات النوافذ لمرشح sinc. زعم توركوفسكي وغابرييل أن مرشح لانكزوس (مع a = 2 ) هو "أفضل حل وسط من حيث تقليل التداخل، والحدة، والحد الأدنى من الرنين"، مقارنةً بمرشح sinc المقتطع ومرشحات sinc ذات نوافذ بارتليت ، وجيب التمام ، وهان ، وذلك لتقليل عدد العينات واستيفاء بيانات الصور ثنائية الأبعاد. [ 1 ] ووفقًا لجيم بلين ، فإن نواة لانكزوس (مع a = 3 ) "تحافظ على الترددات المنخفضة وترفض الترددات العالية بشكل أفضل من أي مرشح (متاح) رأيناه حتى الآن". [ 6 ]

يعد استيفاء لانكزوس مرشحًا شائعًا لـ "ترقية" مقاطع الفيديو في العديد من أدوات الوسائط، مثل AviSynth [ 7 ] و FFmpeg . [ 8 ]

القيود

بما أن النواة تفترض قيمًا سالبة عندما تكون قيمة a أكبر من 1 ، فقد تكون الإشارة المُستكملة سالبة حتى لو كانت جميع العينات موجبة. وبشكل عام، قد يكون نطاق قيم الإشارة المُستكملة أوسع من نطاق قيم العينات المنفصلة. وعلى وجه الخصوص، قد تظهر تشوهات رنينية قبل وبعد التغيرات المفاجئة في قيم العينات، مما قد يؤدي إلى تشوهات قص . ومع ذلك، فإن هذه التأثيرات أقل مقارنةً بمرشح sinc (غير المُقسّم إلى نافذة). فعندما تكون a  =  2 (نواة ثلاثية الفصوص)، تكون نسبة الرنين  أقل من  1%.

عند استخدام مرشح لانكزوس لإعادة تشكيل الصور، يُنتج تأثير التموج هالات فاتحة وداكنة على طول الحواف البارزة. ورغم أن هذه الهالات قد تكون مزعجة بصريًا، إلا أنها تُسهم في زيادة حدة الصورة المُدركة ، وبالتالي تُحسّن من وضوح الحواف . وهذا بدوره قد يُحسّن من جودة الصورة، نظرًا للدور المحوري الذي يلعبه وضوح الحواف في الرؤية . [ 9 ]

في بعض التطبيقات، يمكن تحسين تشوهات القطع في الترددات المنخفضة عن طريق تحويل البيانات إلى نطاق لوغاريتمي قبل التصفية. في هذه الحالة، ستكون القيم المُستكملة عبارة عن متوسط ​​هندسي مرجح، وليس متوسطًا حسابيًا، لعينات الإدخال.

لا تمتلك نواة لانكزوس خاصية تجزئة الوحدة . أي أن المجموعيو(x)=أناZل(x-أنا){\textstyle U(x)=\sum _{i\in \mathbb {Z} }L(x-i)}لا يساوي عدد جميع النسخ المترجمة من النواة دائمًا 1. لذلك، فإن استيفاء لانكزوس لإشارة منفصلة ذات عينات ثابتة لا ينتج عنه دالة ثابتة. يظهر هذا العيب بوضوح عندما تكون قيمة a تساوي 1. كذلك، عند a = 1، تكون مشتقة الإشارة المستوفاة صفرًا عند كل وسيط صحيح. هذا الأمر نظريٌّ إلى حدٍّ ما، لأن استخدام نواة أحادية الفص ( a = 1) يُفقد جميع مزايا منهج لانكزوس ويُنتج مرشحًا ضعيفًا. توجد العديد من دوال التقطيع الأفضل أحادية الفص ذات الشكل الجرسيّ.   

يمكن إدخال تجزئة الوحدة عن طريق التطبيع،

ل(x-أنا)=ل(x-أنا)ج=1-أأل(x-ج){\displaystyle L'(x-i)={\frac {L(x-i)}{\sum _{j=1-a}^{a}L(x-j)}}}

ل0x<1{\displaystyle 0\leq x<1}.

انظر أيضاً

مراجع

  1. 1 2 توركوفسكي، كين؛ غابرييل، ستيف (1990). "مرشحات لمهام إعادة التشكيل الشائعة". في غلاسنر، أندرو س. (محرر). جواهر الرسومات 1. مطبعة أكاديمية. ص 147-165 . CiteSeerX 10.1.1.116.7898 . ISBN   978-0-12-286165-9.
  2. كلود، دوشون (1979-08-01). "ترشيح لانكزوس في بُعد واحد وبُعدين" . مجلة الأرصاد الجوية التطبيقية . 18 (8): 1016-1022 . Bibcode : 1979JApMe..18.1016D . doi : 10.1175/1520-0450(1979)018 < 1016:LFIOAT > 2.0.CO ; 2 .
  3. برجر، فيلهلم؛ بورج، مارك جيه. (2009). مبادئ معالجة الصور الرقمية: الخوارزميات الأساسية . سبرينغر. ص 231-232 . ISBN  978-1-84800-194-7.
  4. لانكزوس، كورنيليوس (1988). التحليل التطبيقي . نيويورك: منشورات دوفر. ص 219-221 . ISBN  0-486-65656-X. OCLC 17650089 . 
  5. مازولي، فرانشيسكو. "شرح استيفاء لانكزوس" . mazzo.li . تم الاسترجاع في 21 أبريل 2025 .
  6. ^ بلين، جيم (1998). ركن جيم بلين: البكسلات القذرة . مورجان كوفمان. ص 26 – 27. ISBN  978-1-55860-455-1.
  7. "تغيير الحجم" . Avisynth. 2015-01-01 . تم الاسترجاع في 2015-07-27 .
  8. "دليل إرشادي: تحويل الفيديو إلى صيغة أعلى باستخدام FFDShow - منتديات Neowin" . Neowin.net. 18-04-2006 . تم الاطلاع عليه بتاريخ 31-07-2012 .
  9. "IPOL: طرق خطية لاستيفاء الصور" . Ipol.im. 2011-09-27 . تم الاطلاع عليه بتاريخ 2012-07-31 .
  • أمثلة على هندسة مكافحة الحبيبات : image_filters.cppتوضح مقارنات إعادة أخذ عينات متكررة من صورة باستخدام نوى مختلفة.
  • imageresampler : فئة إعادة تشكيل الصور في المجال العام بلغة C++ مع دعم لعدة نوى مرشح Lanczos ذات النوافذ.