matlab 素描 code

clear all
clc
I = imread('IMG_4374.JPG');  
imggray = rgb2gray(I);  
N = 255-imggray;     
h1=fspecial('gaussian',3,3);  
g=imfilter(N,h1,'same');  
height=size(I,1);
width=size(I,2);
for i=1:height  
    for j=1:width  
        b = double(g(i,j));  
        a = double(imggray(i,j));  
        temp = a+a*b/(255-b);  
        out(i,j)=min(temp,255);  
    end  
end  

imshow(out/255);  

                        

转自：http://blog.csdn.net/abcd1992719g/article/details/24200551