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

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

ただいまの
回答率

88.61%

配列同士の足し算ができない。要素が増えてしまう結果となる

解決済

回答 1

投稿

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

santamarianeed

score 21

時系列で動画の強度が分かるプログラムを作りました。(intensity1とintensity2)

ここで同じ時間におけるintensity1とintensity2の強度の和をとってplotしたいのですができません。(コードをみたら分かると思います)
どうしたら良いのでしょうか?

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

video_path = "water50_2.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
print(numE)

# 場所を決める
#点滅の場所
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)

# フレーム [numS,numE] の範囲で各フレームの [ymin, ymax]x[xmin, xmax] の画素の平均を計算する。
frame_no = arange(numS,numE)
intensity1 = []
intensity2=[]
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())
    print "frame number ",i


A=array(intensity1)
#print(D)
B=array(intensity2)
C=intensity1+intensity2
print(C)
kazu=len(C)
print(kazu)
#print(intensity1)
#X=stack([intensity1,intensitynoise], axis=1)
#print(X)

#Y=sum(X)
#print(Y)
#Z=stack([Y,tPoints], axis=1)

tPoints=arange(len(intensity1))*(1/fps)


#plt.plot(tPoints,intensitynoise,linewidth=1,color="k")
#Z=stack([Y,tPoints], axis=1)
plt.grid(True)
plt.xlabel("time(s)")#横軸を時間に
plt.ylabel("intensity of PL")
plt.legend(fontsize=17)
#plt.title("1")
plt.xlim(numS/fps,numE/fps)#開始時間、終了時間の設定
plt.legend(loc="lower right")
plt.title("Intensity")
set_printoptions(threshold=inf)
#print(vstack((tPoints,intensity1)).T)
plt.show()

plt.grid(True)
#plt.title(Y1mintoY1max,X1mintoX1max)
plt.title("noise:%d,%d"%(xnoisemin,ynoisemin))
  • 気になる質問をクリップする

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 1

checkベストアンサー

0

せっかくndarrayにしているのですから、C = A + B で良いでしょう。

>>> lst1 = [1, 2, 3]
>>> lst2 = [10, 20, 30]
>>>
>>> lst1 + lst2
[1, 2, 3, 10, 20, 30]
>>>
>>> arr1 = np.array(lst1)
>>> arr2 = np.array(lst2)
>>>
>>> arr1 + arr2
array([11, 22, 33])

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

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

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

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