編集履歴

質問編集履歴

修正

2017/07/13 08:04

投稿

Penefax

スコア38

test CHANGED Viewed

File without changes

test CHANGED Viewed

File without changes

修正

2017/07/13 08:04

投稿

Penefax

スコア38

test CHANGED Viewed

	@@ -1 +1 @@
1	- ~~動的に得られる関数のプロット(3次元、2次元)~~
1	+ real time plot

test CHANGED Viewed

File without changes

追記

2017/07/13 07:41

投稿

Penefax

スコア38

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -24,63 +24,105 @@
 #グラフ
-d= genfromtxt("Planned.csv", delimiter=",")
-fig = pyplot.figure(figsize=(16,9))
-ax1 = fig.add_subplot(121, projection='3d')
-ax1.set_xlabel("X-axis")
-ax1.set_ylabel("Y-axis")
-ax1.set_zlabel("Z-axis")
-ax1.set_xlim(-30, 30)
-ax1.set_ylim(-40, 20)
-ax1.set_zlim(100, 200)
-ax1.plot(d[:,0], d[:,1], d[:,2], "o", color="g", ms=16, mew=0.5)
-ax2 = fig.add_subplot(222)
-ax2.set_xlabel('X-axis')
-ax2.set_ylabel('Y-axis')
-ax2.set_xlim(-30, 30)
-ax2.set_ylim(-40, 20)
-ax2.plot(d[:,0], d[:,1], "o", color="g", ms=24, mew=0.5)
-ax2.grid(True)
-ax3 = fig.add_subplot(224)
-ax3.set_xlabel('Z-axis')
-ax3.set_ylabel('Y-axis')
-ax3.set_xlim(100, 200)
-ax3.set_ylim(-40,20)
-ax3.plot(d[:,2], d[:,1], "o", color="g", ms=24, mew=0.5)
-ax3.grid(True)
-再追記
+# グラフの初期化
+def init_plot(d):
+    fig = plt.figure()
+    ax = fig.add_subplot(111, projection='3d')
+    ax.set_xlabel("X-axis")
+    ax.set_ylabel("Y-axis")
+    ax.set_zlabel("Z-axis")
+    ax.set_xlim(-30, 30)
+    ax.set_ylim(-40, 20)
+    ax.set_zlim(100, 200)
+    # 最初に描画データを与えておく必要あり
+    l1, = ax.plot(d[:,0], d[:,1], d[:,2], "o", color="g", ms=16, mew=0.5) # 初期データ用
+    l2, = ax.plot([0], [0], [0], "x", color="r", ms=8, mew=0.5)           # 追加データ用 初回はダミーをセット
+    return l1,l2 # 呼出元で利用
+# グラフの描画更新
+def update_plot(lines,d):
+    lines.set_data(d[:,0],d[:,1])   # x,y
+    lines.set_3d_properties(d[:,2]) # z
+    plt.pause(.01)
+if __name__ == "__main__":
+    # 最初から表示するデータ
+    d = genfromtxt("Planned.csv", delimiter=",")
+    # グラフの初期化
+    l1,l2 = init_plot(d) # l1=初期データ, l2=追加データ用
+    update_plot(l1,d)
+    time.sleep(1)
+    # データ追加し描画更新
+    a = np.array([[X,Y,Z]]) # 追加データ
+    update_plot(l2,a)
+    time.sleep(1)
+    # ループで追加
+    for i in range(113):
+        a = np.vstack((a,[X,Y,Z]))
+        update_plot(l2,a)
+        time.sleep(0.5)
+    # 最後に確認表示できるように
+    plt.show()
 #	テンプレート画像読み込み
@@ -248,51 +290,7 @@
 追記
-def pause_plot():
-    fig, ax = plt.subplots(1, 1)
-    x4 = X
-    y4 = Y
-    scatter, = ax.plot(x4, y4, "x", color="r", ms=16, mew=0.5)
-    while True:
-        x4 = X
-        y4 = Y
-        ．
-        scatter.set_data(x4, y4)
-        ax.set_xlim(-30, 30)
-        ax.set_ylim(-40, 20)
-        plt.pause(.01)
-        ax.grid(True)
-if __name__ == "__main__":
-    pause_plot()
 ```

再追記

2017/07/11 05:34

投稿

Penefax

スコア38

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -16,7 +16,7 @@
 ```
-というエラーコードが出てしまいます。Planned.csvの乗った散布図に描画されるようにしたい。
+というエラーコードが出てしまいます。
 ###該当のソースコード
@@ -44,8 +44,6 @@
 ax1.set_zlim(100, 200)
 ax1.plot(d[:,0], d[:,1], d[:,2], "o", color="g", ms=16, mew=0.5)
@@ -60,8 +58,6 @@
 ax2.set_ylim(-40, 20)
 ax2.plot(d[:,0], d[:,1], "o", color="g", ms=24, mew=0.5)
 ax2.grid(True)
@@ -78,15 +74,107 @@
 ax3.set_ylim(-40,20)
 ax3.plot(d[:,2], d[:,1], "o", color="g", ms=24, mew=0.5)
 ax3.grid(True)
+再追記
+#	テンプレート画像読み込み
+img_obj = cv2.imread('1027.jpg', cv2.IMREAD_COLOR)
+if img_obj is None:								#	テンプレート画像が取得できない場合
+	print u'画像が取得できません。'
+	sys.exit()
+#	映像取得（カメラ映像）
+src = cv2.VideoCapture('0deg_10fps_x20 ver3.avi')
+if not src.isOpened():								#	カメラ映像が取得できない場合
+	print u'映像が取得できません。'
+	sys.exit()
+retval, frame = src.read()
+height1, width1, channels1 = img_obj.shape
+height2, width2, channels2 = frame.shape
+rec = cv2.VideoWriter('0deg_10fps_x20 改.avi', cv.CV_FOURCC('X','V','I','D'), 6, (width2, height2))
+f = open('result.csv', 'ab') #csvファイルが無ければ作る、の'a'を指定します
+csvWriter = csv.writer(f)
+val = 0
+while True:
+	retval, frame = src.read()        #	1フレーム取得
+        if frame is None:
+		break
+	#	テンプレート・マッチングにより相互相関係数を計算
+	img_ccoeff1 = cv2.matchTemplate(frame, img_obj, cv2.TM_CCOEFF_NORMED)
+	#	[-1, 1] を [0, 1] へ
+	cv2.normalize(img_ccoeff1,img_ccoeff1, 0, 1, cv2.NORM_MINMAX)
+	#	相互相関係数の最小値・最大値とその座標を抽出
+	cMin, cMax, pMin, pMax1 = cv2.minMaxLoc(img_ccoeff1)
+	#	検出領域の中心座標
+	detect = (pMax1[0] + width1/2, pMax1[1] + height1/2)
+再追記終わり
+        x = pMax1[0] + width1/2, pMax1[1] + height1/2
         x1 = (((round(((pMax1[0] + width1/2)-385)*0.327,2))), (round(((pMax1[1] + height1/2)-290)*0.327,2)))
@@ -94,6 +182,52 @@
         y2 = (round(((pMax1[1] + height1/2)-290)*0.327,2))
+        a = np.array([x2,y2])
+        text1 = "" + str(x1) + ""
+        text2 = "" + str(x2) + ""
+        text3 = "" + str(y2) + ""
+        cv2.putText(frame, text1, ((pMax1[0] + width1/2)-65, (pMax1[1] + height1/2)-30), cv2.FONT_HERSHEY_COMPLEX, 0.7, (0, 255, 0), 2, cv2.CV_AA)
+        cv2.putText(frame, "L1", ((pMax1[0] + width1/2)-70, (pMax1[1] + height1/2)-60), cv2.FONT_HERSHEY_COMPLEX, 0.7, (0, 255, 0), 2, cv2.CV_AA)
+        #	探索画像から検出領域を抽出
+	img_crop = frame[pMax1[1]:pMax1[1]+height1, pMax1[0]:pMax1[0]+width1].copy()
+	#	検出領域に赤色の円と十字を描画
+	cv2.circle(frame, detect, width1/2, (0, 0, 255), 1)
+        cv2.line(frame,((pMax1[0] + width1/2) - 25, pMax1[1] + height1/2), ((pMax1[0] + width1/2) + 25 , pMax1[1] + height1/2), (0, 0, 255), 1)
+        cv2.line(frame,((pMax1[0] + width1/2), (pMax1[1] + height1/2 )-25), ((pMax1[0] + width1/2), (pMax1[1] + height1/2) + 25), (0, 0, 255), 1)
+        x1 = (((round(((pMax1[0] + width1/2)-385)*0.327,2))), (round(((pMax1[1] + height1/2)-290)*0.327,2)))
+        x2 = (round(((pMax1[0] + width1/2)-385)*0.327,2))
+        y2 = (round(((pMax1[1] + height1/2)-290)*0.327,2))
         y1 = (round(((pMax2[0] + width1/2)-385)*0.327,2), round(((pMax2[1] + height1/2)-290)*0.327,2))
@@ -124,60 +258,30 @@
     y4 = Y
-    # 初期化的に一度plotしなければならない
-    # そのときplotしたオブジェクトを受け取る必要がある．
-    # listが返ってくるので，注意
     scatter, = ax.plot(x4, y4, "x", color="r", ms=16, mew=0.5)
-    # ここから無限にplotする
     while True:
-        # plotデータの更新
         x4 = X
         y4 = Y
-        # 描画データを更新するときにplot関数を使うと
-        # lineオブジェクトが都度増えてしまうので，注意．
+        ．
-        # 一番楽なのは上記で受け取ったlinesに対して
-        # set_data()メソッドで描画データを更新する方法．
         scatter.set_data(x4, y4)
-        # set_data()を使うと軸とかは自動設定されないっぽいので，
-        # 今回の例だとあっという間にsinカーブが描画範囲からいなくなる．
-        # そのためx軸の範囲は適宜修正してやる必要がある．
         ax.set_xlim(-30, 30)
         ax.set_ylim(-40, 20)
-        # 一番のポイント
-        # - plt.show() ブロッキングされてリアルタイムに描写できない
-        # - plt.ion() + plt.draw() グラフウインドウが固まってプログラムが止まるから使えない
-        # ----> plt.pause(interval) これを使う!!! 引数はsleep時間
         plt.pause(.01)
@@ -206,10 +310,4 @@
 ###補足情報(言語/FW/ツール等のバージョンなど)
-windows10
-python2.7
-matplotlib
+windows10,python2.7,matplotlib,spyder
-spyder

追記

2017/07/10 23:55

投稿

Penefax

スコア38

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -2,7 +2,7 @@
 動画からテンプレートマッチングにより入力画像の中心座標を求めてそれを三角測量の式を用いてX,Y,Zを出しています。
-X,Y,Zを求める計算とほぼ同時のタイミングでグラフにプロットされるようにしたい。
+X,Y,Zを求める計算と同じタイミングでグラフにプロットされるようにしたい。
@@ -10,10 +10,14 @@
 ```
+Using default event loop until function specific to this GUI is implemented
-プロットされない。
+  warnings.warn(str, mplDeprecation)
 ```
+というエラーコードが出てしまいます。Planned.csvの乗った散布図に描画されるようにしたい。
 ###該当のソースコード
 ```Python
@@ -82,7 +86,7 @@
-   #三角測量
         x1 = (((round(((pMax1[0] + width1/2)-385)*0.327,2))), (round(((pMax1[1] + height1/2)-290)*0.327,2)))
@@ -106,6 +110,86 @@
         Z = ((round(50/(l-50)*213,2)))
+追記
+def pause_plot():
+    fig, ax = plt.subplots(1, 1)
+    x4 = X
+    y4 = Y
+    # 初期化的に一度plotしなければならない
+    # そのときplotしたオブジェクトを受け取る必要がある．
+    # listが返ってくるので，注意
+    scatter, = ax.plot(x4, y4, "x", color="r", ms=16, mew=0.5)
+    # ここから無限にplotする
+    while True:
+        # plotデータの更新
+        x4 = X
+        y4 = Y
+        # 描画データを更新するときにplot関数を使うと
+        # lineオブジェクトが都度増えてしまうので，注意．
+        # 一番楽なのは上記で受け取ったlinesに対して
+        # set_data()メソッドで描画データを更新する方法．
+        scatter.set_data(x4, y4)
+        # set_data()を使うと軸とかは自動設定されないっぽいので，
+        # 今回の例だとあっという間にsinカーブが描画範囲からいなくなる．
+        # そのためx軸の範囲は適宜修正してやる必要がある．
+        ax.set_xlim(-30, 30)
+        ax.set_ylim(-40, 20)
+        # 一番のポイント
+        # - plt.show() ブロッキングされてリアルタイムに描写できない
+        # - plt.ion() + plt.draw() グラフウインドウが固まってプログラムが止まるから使えない
+        # ----> plt.pause(interval) これを使う!!! 引数はsleep時間
+        plt.pause(.01)
+        ax.grid(True)
+if __name__ == "__main__":
+    pause_plot()
 ```
 このX,Y,Zをax1(X,Y,Z),ax2(X,Y),ax3(Y,Z)にプロットしたいです。
@@ -114,7 +198,9 @@
 ###試したこと
+[matplotlibでリアルタイム描画](http://qiita.com/hausen6/items/b1b54f7325745ae43e47)を参考にしてコードを考えてみましたが上手くいきませんでした。
-色々、定義したりそれを変えたりもしましたが上手くいきませんでした。
+URL先ではsin関数を描画していますが、散布図でプロットしたいです。