FUNCTION resize,image,img_resize,head

  IF N_ELEMENTS(head) EQ 0 THEN MKHDR,head,image
  sz = size(image) ; *** updating size of image *** 

  nm = sz(1) < sz(2)  ; *** circle must stay circle ***

  CASE sz(3) OF
   1: tmp = BYTARR(nm,nm)
   2: tmp = INTARR(nm,nm)
   3: tmp = LONARR(nm,nm)
   4: tmp = FLTARR(nm,nm)
   ELSE: donothing = 0
  ENDCASE
  tmp = image(0:nm-1,0:nm-1)
  image = REBIN ( tmp, img_resize, img_resize )

  sz = size(image)
RETURN,image
END