Robert算子原理（matlab）

%Roberts
 N1 = 10; 
 A = zeros(1,N1); 
 X = zeros(1,N1);
 tic
 for L=1: N1 
 I=imread([int2str(L),'.bmp']); 
 I=double(I); 
 [M N]=size(I); 
 FI=0; 
 %Robert算子原理，对角方向相邻的两像素之差 
 for x=1:M-1 
     for y=1:N-1 
         FI= FI   (abs(I(x,y)-I(x 1,y 1)) abs(I(x 1,y)-I(x,y 1))); 
     end 
 end 
 A(1,L) = FI;  
 end 
 time=toc
 
  for W = 1:N1 
     C = max(A); 
     D = min(A); 
     E = C-D; 
     R = (A(1,W) - D)/(E); 
     X(1,W) = R; 
 end 
x1=[-50 -40 -30 -20 -10 0 10 20 30 40 ]; 
y1 = [X(1,1) X(1,2) X(1,3) X(1,4) X(1,5) X(1,6) X(1,7) X(1,8) X(1,9) X(1,10)]; 
[p,S]=polyfit(x1,y1,2); 
Y=polyconf(p,x1,y1); 
plot(x1,y1,'c'); 
hold on;

资源下载

此资源下载价格为5D币（VIP免费），请先登录

资源文件列表

luowenqingxidupinjia1.m , 701