pro read_hr_ha_new,fname,data,header

data=readfits(fname,header)
xdim=sxpar(header,'naxis1')
ydim=sxpar(header,'naxis2')
dxx=8
xdim=xdim-dxx
dyy=36
ydim=ydim-dyy
data=data(0:xdim-1,dyy:ydim+dyy-1)
data=data-1.0
header=['The dimension of the image  is changed                    ',$
	'NAXIXX  = 2008                                            ',$
	'NAXIYY  = 2008                                            ',$
	header]

end







