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

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

ただいまの
回答率

88.34%

pythonでのテキストファイルの指定範囲行の読み込み

解決済

回答 1

投稿 編集

  • 評価
  • クリップ 0
  • VIEW 1,849

kamikura

score 16

テキストファイルの指定行から指定行までの読み込み

3次元グラフを製作するためのテキストファイルの読み込みを行っています。
やりたいこととしては、テキストファイルを指定行数(ある行からある行まで)を使用し、ウェーブレット変換に用いることです。
readlineを用いてテキストファイルのリスト化を行っているのですが、指定した行から別の指定した行までのデータを取り出すことが出来ません。
readlineだとそもそも全て読み込みになってしまうからなのでしょうか?
また使用するコード自体が間違っているのでしょうか。
現時点でのコードだと非常に時間がかかりすぎるため、可能であれば時間が短くなればと思います。
現在は5000データから、2000‐4000間の2000データを取り出したい、という状況です。
ご教授お願い致します。

該当のソースコード

-----------------------------------------------------------------------------------------------------------------
def MaxMinNormalization(x):
    x = (x - np.min(x)) / (np.max(x) - np.min(x));
    return x

def morlet(t):
    f0=0.849
    jisuans=np.power(np.pi,-1/4)*(np.cos(-2*np.pi*f0*t)-np.exp(-np.power(2*np.pi*f0,2)/2))*np.exp(-np.power(t,2)/2)
    qyosuans=np.power(np.pi,-1/4)*np.sin(-2*np.pi*f0*t)*np.exp(-np.power(t,2)/2)
    return jisuans,qyosuans

f0=0.849
----------------------------------------------------------------------------------------------------------------
x = np.arange(0.00,2,0.001)
y = open("EMGwave_test.txt")
data=[]
for i in range(5000):
    lines1 = y.readline())
    data.append(float(lines1))
y.close()
y=data[2000:4000]
for line in lines1:

   hz=np.arange(1,100,1)#等差数列
   aa=f0/hz
   summ=[]
   for a in aa:
       wbltb=[]
       for b in x:
            t=(x-b)/a
            wblt1=1/math.sqrt(a)*sum(morlet(t)[0]*y)
            wblt2=1/math.sqrt(a)*sum(morlet(t)[1]*y)
            wblt=np.sqrt(wblt1*wblt1+wblt2*wblt2)
            wbltb.append(wblt)
       summ.append(wbltb)
summ=np.array(summ)

試したこと

readの他のコードや数値の入力等も試してみましたが、引数が~等のエラーによりrun中に止まってしまいました。

補足情報(FW/ツールのバージョンなど)

python3.6.6、Axes3D

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 1

checkベストアンサー

+1

指定行数のみ読み出す、とは考えずに、
最初から順番に読み出していって、指定行になったら取り込む、という考え方で行きましょう
指定行になるまで読み飛ばせばいいし、取り込み分の読み出しが終わったらそこで終わればいいのです

for i in range(開始行数):
    y.readline()
for i in range(読み出し行数):
    lines1 = y.readline()
    data.append(float(lines1))
y.close()

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2018/12/19 13:40

    表示した結果、延々とmorlet(t)[0]が改行込みで出力され続けました。対応する値になっていない、という考え方で良いのでしょうか?丁寧な回答頂いているのに申し訳ありません…。

    キャンセル

  • 2018/12/19 14:06

    morletという関数?がなにか、ということですね。提示されたソースに記述がないのでなんともいえませんが

    キャンセル

  • 2018/12/19 14:34

    morlet wavelet式から導いた関数になります…。一応質問のソースに追記する形で記載しておきました。ここがそもそもソース間違っている可能性もありそうですね。

    キャンセル

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

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

関連した質問

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