2次元配列データの加工方法

Question

２次元の配列データを加工したいのですが、うまくいきません。
素人質問で恐縮ですが、ご回答頂けると幸いです。

■やりたいこと
・時系列のcsvファイルのデータ(3行×1000列)がある
・3行×3列を１セットとする
・1列ずつずらしてセットを作成する
⇒最初は13列目、２番目は24列目
・それらを１つのシートにまとめる

※イメージ
・csvファイルのデータ
t a b c
0.01 1.00 1.02 1.07
0.02 1.20 1.03 1.08
0.03 1.32 1.04 1.07
0.04 1.22 1.08 1.08
0.05 1.05 1.01 1.07
：
↓
・加工後のデータ
0.01 1.00 1.02 1.07
0.02 1.20 1.03 1.08
0.03 1.32 1.04 1.07
0.02 1.20 1.03 1.08
0.03 1.32 1.04 1.07
0.04 1.22 1.08 1.08
0.03 1.32 1.04 1.07
0.04 1.22 1.08 1.08
0.05 1.05 1.01 1.07

下記のコードで実装しましたが、
エラーが出ました。

i = 0
X0 = []
X1 = []
X2 = []
X0 = np.genfromtxt("a.csv", delimiter=',', dtype=np.float32, usecols=(0, 1, 2))

while i < 1000-2:
X1 =X_test0[i:i+2, :]
X2 = np.append(X1)
i+1
print("X_test")

X2 = np.append(X1)
TypeError: append() missing 1 required positional argument: 'values'

説明が分かりにくいかもしれませんが、ご回答頂けると幸いです。
宜しくお願い致します。

Answer 1

こんな感じですか？

データのサンプル：

a=np.arange(30).reshape(3,10)
a=a.transpose()#  10 x 3 配列

新しい配列 b の作り方

i=np.array([0,1,2]).reshape(3,1)
j=np.arange(8)
k= (i + j).flatten('F')
b=a[k,:]

一つ一つprintで内容見れば、何をやっているのかよくわかるはずです。

Comments

[arumon]

ご丁寧にありがとうございます。
printで中身を見ると何をやっているのかわかりやすかったです。
実装できました。

１点、追加質問で恐縮です。
今回はSETの塊を１列ずつずらしていましたが、ずらす列数を指定する事はできるのでしょうか？

今回ご提示いただいたflatten関数を調べましたが、どう設定すればいいのか分からず詰まってしまいました。。
何度も申し訳ありませんが、ご回答頂けると幸いです。
宜しくお願い致します。

[WathMorks]

できますよ。
jを等差数列にすればいいと思います。
2つ飛ばしなら
j=np.arange(0,8,2)
とすればO.K.です。質問の真意をこちらがちゃんとくみ取れているか自信ありませんが。

[arumon]

ご回答ありがとうございます。
実装できました。
私の説明不足で恐縮でしたが、質問の意図はくみ取って頂き、ありがとうございます。

おかげさまで解決しました。
ありがとうございます。

Answer 2

ダブルループの内包表記でもさくっと書けます。

x0[[i + j for i in range(x0.shape[0] - 2) for j in range(3)], :]

これと等価です。

k=[]
for i in range(x0.shape[0] - 2):
    for j in range(3):
        k.append(i+j)
x0[k,:]