Медіанний фільтр (Median Filter)

[radioman]

... Так:...

Прикро, але код модифікує масив. На малюнку у першому стовпчику вміст "динамічного" масиву при включеному сортуванні; у другому - лише циклічне заповнення масиву (симуляція вхідних даних):

Симуляція надходження і сортування даних

PS. Я допустився помилки - Ваш код теж потребує копіювання масиву, адже це модифікація сортування включенням (Insertion sort); на даний час Ваша модифікація найшвидша - 6 мікросекунд з копіюванням масиву.

PPS. Як буде час порівняю час виконання в залежності від розташування елементів в масиві.

[radioman]

Допиляв код. Так як є бажання динамічно змінювати довжину масиву, зіткнувся з перепоною - розмір масиву не може бути заданий змінною; це пов'язано з тим, що розмір глобального масиву повинен бути відомим до початку роботи програми.
Але розмір локального масиву (у функції) можна задати змінною, тому що локальний масив виділяється з динамічної пам'яті та може бути будь якого розміру.
Тому масив для сортування розміщено у функції:

Код: Виділити все

/* функція розрахунку медіани*/
int MedianFilter(int *window, int intArrSize)      	// аргументами функції є посилання на масив та довжина масиву
{
  int newArray[intArrSize];                             // масив для функції

Копіювання масиву за топомогою memcpy(newArray, window, intArrSize); у функції чомусь не працює, тому повернувся до "перезапису" (можна попробувати зробити "вхідне" сортування та подивитись чи буде виграш у часі виконання):

Код: Виділити все

 /* копіювання масиву */
  //memcpy(newArray, window, intArrSize);
  /**/
  byte c;
  for (c = 0; c < intArrSize; c++)
  {
    newArray[c] = window[c];
  }

У такій реалізації найкращі результати показала модифікація алгоритму сортування включенням (Insertion sort) від pawa:

Код: Виділити все

 /* сортування включенням (Insertion sort) */
  int counter = 0;
  for (int i = 1; i < intArrSize; i++) 
  {
    for (int j = i; j > 0 && *(newArray + j - 1) > *(newArray + j); j--) 
    {
      counter++;
      int tmp = *(newArray + j - 1);
      *(newArray + j - 1) = *(newArray + j);
      *(newArray + j) = tmp;
    }
  }

Інший код, практично, без змін.

PS. Можна відмовитись від посилання на масив вхідних даних в аргументі функції. Варто у функцію передавати лише потік даних, а записувати їх у "материнський" масив вже у самій функції.
PPS. У додадому архіві тестовий скетч сортування включенням (Insertion sort) на довільний розмір вікна:

MedianFilter_pawa_210128-1x1.7z

Сортування включенням (Insertion sort)

(1.39 Кіб) Завантажено 2893 разів

[pawa]

Можна зробити масив максимальної довжини (буфер), а використовувати стільки елементів скільки буде потрібно.

[radioman]

Так і є:

Код: Виділити все

/* медіана pawa_radioman */
/* https://radioman.com.ua */
int intFilterLengthWdi = 6;
int window[] = {};

Але, для отримання медіани, потрібно задавати довжину масиву. Вхідний потік записується по лічильнику у "материнський" (буферний) масив (в симуляції - з тестового масиву):

Код: Виділити все

window[interruptCtr] = window0[interruptCtr]; 
 if ( interruptCtr++ >= intFilterLengthWdi - 1 )                                 // лічильник для симуляції
  {
    interruptCtr = 0;
  }

PS. Ви можете переглянути весь скетч за допомогою блокноту.