这个算法原理很简单，对彩色图像的R,G,B 三个通道，分别求梯度，然后将梯度值作为三个通道的值。
clc;
clear all;Image=imread('4.jpg');
Image=double(Image);
R=Image(:,:,1);
G=Image(:,:,2);
B=Image(:,:,3);
R_Gradient=Find_Gradient(R);
G_Gradient=Find_Gradient(G);
B_Gradient=Find_Gradient(B);
Image_new(:,:,1)=R_Gradient;
Image_new(:,:,2)=G_Gradient;
Image_new(:,:,3)=B_Gradient;
imshow(Image_new/255);
function I_Gradient=Find_Gradient(I)
p=3;
q=0;
Gradient_Y=[-1 -p -1
0 q 0
1 p 1];
Gradient_X=[-1 0 1
-p q p
-1 0 1];
Image_y=imfilter(I, Gradient_Y, 'conv');
Image_x=imfilter(I, Gradient_X, 'conv');
alpha=0.5;
I_edge=abs(Image_x)*alpha+abs(Image_y)*(1-alpha);
I_Gradient=I_edge;
原图：
效果图：
