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

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

ただいまの
回答率

90.61%

  • Python

    7539questions

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

  • Python 3.x

    5949questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

2軸グラフ 2つのファイル 異なる時間

解決済

回答 4

投稿

  • 評価
  • クリップ 0
  • VIEW 269

shiroikuma

score 2

 python3.6 異なる時間 2軸グラフの書き方

2つのファイルから作成する2つデータフレームから、
同じグラフ上に2軸でデータをプロットをしたいのですがうまくいきません。

2つのファイルでは、DateTime の解像度が異なり、
一つは秒単位、一つは日単位となっています。

うまく2軸のグラフを描画するにはどうしたらいいのでしょうか?
2軸失敗図
イメージ説明

別々に描画した際
df
イメージ説明

df2
イメージ説明

 該当のソースコード

#!/usr/bin/env python                                                                                                                  

import pandas as pd
import matplotlib.pyplot as plt
import matplotlib.dates as dates
from datetime import datetime,timedelta

font = {'family' : 'meiryo'}

fig, ax1 = plt.subplots()

df = pd.read_csv("401190663.csv",header=None, names=['date','mjd','igrf','lat','lon','atl'])
df2= pd.read_csv("kpap.csv")


#print(df)                                                                                                                             
df['date'] = pd.to_datetime(df['date'])
df2['date'] = pd.to_datetime(df2['date'],format='%Y/%m/%d')
df.index = df['date']
df2.index = df2['date']

#1                                                                                                                                     
df.plot(y='igrf',ax=ax1,figsize=(16,4))
ax2= ax1.twinx()
#2                                                                                                                                     
df2.plot(y='ave_ap',ax=ax2,color='r')
plt.savefig("igrf_index.png")
  • 気になる質問をクリップする

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 4

+3

別々に描画したグラフを見ると、2つのグラフの年月日範囲は概ね一致しています。

秒単位のほうのデータを日単位に変換(秒単位の時間を1日の秒数(86400)で割れば、日単位に変換できる)してから両データをグラフ化すれば、別々に描画したグラフを重ねたようなグラフになると思います。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/30 09:26

    回答ありがとうございます。できました!

    キャンセル

checkベストアンサー

+1

DataFrameにresampleというメソッドがあります。これはデータを指定した単位でまとめるものなので、日次や週次で平均または合計してあげればいいと思います。使い方は以下の通りです。

df.resample('D').mean() #日次平均
df.resample('W').sum()  #週次合計


resampleの引数にまとめたい単位を設定し、その後に行いたい処理を記述します。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/30 09:25

    できました!
    s = df.resample('D').mean() でプロットしたら、意図したグラフになりました。ありがとうございました。

    キャンセル

+1

pandas.merge_asof() を使用して、1日単位のデータdf2に対して 1秒単位のデータdf を DatetimeIndexを揃えた状態で結合してから plot() するというのはどうでしょうか。

以下はサンプル

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt

# df は1秒単位のデータ1年分を適当に生成
N = 60*60*24*365
df = pd.DataFrame({'date':pd.date_range('2017/01/01', freq='1s', periods=N), 'igrf': np.random.randint(0,100,N), 'lat': np.random.randint(0,100,N)})

# df2 は1日単位のデータ1年分を適当に生成
N2 = 365
df2 = pd.DataFrame({'date':pd.date_range('2017/01/01', freq='1d', periods=N2), 'ave_ap': np.random.randint(0,100,N2)})

# それぞれの'date'カラムをIndexに設定
df  = df.set_index('date')
df2 = df2.set_index('date')

# df2 に df を時間軸をそろえてマージ
merged_df = pd.merge_asof(df2,df,left_index=True,right_index=True, direction='nearest')


merged_df.plot(y=['igrf', 'ave_ap'])
plt.show()

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/30 09:29

    回答ありがとうございます。ソースを参考に描き直してみたところ、
    raise ValueError(msg.format(side='right'))
    ValueError: right keys must be sorted
    となり、エラーが出ました。おそらく、データに欠損があり、日単位と秒単位の日付数が微妙に異なるせいかな?と想像します。別の方の回答も見まして、日時平均をとりましたら、解決しました。ありがとうございました。

    キャンセル

0

イメージ説明
回答いただき、日時平均をとりましたらできました。皆様ありがとうございました。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

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

  • ただいまの回答率 90.61%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

  • 受付中

    Matplotlib(Python)のフォントに関して

    前提・実現したいこと Pythonのmatplotlibについて、フォントに関する質問があります。x軸とy軸のラベルフォントを「sans-serif」に設定した上で一部の文字を斜

  • 解決済

    matplotlibでx軸の範囲を変更したい

    x軸の値を変更できずに困っています from matplotlib import pyplot as plt def main(): data = [120,334,44

  • 解決済

    pythonの文字列について

    お世話になっております。 現在、pythonの勉強をしております。 グラフを作成するpythonを書きました。 import sys import os impo

  • 解決済

    pandasで散布図を作成したい

    したいこと 下記のcsvファイルがあります X軸に乗客、y軸にうどんをセットして散布図を作成したいです。 駅,うどん,乗客 青井駅,9,1000 新橋駅,100,1000

  • 解決済

    三次元グラフ同時に表示できるようにしたい

    前提・実現したいこと 同時に2枚(またはそれ以上)表示されるようにしたい 発生している問題・エラーメッセージ 2枚表示はされるが完全に重なった位置に表示されてしまうため、

  • 解決済

    pythonのmatplotlibで3次元の直方体を描きたい

    pythonのmatplotlibのグラフ上に3次元の直方体を描きたいです.調べてみたらmplot3dなどを使えば出来そうな気がしましたが,使い方がよくわかりませんでした.pyth

  • 解決済

    pythonの図示 matpotlib

    pythonで四角形が格子状に区切られそれぞれの点に数字が振られています。それをコンターマップのように(ただし、等高線はいらない)図示したいのですがどのように書けばよいでしょうか

  • 解決済

    Pythonのmatplotlibを使って三次元散布図をシンプルに描写したい

    前提・実現したいこと Python初心者です。 1人3属性からなる1000件ほどのデータを、3次元の散布図として描写したいです。 属性一つ一つがそのままX,Y,Zに対応するイメージ

同じタグがついた質問を見る

  • Python

    7539questions

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

  • Python 3.x

    5949questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。