Блог

Должен быть какой-то алгоритм, который сделает это проще, чем то, что я делаю…

У меня есть два массива, каждый с двумя столбцами. Один столбец в обоих — это временная метка, а другой в обоих — измерение.

Что должно произойти, так это превратить это в единый массив: timestamp, measurement1, measurement2

Проблема в том, что временные метки часто не совпадают точно. В одном массиве может полностью отсутствовать значение за период времени, или временные метки могут быть отключены на незначительную величину (достаточно незначительную, чтобы было нормально назначить оба измерения одной и той же временной метке).

Существует ли какой-нибудь хорошо известный способ выполнения этой операции нечеткого слияния? Простая функция общественного достояния??

Начните с того, что задайте себе следующие вопросы: имеют ли массивы одинаковое количество элементов? Как вы хотите объединить два элемента с одинаковой временной меткой? Как вы хотите объединить два элемента с разными временными метками?

Вероятно, вам придется написать алгоритм самостоятельно. Что-то подобное было бы легко реализовать:

1. Начните с сортировки каждого из массивов по отдельности в порядке временных меток.
2. Объявляйте два итератора в начале каждого входного массива соответственно и пустой выходной массив.
3. Затем вы проверяете, какой из массивов имеет самую раннюю временную метку. Вызовите его РАНО, а другой ПОЗДНО.
   • Если РАННЕЕ значение близко к ПОЗДНЕМУ (меньше некоторой константы), вы применяете операцию слияния и вставляете результат в конец выходного массива. Увеличьте оба итератора и вернитесь к 3.
   • В остальном, РАНО — это далеко не ПОЗДНО. Вам нужно обработать пропущенное значение в ПОСЛЕДНЕМ массиве, возможно, повторив предыдущее значение или интерполировав его с помощью какой-либо функции. Решите, вставлять или нет значение в выходной массив. В этом случае вы только увеличиваете НАЧАЛЬНЫЙ итератор массива и возвращаетесь к 3.
4. Если вы достигли конца любого из массивов, остальная часть другого массива ЗАПАЗДЫВАЕТ. Возможно, вы захотите интерпретировать это как пропущенные значения, а также повторить или интерполировать измерения.
5. Возвращает выходной массив.

~

 OutputArray merge(InputArrayamp; a, InputArrayamp; b) {
    InputArray::iterator a_it = a.begin();
    InputArray::iterator b_it = b.begin();
    while(a_it != a.end() amp;amp; b_it != b.end()) {
        InputArray::iteratoramp; early = *a_it.timestamp < *b_it.timestamp ? a_it : b_it;
        InputArray::iteratoramp; late = *a_it.timestamp < *b_it.timestamp ? b_it : a_it;
        if(*late.timestamp - *early.timestamp < TIMESTAMP_CLOSE_ENOUGH) {
            output.timestamp = (*late.timestamp   *early.timestamp) / 2; // mean value
            output.measure1 = *a_it.measure;
            output.measure2 = *b_it.measure;
            outputArray.push_back(output);
            a_it  ; b_it  ;
        }
        else {
            output.timestamp = *early.timestamp;
            output.measure1 = *a_it.timestamp < *b_it.timestamp ? *a_it.measure : outputArray.back.measure1; // previous value if missing
            output.measure2 = *a_it.timestamp < *b_it.timestamp ? outputArray.back.measure2 : *b_it.measure;
            outputArray.push_back(output);
            early  ;
        }
    }

    InputArray::iteratoramp; late = a_it != a.end() ? a_it : b_it;
    InputArray::iterator late_end = a_it != a.end() ? a.end() : b.end();
    while(late != late_end) {
            output.timestamp = *late.timestamp;
            output.measure1 = a_it != a.end() ? *a_it.measure : outputArray.back.measure1; // previous value if missing
            output.measure2 = a_it != a.end() ? outputArray.back.measure2 : *b_it.measure;
            outputArray.push_back(output);
            late  ;
    }
    return outputArray;
}