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

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

ただいまの
回答率

88.80%

python heatmapのメモリ位置をずらす方法、右端のデータを表示する方法

解決済

回答 2

投稿 編集

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

mor1mi

score 1

前提・実現したいこと

Pythonで実験で得たデータを解析しています。
例えばある物理量の時間依存の温度依存性をとったとして、
色で強度を表した二軸のカラーマップを作成したいです。
plt.pcolor()を使用してかなり近いことはできているんですが、
軸が意図したものと違っていて困っています。

x=np.arange(1,1001,1)*1e-3 #time
y=np.arange(20,151,10) #temp
(i,j)=(0,0)
Z=np.zeros((len(x),len(y)))
for X in x:
    for Y in y:
        Z[i][j]=10*(1-np.exp(-X/10*Y))
        if Y==20:
            Z[i][j]=1
        if Y==150:
            Z[i][j]=1
        j+=1
    i+=1
    j=0
plt.pcolor(y,x,np.log(Z))
plt.yscale('log')
plt.show()
plt.close()


出力画像
出力画像
20Cなどのメモリ位置がデータの左側に
表示されるため、25Cのデータのように見える。
X軸をずらしてメモリ位置を各データの真下にできますか?
また150Cのデータが表記されません。
この問題を解決する方法を教えてください。

試したこと

plt.pcolor()ではlog表記できますが、X軸の位置がおかしいです。seaborn.heatmap()は、X軸の位置は良く全てのデータを表示できるのですが、
Y軸のlog表記を達成する方法が分かりませんでした。

df = pd.DataFrame(data=Z, columns=y,index=x)
ax = seaborn.heatmap(df)


出力画像
イメージ説明

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

質問への追記・修正、ベストアンサー選択の依頼

  • meg_

    2020/10/18 11:32

    > また150Cなど、右端に来るべきデータが表記されません。
    これについてはグラフにデータは表示されています。軸ラベルが表示されていないだけですね。

    キャンセル

  • mor1mi

    2020/10/18 11:40 編集

    meg_さん:
    ありがとうございます。
    すみません、分かりにくくて申し訳ありませんでした。
    追加しました。これでどうでしょうか。
    >これについてはグラフにデータは表示されています。軸ラベルが表示されていないだけですね。
    本当ですか?

    キャンセル

  • toast-uz

    2020/10/18 12:16

    print(plt.xlim())で表示範囲が出ますので確認してみてください

    キャンセル

回答 2

checkベストアンサー

+1

seabornでやってみました。
※手動でラベル付けしているため正確な値からは多少ずれることがあります。

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.ticker as ticker
import pandas as pd
import seaborn as sns


x=np.arange(1,1001,1)*1e-3 #time
y=np.arange(20,151,10) #temp


(i,j)=(0,0)
Z=np.zeros((len(x),len(y)))
for X in x:
    for Y in y:
        Z[i][j]=10*(1-np.exp(-X/10*Y))
        j+=1
    i+=1
    j=0

df = pd.DataFrame(Z, columns=y,index=x)

tick = ticker.ScalarFormatter(useOffset=False, useMathText=True)
tick.set_powerlimits((0,0))

tg = ["${}$".format(tick.format_data(10**i)) for i in np.arange(np.floor(np.min(np.log10(x))), np.ceil(np.max(np.log10(x))) + 1, 1)]

ax = sns.heatmap(df, yticklabels=len(df)//len(tg))
ax.set_yticklabels(tg)


plt.show()

イメージ説明

参考:How to use scientific notation in seaborn heatmap labels?

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2020/10/18 14:23

    x軸を離散値と捉えて、飛び飛びの値を中心としたヒストグラムにしたのですね。温度という連続量なので、ヒストグラムという発想には気付きませんでした。確かに、質問者様とののやりとりを振り返ると、ヒストグラムを要望されているように感じます。

    キャンセル

  • 2020/10/19 03:38

    非常に助かりました!meg_さん、toast-uzさん、ありがとうございました!

    キャンセル

+1

質問者様の要望を以下のように解釈しました。

  • x軸は20〜150のデータがあり、150まで描画されているが、データ範囲の外側もマージンをとって描画したい。
  • y軸は、既にlog表記はできているが、さらに目盛数値の指数部だけを表示したい。

plt.show()の前に以下2行を入れてください。

plt.yticks([0.001, 0.01, 0.1, 1], ['-3', '-2', '-1', '0'])
plt.xlim(10, 160)

なお、yticksを自動描画する方法も無きにしも非ずですが、このやり方が手っ取り早いようです。きっと、どなたかがチャレンジしてくれるでしょう。subplotを使う例しか見当たらず、調べきれませんでした。

描画結果は以下です。

イメージ説明

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2020/10/18 21:49

    ありがとうございます。申し訳ないのですが、やはりplt.pcolor()では、150のデータは表示されていません。データを見て頂けますと、150には、確認のために、20Cと同じデータを入れております。しかしながら、表示ではそれが確認できません。

    キャンセル

  • 2020/10/18 22:06

    20, 30, 40, ・・・, 150 を離散値としてみて、ヒストグラムを表示させたいのでしょうか?
    ヒストグラムであれば、150が「一定の幅を持って」表示されるべきなのは理解します。
    私は、温度なので連続値であると解釈しまたした。連続値だとすると150は「右端」ですので、150のところの値は縦線未満として見えなくなってしまいます。

    キャンセル

  • 2020/10/19 03:32

    なるほど。離散値とみて、ヒストグラムを表示させたいというのが、toast-uzさんの言葉での私の要望のようです。

    キャンセル

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

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

関連した質問

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