帆海止水-优快云博客

原创自动阈值法.m

%自动阈值法：Otsu法用MATLAB实现Otsu算法：clc;clear;close;I=imread('e:\role0\003i.bmp'); subplot(1,2,1),imshow(I);title('原始图像')grid on; %显示网格线axis on; %显示坐标系level=graythres

2015-05-12 15:29:24 1389

原创求线性滤波器的频率响应.m

h=fspecial('gaussian');freqz2(h);

2015-05-12 15:10:09 746

原创膨胀操作.m

%膨胀操作clc;clear;close;I=imread('e:\role0\003i.bmp'); %载入图像I1=rgb2gray(I);subplot(1,2,1);imshow(I1);title('灰度图像')grid on; %显示网格线axis on; %显示坐标系se=strel

2015-05-12 15:09:27 792

原创开启和闭合组合操作.m

%开启和闭合组合操作clc;clear;close;%clc清除命令行，clear清除存在内存里的数据，close关闭打开了的文件，I=imread('e:\role0\003i.bmp'); %载入图像subplot(3,2,1),imshow(I);title('原始图像');axis on; %显示坐标系I1=rgb2gray(I

2015-05-12 15:09:04 620

原创双曲线例子.m

plot(x,sin(x),'co',x,cos(x),'g*');;xlabel('input value');ylabel('function value');legend('y=sin(x)','y=cos(x)');grid on

2015-05-12 15:08:40 758

原创生成跳动的红心gif图片.m

clc ;clearx = linspace(-2,2,100);[X,Y,Z] = meshgrid(x,x,x);I1 = (X.^2+9/4*Y.^2+Z.^2-1).^3-X.^2.*Z.^3-9/80*Y.^2.*Z.^3;p = patch(isosurface(X,Y,Z,I1,0));set(p, 'FaceColor', 'red', 'EdgeColor', 'non

2015-05-12 15:06:32 3318

原创 %开启和闭合操作用 MATLAB实现开启和闭合操作

%开启和闭合操作用 MATLAB实现开启和闭合操作clc;clear;close;%clc清除命令行，clear清除存在内存里的数据，close关闭打开了的文件，I=imread('e:\role0\003i.bmp'); %载入图像subplot(2,2,1),imshow(I);title('原始图像');axis on; %显示坐

2015-05-12 15:03:41 3365

原创特殊的二维图形的绘制.m

bar(x,y)comet(x,y)errorbar(x,y,1,u)fplot(x,y)polar(x,y)stairs(x,y)hist(x,y)stem(x,y)fill(x,yfill(x,y)feather(x,y)bar(x,y)compass(x,y)quiver(x,y)

2015-05-12 14:58:44 417

原创梯度算子检测边缘用.m

%梯度算子检测边缘用 MATLAB实现如下： clc;clear;close;I=imread('e:\role0\003i.bmp'); subplot(2,3,1); imshow(I); title('原始图像'); grid on; %显示网格线 axis on; %显示坐标系 I1=im2bw(I);

2015-05-12 14:57:37 910

原创图片的剪切.m

img_ii=imread('apple.jpg');i2=imcrop(img_ii,[30 60 120 160]);%没参数时动鼠标剪切imshow(img_ii);figure,imshow(i2);

2015-05-12 14:56:55 484

原创图像的乘法操作immultiply.m

img_ii=imread('apple.jpg');j=immultiply(img_ii,0.5);subplot(1,2,1);imshow(img_ii);subplot(1,2,2);imshow(j);

2015-05-12 14:56:29 2459

原创添加颜色条.m

%添加颜色条RGB = imread('e:\role0\003i.bmp'); I = rgb2gray(RGB);h = [1 2 1; 0 0 0; -1 -2 -1];I2 = filter2(h,I);imshow(I2,[]), colorbar('vert')

2015-05-12 14:55:38 483

原创计算方框图像在0°到45°方向上的radon变换.m

%计算方框图像在0°到45°方向上的radon变换i=zeros(100,100);i(25:75,25:75)=1;%产生一个正方形的黑框imshow(i);[r,xp]=radon(i,[0 45]);figure;%显示图像在0°方向上的radon变换plot(xp,r(:,1));title('r_{0^o}(x\prime)');figure;%显示图像在45°方向上

2015-05-11 07:25:22 1363

原创计算方框图像从0°~180°每隔1°增加多个radon变换.m

i=zeros(100,100);i(25:85,25:85)=1;theta=0:180;[r,xp]=radon(i,theta);imagesc(theta,xp,r);title('r_{\theta}(x\prime)');xlabel('\theta(degrees)');ylabel('x\prime');set(gca,'xtick',0:20:180);colo

2015-05-11 07:23:50 744

原创灰度线性变换.m

%灰度线性变换 MATLAB 程序实现如下：clc;clear;close;I=imread('e:\role0\003i.bmp'); subplot(2,2,1),imshow(I); title('原始图像'); axis on; %显示坐标系 I1=rgb2gray(I); subplot(2,2,2),imshow(I1); title(

2015-05-11 07:21:34 694

原创基本三维图形的绘制.m

t=0:pi/50:10*pi;x=sin(t);y=cos(t);z=t;h=plot3(x,y,z);set(h,'linewidth',4*get(h,'linewidth'));grid on

2015-05-11 07:20:59 557

原创行程编码解码.m

close all;clear all;clc; %关闭所有图形窗口，清除工作空间所有的变量，清除命令行I1=imread('e:\role0\lena.jpg'); %读入图像I2=I1(:); %将原始图像写成一维的数据并设为I2I2length=length(I2); %计算I2的长度I3=im2bw(I1,0.5); %将原图转换成二值图像阀值为0.5%以下程序为对原图像进行行

2015-05-11 07:18:37 2833 2

原创函数分布的快速绘图.m

fplot('sin(x)/x',[-20 20 -0.4 1.2]);title('fplot of f(x)=sin(x)/x');xlabel('x'),ylabel('f(x)');

2015-05-11 07:18:34 552

原创哈夫曼图片编码.m

clc;clear;close;[Filename,Pathname]=uigetfile(... {'*.tif';'*.bmp';'*.png';'*.gif';'*.jpg'},'选择图片');str=[Pathname Filename];if (Filename~=0) z=imread(str); imshow(z); %else %cle

2015-05-11 07:14:06 675

原创哈夫曼编码.m

clc;close;clear;p=input('please input a number:') %提示输入界面eg:[0.32 0.18 0.05 0.10 0.15 0.2]n=length(p);for i=1:n if p(i)<0 fprintf('\n The probabilities in huffman can not less than 0!

2015-05-11 07:13:00 626

原创哈夫曼编码的MATLAB实现（基于0、1编码）：

%哈夫曼编码的MATLAB实现（基于0、1编码）：clc;clear;A=[0.5,0.19,0.19,0.12];%信源消息的概率序列A=fliplr(sort(A));%按降序排列T=A;[m,n]=size(A);B=zeros(n,n-1);%空的编码表（矩阵）for i=1:n B(i,1)=T(i);%生成编码表的第一列endr=B(i,1)+B(i-1,1

2015-05-11 07:12:08 10776

原创腐蚀操作.m

% 16.腐蚀操作 MATLAB 实现腐蚀操作clc;clear;close;I=imread('e:\role0\003i.bmp'); %载入图像I1=rgb2gray(I);subplot(1,2,1);imshow(I1);title('灰度图像')grid on; %显示网格线axis on;

2015-05-11 07:11:53 812

原创二进制图像求补.m

img_i2=imcomplement(i);subplot(1,2,1);imshow(i);subplot(1,2,2);imshow(img_i2);

2015-05-11 07:06:00 492

原创非线性变换.m

%非线性变换 MATLAB 程序实现如下： I=imread('e:\role0\003i.bmp'); I1=rgb2gray(I); subplot(1,2,1),imshow(I1); title(' 灰度图像'); grid on; %显示网格线 axis on; %显示坐标系J=double(I1); J

2015-05-11 07:05:38 539

原创放大图片.m

img_ii=imread('apple.jpg');img_jj=imresize(img_ii,1.5);imshow(img_ii);figure,imshow(img_jj);

2015-05-11 07:05:01 371

原创对灰度图像进行互补运算.m

j=imread('ice.jpg');m=imcomplement(i);imshow(i),figure,imshow(m);

2015-05-11 07:03:08 559

原创打印和其他选项.m

x=linspace(-2*pi,2*pi,60)y=sin(x).^2./(x+eps)plot(x,y)[a,b]=ginput(8)hold onplot(a,b,'co')hold off

2015-05-11 06:53:59 392

原创创建avi视频sin函数的动画.m

clc;clear;fig=figure;aviobj=avifile('e:\role0\1\21.avi');n=50;x=0:pi/n:2*pi;y=sin(x);k=0;for t=0:pi/n:2*pik=k+1;x(k)=t;y(k)=sin(t);H=plot(x,y,x(k),y(k),'or');gridMOV=getframe(fig);aviob

2015-05-11 06:52:29 648

原创饼图指令pie.m

a=[1,1.6,1.2,0.8,2.1];subplot(1,2,1),pie(a,[1 0 1 0 0]),legend({'1','2','3','4','5'})subplot(1,2,2), b=int8(a==min(a)) ;pie3(a,b)colormap(cool)

2015-05-11 06:49:50 844

原创查找矩阵A中最大元素及下表.m

close all;clear all; clc;%关闭所有图形窗口，清除工作空间所有变量，清除命令行A=magic(5); %查找矩阵A中最大元素及下表a=A(1);for i=2:25 if(A(i)>a) a=A(i); n=i; endend

2015-05-11 06:49:18 698

原创保存gif文件.m

filename=strcat('C:\Users\Public\Videos\Sample Videos\imw\imwrite\',int2str(i),'.','jpg');%将每祯转成jpg的图片搜索clear;clc;figure %新建一张图axis([0 5 0 2])%定义x轴（从0到5）和y轴的范围（从0到2）for i=1:4 if i==1 t

2015-05-11 06:44:43 1009

原创边界跟踪（bwtraceboundary函数）.m

%11.边界跟踪（bwtraceboundary函数）clc; clear; I=imread('e:\role0\003i.bmp');figure;imshow(I);title('原始图像');I1=rgb2gray(I); %将彩色图像转化灰度图像threshold=graythresh(I1); %计算将灰度图像转化为二值图像所

2015-05-11 06:42:48 6864 5

原创 n个图片生成gif图片.m

save_file='c:\Users\Public\Videos\Sample Videos\imw\imwrite\123out.gif';path='c:\Users\Public\Videos\Sample Videos\imw\imwrite\';path2='*.jpg';s=strcat(path,path2);files=dir(fullfile(s));len=leng

2015-05-10 20:43:30 738 1

原创 ch10_2_4.m

[f1,f2] = freqspace(25,'meshgrid');Hd = zeros(25,25); d = sqrt(f1.^2 + f2.^2) < 0.5;Hd(d) = 1;mesh(f1,f2,Hd)

2015-05-10 20:32:49 367

原创 ch10_2_3窗口法.m

%窗口法%利用函数find1Hd = zeros(11,11); Hd(4:8,4:8) = 1;[f1,f2] = freqspace(11,'meshgrid');mesh(f1,f2,Hd);axis([-1 1 -1 1 0 1.2]), colormap(jet(64))h = fwind1(Hd,hamming(11));figure, freqz2(h,[

2015-05-10 20:32:20 489

原创 ch10_2_2频率采样法.m

%频率采样法Hd = zeros(11,11); Hd(4:8,4:8) = 1;[f1,f2] = freqspace(11,'meshgrid');mesh(f1,f2,Hd), axis([-1 1 -1 1 0 1.2]), colormap(jet(64))h = fsamp2(Hd);figure, freqz2(h,[32 32]), axis([-1 1 -1 1 0 1

2015-05-10 20:25:01 496

原创 ch10_2_1利用频率变换法生成一个2-D滤波器.m

%利用频率变换法生成一个2-D滤波器b = firpm(10,[0 0.4 0.6 1],[1 1 0 0]);h = ftrans2(b);[H,w] = freqz(b,1,64,'whole');%得到b的1-D频率响应colormap(jet(64))%图形颜色映射plot(w/pi-1,fftshift(abs(H)))figure, freqz2(h,[32 32])

2015-05-10 20:23:32 568

原创 5－1－1－3 采用对数形式的变换函数进行动态范围压缩源程序：

%%%% 5－1－1－3 采用对数形式的变换函数进行动态范围压缩源程序：clear all;X1=imread('E:\role0\001i.bmp');figure,imshow(X1)c=255/log(256);x=0:1:255;y=c*log(1+x);figure,plot(x,y)axis tight,xlabel('f'),ylabel('g')title(

2015-05-10 20:18:44 1218

原创 ch5_1_1 线性变换的matlab源程序：.m

%%%%% 5－1－1－1 线性变换的matlab源程序：X1=imread('E:\role0\001i.bmp');figure,imshow(X1);f0=0;g0=0;f1=10;g1=10f2=180;g2=1800;f3=255;g3=255;figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3])axis tight,xlabel('f'),

2015-05-10 20:17:23 1728

原创 5－1－1－2 分段线性变换的源程序：

%%%% 5－1－1－2 分段线性变换的源程序：clear all;X1=imread('E:\role0\001i.bmp');figure,imshow(X1)f0=0;g0=0;f1=20;g1=10;f2=180;g2=230;f3=255;g3=255;figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3])axis tight,xlabe

2015-05-10 20:17:12 1043