北大心理学系.Psychtoolbox1重点分析.pptx

Psychtoobox Lihan Chen 2015.4.28 MATLAB 图像刺激 图像处理 图像基础 图像的变换与操纵 计算生成图像 图像基础 像素、颜色、索引图像 图像基础——像素 数字图像 离散的坐标：像素 离散的强度：0~255 (黑~白) 图像基础——颜色 RGB颜色系统 每个像素的红、绿、蓝均由uint8表示 (0~255) alpha通道：透明度 24位/32位真彩色 imtool( photo.jpg); 用MATLAB存取并显示图片 读取并显示图片 im = imread(tmp.png); size(im) imshow(im) 写入图片 imwrite(im, tmp.jpg); a=[ones(400,40),zeros(400,40)]; %先生成一根白条和一根黑条 a1=repmat(a,[1,5]); % 将这个组合复制5遍 a2=xor(a1,a1’);%将这个黑白条纹与自己的转置做异或逻辑运算 imshow(a2) %棋盘格就生成了 实例：红蓝3D与双眼竞争 im1 = imread(Koala.jpg); im2 = imread(Penguins.jpg); % 用第二幅图替代第一幅图的红色通道 im1(:,:,1) = im2(:,:,1); % 把绿色通道清零 im1(:,:,2) = 0; imshow(im1) image与colormap image函数是MATLAB最基本的图像显示函数，可以绘制索引图像，即每个像素的值对应颜色查找表中的索引 colormap：定义图像显示用的颜色查找表 imagesc将数据scale后绘制成图（例如绘制相关矩阵） image与imshow Imshow仅用于显示由RGB或灰度值定义的图像（image也可以） 无论是哪个函数，若图像是以uint8表示的，则取值范围为0~255，若以double表示，则取值范围是0~1 image 绘制数据 imagesc colormap 绘制相关矩阵 a = rand(10); for i = 2:10 a(:,i) = a(:,i-1)+rand(10,1)*0.5; end imagesc(corr(a)) colorbar colormap autumn 图像的变换与处理 调整曲线、窗口化 像素亮度的变换操作 (PS中的“曲线”) 实例：去除扫描时背面透过来的影 im = imread(book-000123.png); g = rgb2gray(im); d = double(g); y = uint8(d*1.7-85); imshow(y); imwrite(y, out.png); 图像的“窗口” 操纵每个像素的亮度：乘以一个系数 生成高斯窗口 窗口化图像 imagesc(GaWindow) tmp = double(im).*GaWindow; Windowed = uint8(tmp); imshow(Windowed) 邻居处理 高斯模糊、边缘检测 空间滤波 定义一个矩阵(称为滤波器filter) 将该矩阵覆盖在每个像素点上 对所有被覆盖的点 将原图像矩阵和该矩阵的对应点相乘 求和变成新图像的点 由邻居决定新的像素点的强度 高斯模糊 改变窗口大小，令x = y = 10, stdev = 5，生成一个“高斯窗口” 让窗口矩阵的和等于1 GaFilter = GaWindow/sum(sum(GaWindow)); 用imfilter或filter2或conv2对原图进行滤波 imf = imfilter(im, GaFilter); imshow(imf) 高斯模糊的图片 边缘检测 原始图像 gray edge(gray) imfilter(gray, fspecial(sobel)) 更多的滤波器 用fspecial生成二维的特殊滤波器 average 简单平均 gaussian 高斯滤波 motion 模拟摄像机抖动 sobel 检测边缘算子 doc fspecial 计算生成图像 正弦光栅, Gabor Patch, … 生成正弦光栅 [X, Y]=meshgrid(1:50,1:100); imshow(sin(X)) imshow(sin(X+Y)) imshow(sin((X+Y)/2)) imshow(sin((X+2*Y)/2)) imshow(0.5+0.5*sin((3*X+4*Y)/6)) 完整代码 [x,y]=meshgrid(1:Boundx, 1:Boundy); % define the costant for the cosine wo=SAngle/180*pi; % conversion from deg to rad