%函數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(:));                %將此模板的內容依照大小排序