Laplacian filter in IDL

;####################################################
;                       Laplacian filter
;                  Author : S.Ganesh Babu
;####################################################

PRO EX13

file = DIALOG_PICKFILE(PATH='E:\Image')
im = READ_IMAGE(file)
sz=SIZE(im,/DIMENSIONS)
nv=SIZE(sz,/DIMENSIONS)
nim=im
ker=[[0,1,0],[1,-4,1],[0,1,0]]



IF nv EQ 2 THEN BEGIN
nim=CONVOL(im,ker)
WINDOW,0,XSIZE=sz[0],YSIZE=sz[1],TITLE='INPUT IMAGE'
TVSCL,im
WINDOW,1,XSIZE=sz[0],YSIZE=sz[1],TITLE='Laplacian filtered'
TVSCL,nim
ENDIF ELSE BEGIN

WINDOW,0,XSIZE=sz[1],YSIZE=sz[2],TITLE='INPUT IMAGE'
TVSCL,im,TRUE=1
WINDOW,1,XSIZE=sz[1],YSIZE=sz[2],TITLE='Laplacian filtered'
TVSCL,nim,TRUE=1
ENDELSE


END