clear,
im=imread('1.png');
img=rgb2gray(im);
img=im;
size_im=size(im);
width=size_im(1);
height=size_im(2);
oilpaint=uint8(zeros(width,height)); %生成一个和原图同样大小的矩阵用0填充
radius=5;
intensity_level=150;
img=floor(img*(intensity_level/255));
for x=1:width
for y=1:height
intensity_counter=uint8(zeros(1,intensity_level+1));
r=uint8(zeros(intensity_level+1));
g=uint8(zeros(intensity_level+1));
b=uint8(zeros(intensity_level+1)); %生成一个一维矩阵intensity_counter
for i=(x-radius):(x+radius)
for j=(y-radius):(y+radius)
if i>0 && i<=width &&j>0 && j<=height
intensity =img(i,j)+1;
intensity_counter(intensity)=intensity_counter(intensity)+1;
r(intensity)=im(i,j,1);
g(intensity)=im(i,j,2);
b(intensity)=im(i,j,3); %将原图每个像素点的rgb三个通道的分量值赋给intensity，统计每个颜色出现的次数
end
end
end
[max_ index] =max(intensity_counter); %intensity_counter取最大值赋给index，即重复率最高的颜色
oilpaint(x,y,1)=r(index);
oilpaint(x,y,2)=g(index);
oilpaint(x,y,3)=b(index); %index的rgb通道值赋给最终输出的油画效果图像
end
end
figure,imshow(im,[]);
figure,imshow(oilpaint,[]);
