FUNCTION it_square,array

  sz=SIZE(array)
  xcen=sz(1)/2
  ycen=sz(2)/2
  image=array(xcen-ycen:xcen+ycen-1,0:sz(2)-1)

RETURN,image
END