Roberts滤波器

Halcom

1. function im1 = Roberts_fspecial(im,type)
   
2.     if nargin < 2
   
3.         type = 'Roberts';
   
4.     end
   
5.     if ~isa(im,'double')
   
6.         im = double(im)/255;
   
7.     end
   
8.     [a, b] = size(im(:,:,1));  %行＼列
   
9.     im1(:,:,1) = zeros(a,b);  % R1
   
10.    switch type
    
11.        case 'Roberts'  % Roberts filter
    
12.           for i = 1:size(im,3)  % 矩阵的维数，3D（RGB图像） 或者 2D（灰度图像）
    
13.               for j=1:a-1  %　行数
    
14.                   for k=1:b-1 % 列数
    
15.                       im1(j,k,i) = abs( im(j+1,k+1,i)-im(j,k,i) ) + abs( im(j+1,k,i)-im(j,k+1,i) );
    
16.                   end
    
17.               end
    
18.           end
    
19.    end
    
20. end

复制代码