北大心理学系.Psychtoolbox2分析.pptx

Psychtoobox;MATLAB 图像刺激; ;图像基础;图像基础——像素;图像基础——颜色;用MATLAB存取并显示图片;a=[ones(400,40),zeros(400,40)]; %先生成一根白条和一根黑条 a1=repmat(a,[1,5]); % 将这个组合复制5遍 a2=xor(a1,a1’);%将这个黑白条纹与自己的转置做异或逻辑运算 imshow(a2) %棋盘格就生成了;实例：红蓝3D与双眼竞争;image与colormap;image与imshow;;clear all face = imread (face.jpg); house = imread (house.jpg); mix (:,:,1) = face(1:680,1:550,1); mix (:,:,2) = house(1:680,221:770,2); mix (:,:,3) = 0; imshow (mix);;image 绘制数据;imagesc colormap;data.mat 中保存了一组虚拟的fMRI数据，数据里记录了V1,V2,V3,V4,V5,ips 6个脑区 100个时间点的值，请画出其6x6 的相关矩阵。; load data.mat fdata=[V1,V2,V3,V4,V5,ips]; imagesc(corr(fdata));马赫带效应;clear all close all colormap(gray(256)) %将颜色查找表设置为灰度图 img = reshape(1:256,16,16); image(img); axis square %将长宽设置为等长 axis off pause for i=1:200 paintpots = rand(256,3); %将颜色查找表设为随机，图片的像素颜色也变成随机 colormap(paintpots); drawnow % 立刻呈现 end;图像的变换与处理;像素亮度的变换操作 (PS中的“曲线”);实例：去除扫描时背面透过来的影;图像的“窗口”;生成高斯窗口;窗口化图像;先生成一个200x200的高斯滤波器，sigma=50，读取Simpson.jpg, 使得二维高斯覆盖在图像上，高斯覆盖范围外值为0，（如3.jpg）。并且将这个过程做成动态。通过鼠标控制高斯覆盖的位置，刷新的频率定为25Hz （既每次呈现一帧图像 要pause 0.04s），点击鼠标右键退出循环。;clear all cdata=imread(Simpson.jpg); [y,x]=size(cdata); [X,Y] = meshgrid(1:x, 1:y); while 1 [x0,y0,button] = GetMouse ; if button(3)~=0 close all break end ga = exp((-1/2/50^2)*((X-x0).^2+(Y-y0).^2)); %ga = ga./max(max(ga)); tmp = double(cdata).* ga; imshow(tmp/255) pause (0.04) end;邻居处理;空间滤波;我们对下幅图片用3x3的平均滤波，即滤波的算子是一个3x3的平均矩阵。设原矩阵名为x（图），x(3,3)的值为74，则经过平均滤波之后，这个点的值变为（20+65+1+102+74+54+58+98+50）/9=58。滤波的过程就是对所有点都做这个计算，得到一个新的值。 ?;高斯模糊;高斯模糊的图片;边缘检测;更多的滤波器;（1）将face进行高斯滤波，高斯算子的边长是100个pixel， sigma为 30。 clear all face = imread (face.jpg); op=fspecial(Gaussian,100,30); filface=imfilter(face,op); imshow(filface);? ? （2）对house进行边缘检测。 clear all house=imread(house.jpg); ghouse=rgb2gray(house); imshow(edge(ghouse,canny))? ? （3） 制作一个200x200 竖直朝向顺时针偏转15度的gabor，其他参数可自己定,并解释各个参数的含义。 ?;clear all [X,Y]=meshgrid(1:200,1:200); grating=sin((tan(75/180*pi)*X+Y)/25); % (tan(75/180*pi)表示竖直朝向顺时偏转15度， % gabor一个周期是25*2*pi/sqrt(tan(75/180*pi)^2+1)个pixel，