Matlab图像处理教程讲义.pptVIP

X为原始样本组成n*d的矩阵，其每一行是一个样本特征向量，每一列表示样本特征向量的一维。 COEFF：主成分分量，也是样本协方差矩阵的本征向量。 SCORE：主成分，X的低维表示。 latent：一个包含着样本协方差矩阵本征值的向量。 X=[1,2;3,3;3,5;4,5;5,6;6,5;8,7;9,3;7,6]; X [COEFF,SCORE,latent]=princomp(X);%主成分分析 COEFF%主成分分量 SCORE%主成分 latent%X样本协方差矩阵的本征值 Sobel交叉梯度 对水平边缘有较大响应 对垂直边缘有较大响应 拉普拉斯算子 I=imread(nir.bmp); I=double(I);%双精度化 w1=[-1 0;0 1]; w2=[0 -1;1 0]; G1=imfilter(I,w1,‘corr’,‘replicate’);%正45°梯度 G2=imfilter(I,w2,‘corr’,‘replicate’);%负45°梯度 G=abs(G1)+abs(G2);%计算Robert梯度 figure;imshow(G,[]); figure;imshow(abs(G1),[]); figure;imshow(abs(G2),[]); Robert交叉梯度 I=imread(nir.bmp); Id=double(I);%双精度化 h_1=fspecial(log,5,0.5);%大小为5，sigma=0.5的LOG算子 I_1=imfilter(Id,h_1,corr,replicate); figure;imshow(uint8(abs(I_1)),[]); h_2=fspecial(log,5,2);%大小为5，sigma=2的LOG算子 I_2=imfilter(Id,h_2,corr,replicate); figure;imshow(uint8(abs(I_2)),[]); 高斯-拉普拉斯锐化 /5、频率域图像增强 用傅里叶变换表示的函数特征可以完全通过傅里叶反变换进行重建而不丢失任何信息 吉布斯现象Gibbs phenomenon（又叫吉布斯效应）： 将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯现象。 幅度谱 相位谱 功率谱 一、傅里叶变换 /5、频率域图像增强 I=fft2(x);%快速傅里叶变换 I=fft2(x,m,n); x为输入图像；m和n分别用于将x的第一和第二维规整到指定的长度。 当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。 I1=abs(I);%计算I的幅度谱 I2=angle(I);%计算I的相位谱 Y=fftshift(I);%频谱平移 I=ifft2(x);%快速傅里叶逆变换 I=ifft2(x,m,n); 幅度谱 相位谱 平移后的相位谱 function out = imidealflpf(I,freq) %imidealflpf函数 构造高斯频域低通滤波器 [M,N] = size(I); out = ones(M,N); for i=1:M for j=1:N if (sqrt(((i-M/2)^2+(j-N/2)^2))freq) out(i,j)=0; end end end 理想低通 function out = imgaussflpf(I,sigma) %imgaussflpf函数 构造高斯频域低通滤波器 [M,N] = size(I); out = ones(M,N); for i=1:M forj=1:N out(i,j) = exp(-((i-M/2)^2+(j-N/2)^2)/2/sigma^2); end end 高斯低通 /6、彩色图像处理 色彩模型：RGB模型、CMY模型、CMYK模型、HIS模型、 HSV模型、YUV模型、YIQ模型。 RGB模型 国际照明委员会（CIE）规定以蓝（435.8nm）、绿（546.1nm）和红（700nm）作为主原色。 RGB_image=cat(3,PR,PG,PB);%将PR、PG、PB三个矩阵在第3个维度上进行级联，进行图像合成 PR=RGB_image(:,:,1);%提取红色分量 PG=RGB_image(:,:,2);%提取绿色分量 PB=RGB_image(:,:,3);%提取蓝色分量 Matlab中一幅