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

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

ただいまの
回答率

90.50%

  • Python 2.7

    1266questions

    Python 2.7は2.xシリーズでは最後のメジャーバージョンです。Python3.1にある機能の多くが含まれています。

  • NumPy

    448questions

    NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

pythonで複素数のデータを読み書き

解決済

回答 1

投稿

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

 pythonで複素数のデータを読み書き

python2.7で複素数のデータを扱いたいのですが、()がついたり、読み込みがうまく行かなかったりでうまく動作しません。

 該当のソースコード

#!/usr/bin/env python2
# -*- coding: utf-8 -*-
import numpy as np

(a1) = np.genfromtxt("○○○.csv",unpack=True, dtype=np.complex128,delimiter=',')

print a1
ファイルの書き込みは
np.savetxt("○○○.csv",data , fmt='%.16f',delimiter=',')
で行いました。

dataには(1.9464469341504152e-12-3.1456638809610314e-12j), (8.78883381959785e-13-3.852325732580186e-13j), (-9.583868301245934e-13-1.873118734890551e-12j),こんな感じで複素数が入ってます。

 試したこと

結果として

[-2.021478e-10+9.414360e-10j -1.836750e-11+8.232035e-10j
6.789447e-10+3.419484e-10j ...           nan+0.000000e+00j
nan+0.000000e+00j -1.839600e-12+1.562000e-13j]

のように途中でnanと欠損していしまいます。csvにはしっかりと書き込まれいています。
読み込みが悪いのか、書き込みが悪いのかもよくわかりません。
よろしくお願いします。

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 1

checkベストアンサー

+2

numpyのバグで、虚部の符号がマイナスの場合、

出力csvファイル

 (1.9464469341504152e-12+-3.1456638809610314e-12j)
 (8.7888338195978500e-13+-3.8523257325801861e-13j)
 (-9.5838683012459342e-13+-1.8731187348905510e-12j)

のように+-というフォーマットで出力されてしまうという問題です。そのため読み込むことができません。

リンク:
python - numpy.genfromtxt give nan for complex numbers with minus signs? - Stack Overflow
savetxt complex number bug · Issue #7895 · numpy/numpy · GitHub
loadtxt does not accurately save float128 values · Issue #10436 · numpy/numpy · GitHub

悲しいことに、2016年だかに報告されて、まだ直る気配がないようです。

回避策は上のリンクの1つめのstackoverflowに載っていますが、要するに自分でコンバータ作れば良いので、やり方はそこの方法に限らず色々あると思います。

面倒なら最初からcsv出力はやめて、np.save()np.load()でやった方が正直良いと思います。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2018/07/28 14:41

    返信遅れて申し訳ありません。
    なるほどnumpyのバグなのですね。

    np.save()とnp.loadで複素数の部分の計算はやることにしました。
    ありがとうございました。

    キャンセル

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

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

関連した質問

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

  • Python 2.7

    1266questions

    Python 2.7は2.xシリーズでは最後のメジャーバージョンです。Python3.1にある機能の多くが含まれています。

  • NumPy

    448questions

    NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。