トップページ⇒tech
152コメント50KB

Matplotlib でグラフプログラミング【Python】

■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん2009/10/24(土) 21:06:22

高機能で柔軟性があり、使い易く、数式表示がデフォルトで可能で、
美しいグラフが作成できる人気のグラフライブラリである Matplotlib
のスレッドです。


ギャラリー
http://matplotlib.sourceforge.net/gallery.html
http://matplotlib.sourceforge.net/users/screenshots.html

ユーザーズガイド
http://matplotlib.sourceforge.net/contents.html

ダウンロード
http://sourceforge.net/projects/matplotlib/files/matplotlib/matplotlib-0.99.1/
0002デフォルトの名無しさん2009/10/24(土) 21:07:25
このスレッドは天才チンパンジー「アイちゃん」が
言語訓練のために立てたものです。

アイと研究員とのやり取りに利用するスレッドなので、
関係者以外は書きこまないで下さい。

                  京都大学霊長類研究所
0003デフォルトの名無しさん2009/10/24(土) 21:07:30
http://pc12.2ch.net/test/read.cgi/unix/1164719425/708 より:


gnuplot vs matplotlib

2Dプロット系: "できること" で言えば、パイチャートなどがデフォルト
でプロットできる matplotlib の方が上。ちょっと混み入った設定や
アノテーションの加え方も matplotlib の方が楽。

3Dプロット系: 完成度と速度の面で gnuplot の勝ち。

プロットプログラムの再利用性: これは matplotlib の余裕勝ち。

見た目の綺麗さ: matplotlibではデフォでグラフにアンチエイリアスや
ヒンティングがかかっている。gnuplot だと wx 末端を使えば
matplotlibのそれとほとんど同様になるが、他の末端も
含めて評価するとmatplotlibに軍配が上がる。

数式: TeXのフォントがデフォルトで入っていて、インストール直後
LaTeXのインストール無しでも即効TeXと同じように書けるので matplotlib の勝ち。

結論としては、多種多様のプロットを生成したいのであれば、
matplotlibを使った方が長期的に見て効率が良さそう。特に python 使い
であれば、断然 matplotlib の方が有利。それと、gnuplotはあまり融通
が効かず、ちょっとした変更が時に面倒だったり、"できそうでできない"
ことが何気に結構あったりしたが、matplotlibなら、その心配は無さそう。

しかし、「簡単にsin(x)のプロットが作れればいいやー」程度のユーザーなら、
gnuplotでいいと思うし、そういった単純なグラフなら gnuplot の方が簡単。
0004デフォルトの名無しさん2009/10/24(土) 21:14:56
>>1乙。ありがとう!
0005デフォルトの名無しさん2009/10/24(土) 21:18:52
自分でスレ立てといて言うのもなんだが、
google、インデックス早っ!
もう既にこのスレが検索結果(結構上位)に引掛ってるし。

ttp://www.google.com/search?q=matplotlib&hl=ja&lr=lang_ja
0006デフォルトの名無しさん2009/10/24(土) 21:29:06
Debianでmatplotlibインストールしたら、バージョンが0.98.1だったぁ!

3次元プロットがギリできない〜
0007デフォルトの名無しさん2009/10/24(土) 21:29:49
三次元空間において、
曲面 z = f(x,y)と
曲面 z = g(x,y)
が交わる曲線をプロットする良い方法はあるかな?
0008デフォルトの名無しさん2009/10/24(土) 21:34:20
>>7
普通に同じaxesにfとgを加えればいいでしょ。
使ったことないけど、 ax.add_collection3d 見てみれば?
ttp://matplotlib.sourceforge.net/examples/mplot3d/polys3d_demo.html
0009デフォルトの名無しさん2009/10/25(日) 01:16:42
これ良くできてるけど、意外と知られてないよね。

gnuplotだけで頑張ってる人は一見の価値あり。
0010デフォルトの名無しさん2009/10/25(日) 09:58:48
gnuplotだと付属のヘルプで大体マスターできるけど、
matplotlibって何か簡潔で分かり易いのない?
0011デフォルトの名無しさん2009/10/25(日) 11:25:20
gnuplotで頑張ってきた人だけど、おもしろそうなので
ちょっと遊んでみます。
0012デフォルトの名無しさん2009/10/25(日) 21:55:41
gnuplotだと『splot 'data.dat' using 5:6:18』とすればプロットしてくれたけど、
こっちでは自分で全部書かないといけないのかな?
自分は↓こうしているんだけど、もっと簡単な方法ってある?

#!/usr/bin/env python
import re
myRegExp = re.compile('^\s*#')

column1 = 5
column2 = 6
column3 = 18

X = []
Y = []
Y = []

fname = open('data.dat', 'r')
for lines in fname:
  if myRegExp.match(lines):
    print '1 line ignored.'
  else:
    itemList = lines.split()
    X.append( float( itemList[ column1 ] ) )
    Y.append( float( itemList[ column2 ] ) )
    Z.append( float( itemList[ column3 ] ) )
0013デフォルトの名無しさん2009/10/26(月) 03:49:21
基本的にプロット用のライブラリだから、データの用意は自分で...って感じだけど、
テキスト形式を読みこむ関数はある(実際はmatplotlibのバックエンドのnumpyの機能)。

import pylab
data = pylab.loadtxt("data.dat", comments="#")
X, Y, Z = data[:,(5, 6, 8)]
0014デフォルトの名無しさん2009/10/26(月) 04:00:13
>>12
即席で書いたやつだが、こうすれば良し。

LL = [ line.strip().split() for line in file('data.dat') ]
columns = map(list, zip(*LL))

columns[4], columns[5], column[18]
0015デフォルトの名無しさん2009/10/27(火) 05:16:57

matplotlibはpythonのライブラリだけど、
コアな計算部分は numpy や PIL だから、
中身は大体Cの部分が実行されてるよね?

スピードはgnuplotと比べてどのぐらい違うんだろうか?
0016122009/10/27(火) 12:12:57
>>13-14
ありがとう!だいぶ短くなるんですね。

pythonを使うのが初めてなので、まず言語の勉強から始めなければ。
0017デフォルトの名無しさん2009/10/27(火) 17:19:13
>>13
NumPyでこんなスライスができるのは知らなかった
0018デフォルトの名無しさん2009/10/27(火) 22:45:51
これはまれにみる糞スレage
0019デフォルトの名無しさん2009/10/28(水) 23:16:55
python 使いで gnuplot の限界を感じた人が
使うってだけじゃなくて matplotlib 自体がお目あてで
来る人もいるんだねぇ

そういう人むけのチュートリアルとかが出来るといいね
0020デフォルトの名無しさん2009/11/01(日) 02:19:36
Matplotlibでgriddataを使っているヒトいますか?ちゃんと動きますか?

Debianにgriddataをソースからインストールしたのですが、
セグメントエラーを吐くようになってしまいました。
0021デフォルトの名無しさん2009/11/06(金) 09:06:05
>>1
ギャラリーにSICPで見かけたおねえさんがw
0022デフォルトの名無しさん2009/11/06(金) 12:15:54
lenaだよ
0023デフォルトの名無しさん2009/11/06(金) 18:15:36
pngファイルを読み込んでそれをそのままmatplotで表示って、
皆さんどうやってますか?
0024デフォルトの名無しさん2009/11/06(金) 18:20:08
>>23
うーん・・・。
背景にしたいの?
0025デフォルトの名無しさん2009/11/06(金) 18:29:56
>>24
たとえば入力した画像の明るさを読み取って、ある値以上の明るさの画素は白、他は黒にして出力
などといった画像処理につなげたいのですが…
0026デフォルトの名無しさん2009/11/06(金) 19:34:28
import pylab
import scipy.misc
m = scipy.misc.imread("grayscale.png")
pylab.imshow(m < 127)
pylab.savefig("output.png")

こんなんでどう?
0027デフォルトの名無しさん2009/11/06(金) 22:28:28
>>26
おお!ありがとうございます!
今ちょっと手元に無いので明日試してみます
0028デフォルトの名無しさん2009/11/08(日) 02:46:26
研究室の友人(C言語しか知らない)にPython勧めたら、その場に居合わせた他の研究室の奴に、
「なんでそんな聞いたこともない言語使うの?VBでいいじゃん」
って言われてしまった。理系だったらPythonのがいいと思うんだけどなぁ。。。
0029デフォルトの名無しさん2009/11/08(日) 02:56:20
聞いたことも無いのはお前が不勉強だからだってすぐに言えば良かったのに
0030デフォルトの名無しさん2009/11/08(日) 02:57:45
python しらねーの?w
情弱 プギャー w
でよかったんじゃね?
0031デフォルトの名無しさん2009/11/08(日) 03:27:41
>>28
両方使えるのがいい。
VBはExcel操作できるから、知っていれば意外と使える。
0032デフォルトの名無しさん2009/11/08(日) 10:57:21
excel操作なんか大抵の言語でできるよ(当然pythonも)
0033デフォルトの名無しさん2009/11/08(日) 11:30:12
VBで出来ないことまで出来るよねw
0034デフォルトの名無しさん2009/11/08(日) 11:50:40
Fランク大学は放置しる
0035デフォルトの名無しさん2009/11/08(日) 11:51:53
Excelに10000行くらいのデータを貼り付けてプロットしてる人を見るとアホかと思う
0036デフォルトの名無しさん2009/11/12(木) 18:17:20
>>35
瞬時に結果表示できるなら、それで問題ない。
0037デフォルトの名無しさん2009/11/12(木) 19:11:15
データをとるたびにそれを繰り返してるってことでしょ。
0038デフォルトの名無しさん2009/11/12(木) 19:23:59
Excelって小学生並みのアホなグラフしか描けないのに何で使うのかな
0039デフォルトの名無しさん2009/11/12(木) 20:27:31
商品として優れてるからだろ。
いい製品なんて大抵ゴミだし。
0040デフォルトの名無しさん2009/11/12(木) 21:44:39
>>37
まぁ、コピペの部分ぐらいはマクロ使うかな。
0041デフォルトの名無しさん2009/11/16(月) 13:26:05
あんま、話題ないんね。
なんか良い入門ないかな。
28だけど、研究室に布教しようと努力中。。
0042デフォルトの名無しさん2009/11/16(月) 13:30:01
最近使いはじめたよ。
便利でキレイでいいね。
gnuplot みたいな手軽さは無いけど。

サクっとグラフを見たいときは gnuplot で
貼り付けたい時なんかは、もう Matplotlib だね。
0043デフォルトの名無しさん2009/11/17(火) 07:53:32
入門書と言えばBeginning Python Visualization位しか持ってないんだけど、
他にある?
この本でもChapter 6で扱ってるだけ。
あんまり詳しくない。
もっと網羅的な奴ないかな。
0044デフォルトの名無しさん2009/12/01(火) 11:22:44
matplotlib でplot 出力するのは
たいてい png な説明がでてるけど
ps/eps にしてもなんか"画像(ラスタイメージ)"になってるね?

細かいデータを打って拡大すると
ticとかが同じ線幅のはずなのに太かったり細かったりって
なっちゃうのはどうにかならないのかな?
(savefig で dpi=200 とかしてもダメみたい…)

あと,"真ん中あたりに小さくグラフが描画される" けど
もっと出力サイズ一杯に描画する方法ないんでしょうか?
(.ps だと a4目一杯)

savefig にはそういうオプションなさそうだけど…
0045デフォルトの名無しさん2009/12/01(火) 12:06:02
手元の環境では普通にベクトルデータで生成されるよ。svgはどう?
0046デフォルトの名無しさん2009/12/01(火) 12:59:46
あれ,そうですか.
確かに ps データのなかを改めて見ると
そんな感じに見えますね… orz

うーんなんでだろう.
gv で見るときの問題なのかなぁ...
0047デフォルトの名無しさん2010/01/14(木) 09:27:49
matplotlib で line なグラフを描くとき
線種って結構限られていますよね?

実線('-') 破線('-.') 点線(':')

gnuplot だと(指定方法が色とごっちゃになってる問題はありますが)
10種類以上くらいの破線種類なんかがあったりしますが,
破線の種類なんかをもっと沢山用意してプロットするのは
できないんでしょうか?
-・-・- とか −−・−−・−−・−− とか

line の属性をほいほいっと弄れば出来るとか…
0048デフォルトの名無しさん2010/01/14(木) 16:27:03
(line1, ) = plot([1, 2, 3], [3, 4, 5])
line1.set_dashes((2, 1, 2, 1, 2, 5)) # インクのon/offを交互に指定
savefig("asdf.png")

もっといい方法があるかもしれないですが
0049デフォルトの名無しさん2010/01/15(金) 17:07:16
くだすれPython(超初心者用)から来ました。

オライリーの「集合知プログラミング」にnumpyを使った非負値行列因子分解と言うのがあるのですが
(ソースは英語版ですが http://examples.oreilly.com/9780596529321/ から落とせます)、
これを乗法的更新ルールというので行うと
http://pc12.2ch.net/test/read.cgi/tech/1256350845/785 のようにすぐにNaNになってしまいます
(10回に一回のprint costを毎回print costするようにすると3周目でnanになるようです)。

このスレの方ならnumpyのことをご存知で無いかと思い質問したしだいです。よろしくお願いします。
0050デフォルトの名無しさん2010/01/15(金) 17:15:54
ソースコードは10章のnmf.py(データ取得部はnewsfeatures.py)を見てください。
0051デフォルトの名無しさん2010/01/15(金) 18:29:38
numpyのバージョンは?
0052デフォルトの名無しさん2010/01/16(土) 00:53:45
>>51
numpy-1.4.0-win32-superpack-python2.6
をインストールしました。
0053デフォルトの名無しさん2010/01/16(土) 13:15:11
>>52
nnmf.pyを落としてきて,ちょっと見てみたけど
29行目のループ終了条件

if cost == 0: break

を

if cost < 1e-8: break

とかにすればいいんじゃない?
0054デフォルトの名無しさん2010/01/17(日) 06:13:27
>>53
レスありがとうございます。ただcosr==0も終了条件なのですが、
繰り返し回数の最大はiterで引数で与えているので、あまり関係ないかと。
3回繰り返すとnanが出てくるようです。
0055デフォルトの名無しさん2010/01/17(日) 08:53:54
いやこれ繰り返し回数じゃなくて
二乗距離でしょ?
かえてやってみたのかと
まぁいいや
0056デフォルトの名無しさん2010/01/17(日) 10:54:49
>>55
すみませんおっしゃる事が良くわかりません。
これってcostのことですか?iterのことですか?それ以外ですか?
>>53 の変更はやってみたのですが、結果は変わりませんでしたが、
そういう意味ではないのでしょうか?
0057デフォルトの名無しさん2010/01/17(日) 13:51:06
>>56
そっか、変えてやってみてもダメでしたか。すまんかった。
書籍持ってなくてnnmf.pyをちら見しただけでコメントしたんだけど
ウェイト更新とかかいてあるんで、ウェイト自体が微小になってしまうことを避ける
条件文なのかと思ってた。

いま改めてみると、ランダムに生成したWとHの積と
与えられたVが等しくなるかどうかのチェックなんだね。

ちょっとググってみて、非不値行列行列因子分解の乗法更新式に関する文献↓を流し読みしてみたけど
http://www.brl.ntt.co.jp/people/kameoka/publications/Kameoka2009ASJ09.pdf
ソースのアルゴリズムは文献のLS法どおりなので問題なさそう。

Pythonには組み込み関数の"iter"があるのに、引数名に"iter"を使っているとか
最初に"h"を生成しているリスト内包表記が、内側と外側で同じループ変数名を使っているとか
スタイルで気になるところはあるけど。

"h"と"w"を更新している箇所は要素同士の計算なので
h, hn, hdやw, wn, wdのある要素[i,j]を更新前後で表示してみれば
NaNになる理由がわかると思うよ。
0058デフォルトの名無しさん2010/01/18(月) 17:33:11
>>57
print デバッグですね。分りました。
0059デフォルトの名無しさん2010/01/19(火) 18:16:55
>>51 factorize(..) の引数にはどんなデータを与えているんだ。

俺の環境:Enthought python 2.5 では下のように小さな値の行列でも NaN は出てこない

python temp.py
from numpy import *

def difcost(a,b):
  dif=0
  for i in range(shape(a)[0]):
    for j in range(shape(a)[1]):
      # Euclidean Distance
      dif+=pow(a[i,j]-b[i,j],2)
  return dif

def factorize(v,pc=10,iter=50):
  ic=shape(v)[0]
  fc=shape(v)[1]

  # Initialize the weight and feature matrices with random values
  w=matrix([[random.random() for j in range(pc)] for i in range(ic)])
  h=matrix([[random.random() for i in range(fc)] for i in range(pc)])

  # Perform operation a maximum of iter times

続く
0060デフォルトの名無しさん2010/01/19(火) 18:17:35
  for i in range(iter):
    wh=w*h
   
    # Calculate the current difference
    cost=difcost(v,wh)
   
    if i%10==0: print cost
   
    # Terminate if the matrix has been fully factorized
    if cost==0: break
   
    # Update feature matrix
    hn=(transpose(w)*v)
    hd=(transpose(w)*w*h)
 
    h=matrix(array(h)*array(hn)/array(hd))

    # Update weights matrix
    wn=(v*transpose(h))
    wd=(w*h*transpose(h))

    w=matrix(array(w)*array(wn)/array(wd))
   
  return w,h


import scipy as sc
print factorize(sc.randn(10,10)/1000)
0061デフォルトの名無しさん2010/01/29(金) 09:36:08
& nbsp;には最後にセミコロン付ける必要があるぞw
0062デフォルトの名無しさん2010/02/12(金) 09:12:09
これはひどい
0063デフォルトの名無しさん2010/02/17(水) 12:53:49
そういえば統計方面の R言語のスレってどこかの板にある?

あるなら関連スレとして参照しておくくらいは欲しい
0064デフォルトの名無しさん2010/02/17(水) 14:00:06
ここにあったのを見た覚えがある
0065デフォルトの名無しさん2010/02/17(水) 14:03:02
ここにRみたいだ
http://science6.2ch.net/test/read.cgi/math/1224142396/
0066デフォルトの名無しさん2010/02/22(月) 15:15:45
wxバックエンドの場合、
どうやってバックグラウンドカラーを変更できるの
(デフォはグレイ)?
0067デフォルトの名無しさん2010/02/27(土) 10:59:30
質問させてください。

http://d.hatena.ne.jp/noazoh/20091023/1256293199
を参考にローソク足チャートを描きました。また
http://matplotlib.sourcearchive.com/documentation/0.90.1/namespacematplotlib_1_1finance_828611c3468bfc1d1b336155704f4b12.html#828611c3468bfc1d1b336155704f4b12
を見て、quotesにvolumeの値を加えれば、volume_overlay3 がそのまま使えると思ったのですが、
実際に描画してみると、全部同じ高さの棒になってしまいます(上に突き抜けてる?)。

どなたかこのvolume_overlay3の使い方が分かる方いらっしゃいませんでしょうか。
0068672010/02/28(日) 02:31:01
コード書きましたようはこんな感じでやりたいんですが

def downyahoo():
 date1 = datetranslate("2009/02/18")
 date2 = datetranslate(time.strftime("%Y/%m/%d"))

 quotes=[]
 quotes = matplotlib.finance.quotes_historical_yahoo(^GSPC, date1, date2, False, True, None)
 
 fig = figure()
 ax = fig.add_subplot(111)
 matplotlib.finance.candlestick(ax, quotes, 0.2, 'r', 'k')
 matplotlib.finance.volume_overlay3(ax, quotes)
 ax.xaxis_date()
 fig.autofmt_xdate()
 show()

def datetranslate(date):
 d = date.split('/')
 startDate = datetime.date(int(d[0]),int(d[1]),int(d[2]))
 return startDate

def main():
 downyahoo()
0069デフォルトの名無しさん2010/02/28(日) 17:50:19
>>68

volume_overlay3はメンテされてないみたい
ほとんど同じ内容の "candlestick2" を参考にして
volume_overlay3を書き直してみたよ

http://codepad.org/vvKwA6wO
0070デフォルトの名無しさん2010/02/28(日) 18:10:55
users-mlで同じ症状を訴えている人がいるんだけど
思いっきり無視されてるんで、メンテされていないのかなと

http://www.mail-archive.com/[email protected]/msg12590.html
0071デフォルトの名無しさん2010/02/28(日) 23:28:07
>>69
ありがとうございます。本家よりこっちの方が進んでいると言うことですかw
0072デフォルトの名無しさん2010/03/03(水) 17:27:23
date2numした数字をdate2numで元に戻したいのですが、

dates=[matplotlib.dates.num2date(d[i][0]).date() for i in range(len(d)) if d[i][5]==d[i][5]]

とやると、printしたときにdatesの中に

[datetime.date(2010, 2, 25), … datetime.date(2010, 3, 1)]

とか入るんですが、これを

['2010, 2, 25',…'2010, 3, 1']

のように短くしたいのですが、どうやればいいか教えて下さい。
0073デフォルトの名無しさん2010/03/03(水) 17:30:16
num2dateで元に戻したいでした…
0074デフォルトの名無しさん2010/03/03(水) 18:12:13
strftime
0075デフォルトの名無しさん2010/03/03(水) 21:31:10
matplotlibスレがあるなんて思わなかった!
0076デフォルトの名無しさん2010/03/03(水) 22:07:11
この本買ったけど、基礎的な事を学ぶには良かったよ。

Matplotlib for Python Developers
http://www.amazon.co.jp/Matplotlib-Python-Developers-Sandro-Tosi/dp/1847197906

Qt4のWidgetとする方法とか、参考になった。
0077デフォルトの名無しさん2010/03/04(木) 07:52:49
2次元グラフ作成ソフト「Ngraph for Windows」がフリーソフト化
ttp://www.forest.impress.co.jp/docs/news/20100303_352467.html
0078デフォルトの名無しさん2010/03/04(木) 11:37:02
>>74
strftime(%Y/%m/%d)で出来ました。ありがとう。
0079デフォルトの名無しさん2010/03/04(木) 12:35:12
本なんてあるんだ。
0080デフォルトの名無しさん2010/03/04(木) 15:30:33
とりあえずお布施で買っておくか
0081Zeno ◆5nZQbNmQPs 2010/03/04(木) 16:28:01

吾も情報を落とすことにするか。

>>76
早速 "Matplotlib for Python Developers" を手に入れてみた(英語版)。
Web アプリや Matplotlib のグラフィックインターフェースへの組み込み方の
説明に力を入れてはいるものの、肝心の matplotlib の活用/応用系のチャプターは
短絡的すぎる。実際、そもそも作者の言う "advanced" は "basic" レベル止まりで、
全体的に内容が浅くなっている件は否定できない。

もし興味があれば
ttp://www.amazon.com/Beginning-Python-Visualization-Transformation-Professionals/dp/1430218436/
も平行して読んでみると良いかもしれん (こちらも初歩的なプロティングしか触れていない) が、
個人的に、 ttp://matplotlib.sourceforge.net/gallery.html を参照するのが
上達の一番の近道だと思っている。
0082デフォルトの名無しさん2010/03/04(木) 21:05:57
>>81
確かに内容浅いね。広く浅くという感じか。
俺は結局公式ドキュメントやexamples、場合によっては本体のソースを見ながらやってる。
0083デフォルトの名無しさん2010/03/05(金) 20:33:20
意外に欲しいメソッドが無かったりするね。
例えば、tick labelのフォントを変えたりとか。
figureの保持するlegendを取得するメソッドも無い。axesにはあるのに。
まぁ、俺が知らないだけかも知れないけど。
0084デフォルトの名無しさん2010/03/06(土) 00:00:21
結構どんどん追加されていってるからね

伊達にv1.0未満じゃないというか…

オンラインドキュメント調べて
あれこの機能が動作しないと思ったら
最新版じゃないと駄目だったなんてよく起きるし
0085デフォルトの名無しさん2010/03/06(土) 22:54:12
あるねぇ。
リリース版では「class IndexFormatter」も「Formatter」を継承してなくて、動くはずのものがエラーになった。
0086デフォルトの名無しさん2010/05/09(日) 21:12:20
gnuplotとGhostscript と GSviewあたりを最新版にしたら、うまくPNGとか作成されなくなった
色々と設定し直せばいいんだが、面倒になったきた
そこで、pythonで別のグラフ作成方法ないかと思って、これにたどり着いたが・・・
python3使ってるわ・・・残念


>>25
これ他方法でもできないかな
0087デフォルトの名無しさん2010/05/12(水) 11:06:47
試していないから間違っていたら聞き流して欲しいんだけど,awkで
http://gauc.no-ip.org/awk-users-jp/blis.cgi/DoukakuAWK_114
http://gauc.no-ip.org/awk-users-jp/blis.cgi/DoukakuAWK_254
を参考にしたら出来そうな気がする.
0088デフォルトの名無しさん2010/05/16(日) 14:43:45
Pythonで学ぶ萌え萌え統計学 ー 正規分布からマルコフ連鎖モンテカルロまで 2300円
0089デフォルトの名無しさん2010/12/24(金) 15:50:18
そういやこのスレ忘れ去られているんでネタ投下

scipy conference japan やろうぜって話が
python workshop で出てましたねー
この界隈に興味ある人はいるのかな?
0090デフォルトの名無しさん2010/12/27(月) 04:25:56
hage
0091デフォルトの名無しさん2011/01/30(日) 02:38:10
numpyとmatplotlibがPython3に対応したようですね。
0092デフォルトの名無しさん2011/01/30(日) 02:39:15
まじで
0093デフォルトの名無しさん2011/01/30(日) 05:36:09
PILはまだか
0094デフォルトの名無しさん2011/01/30(日) 18:03:55
Python(x,y)とMacPortsでも対応したら本気出す>3
0095デフォルトの名無しさん2011/02/01(火) 00:10:14
>>91
対応してねーじゃねぇか、ぷんすか。
0096デフォルトの名無しさん2011/02/24(木) 18:01:22.07
質問させてください。
matplotlib で グラフのフレームの色を変更することは可能でしょうか?

・やりたいことは背景色を黒、メモリなどを白としたい。
1.Figureに facecolor='black'を指定し、set_axis_bgcolorで黒を
  指定することで背景を黒にできました。
2.set_xticklabels, set_yticklabelsで白を指定しメモリの数字は
 白くできました。

しかし、枠およびメモリ自体を白くする方法がわかりません。
わかる方おられましたら、お願いします。
0097デフォルトの名無しさん2011/03/03(木) 22:38:15.74
>>91
numpyは1.5で対応したみたいですね。
matplotlibは、2.7までのようです。
0098デフォルトの名無しさん2011/03/04(金) 16:36:25.26
matplotlibで関数の動きって表現できる?
例えばy = sinx * nとして
nが1から10まで動くとどうなるかっていうのを動的に表現したいんだけど
0099デフォルトの名無しさん2011/03/04(金) 19:24:19.40
できますが
nが1から10まで動くと
振幅が変わるだけではないでしょうか
0100デフォルトの名無しさん2011/03/05(土) 00:47:43.99
ごめん 例えが悪かったわ
結局は関数を更新していってその動きを見たいんだけど
そういうのはmatplotlibに組み込まれてる機能でできるの?
よければ参考になるページとか教えて貰いたいんだが
0101デフォルトの名無しさん2011/03/05(土) 13:52:37.64
ふう
0102デフォルトの名無しさん2011/03/06(日) 19:02:55.04
なんのために汎用言語であるpythonを使っているんだ?
for i in range(10):
draw_hoge()
time.sleep(1)
0103962011/03/16(水) 16:51:06.66
自己解決
#個々のグラフ領域の枠を白くする
ax = fig.add_subplot(1, 1, 1)
for child in ax.get_children():
if isinstance(child, matplotlib.spines.Spine):
child.set_color('white')

#X軸Y軸のメモリの色を白くする
#メジャーとマイナーは一度には変更できないみたい
#なので2回に分けて行う
for liney in ax.yaxis.get_ticklines(minor=True):
liney.set_color('white')
for liney in ax.yaxis.get_ticklines(minor=False):
liney.set_color('white')

0104デフォルトの名無しさん2011/03/22(火) 13:24:03.28
線の色指定をするときに、任意のカラーマップから色を指定することって可能?
Gray Scaleの場合は、color = '0.75'のようにすればいいけど、
そこを別のカラーマップで0-1の間の数値を入れて指定できるようにしたいです。
0105デフォルトの名無しさん2011/03/27(日) 13:27:24.41
二項分布とか正規分布とかを正確に表示したいときは自分で関数つくんないとダメ?
0106デフォルトの名無しさん2011/03/27(日) 14:11:39.66
>>105
scipy.stats以下にたくさん分布関数がある
0107デフォルトの名無しさん2011/04/25(月) 15:02:35.40
Python3に対応したから、初めて使ってみたらmatplotlibすげー
軽くgnuplotみたいに画面で直接拡大とかも出来るんだな
0108デフォルトの名無しさん2011/04/25(月) 15:09:27.89
お、 Python3 対応きたんだあ。
0109デフォルトの名無しさん2011/04/25(月) 15:35:17.49
>>100,102
同じところで詰まった…
スレッドとか必要なのかな
Tkでやってた時もよくわからなくてpygameに逃げた
0110デフォルトの名無しさん2011/04/25(月) 16:39:03.66
出来た出来た
適当にバラまいたノードがランダムに動くプログラム
でも重すぎるのか途中で処理が止まる…
本当はノード数を1000以上にしたいのに
それと、きちんと再描画されてなく、毎回、点が増えていく…

pylab.ion()

NodeN=int(input("ノード数入力:"))
lisNode =[0]*NodeN
for n in range(NodeN):
 lisNode[n]=random.randint(-200,200),random.randint(-200,200)

for i in range(100):
 pylab.plot([x for x in lisNode] ,[y for y in lisNode],"*")
 pylab.draw()
 lisNode=[[lisNode[n][0]+random.randint(-1,1),lisNode[n][1]+random.randint(-1,1)] for i in range(NodeN)]
0111デフォルトの名無しさん2011/04/25(月) 18:39:48.95
上書きがデフォルトなので、毎回図をクリアする必要がある。
あと、そのプログラムが望む動作をするかどうか、落ち着いてコードを読み直してみてね
0112デフォルトの名無しさん2011/04/25(月) 20:06:07.41
pylab.plot([x[0] for x in ………
ですね。
手打ちで書いてたので間違えました。
とりあえず、処理が止まらなければ望む動きはしています。
上書きを除いて

画面の図、クリアを調べてみます
0113デフォルトの名無しさん2011/05/12(木) 11:49:54.78
外枠を太くしてgnuplotと体裁が同じになるようにしたい
何か方法ないのかな
0114デフォルトの名無しさん2011/05/12(木) 22:55:32.05
>>110
こんな感じでやると、描写中の図を動かそうとすると止まっちゃうんだけどPCスペック低いからかな?
あと描写スピードが遅いのと最後に自動で図が消えちゃうのってどうすればいい?
0115デフォルトの名無しさん2011/05/12(木) 23:37:13.09
>>114
>描写中の図を動かそうとすると止まっちゃう
知らんがな
>描写スピードが遅い
プログラムが悪いか、データ点が多すぎるか
>最後に自動で図が消えちゃう
たぶんあなたが消してる

>>110は真似しない方が・・・
0116デフォルトの名無しさん2011/05/13(金) 00:10:12.82
>>114
>110のプログラムは試してないからよくは分からないけど。。

windowを動かすと描画が止まる、ってことなら、うちでもそうなる。おそらくバックエンド部分のバグじゃないかと疑ってる。最期に消えちゃうのは、プログラム最後にshow()を入れると止められる。
0117デフォルトの名無しさん2011/05/13(金) 01:50:14.28
>>115
>>116
ありがと やっぱ動かすと止まっちゃうのね

import matplotlib.pyplot as plt
import random

plt.ion()

x = [0]*100
y = [0]*100

random.seed();
for i in range(100):
x[i] = random.randint(-100, 100)
y[i] = random.randint(-100, 100)

plt.axis([-110, 110, -110, 110])
plt.grid(True)

for i in range(len(x)):
plt.plot(x[i], y[i] , 'kD')
plt.draw()

plt.show()

これで100個プロット打つのに9秒くらいかかるんだけどこんなものかな
あとshow入れても最後にはやっぱ消えちゃうわ
0118デフォルトの名無しさん2011/05/15(日) 11:21:02.45
>>117
show()の前にioff()を入れると画面は消えない。
0119デフォルトの名無しさん2011/05/16(月) 15:50:46.45
>>117
plt.ion()はコメントアウトして、
importの下からplt.show()まで時間計ったら0.15秒

for i in range(len(x)):
 plt.plot(x[i] (ry
の部分を、for文やめて
plt.plt(x,y,'kD')
に変更して、0.1秒


Windows7 SP1 32bit
i3 550 3.20GHz
メモリ 2G
Python3.1
0120デフォルトの名無しさん2011/05/16(月) 21:24:39.24
drawが時間喰ってる@Python2.7
0121デフォルトの名無しさん2011/05/16(月) 21:29:50.63
>>119
draw()もコメントアウトしてある
0122デフォルトの名無しさん2011/05/16(月) 22:57:10.63
>>121
元の動機が図をダイナミックに更新したい、ということだからそれだと意味なくない?
0123デフォルトの名無しさん2011/05/17(火) 17:03:36.09
Gridのmajorとminorの線種をそれぞれ別で指定する方法ありますか?
01241232011/05/17(火) 19:46:58.24
解決しました
0125デフォルトの名無しさん2011/05/23(月) 21:07:48.89
windowsだとtkですら、途中で応答無くなるのなんとかならないかな
0126デフォルトの名無しさん2011/05/24(火) 21:56:18.64
横軸の線の先を矢印にすることって出来ませんか?
0127デフォルトの名無しさん2011/06/13(月) 00:05:31.79
eps化してtexに貼りたいんですが、texのフォントを(Computer Modern)を使うことは可能ですか?
またpdf化の際に埋め込めますか?
0128デフォルトの名無しさん2011/06/28(火) 17:44:57.28
ttp://matplotlib.sourceforge.net/users/image_tutorial.html
の例(をsubplotsを使って書き直そうとしたもの)なんですが、カラーバーの位置がおかしいです。
どうしてax[0].figureで欲しいものが取れないのでしょうか?

=========================
import matplotlib.pyplot as plt
import matplotlib.image as mpimg

img = mpimg.imread('stinkbug.png')

fig,ax = plt.subplots(ncols=2)
ax[0].set_aspect(1)
img1 = ax[0].imshow(img[:,:,0])
ax[0].figure.colorbar(img1,orientation ='horizontal')
plt.show()
01291282011/06/28(火) 18:03:15.11
ごり押しでいけた。
plt.subplot(ax[0]).figure.colorbar(img1,orientation ='horizontal')

・・・明らかに設計ミスだと思うんだけど。
0130デフォルトの名無しさん2011/06/29(水) 21:34:04.98
subplot 121, 122 で分けてるのですが、
左右の余白大きすて調整したいです。
subplots_adjust では埋められません。
ほかに何かいい方法はありませんか。
0131デフォルトの名無しさん2011/06/30(木) 12:08:40.08
type 1フォントをベクトルデータでepsに埋め込めませんか?
0132天使 ◆uL5esZLBSE 2011/07/02(土) 18:22:05.67
> 元の動機が図をダイナミックに更新したい、ということだからそれだと意味なくない?
ハッアァァァァァァァァァァァァァァアァ??????????????
なるほどゴミなわけだ
ゴミじゃねーか
死ねよゴミ
0133デフォルトの名無しさん2011/07/06(水) 01:06:18.37
>>130
matplotlibrcをいじってみるとかはどう?
0134デフォルトの名無しさん2011/09/04(日) 21:06:24.11
これCUIオンリーの環境じゃ動かないの?
/var/lib/python-support/python2.5/gtk-2.0/gtk/__init__.py:72: GtkWarning: could not open display
からずらずらエラーで動かないよ
plt.savefig("a.png")だけでもダメなのかな
show()であかんのは納得するんだけど。。
0135デフォルトの名無しさん2011/09/04(日) 21:57:56.62
うちでは動くよ
pngに変換するためのバックエンド(Agg)がインスコされてないんじゃない?
http://matplotlib.sourceforge.net/faq/installing_faq.html#backends
0136デフォルトの名無しさん2011/09/05(月) 09:48:22.44
>134
同じく…
あきらめて X11 環境下でのみ作業してる

たまに ssh login してグラフの作り直しだけしたい
みたいなときに不便なので対策があったら知りたい
0137デフォルトの名無しさん2011/09/06(火) 19:10:52.18
>>134
>>136
だめだね
import matplotlib.pyplot as plt
これ一行だけで落ちるからX必須みたいだ
どうにかならんか調べてみます
0138デフォルトの名無しさん2011/09/06(火) 19:33:23.38
いやあの>>135をですね

import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
plt.plot([1,2,3,4])
plt.ylabel('some numbers')
plt.savefig('somenum.png')
0139デフォルトの名無しさん2011/09/07(水) 09:18:37.63
>>138
すいません135の意味がよくわからなかったのです。
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
ありがとうございます。 解決しました
0140デフォルトの名無しさん2011/09/08(木) 23:15:48.01
Numpy/Scipy,Matplotlibのチュートリアル的なサイトって
http://www.ike-dyn.ritsumei.ac.jp/~uchida/scipy-lecture-notes/index.html
実際ここくらいですか?
0141デフォルトの名無しさん2011/09/09(金) 16:08:40.70
うへ,なにこれ すげーボリューム
0142デフォルトの名無しさん2011/09/13(火) 08:09:30.74
numplotlib初心者です。
CentOS5のPython2.6上で使いたいのですが、
numpyインストール後matplotlib入れようと
するとド頭の環境チェックでnumpy: noとか
出ちゃってインストールできません。
numpyは1.6.1、matplotlibは1.0.1のtar拾ってきて
どちらもsetup.py installで入れようとしています。
numpy.arangeとかの結果は返ってくるので
numpyは一応入ったようなのですが、
こういう場合の対処方法ありませんか?
01431422011/09/13(火) 11:17:01.71
自己解決しました。
setup.py実行時に実行フォルダ内のゴミを
拾ってしまっていたようで、一旦cleanして
再度installしたところ無事に入りました。
01441422011/09/14(水) 23:24:10.47
白黒印刷用の棒グラフを作りたいのですが、
ハッチング等の白黒パターンで棒を塗り潰す
ことって可能でしょうか?
0145デフォルトの名無しさん2011/10/28(金) 12:52:21.33
matplotlibはデフォルトで2つのデータ点を直線で結びますが、
一部分だけ線を引かないようにすることは出来ますか。

gnuplotなら、渡すデータを
1 1
2 2

3 3
4 4
のように空行を挟んどけば、(2,2)と(3,3)はwith lineの対象外になるのですが
0146デフォルトの名無しさん2011/10/30(日) 09:51:54.81
>>145
matlabだとnanなら線が途切れるな
0147デフォルトの名無しさん2011/11/08(火) 11:52:07.56
去年辺りから、海外ではNumpy,Scipy,matplotlibを既存の科学用DSLの代替にしよう、
って流れが盛んですね。
米アマゾンでも、けっこうな勢いで解説書が出てるし。
既に全部、非公式ながらPython3にも対応して、前途も有望そう。

このスレも、この辺の統合スレにしたほうが盛り上がりそうですけど、
2年間で140レスか・・・
Rのスレは数学板にしかないし、ム板では盛り上がらないのかな?
0148デフォルトの名無しさん2011/11/09(水) 13:53:06.66
Rは昔ここにもあった
0149デフォルトの名無しさん2011/11/09(水) 14:35:55.50
ム板なんてITドカタしかいないからな
まともに数学の話出来る訳がない
0150デフォルトの名無しさん2011/11/11(金) 21:13:14.56
Rやmaximaあってもいいけどな。
>>147
そういった本ほんと盛んだよな。でももっとやすうなってくれんかな
と思ってしまう。Pythonは科学関連のライブラリはかなりほうふだから
基盤として使えるようにしておくといいよね。
p3 に対応したのか。
0151デフォルトの名無しさん2011/11/19(土) 12:41:07.85
wekaとscipyどっちがオヌヌメ?
0152デフォルトの名無しさん2011/12/06(火) 22:27:37.37
scipyかな
■ このスレッドは過去ログ倉庫に格納されています