lum.m

DigitalImageProcessing_base_on_Matlab 基于Matlab的数字图像处理

%函数LUM.m: LUM 滤波器%
function [O_data]=LUM(I_data,k,l,w)
%w=input('please input window size (5 or 7):') 模板大小%
%k=input('please input k parameter:') 平滑参数%
%l=input('please input l parameter:') 锐化参数%

I_data1=I_data; 
w1=(w-1)/2;
%将原始图像依照模板的尺寸扩大%
[M,N]=size(I_data);
A=[];
A=[I_data(:,1:w1) I_data I_data(:,N-w1+1:N)];
A=[A(1:w1,:);A;A(M-w1+1:M,:)];
I_data=A;

%进行LUM滤波%
for i=1:1:M
   for j=1:1:N
      block=I_data(i:i+w-1,j:j+w-1);   %取一个W×W图像模板进行滤波
      b=sort(block(:));                %将此模板的内容依照大小排序