Как сгладить трехмерный массив скользящего окна в двухмерный массив? - PullRequest
Новая
       25

Как сгладить трехмерный массив скользящего окна в двухмерный массив?

3 голосов
/ 26 июня 2019

Я использую Autoencoder LSTM в Python (Keras).У меня есть многовариантный ввод, и я использую метод скользящего окна, чтобы преобразовать его в правильный формат ввода LSTM.В конце я получаю вывод той же формы, что и окно.Затем я хочу преобразовать этот массив в исходную форму ввода.Может кто-нибудь помочь мне, как я должен сделать это?

Это мой код, чтобы поставить скользящее окно на многомерный сигнал: 


def window(samples, windows_size, step):
    m, n = samples.shape 
    print("\nold shape: ", m, "*", n)
    num_signals = n    
    num_samples = (samples.shape[0] - windows_size) // step + 1
    aa = np.empty([num_samples, windows_size, num_signals])

    for j in range(num_samples):
        for i in range(num_signals):
            aa[j, :, i] = samples[(j * step):(j * step + windows_size), i]
    samples = aa
    m ,n, k = samples.shape
    print("new shape: ", m, "*", n, "*", k)
    return samples

x = np.asarray([[1,0.1,0.1],[2,0.2,0.2],[3,0.3,0.3],[4,0.4,0.4],
                [5,0.5,0.5],[6,0.6,0.6],[7,0.7,0.7],[8,0.8,0.8]])

window(x, 3, 2)

old shape:  8 * 3
new shape:  3 * 3 * 3
Out[65]: 
array([[[1. , 0.1, 0.1],
        [2. , 0.2, 0.2],
        [3. , 0.3, 0.3]],

       [[3. , 0.3, 0.3],
        [4. , 0.4, 0.4],
        [5. , 0.5, 0.5]],

       [[5. , 0.5, 0.5],
        [6. , 0.6, 0.6],
        [7. , 0.7, 0.7]]])

1 Ответ

0 голосов
/ 26 июня 2019

Вы можете использовать это:

Примечание: шаг - это та же концепция, что и в CNN, количество элементов, которые вы пропускаете, чтобы получить следующее окно. 

inp = np.array([[[1. , 0.1, 0.1],
        [2. , 0.2, 0.2],
        [3. , 0.3, 0.3]],

       [[3. , 0.3, 0.3],
        [4. , 0.4, 0.4],
        [5. , 0.5, 0.5]],

       [[5. , 0.5, 0.5],
        [6. , 0.6, 0.6],
        [7. , 0.7, 0.7]]])

def restitch(array, stride):
    flat = array.flatten().reshape(-1,array.shape[2])
    keep = [i for i in range(len(flat)) if not(i%(stride+1)==0 and i>0)]
    return flat[keep]

restitch(inp, 2)

array([[1. , 0.1, 0.1],
       [2. , 0.2, 0.2],
       [3. , 0.3, 0.3],
       [4. , 0.4, 0.4],
       [5. , 0.5, 0.5],
       [6. , 0.6, 0.6],
       [7. , 0.7, 0.7]])
...