過去記事 coffee-guhaw.hateblo.jp loadtxtでcsvから指定した列のデータを読み込む サンプルコード import numpy as np fname="hoge.csv" a = np.loadtxt(fname, delimiter=",", dtype="unicode", skiprows=1, usecols=[0, 3, 4]) 説明 今回は、いつもやって…

pandas.DataFrameの作り方あれこれ #Python - Qiitaここに色々と載ってる個人的にやりやすいな・・・と思ったのは、 import pandas as pd d = {"x": [1, 3], "y": [2, 4]} pd.DataFrame(d) の方法pandas.DataFrame, Seriesを辞書に変換（to_dict） | note.nk…

こんな感じのテキストデータがあって、これを読み込む 1 21 2 22 3 23 4 24 5 25 2列目のみ解析に使う、2列目のみ使うので1列目は捨てる import numpy as np fname="hoge.txt" a = (np.loadtxt(fname, unpack=True))[1] unpackのオプションがなんで必要かは…

python - How to convert an array of strings to an array of floats in numpy? - Stack Overflow サンプルコード import numpy as np fname="hoge.txt" a = np.loadtxt(fname, unpack=True, dtype=str) y = a["x"].astype(np.float) >|python| hoge.txtと…

import numpy as np fname="hoge.txt" a = np.loadtxt(fname, delimiter=",", unpack=True, comments="%") comments="%"でOK ランキング参加中です ↓クリックしていただけると嬉しいです〜ランキング参加中Python ランキング参加中HTML ランキング参加中Mac

hogeというarrayがあるとする このarrayの要素が1以上かつ10以下のものだけを解析に使いたいとすると・・・hoge[(1 <= hoge) & (hoge <= 10)]&の代わりにandではダメ 条件式を()で囲わないとダメ ランキング参加中です ↓クリックしていただけると嬉しいです…

実際使うことはなかったけど、今後必要な気もするのでメモ2つのarrayがある、例えばxとyとするxはデータのindex相当のもので、順番がぐちゃぐちゃになってる yはデータの値 xをソートするときに、yの順番も同時に同じようにソートしてほしい■ 参考 : argsort…

■ 参考 : Numpyで特定のindex以外を抽出する方法■ 参考 : Numpy で特定のインデックス以外の要素を取り出すarrというndarrayに条件を掛けて、条件を満たす成分をarrから抜いて、別のarrayを作成する 例えば、arrの中身がx0より大きいものだけ取り除きたいと…

ndarrayに文字列のリストが詰まっている その文字列の中にあるキーフレーズが含まれているかどうかを検索して、そのindexを取得して、 その文字列を取り出し、必要な部分のみ切り出したい最初は1度リストに変換してから、inを使って判定してたけど ndarrayの…

■ 過去記事 : matplotlibの凡例（legend）レイアウト関連メモmatplotlibでプロットするときにlabel="aaa"と書いて、ax.legend() みたいにlegend機能をオンにすると 最初から図にlegendを付けられる# ↑のリンクからお借りしたimport matplotlib.pyplot as plt…

ndarrayなことを活用しつつ、最大限高速で処理を終わらせたいndarrayを区間ごとに分割するのはnp.splitを使えばいい 割り切れない場合は、numpy.array_split()を使うとndarrayをできるだけ等分割で分割できる numpy.split()は等分割■ 参考 : NumPy配列ndarra…

IIRフィルターについてはscipyのドキュメントが詳しいのでそちらを一読がおすすめ■ 参考 : Signal Processing with SciPy: Linear Filters(Warren Weckesser)IIRフィルターには3つの表し方の形式がある ・伝達関数 ・zero, pole, gain(ZPK方式) ・second ord…

import random import numpy as np import matplotlib.pyplot as plt a = [random.random() for i in range(100)] val1, base1 = np.histogram(a, bins=100, range=(0.4, 0.8)) みたいな感じでヒストグラムを作る val1とbase1にはヒストグラムのあるビンのサ…

今のところやり方が見つかっていないarrayじゃなくてリストだと部分内包表記で書くことができる[c_match for c_match in c_list if "hoge" in c_match] [c_match for c_match in c_list if "hoge" == c_match]とかもしリストの要素数が多い場合は、numpyのar…

普通の累積分布はy軸は0から1に向かって増えていく 今回は1から0に向かって減っていく累積分布を書く あと、縦軸を0~1の範囲に収まるように規格化する必要がある import matplotlib.pyplot as plt import numpy as np val1, base1 = np.histogram(input1, bi…

sin波の位相と振幅を調べる必要があった 振幅はsin波のmaxを調べればいい →np.argmaxを使えばOK位相はzero crossする時刻から割り出せばよいのかな？と思った ひとまず、zero crossするindexを割り出す方法を調べてみた■ 参考 : Efficiently detect sign-cha…

numpy.argmax numpyのarrayの中から最大値になってるindexを返してくれる使い方例 : datasetというarrayのうち、i_beg からi_endの要素の中から最大値を見つけたい 返ってくるindexはi_begから何番目かなので、最後にそれにi_begを足す i_max = numpy.argmax…

dataというnumpy arrayがあったとして、それが2次元リストか1次元リストか知りたい # 1次元だけの空のarray data = [] # 1次元のarray data = [1, 2, 3] # 2次元のarray data = [[1, 2, 3], [4, 5, 6]] これらを区別したい - そのために len(data) data.shap…

まずはcsvから日付の部分を読み込む 日付データは 2022/03/28 というフォーマットになっている これをdatetime型で読み込めば良さそうだけど、他の部分が実数で読み込まないといけないのでこういう実装にしてるデータを読み込んで文字列にする ↓ / を - に置…