基于偏微分方的图像处.docVIP

clear; image_I=imread(C22.bmp); subplot(2,2,1); imshow(image_I); matrix_R(:,:,1)=image_I(:,:,1); matrix_R(:,:,2)=0; matrix_R(:,:,3)=0; subplot(2,2,2); imshow(matrix_R); title(R分量); matrix_G(:,:,2)=image_I(:,:,2); matrix_G(:,:,1)=0; matrix_G(:,:,3)=0; subplot(2,2,3); imshow(matrix_G); title(G分量); matrix_B(:,:,3)=image_I(:,:,3); matrix_B(:,:,1)=0; matrix_B(:,:,2)=0; subplot(2,2,4); imshow(matrix_B); title(B分量)  clear; Image=imread(I22_256.bmp); imshow(Image); [m,n]=size(Image); image11=Image; for i=1:m for j=1:n image12(i,j)=image11(i,n-j+1); end end for j=1:n for i=1:m image21(i,j)=image11(m-i+1,j); image22(i,j)=image12(m-i+1,j); end end image1=[image11,image12]; image2=[image21,image22]; image=[image1;image2]; figure imshow(image) clear; Image=imread(I22_256.bmp); Image1=0.5*double(Image); subplot(3,2,1); imshow(Image); subplot(3,2,2); imhist(Image); subplot(3,2,3); imshow(uint8(Image1)); image1=uint8(Image1); imwrite(image1,I22_256_1.bmp) subplot(3,2,4); imhist(image1); Image_hist=histeq(image1); subplot(3,2,5); imshow(Image_hist); imwrite(Image_hist,I22_256_2.bmp) subplot(3,2,6); imhist(Image_hist); function level(image_name,N) %求一图像的水平集、线 %输入参数：image_name---图像文件名 % ：N---阀值(0-255) if nargin1 disp(请输入文件名!!!); return; elseif nargin2 N=128; end image_input=imread(image_name); subplot(2,2,1); imshow(image_input); title(原图像); [size_m,size_n]=size(image_input); matrix_temp=zeros(size_m,size_n); % ?ó?????ˉ for row=1:size_m for col=1:size_n if image_input(row,col) N matrix_temp(row,col)=1; end end end subplot(2,2,2); imshow(matrix_temp,[]); title(图像的水平集); imwrite(matrix_temp,level_setzhan.bmp); %图像矩阵扩展 赋值 便于处理边界 matrix_ex=zeros(size_m+2,size_n+2); for row=1:size_m for col=1:size_n matrix_ex(row+1,col+1)=matrix_temp(row,col); end end %四邻域反填充 得水平线 matrix_new=matrix_temp; for row=2:size_m+1 for col=2:size_n+1 if mat