8.6.2 reimage(image,data,xy,dxy,zero_offset,xy_scale,ii,jj,kk)

reimage $<$ f1 $<$ BoaF1.f90

! NAM: reimage (subroutine)
! HIS: FB040510 created
! DES: fill an image array with data in pixel boxes of given size
!        image = reimage(image,data,xy,dxy,zero_offset,xy_scale)
! INP: 
!        image       : input rank-2 array('f') with bounds (ii,jj)
!        data        : input rank-1 array('f') with bounds (kk)
!        xy          : input rank-2 array('f') with bounds (2,kk)
!        dxy         : input rank-2 array('f') with bounds (2,kk)  -- half width!
!        zero_offset : input rank-1 array('f') with bounds (2)
!                      bottom left corner user coordinate (eg arcsec) offset.
!        xy_scale    : input rank-1 array('f') with bounds (2)
!                      image scale, pixel per user coordinate (eg pix/arcsec)
! OUT:  
!        image       : rank-2 array('f') with bounds (ii,jj)
! USE: Example:
!        from Numeric import *
!        from fortran import f90
!        im = zeros((11,11),'f')  # odd size since coordinates refer to pix center
!        data    = array([10,20,30],'f')
!        xyscale = array([0.1,0.1],'f')
!        xy      = array([ [-50,0,20],[-30,0,20] ],'f')
!        dxy     = array([ [11,6,11] , [11,6,11] ],'f')
!        zerooffset = array([-50,-50],'f') # blc pix refers to -50,-50
!        im = f90.f1.reimage(im,data,xy,dxy,zerooffset,xyscale)
!        print im