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

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

ただいまの
回答率

88.32%

pythonを使用した行列をテキストファイルに出力

解決済

回答 1

投稿

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

kamikura

score 16

pythonを使用して、指定のテキストファイルからデータを読み込み、ソースコード中の計算式で計算した結果を出力しようとしています。
現在、出力という行為自体はうまくいっておりテキストファイル、またpython shellへの出力は出来ているのですが…。
行列として出力を行うと以下のような結果になってしまいます。

[1 2 3 4 5 6
7 8 9 10 11 12
13 14 ...]
[1a 2a 3a 4a 5a 6a
7a 8a 9a 10a 11a 12a
13a 14a...]

上記の様に6列複数行の形で出力されてしまいます()。
データ量が膨大なため、何十列となるのは分かっていますがこれを1行ずつの行列に変換できないでしょうか?
理想としては上記の出力結果を

[1 2 3 4 5 6 7 8 9 10 11 12 13 14 ...]
[1a 2a 3a 4a 5a 6a 7a 8a 9a 10a 11a 12a 13a 14a ...]

のように出来ればと考えております。
質問内容が分かりにくければコメントでご指摘お願い致します。
宜しくお願いします。

以下、使用しているソースコードになります(3次元グラフも出力しますが、ソースコードは省いてあります)。

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
import math
import numpy as np

np.set_printoptions(threshold=np.inf)

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,1.00,0.001)
y = open("EMGwave_test.txt")

lines1 = y.read()
y.close()

lines2 = lines1.split("\n")
data=[]
for line in lines2:
    data.append(float(line))
y=data[0:1000]

#Wavelet
Hz=np.arange(50,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)


file = open("writing.txt","w")
file.write(str(summ))
file.close()
  • 気になる質問をクリップする

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

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

  • can110

    2018/12/21 11:05

    出力結果は「[~]」で囲む、および区切り文字は空白である必要はありますか?
    すなわち「[1 2 3 4~]」ではなく「1,2,3,4,~」という出力でもかまいませんか?

    キャンセル

  • kamikura

    2018/12/21 11:09

    数値ごとにカンマで区切る、ということでしょうか?
    出力されたデータをmatlabで流用できる形であれば、どのような形でも構いません。
    また流用出来ない、という場合でも現状では問題ない為どちらの形でも大丈夫です。

    キャンセル

回答 1

checkベストアンサー

+1

これは "略式の" つまり読み易い object の文字列表現です。

https://docs.python.jp/3/library/stdtypes.html#str

str(summ)は人間が見て読みやすい文字列表現を返します。適当に略されたり丸められたりしますし、一般的な入力フォーマットとの親和性も低いので、データをdumpするような目的にはまったく向いていません。

そういう用途であればnumpy.savetxtあたりを使ってください。

numpy.savetxt — NumPy v1.15 Manual

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2018/12/21 11:47

    頂いたコードで、上記のような形で出力されました!
    ありがとうございます。
    こちら、最終列はどうしても改行されてしまうのでしょうか…?
    その形でもdumpするような目的で使用可であれば、この形でも大丈夫です。

    キャンセル

  • 2018/12/21 11:50

    >こちら、最終列はどうしても改行されてしまうのでしょうか…?
    ちょっと言いたいことがよくわからなかったので、詳しく説明してもらえると助かります。

    np.savetxtで書き出したものは基本的にnp.loadtxtで読み込めます。使い方とかは検索すれば出てくるのでここでは書きませんが。

    キャンセル

  • 2018/12/21 12:29

    ごめんなさい、こちらの表示上の問題でした。
    2行3列の形で出力出来ていることを確認できました。
    ご丁寧にありがとうございました!
    np.loadtxtについてはこちらで検索してみます。

    キャンセル

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

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

関連した質問

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