利用另外一种算法完成素描特效的生成。

%%% Sketch

clc;

clear all;

Image=imread('4.jpg');

Image=double(Image);

[row,col,layer]=size(Image);



Filter_size=30;

sigma=Filter_size/6;

F_gaussian=fspecial('gaussian', Filter_size, sigma);



for i=1:row

    for j=1:col

        I(i,j)=max(Image(i,j));

    end

end

I=I/255;

p=1;

Gradient_X=[-1 -p -1

            0   0  0

            1   p  1];

        

Gradient_Y=[-1 0 1

            -p 0 p

            -1 0 1];

Image_y=imfilter(I, Gradient_Y, 'conv');        

Image_x=imfilter(I, Gradient_X, 'conv');



alpha=0.5;

% Image_edge=(abs(Image_x)*alpha+(1-alpha)*abs(Image_y));

Image_edge=(abs(Image_x)+abs(Image_y));



F_image=imfilter(Image,F_gaussian,'conv');

[L,a,b]=Rgb2Lab(F_image);

I=sqrt(L.*L+a.*a+b.*b);

I_mean=mean(mean(I));

F_S=I;

for i=1:row

    for j=1:col

        if(I(i,j)<I_mean)

            F_S(i,j)=(1-I(i,j)/I_mean).^2;

        else

            F_S(i,j)=0;

        end

    end

end



Image_out=(1-Image_edge).*(0.8-F_S);



%% 设置底色层

    Image_Dodge(:,:,1)=Image_out;

    Image_Dodge(:,:,2)=Image_out;

    Image_Dodge(:,:,3)=Image_out;

    Base_layer=Image_out;

    Base_layer(:,:,1)=210/255;

    Base_layer(:,:,2)=225/255;

    Base_layer(:,:,3)=105/255;

   alpha=0.7;

   Image1=alpha*Image_Dodge+(1-alpha)*Base_layer;

   Image3=Image_Dodge.*Image1;

   figure, imshow(Image3);

 



%%%  Rgb2Lab

%%%  

function [L,a,b]=Rgb2Lab(Image)

Image=double(Image)/255;

R=Image(:,:,1);

G=Image(:,:,2);

B=Image(:,:,3);



X=0.5767309*R + 0.1855540*G + 0.1881852*B;

Y=0.2973769*R + 0.6273491*G + 0.0752741*B;

Z=0.0270343*R + 0.0706872*G + 0.9911085*B;

[row, col]=size(R);



L=R;

a=R;

b=R;



for i=1:row

    for j=1:col        

        L(i,j)=116*F_Correction(Y(i,j))-16;

        a(i,j)=500*(F_Correction(X(i,j))-F_Correction(Y(i,j)));

        b(i,j)=200*(F_Correction(Y(i,j))-F_Correction(Z(i,j)));

    end

end

原图

效果图