function Out=unsharp_mask(Ima,k,sigma)

G=fspecial('gaussian',6*sigma,sigma);
I_edges=Ima-imfilter(Ima,G);
Out = Ima + k*I_edges;
