数字图像处理实验报告(4)

4. RGB 彩色图像分割

function [C,m]=covmatrix(X) [K,n]=size(X); X=double(X); if n==1 C=0; m=x; else

m=sum(X,1)/K;

X=X-m(ones(K,1),:); C=(X'*X)/(K-1); m=m'; end

function I=colorseg(varargin) f=varargin{2};

if(ndims(f)~=3)|(size(f,3)~=3) error('Input image must be RGB!'); end

M=size(f,1); N=size(f,2); [f,L]=imstack2vectors(f); f=double(f); I=zeros(M*N,1); T=varargin{3}; m=varargin{4};

m=m(:)';%make sure that m is a row vetor if length(varargin)==4 method='euclidean';

elseif length(varargin)==5 method='mahalanobis'; else

error('wrong number of inputs!'); end

switch method case 'euclidean' p=length(f);

D=sqrt(sum(abs(f-repmat(m,p,1)).^2,2)); case 'mahalanobis' C=varargin{5};

D=mahalanobis(f,C,m); otherwise

error('Unknown segmentation method!'); end

J=find(D<=T); I(J)=1;

I=reshape(I,M,N);

function [X,R]=imstack2vectors(S,MASK) [M,N,n]=size(S); if nargin==1

MASK=true(M,N); else

MASK=MASK~=0; end

[I,J]=find(MASK); R=[I,J]; Q=M*N;

X=reshape(S,Q,n);

MASK=reshape(MASK,Q,1); X=X(MASK,:);

%%%%%%RGB 图像分割主程序

clear all

clc

f=imread('D:\\Fig0630(01)(strawberries_fullcolor).tif'); mask=roipoly(f);

subplot(2,3,1),imshow(mask)%mask 是一幅2 值图像（大小与f 相同） red=immultiply(mask,f(:,:,1)); subplot(2,3,2),imshow(red) title('red 分量');

green=immultiply(mask,f(:,:,2)); subplot(2,3,3),imshow(green) title('green 分量');

blue=immultiply(mask,f(:,:,3)); subplot(2,3,4),imshow(blue) title('blue 分量');

g=cat(3,red,green,blue); subplot(2,3,5),imshow(g)

%计算ROI 中的点的均值向量和协方差矩阵 [M,N,K]=size(g); I=reshape(g,M*N,3); idx=find(mask);

I=double(I(idx,1:3)); [C,m]=covmatrix(I); d=diag(C); sd=sqrt(d)'

E25=colorseg('euclidean',f,25,m);

%colorseg('m',f,T,m)T 值可取的值为25 的倍数时得到的不同结果

实验六小波变换

一、 实验目的

1 Haar、尺度和小波函数；

2 比较函数wavefast 和函数wavedec2 的执行时间； 3 小波的方向性和边缘检测。

二、 实验步骤

1.1. Haar、尺度和小波函数

[Lo_D,Hi_D,Lo_R,Hi_R]=wfilters('haar')%Haar 变换的分解和重构滤波器的长度 为2

waveinfo('haar');

[phi,psi,xval]=wavefun('haar',10); xaxis=zeros(size(xval));

subplot(121);plot(xval,phi,'k',xval,xaxis,'--k'); axis([0 1 -1.5 1.5]);axis square; title('haar scaling function');

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库数字图像处理实验报告(4)在线全文阅读。