Как получить сигнал сердцебиения из этих данных?

Question

Как получить сигнал сердцебиения из этих данных?

Я делаю свои первые шаги в анализе данных, сбора / очистки.

Чтобы узнать, я пытаюсь создать простой код, который может обнаружить сердцебиение от цветовых вариаций от изображения, поступающего с камеры на iPhone.

Это то, что я сделал до сих пор:

Я создал код, который распознает лица и вырезает их из изображения.
тогда я получаю средний цвет от этого изображения лица.
Поскольку значения захватываются с нерегулярными временными интервалами, я интерполирую данные. Теперь у меня есть ряд значений, "взятых" через равные промежутки времени.

Графики данных у меня есть это:

Я думаю, что вижу что-то похожее на сердцебиение.

Чтобы подтвердить, что я не вижу артефактов с камеры, я выбираю прямоугольник на заднем плане и получаю тот же цвет с очень небольшим случайным отклонением. Цвет фона стабильный и практически не изменяется.

Теперь я применяю прямое косинусное преобразование к данным и получаю это:

У DCT есть отрицательные долины.

Двумя первыми элементами DCT являются пики при t=0 и t=0,10482311 секунд.

Я прочитал статью, где исследователь говорит, что сердцебиение возникает в диапазоне от 0,4 до 4 Гц. Я предполагаю, что это имеет отношение, почему первые два элемента DCT являются пиками...

Что мне теперь делать?

Я попытался применить пороговое значение к нулю каждый член из DCT ниже 150 и сделать обратный DCT, чтобы восстановить сигнал и получил это

Я уменьшил порог до 85 и получил это

Я не уверен, вижу ли я сердцебиение там. Я имею в виду пики, где начинаются ритмы.

Это данные через регулярные промежутки времени (256 выборок):

и это регулярные интервалы времени

Это код, который я до сих пор...

let numberOfSamples = interpolatedValues.count


// perform the Discreet Cosine Transform
let data = Array(UnsafeBufferPointer<Float>(start:interpolatedValues, count: numberOfSamples))

let DCTSetupForward = vDSP_DCT_CreateSetup(nil,
                                           vDSP_Length(numberOfSamples),
                                           vDSP_DCT_Type.II)
var dctResult = Array<Float>(repeating:0.0,
                       count: numberOfSamples)

vDSP_DCT_Execute(DCTSetupForward!, data, &dctResult)

// apply a threshold    
let stride = vDSP_Stride(1)
var threshold : Float = 85

vDSP_vthres(dctResult,
            stride,
            &threshold,
            &dctResult,
            stride,
            vDSP_Length(numberOfSamples))

// recreate the signal
let DCTSetupInverse = vDSP_DCT_CreateSetup(nil, vDSP_Length(numeroDados), vDSP_DCT_Type.III)

var dctSetupInverse = Array<Float>(repeating:0.0, count: numberOfSamples)
vDSP_DCT_Execute(DCTSetupInverse!,
                 dctResult,
                 &dctSetupInverse)

0

swift signal-processing dct accelerate

Источник

user316469 11 фев '19 в 16:51

0 ответов

Другие вопросы по тегам swift signal-processing dct accelerate