半年ぐらい棚ざらしにしてあったのですが、ちょっと時間があったのでいじることに。（でも、最終成果物が出来るのは、さらに一年後ぐらい、かなぁ？orz）
＜＜追記＞＞
こちらの記事に続く
とりあえずcvCalcOpticalFlowBM関数のサンプルコードが動くようになったので無くさないようにblog保存。
ググる:オプティカルフローとは
ビデオカメラ、YouTubeやニコニコ動画などのFLVデータ、などなど
地デジが映るのもオプティカルフローのおかげ。
以下はコンピューターのプログラム＆チラシの裏。
最終的に作ろうと思ってるものに一番近いのがコレ
既に有るんやから作らなくても。。。と言いたいところですが、いろんな都合もあって、自分でもチャレンジしてみることに。
っていうか、そう思ったのが去年の頭だから、とりかかるまでに時間かかりすぎですが（自分）
（実は、Pythonでは、オプティカルフローの量を検出するところまでは書けても、それを再利用する部分のライブラリが無い。 Distortionマッピングとか、 Displaceマッピングと呼ばれる機能はFlashにも有るので、将来、PILなどにも実装されると嬉しいなぁ。Cで書かれたライブラリを探したけど見つからず。なので、一部、GIMPの機能を使うことに）
次は、これをGIMPプラグイン用に書き直さないと （の前にGIMPをアンインストールしたので、もう一度入れなおさないと。。。）
プログラム言語という魔法の呪文
ちょっと休憩。
私自身も、使ってないプログラム言語はサッパリサッパリで魔法の呪文に見えます。
（とくに人の書いたプログラムコードは「進んだ科学は魔法と見分けがつかない」ってのを、実感したりするわけですけれども）
だから、判らないもの＝魔法、的な視線をなげかけるのは（私も含めて）しょうがないことではあるのですけれども
ただ、魔法との一番の違いって。。。
- 杖を一振りして呪文となえるだけ、ってわけにはいかないよなぁとか。
- そのかわり、手順を公開して、ステップ・バイ・ステップで行えば、ある程度、再現性があるよね、という
（イチローとか、卓球の愛ちゃんみたく、特殊な技能、才能を必要とすることもありますけれども。）
文明とか、科学技術ってのは、そういう意味で（文字通り、魔法をタネあかしすると、手品になってしまう、というのと同じ文脈で） 『知の特権性』を奪ってしまう部分があるのかなぁ、と。
（いいか悪いかは別として、コンピューターの根本的な原理を知らなくてもプログラムは書ける。また、誰でも出来るようになると陳腐化してしまう。）
それをつまらない、と思うか、誰でも（特に自分に）出来るようになるのはすばらしい、と思うかは人それぞれ、だとは思うのですが。
使ってるモジュールはswig版です。
_import opencv
_from opencv import cvCalcOpticalFlowBM, cvSize,\
_ cvPoint,cvLine,cvScalar,cvCreateMat,cvSetZero,\
_ CV_32FC1,cvGetReal2D,CV_RGB,CV_AA
_
_import opencv.highgui as highgui
_from opencv.highgui import cvLoadImage,cvNamedWindow,\
_ cvLoadImageM,cvShowImage,\
_ CV_LOAD_IMAGE_GRAYSCALE,CV_LOAD_IMAGE_COLOR
_
_import math
_from math import ceil
_import Image
_
_"""
_cvCalcOpticalFlowBM(CvArr prev, CvArr curr, CvSize block_size, CvSize shift_size,
_ CvSize max_range, int use_previous,
_ CvArr velx, CvArr vely)
_"""
_fnames=["01.png","02.png","03.png","04.png"]
_img1=cvLoadImage (fnames[0],CV_LOAD_IMAGE_GRAYSCALE)
_img2=cvLoadImage (fnames[1],CV_LOAD_IMAGE_GRAYSCALE)
_dst_img = cvLoadImage (fnames[1], CV_LOAD_IMAGE_COLOR)
_
_block_size = 10;
_shift_size = 1;
_rows = int (ceil (img1.height / float(block_size)))
_cols = int (ceil (img1.width/ float(block_size)))
_velx = cvCreateMat (rows, cols, CV_32FC1)
_vely = cvCreateMat (rows, cols, CV_32FC1)
_cvSetZero (velx)
_cvSetZero (vely)
_
_
_blocksize = cvSize (block_size ,block_size )
_shiftsize = cvSize (shift_size,shift_size)
_max_range = cvSize (50, 50);
_cvNamedWindow("opencv_test")
_cvShowImage("opencv_test",dst_img)
_
_cvCalcOpticalFlowBM (img1, img2, blocksize, shiftsize,\
_ max_range,0,velx,vely)
_
_# (3)計算されたフローを描画
_for i in xrange(velx.width):
_ for j in xrange (vely.height) :
_ dx = int(cvGetReal2D (velx, j, i))
_ dy = int(cvGetReal2D (vely, j, i))
_ cvLine (dst_img, cvPoint (i * block_size, j * block_size),
_ cvPoint (i * block_size + dx, j * block_size + dy), CV_RGB (255, 0, 0), 1, CV_AA, 0);
_
_
_cvShowImage("opencv_test",dst_img)