« Nature Neurosceince 2月号 | 最新のページに戻る | Plaid (Movshon vs. Born) »
■ plaid (Movshon vs. Born)
図1
いつも先走りぎみなのですこし丁寧にいってみます。
"Dynamics of motion signaling by neurons in macaque area MT." J Anthony Movshon。格子模様(plaid:発音注意。pl[ae]dです)の話。Gifアニメーションで作ってみました(図1)。縞模様が動いているところ(右上方向)にもうひとつ角度の違う縞模様が動いているの(右下方向)を重ねると格子模様が右真横方向に動いているように見えるわけです。なお､この図はガンマ補正のことをまったく考慮していない簡素なつくりですのでご注意を。以下がmatlabのコード。
clear all;close all; M=[]; angledeg=60;anglerad=angledeg*pi/180; freq=15;f=2*pi/freq; [XW,YW]=meshgrid(-63:63,-63:63); gauss=exp(-(XW/32).^2-(YW/32).^2); a1=cos(anglerad)*f; b1=sin(anglerad)*f; for j=1:3 for theta=5:25:1000 [XI,YI]=meshgrid(1-theta:127-theta,1-theta:127-theta); ZI1=gauss.*sin(a1*XI+b1*YI); ZI2=flipud(ZI1); ZI3=(ZI1+ZI2)/2; if j==1, Z=ZI1; elseif j==2, Z=ZI2; elseif j==3, Z=ZI3; end ZR(:,:,1)=(Z+1)/2; ZR(:,:,2)=(Z+1)/2; ZR(:,:,3)=(Z+1)/2; h=imshow(ZR); M=[M;getframe]; end end
これをmovie2avi.mでaviファイルに変換して､それをフリーソフトでGIFアニメに変換しました。(ところで、いつのまにかUNISYSのGIF特許は失効していたのですね。)
いまだにmeshgridの使い方がマスターできていないのでYahoo groupのPsychophysics ToolboxでのDenis Pelli@NYUのリプライを元ネタに改変しました。なお、Denis Pelliは Psychophysics Toolboxの元となったThe VideoToolboxの作者として有名だと思いますが､わたしはQUEST: a Bayesian adaptive psychometric method (pdf file)の著者としてはじめてその名を知りました。
つづきます。