/* Auth:                                                        minas */
/* Date:                                 Wed Oct 26 09:26:09 EDT 2016 */
/*--Break--*/
im=read_img("dsc01482.pgm");
/*--Break--*/
gshow(im);
/*--Break--*/
im=to_fimg(im);
/*--Break--*/
img=im(*)mk_gauss_tmpl(1.5); 
/*--Break--*/
gshow(img);
/*--Break--*/
imx=D_x(img);
/*--Break--*/
gshow(imx,:rescale=t);
/*--Break--*/
imy=D_y(img);
/*--Break--*/
gshow(imy,:rescale=t);
/*--Break--*/
imsq=imx*imx+imy*imy;
/*--Break--*/
gshow(imsq,:rescale=t);
/*--Break--*/
gshow(glog(imsq+20),:rescale=t);
/*--Break--*/
gshow(imsq>30,:rescale=t);
/*--Break--*/
gg=mk_gauss_tmpl(1.5);
/*--Break--*/
Exx= (imx*imx)(*)gg;
Exy= (imx*imy)(*)gg;
Eyy= (imy*imy)(*)gg;
/*--Break--*/
ExxEyy = Exx+Eyy;
Exx_Eyysq = Exx-Eyy;
Exx_Eyysq *= Exx_Eyysq;
/*--Break--*/
lam1 = (ExxEyy+gsqrt(Exx_Eyysq+4*Exy*Exy))/2.0;
lam2 = (ExxEyy-gsqrt(Exx_Eyysq+4*Exy*Exy))/2.0;
/*--Break--*/
gshow(lam1,:rescale=t);
/*--Break--*/
gshow(lam2,:rescale=t);
/*--Break--*/
