灰度图像处理——纹理分割（Matlab）

主函数：

clc;
clear all;
close all;
%纹理分割测试图像
I=im2double(imread('D:\Gray Files\9-43.tif'));
%获得图像大小
[M,N]=size(I);
%灰度闭操作
r=31;
J_Closing=GrayscaleClosing(I,r);
%灰度开操作
r=60;
J_Opening=GrayscaleOpening(J_Closing,r);
%灰度梯度
n=3;
J_Gradient=GrayscalGradient(J_Opening,n);
J=I+J_Gradient;
imshow(J)

灰度闭操作

GrayscaleClosing函数如下：

%灰度闭操作
% I 输入图像
% n 单一结构元素矩阵的大小
function J_GrayscaleClosing=GrayscaleClosing(I,r) 

    n_B=2*(r-1)+1;
    B=ones(n_B,n_B);
    %将距离大于40的点置零
    for i=1:n_B
        for j=1:n_B
            d=sqrt((i-r)^2+(j-r)^2);
            if d>r
                B(i,j)=0;
            end
        end
    end
    
    ind=find(B==0);
    n_l=floor(n_B/2);
    [M,N]=size(I); 
    J_GrayscaleClosing=zeros(M,N);  
    I_pad=padarray(I,[n_l,n_l],'symmetric');
       
    %膨胀图像
    J_Dilation=zeros(M,N);
    for i=1:M
        for j=1:N
            %获得图像子块区域
            Block=I_pad(i:i+2*n_l,j:j+2*n_l);
            C=Block.*B;
            %删除0值，保留4连通数值
            C=C(:);
            C(ind)=[];
            J_Dilation(i,j)=max(C);
        end
    end
    %对膨胀图像进行扩展
    J_Dilation_pad=padarray(J_Dilation,[n_l,n_l],'symmetric');
    %腐蚀操作
    for i=1:M
        for j=1:N
            %获得图像子块区域
            Block=J_Dilation_pad(i:i+2*n_l,j:j+2*n_l);
            C=Block.*B;
            %删除0值，保留4连通数值
            C=C(:);
            C(ind)=[];
            J_GrayscaleClosing(i,j)=min(C);
        end
    end
end

灰度开操作

GrayscaleOpening函数如下：

%灰度开操作
% I 输入图像
% n 单一结构元素矩阵的大小
function J_GrayscaleOpening=GrayscaleOpening(I,r)
    n_B=2*(r-1)+1;
    B=ones(n_B,n_B);
    %将距离大于40的点置零
    for i=1:n_B
        for j=1:n_B
            d=sqrt((i-r)^2+(j-r)^2);
            if d>r
                B(i,j)=0;
            end
        end
    end
    ind=find(B==0);
    n_l=floor(n_B/2);
    I_pad=padarray(I,[n_l,n_l],'symmetric');
    [M,N]=size(I);

    %-------------------------------灰度开操作---------------------------------
    J_GrayscaleOpening=zeros(M,N);
    %腐蚀操作
    J_Erosion=zeros(M,N);
    for i=1:M
        for j=1:N
            %获得图像子块区域
            Block=I_pad(i:i+2*n_l,j:j+2*n_l);
            C=Block.*B;
            %删除0值，保留4连通数值
            C=C(:);
            C(ind)=[];
            J_Erosion(i,j)=min(C);
        end
    end
    %对腐蚀图像进行扩展
    J_Erosion_pad=padarray(J_Erosion,[n_l,n_l],'symmetric');
    %膨胀图像
    for i=1:M
        for j=1:N
            %获得图像子块区域
            Block=J_Erosion_pad(i:i+2*n_l,j:j+2*n_l);
            C=Block.*B;
            %删除0值，保留4连通数值
            C=C(:);
            C(ind)=[];
            J_GrayscaleOpening(i,j)=max(C);
        end
    end
end

灰度梯度

GrayscalGradient函数如下：

%灰度梯度
% I 输入图像
% n 结构元素大小
function J_Gradient=GrayscalGradient(I,n)
    B=ones(n,n);
    n_l=floor(n/2);
    %对边界图进行扩充，目的是为了处理边界点,这里采用边界镜像扩展
    I_pad=padarray(I,[n_l,n_l],'symmetric');
    [M,N]=size(I);
    J_Erosion=zeros(M,N);
    J_Dilation=zeros(M,N);
    J_Gradient=zeros(M,N);
    for i=1:M
        for j=1:N
            %获得图像子块区域
            Block=I_pad(i:i+2*n_l,j:j+2*n_l);
            C=Block.*B;
            C=C(:);
            %腐蚀操作
            J_Erosion(i,j)=min(C);
            %膨胀操作
            J_Dilation(i,j)=max(C);
            %灰度梯度
            J_Gradient(i,j)=J_Dilation(i,j)-J_Erosion(i,j);
        end
    end
end