matlab 编程_matlab

clear; clear all;close all; %%%%%灰度化、中值滤波、二值化、（去除孤立噪声）、灰度重心法质心标定 %%%%%%%%%%%%%%%%%%%%%实验图处理用%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%当使用CCD采集到的8位图进行处理时，应使用图像二值化最优阈值迭代，其它均相同 % image=imread('6.bmp');%%%%%%%%%直接读取CCD图片%%%%%%%%%%%%%%%%%%%%%%%%%%%% % gray=image; % figure,imshow(image); % title('原始图片'); % gray=rgb2gray(image)-50; % 彩色图像灰度化 % figure,imshow(gray); % title('原始图片灰度图'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%仿真用%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%直接从数值矩阵读取 load i figure imagesc(i); title('原始图片'); i=i./max(max(i));%%归一化[0,1] %%%%%主要噪声：CCD相机的读出噪声(高斯分布)和光子起伏噪声(泊松分布) i=imnoise(i,'gaussian');%%原始图片加入高斯噪声 i=imnoise(i,'poisson');%%原始图片加入泊松噪声 figure,imshow(i); title('原始图片加入高斯噪声和泊松噪声'); B=mat2gray(i);%%%%彩色图像灰度化 gray=B; figure,imshow(B); title('原始图片灰度图'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%中值滤波%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% d=3; med=medfilt2(gray,[d,d]);%中值滤波，d为中值滤波所选取的模板大小%%%中值滤波去除椒盐噪声 figure,imshow(med); title('第一步中值滤波'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%二值化%%%%%%%%%%%%%%%%%%%%%%%%%%%%实验图处理用%%%%%%%%%%%%%%%%%%%%%%%%%%% % bw=optimumthreshold(med); %图像二值化最优阈值迭代(对图像二值化并且使二值图像光滑) % figure,imshow(bw); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%二值化%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%普通二值化处理（仿真用）%%%%%%%%%% bw=im2bw(med); figure,imshow(bw); title('第二步二值化处理'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%去除孤立噪声%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % bwdn=denoise(bw); %数学形态学远算去除孤立噪声（先腐蚀后膨胀） % figure,imshow(bwdn); % title('第三步去除孤立噪声'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% bwdn=bw; figure,imshow(bwdn);%%此处作图供区域标定使用 title('区域标定'); [L,n]=bwlabel(bwdn); %区域标定 hold on for j=1:n [r,c]=find(L==j); x0=mean(r); y0=mean(c); plot(y0,x0,'Marker','o','MarkerEdgeColor','k','MarkerFaceColor','k','MarkerSize',10); plot(y0,x0,'Marker','*','MarkerEdgeColor','w'); Z(1,1,j)=x0; Z(1,2,j)=y0; % Z=[x0,y0] % X=round(x0); % Y=round(y0); % D=[X,Y] end n %%%%由矩阵Z保存每个点的坐标信息，前提是每个光斑焦点均被正确检测到