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

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

ただいまの
回答率

88.79%

信号のフーリエ変換におけるエラーの解決法

受付中

回答 1

投稿

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

santamarianeed

score 21

 動画内における光を強度を信号として得たあとのフーリエ変換したプログラムを作成しましたがどうしてもエラーが生じてしまいます。詳細は以下の通りです。何がいけないのでしょうか?

ソースコード

# -*- coding: utf-8 -*-
import cv2
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
from scipy import signal

video_path = "firstcrop.mkv"
cap = cv2.VideoCapture(video_path)

count = cap.get(cv2.CAP_PROP_FRAME_COUNT)
fps = cap.get(cv2.CAP_PROP_FPS)

numS=raw_input("start frame:  ")
numS=int(numS)

#numE=raw_input("end frame:  ")
#numE=int(numE)
numE=count

# 場所を決める
#点滅の場所
x1min=raw_input("点滅1左上のx座標:")
y1min=raw_input("点滅1左上のy座標:")
x1max=raw_input("点滅1右下のx座標:")
y1max=raw_input("点滅1右下のy座標:")

X1min=int(x1min)
Y1min=int(y1min)
X1max=int(x1max)
Y1max=int(y1max)

#x2min=raw_input("点滅2左上のx座標:")
#y2min=raw_input("点滅2左上のy座標:")
#x2max=raw_input("点滅2右下のx座標:")
#y2max=raw_input("点滅2右下のy座標:")

#X2min=int(x2min)
#Y2min=int(y2min)
#X2max=int(x2max)
#Y2max=int(y2max)

#x3min=raw_input("点滅3左上のx座標:")
#y3min=raw_input("点滅3左上のy座標:")
#x3max=raw_input("点滅3右下のx座標:")
#y3max=raw_input("点滅3右下のy座標:")

#X3min=int(x3min)
#Y3min=int(y3min)
#X3max=int(x3max)
#Y3max=int(y3max)

#ノイズの場所
xnoisemin=raw_input("左上ノイズx座標:")
ynoisemin=raw_input("左上ノイズy座標:")
xnoisemax=raw_input("右下ノイズx座標:")
ynoisemax=raw_input("右下ノイズy座標:")

Xnoisemin=int(xnoisemin)
Ynoisemin=int(ynoisemin)
Xnoisemax=int(xnoisemax)
Ynoisemax=int(ynoisemax)

# フレーム [numS,numE] の範囲で各フレームの [ymin, ymax]x[xmin, xmax] の画素の平均を計算する。
frame_no = np.arange(numS,numE)
intensity1 = []
intensity2=[]
intensity3=[]
intensitynoise=[]
for i in frame_no:
    # フレーム frame_no を取得する。
    cap.set(cv2.CAP_PROP_POS_FRAMES, i)
    ret, frame = cap.read()
    if not ret:
        print('Failed to grab frame.')
        break
    # 平均画素値を計算する。
    intensity1.append(frame[Y1min:Y1max+1, X1min:X1max+1].mean())
    #intensity2.append(frame[Y2min:Y2max+1, X2min:X2max+1].mean())
    #intensity3.append(frame[Y3min:Y3max+1, X3min:X3max+1].mean())
    intensitynoise.append(frame[Ynoisemin:Ynoisemax+1, Xnoisemin:Xnoisemax+1].mean())
    print "frame number ",i


#Fourier transform
A=np.fft.fft(intensity1);
B=np.fft.fft(intensitynoise);
amplitudeA=abs(A)
amplitudeB=abs(B)

# データのパラメータ
N=len(frame_no)#サンプル数
dt = 0.01          # サンプリング間隔
freq = np.linspace(0, 1.0/dt, N) # 周波数軸

plt.subplot(2,1,1)
t1Points=np.arange(len(intensity1))*(1/fps)
tnoisePoints=np.arange(len(intensitynoise))*(1/fps)
plt.title("point6 and noise")
plt.plot(t1Points,intensity1,linewidth=1,color="r")
plt.plot(tnoisePoints,intensitynoise,linewidth=1,color="k")
plt.grid(True)
plt.xlabel("time(s)", fontsize=17)
plt.ylabel("intensity of PL",fontsize=17)
plt.legend(fontsize=17)
plt.xlim(numS/fps,numE/fps)#開始時間、終了時間の設定
plt.ylim(0,255)


plt.subplot(2,1,2)
plt.plot(freq,amplitudeA,color="r")
plt.plot(freq,amplitudeB,color="k")
plt.xlabel('Frequency', fontsize=17)
plt.ylabel('Amplitude', fontsize=17)
plt.ylim(0,10000)
plt.grid()
plt.legend(fontsize=17)

plt.show()

エラーメッセージ

/usr/lib/python2.7/dist-packages/matplotlib/axes/_axes.py:519: UserWarning: No labelled objects found. Use label='...' kwarg on individual plots.
warnings.warn("No labelled objects found. "
Traceback (most recent call last):
File "特定領域の強度とFFT_1.py", line 111, in <module>
plt.plot(freq,amplitudeA,color="r")
File "/usr/lib/python2.7/dist-packages/matplotlib/pyplot.py", line 3154, in plot
ret = ax.plot(*args, kwargs) File "/usr/lib/python2.7/dist-packages/matplotlib/init.py", line 1814, in inner return func(ax, *args, kwargs)
File "/usr/lib/python2.7/dist-packages/matplotlib/axes/_axes.py", line 1424, in plot
for line in self._get_lines(*args, **kwargs):
File "/usr/lib/python2.7/dist-packages/matplotlib/axes/_base.py", line 386, in _grab_next_args
for seg in self._plot_args(remaining, kwargs):
File "/usr/lib/python2.7/dist-packages/matplotlib/axes/_base.py", line 364, in _plot_args
x, y = self._xy_from_xy(x, y)
File "/usr/lib/python2.7/dist-packages/matplotlib/axes/_base.py", line 223, in _xy_from_xy
raise ValueError("x and y must have same first dimension")
ValueError: x and y must have same first dimension


(program exited with code: 1)
Press return to continue

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

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

  • 1T2R3M4

    2019/05/30 20:59

    /usr/lib/python2.7ですか?
    TAGは3.xですが。

    キャンセル

  • santamarianeed

    2019/05/31 08:35

    そうです!python2.7です。
    間違えました

    キャンセル

回答 1

0

エラーメッセージを見る感じ、freqとamplitudeAが1次元かつ長さが同じでなければならないようです。大きさにずれはないでしょうか。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

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

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

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