随笔分类 - Matlab练习程序
摘要：cl;H=1; %索引pix中第一个元素，即高度W=2; %索引pix中第二个元素，即宽度jiaodu=45; %要旋转的角度，旋转方向为顺时针img=imread('Corner.png'); %这里v为原图像的高度，u为原图像的宽度imshow(img); %这里y为变换后图像的高度，x为变换后图像的宽度[v u]=size(img);theta=jiaodu/180*pi;rot=[cos(theta) ...
阅读全文
摘要：cl;delta_x=67; %要求为整数，水平移动的偏移量，正为向右，负为向左delta_y=-89; %要求为整数，竖直移动的偏移量，正为向下，负为向上img=imread('Corner.png'); %这里v为原图像的高度，u为原图像的宽度imshow(img); %这里y为变换后图像的高度，x为变换后图像的宽度[v u]=size(img);imgn=zeros(v+abs(delta_y),u+abs(delta_x));rot=[1 0 0;0 1 0;delta_y delta_x 1]; in...
阅读全文
摘要：cl;w=0.6; %放大或缩小的宽度h=1.4; %放大或缩小的高度img=imread('Corner.png');imshow(img);[m n]=size(img);imgn=zeros(h*m,w*n);rot=[h 0 0;0 w 0;0 0 1]; %变换矩阵x=h*u,y=w*vinv_rot=inv(rot);for x=1:h*m for y=1:w*n pix=[x y 1]*inv_rot; ...
阅读全文
摘要：cl;w=0.7; %放大或缩小的宽度h=1.5; %放大或缩小的高度img=imread('Corner.png');imshow(img);[m n]=size(img);imgn=zeros(h*m,w*n);rot=[h 0 0;0 w 0;0 0 1]; %变换矩阵x=h*u,y=w*vfor u=1:m for v=1:n pix=[u v 1]*rot; if (pix(1)>=0.5 && pix(2)>=0.5) ...
阅读全文
摘要：cl;bw=[eye(3) eye(3);eye(3) eye(3)]; %要处理的矩阵f=eye(3); %模板矩阵,其实是需要使用到的逻辑lut=makelut(@(nbood) isequal(nbood,f),3); %生成查表向量，将f矩阵这样的逻辑应用到3*3的nbood上bw1=applylut(bw,lut); %应用查表向量。方法：遍历bw所有pixel,如果哪个pixel和f完全相同的,则bw1中相应位置的pixel置1，否则置0bw,f,bw1 ...
阅读全文
摘要：cl;img1=imread('gaosi.jpg');[m,n]=size(img1);w=fspecial('gaussian',[3 3]);img2=imresize(imfilter(img1,w),[m/2 n/2]);img3=imresize(imfilter(img2,w),[m/4 n/4]);img4=imresize(imfilter(img3,w),[m/8 n/8]);img5=imresize(imfilter(img4,w),[m/16 n/16]);imshow(img1);figure,imshow(img2);figure,
阅读全文
摘要：cl;img=imread('mask.jpg');imshow(img);[x,y]=size(img);img_man=zeros(x,y);img_com=zeros(x,y);%% 直方图均衡化算法Max=max(max(img));Min=min(min(img));Hist=zeros(1,256);for i=1:x for j=1:y Hist(img(i,j)+1)=Hist(img(i,j)+1)+1; endendfigure,plot(Hist);p=zeros(1,256);for i=1:256 p(i)=Hist(i)/(x*y);...
阅读全文
摘要：cl;img=imread('fupeng.jpg');z=3;[x y]=size(img);imshow(img);img=double(img);img_cai=zeros(x,y,z);Max=max(max(img));Min=min(min(img));img=(255/(Max-Min))*img-(255*Min)/(Max-Min);figure,imshow(uint8(img));r=1;g=2;b=3;for i=1:x for j=1:y temp=(2*pi/(Max-Min))*img(i,j)-(2*pi*Min)/(Max-Min); ...
阅读全文
摘要：cl;img=imread('TestDatabase/1.bmp');imshow(img);img=double(img);figure,imshow(mat2gray(img))ma=double(max(max(img)));mi=double(min(min(img)));img=(255/(ma-mi))*img-(255*mi)/(ma-mi);img=uint8(img);figure,imshow(img)
阅读全文
摘要：cl;m=31;n=31;img=zeros(m+1,n+1);img=double(img);pi=3.1415926;sigma=10;for i=-(m/2):m/2 for j=-(n/2):n/2 img(i+m/2+1,j+n/2+1)=(1/(2*pi*sigma*sigma))*exp(-(i*i+j*j)/(2*sigma*sigma)); endendimg=mat2gray(img);imshow(img);imwrite(img,'pic.bmp');%一个更好的实现，2012.12.18修改h=5;w=2;[x y]=meshgrid(-w...
阅读全文
摘要：cl;raw=zeros(200,256,30);for i=1:30 filename=strcat('F:\算法实验\data\seq3\',int2str(i),'.bmp'); raw(:,:,i)=imread(filename);end方法二：cl;img=cell(1,30);for i=1:30 filename=strcat('F:\算法实验\data\seq3\',int2str(i),'.bmp'); img{i}=imread(filename);end
阅读全文
摘要：cl;img=imread('1.bmp');[m n]=size(img);imshow(img)imgn=img;temp=[];%开操作for i=2:m-1 %腐蚀 for j=2:n-1 temp=img(i-1:i+1,j-1:j+1); imgn(i,j)=min(temp(:)); endendfor i=2:m-1 %膨胀 for j=2:n-1 temp=img(i-1:i+1,j-1:j+1); imgn(i,j)=max(...
阅读全文
摘要：cl;img_gray=imread('fupeng.jpg');img_erzhi=imread('erzhi_fupeng.jpg');imshow(img_gray)figure,imshow(img_erzhi)[m n]=size(img_gray);img_gray_fu=zeros(m,n);img_gray_peng=zeros(m,n);img_erzhi_fu=zeros(m,n);img_erzhi_peng=zeros(m,n);for i=2:m-1 for j=2:n-1 img_gray_fu(i,j)=min(min(img_gr
阅读全文
摘要：cl;img=imread('1.bmp');img=double(img);imshow(mat2gray(img));[m n]=size(img);imgn=zeros(m-3,n-3);temp=[];for i=1:m-3 for j=1:n-3 temp=img(i:i+3,j:j+3); imgn(i,j)=median(temp(:)); endendfigure,imshow(mat2gray(imgn))imgn=img(1:m-3,1:n-3)-imgn;figure,imshow(mat2gray(imgn))imgn=i...
阅读全文
摘要：clear;clc;width=3;xwidth=(width-1)/2;imgn=imread('1.bmp');imshow(imgn,[]);imgn=double(imgn);[m n]=size(imgn);imgn1=imgn;z=zeros(4,width);tem=1;for i=1+xwidth:1:m-xwidth for j=1+xwidth:1:n-xwidth for k=-xwidth:1:xwidth z(1,tem)=imgn(i,j+k); z(2,tem)=imgn(i+k,j)...
阅读全文
摘要：clear;clc;width=3;xwidth=(width-1)/2;imgn=imread('1.bmp');imshow(imgn,[]);imgn=double(imgn);[m n]=size(imgn);imgn1=imgn;z=zeros(4,width);tem=1;for i=1+xwidth:1:m-xwidth for j=1+xwidth:1:n-xwidth for k=-xwidth:1:xwidth z(1,tem)=imgn(i,j+k); z(2,tem)=imgn(i+k,j)...
阅读全文
摘要：cl;img=imread('3.17.tif');imgn=imnoise(img,'salt & pepper',0.02);[m n]=size(img);%h=fspecial('average',[3 3]); % 3*3的均值滤波%imgn1=imfilter(imgn,h,'replicate');%imgn=medfilt2(imgn,[3 3]); % 3*3的中值滤波%%%%自己编写3*3均值滤波%{tem=0;for i=2:1:m-1 for j=2:1:n-1 for p=i-1:1:i+1 ..
阅读全文
摘要：cl;img=imread('15.bmp');%img=double(img);f=fft2(img); %傅里叶变换f=fftshift(f); %使图像对称r=real(f); %图像频域实部i=imag(f); %图像频域虚部margin=log(abs(f)); %图像幅度谱，加log便于显示phase=log(angle(f)*180/pi); %图像相位谱l=log(f); subplot(2,2,1),imshow(img),title('源图像');%subplot(2,2,2)...
阅读全文
摘要：cl;img=imread('15.bmp');%imshow(img);[m n]=size(img);max=0;min=256;avg=0;for i=1:1:m for j=1:1:n if img(i,j)<min min=img(i,j); end if img(i,j)>max max=img(i,j); end avg=avg+double(img(i,j)); endend%avg=mean(mean(img));avg=uint8(avg/(m*n)...
阅读全文