質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
Matplotlib

MatplotlibはPythonのおよび、NumPy用のグラフ描画ライブラリです。多くの場合、IPythonと連携して使われます。

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

pandas

Pandasは、PythonでRにおけるデータフレームに似た型を持たせることができるライブラリです。 行列計算の負担が大幅に軽減されるため、Rで行っていた集計作業をPythonでも比較的簡単に行えます。 データ構造を変更したりデータ分析したりするときにも便利です。

Q&A

解決済

1回答

8813閲覧

時系列データフレームのn行ごとの平均値作成について

wtxa111

総合スコア18

Matplotlib

MatplotlibはPythonのおよび、NumPy用のグラフ描画ライブラリです。多くの場合、IPythonと連携して使われます。

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

pandas

Pandasは、PythonでRにおけるデータフレームに似た型を持たせることができるライブラリです。 行列計算の負担が大幅に軽減されるため、Rで行っていた集計作業をPythonでも比較的簡単に行えます。 データ構造を変更したりデータ分析したりするときにも便利です。

0グッド

1クリップ

投稿2018/10/26 10:37

編集2018/10/28 22:31

いつもお世話になります。プログラミング初心者です。
10分毎に観測したデータのデータフレームより、
n行ごとに平均化したデータフレームを作成してグラフ化まで
もっていきたいと考えています。

datetime a b group

0 2018-05-29 20:00:00 19.3 46.2 0
1 2018-05-29 20:10:00 20.2 46.7 0
2 2018-05-29 20:20:00 21.7 47.4 0
3 2018-05-29 20:30:00 23.7 47.3 0
4 2018-05-29 20:40:00 25.3 45.7 0
5 2018-05-29 20:50:00 26.9 45.7 0
6 2018-05-29 21:00:00 27.0 48.5 1
7 2018-05-29 21:10:00 32.1 40.6 1
8 2018-05-29 21:20:00 32.6 38.0 1
9 2018-05-29 21:30:00 31.8 37.3 1
10 2018-05-29 21:40:00 31.4 42.4 1
11 2018-05-29 21:50:00 32.9 42.4 1

例えば、上記のようなデータセットを読み込み、
他の皆様の質問を頼りに、groupbyを用いて以下のようにしてみましたが、

python

1#モジュールのインポート 2import matplotlib.pyplot as plt 3import pandas as pd 4import numpy as np 5 6#csvの読み込み 7df = pd.read_csv('XXX.csv',encoding='ANSI',parse_dates=[0]) 8 9#datetimeindexの作成 10index = pd.date_range 11 12#グループ番号を付け加える 13n = 6 14df['group'] = [x // n for x in range(0,len(df))] 15 16#n要素ごとの平均値を求める 17dfg = df.groupby("group").mean() 18 19print(df) 20print(dfg)

こうすると、以下のようなデータセットが出来上がります。

a b

group
0 22.85 46.500000
1 31.30 41.533333

見ると2要素の3行毎の平均値ができてはいるのですが、
これですと、datetimeの列が抜けており、時系列グラフに表す時に少し不都合です。

例えば、
group datetime a b
o 2018-5-29 20: 22.85 46.500000
1 2018-5-29 21: 31.30 41.533333

のような形、つまり6行ごとにまとめた日時を挿入したデータフレームを作りたいのです。
何か良い方法はありませんでしょうか?初歩的な質問で申し訳ありません。
よろしくお願い致します。

※追記

#モジュールのインポート import matplotlib.pyplot as plt import pandas as pd import numpy as np #csvの読み込み df = pd.read_csv('ocean_current1.csv',encoding='ANSI',parse_dates=[0],index_col=[0]) #期間毎の平均値を求める(1H) dfg = df.resample('1H').mean() #列データに記号を指定 dt = dfg['datetime'] a = dfg['a'] b = dfg['b'] #流向の度数→radian変換 c = np.radians(b) #ベクトルの成分計算 u = a*np.sin(c) v = a*np.cos(c) KeyError: 'datetime'

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

「6行ずつ」ではなく「1時間単位で」となりますが、今回の場合は DataFrame.resample() を使かってまとめるのが適当ではないでしょうか。

Python

1import pandas as pd 2import io 3 4data=""" 5datetime,a,b 62018-05-29 20:00:00,19.3,46.2 72018-05-29 20:10:00,20.2,46.7 82018-05-29 20:20:00,21.7,47.4 92018-05-29 20:30:00,23.7,47.3 102018-05-29 20:40:00,25.3,45.7 112018-05-29 20:50:00,26.9,45.7 122018-05-29 21:00:00,27.0,48.5 132018-05-29 21:10:00,32.1,40.6 142018-05-29 21:20:00,32.6,38.0 152018-05-29 21:30:00,31.8,37.3 162018-05-29 21:40:00,31.4,42.4 172018-05-29 21:50:00,32.9,42.4 18""" 19 20df = pd.read_csv(io.StringIO(data), parse_dates=['datetime'],index_col='datetime') 21 22print(df.resample('1H').mean()) 23# a b 24#datetime 25#2018-05-29 20:00:00 22.85 46.500000 26#2018-05-29 21:00:00 31.30 41.533333

投稿2018/10/26 11:50

magichan

総合スコア15898

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

wtxa111

2018/10/28 22:27

ご回答ありがとうございます。きれいに1時間毎の平均が得られました。 しかし、これを基にして計算、グラフ化するのに新たに問題が出てきました。 例えば、1時間毎の平均を表す新しいデータフレームを #期間毎の平均値を求める(1H) dfg = df.resample('1H').mean() として指定します。そのあと、各列を記号に代入しますと、 #列データに記号を指定 dt = dfg['datetime'] a = dfg['a'] b = dfg['b'] #流向の度数→radian変換 c = np.radians(b) #ベクトルの成分計算 u = a*np.sin(c) v = a*np.cos(c) KeyError: 'datetime'が出てうまくいきません。 後で時系列ベクトル図を描く場合、これができないと不都合です。 何が原因となっているのでしょうか?
magichan

2018/10/28 23:51

追記されましたコードでは read_csv() のパラメータに index_col=[0]が指定されておりますので datetime列は Indexに設定されております。 datetime列の取得は dt = dfg['datetime'] ではなく dt = dfg.index とするとよいのではないでしょうか。
wtxa111

2018/10/29 01:31

ありがとうございます。 インデックス指定した場合は.index
wtxa111

2018/10/29 01:46

になるのですね。勉強になりました。 またどうぞよろしくお願い致します。 頑張ります。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問