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

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

新規登録して質問してみよう
ただいま回答率
85.48%
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

600閲覧

pythonのimport、from文について

kamikura

総合スコア16

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2018/12/29 01:30

編集2018/12/29 01:47

pythonでimportをしたいのですが、ファイル名を入力しても存在しないと返ってきてしまいます。
from import文でも試しましたが、同じく存在しないという結果でした。
ファイル内を見てみるとどうやらコンパイル自体が出来てない様子です(import出来てないから当たり前と言えばそうですが)

python

1from __future__ import division 2import numpy as np 3import matplotlib.pyplot as plt 4import pycwt as wavelet 5from pycwt.helpers import find

呼び出したいファイルとしては
Python36\Lib\site-packages\pycwt\pycwt\helpers
のfind関数と
Python36\Lib\site-packages\pycwt\pycwt\mothers
のMorlet関数になります。
余りに初歩的な問題かもしれませんが、宜しくお願いします。

質問がありましたので、ソースコード全体とエラー文を追記します。
ソースコードから

python

1#ライブラリ読み込み 2from __future__ import division 3import numpy as np 4import matplotlib.pyplot as plt 5import pycwt as wavelet 6from pycwt.helpers import find 7 8#分析データ読み込み 9url = "http://paos.colorado.edu/research/wavelets/wave_idl/nino3sst.txt" 10dat = np.genfromtxt(url, skip_header=19) 11 12#メタデータ設定 13title = "xion t" 14label = "xion sst" 15unit = "deg" 16 17t0 = 1871.0 18dt = 0.25 19 20#時間データ生成 21N = dat.size 22t = np.arange(0,N) * dt + t0 23 24#分析対象データの整形 25p = np.polyfit(t - t0, dat, 1) 26dat_notrend = dat - np.polyval(p, t - t0) 27std = dat_notrend.std()#標準偏差 28var = std ** 2#分散 29dat_norm = dat_notrend / std#正規化 30 31#ウェーブレット変換のパラメタ設定 32mother = mothers.Morlet(6) 33s0 = 2 * dt 34dj = 1 / 12 35J = 7 / dj 36alpha, _, _ = wavelet.ar1(dat) 37 38#ウェーブレット変換と逆ウェーブレット変換 39wave, scales, freqs, coi, fft, fftfreqs = wavelet.cwt(dat_norm, dt, dj, s0, J, mother) 40iwave = wavelet.icwt(wave, scales, dt, dj, mother) * std 41 42#ウェーブレットとフーリエの各スペクトル算出 43power = (np.abs(wave)) ** 2 44fft_power = np.abs(fft) ** 2 45period = 1 / freqs 46 47#規格化 48power /= scales[:, None] 49 50#パワースペクトル95%信頼区間での優位性検証 51signif, fft_theor = wavelet.significance(1.0, dt, scales, 0, alpha, significance_lavel=0.95, wavelet=mother) 52sig95 = np.ones([1, N]) * signif[:, None] 53sig95 = power / sig95 54 55#グローバルウェーブレットスペクトルとその優位性の算出 56glbl_power = power.mean(axis=1) 57dof = N - scales 58glbl_signif, tmp = wavelet.signifcance(var, dt, scales, 1, alpha, significance_lavel=0.95, dof=dof, wavelet=mother) 59 60#スケールの平均値とその優位性の算出 61sel = find((period >=2) & (period < 8)) 62Cdelta = mother.cdelta 63scale_avg = (scales * np.ones((N, 1))).transpose() 64scale_avg = power / scale_avg 65scale_avg = var * dj * dt / Cdelta * scale_avg[sel, :].sum(axis=0) 66scale_avg_signif, tmp = wavelet.significance(var, dt, scales, 2, alpha, significance_lavel=0.95, dof=[scales[sel[0]],scales[sel[-1]]], wavelet=mother) 67 68#ウェーブレット解析結果の可視化 69plt.close("all") 70plt.ioff() 71figprops = dict(figsize=(11, 8), dpi=72) 72fig = plt.figure(**figprops) 73 74ax = plt.axes([0.1, 0.75, 0.65, 0.2]) 75ax.plot(t, iwave, "-", linewidth=1, color=[0.5, 0.5, 0.5]) 76ax.plot(t, dat, "k", linewidth=1.5) 77ax.set_title("a) {}".format(title)) 78ax.set_ylabel(r"{} [{}]".format(lavel, units)) 79 80bx = plt.axes([0.1, 0.37, 0.65, 0.28], sharex=ax) 81levels = [0.0625, 0.125, 0.25, 0.5, 1, 2, 4, 8, 16] 82bx.contourf(t, np.log2(period), np.log2(power), np.log2(levels), extend="both", cmap=plt.cm.viridis) 83extent = [t.min(), t.max(), 0, max(period)] 84bx.contour(t, np.log2(period), sig95, [-99, 1], colors="k", linewidths=2, extent=extent) 85bx.fill(np.concatenate([t, t[-1:], +dt, t[-1:], + dt, t[:1] - dt, t[:1] - dt]), 86 np.concatenate(np.log2(coi), [1e-9], np.log2(period[-1:]), np.log2(period[-1:]), [1e-9]), "k", alpha=0.3, hatch="x") 87bx.set_title("b){}Wavelet Power Spectrum ({})".format(label, mother.name)) 88bx.set_ylabel("Period (years)") 89 90Yticks = 2 ** np.arange(np.ceil(np.log2(period.min())), np.ciel(period.max())) 91bx.set_yticks(np.log2(Yticks)) 92bx.set_yticklabels(Yticks) 93 94cx = plt.axes([0.77, 0.37, 0.2, 0.28], sharey=bx) 95cx.plot(glbl_signif, np.log2(period), "k--") 96cx.plot(var * fft_theor, np.log2(period), "--", color="#ccccc") 97cx.plot(var * fft_power, np.log2(1./fftfreqs), "-", color="#ccccc", linewidth=1.) 98cx.plot(var * glbl_power, np.log2(period), "k-", linewidth=1.5) 99cx.set_title("c) Global Wavelet Spectrum") 100cx.set_xlabel(r"Power [({})^2]".format(units)) 101cx.set_xlim([0, glbl_power.max() +var]) 102cx.set_ylim(np.log2([period.min(), period.max()])) 103cx.set_yticks(np.log2(Yticks)) 104cx.set_yticklabels(Yticks) 105plt.setp(cx.get_yticklabels(), visible=False) 106 107dx = plt.axes([0.1, 0.07, 0.65, 0.2], sharex=ax) 108dx.axhline(scale_avg_signif, color="k", linestyle="--", linewidth=1.) 109dx.plot(t, scale_avg, "k-", linewidth=1.5) 110dx.set_title("d) {}--{} year scale-averaged power".format(2, 8)) 111dx.set_xlabel("Time (year)") 112dx.set_ylabel(r"Average variance [{}]".format(units)) 113ax.set_xlim([t.min(), i.max()]) 114 115plt.show()

エラー文

python

1Traceback (most recent call last): 2 File "Z:\03作業用フォルダ\学生\テキスト読み込みできるwavelet変換\wavelet_power_spectle.py", line 6, in <module> 3 from pycwt.helpers import find 4ModuleNotFoundError: No module named 'pycwt.helpers'

以上です。よろしくお願いします。

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

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

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

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

can110

2018/12/29 01:37

実際にエラーが発生するソースコードと発生したエラー内容(Traceback全文)をそのまま記載ください。 また、pycwt はインストールしていますか?あるいは、どのような手順でインストールしましたか?
kamikura

2018/12/29 01:51

pipを使ったものではなく、ファイルを直接DLしてそれをsite_packages内に入れました。 pythonのDLが悪かったのか、pipが使用できない(バージョン確認も出来ない)状態の為、上記の様にしています。
can110

2018/12/29 01:54

どこからどれをDLしましたか?
can110

2018/12/29 01:59

回答しました。インストール作業が必要です。
guest

回答1

0

ベストアンサー

Installationに従い、コードをダウンロードしpython setup.py installする必要があります。

投稿2018/12/29 01:58

can110

総合スコア38266

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

kamikura

2018/12/29 02:10

そもそもインストールが出来ていなかったんですね…。 ありがとうございました、試してみます。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問