PRO dostack_header,mapHeader,imgStack,stackHeader

extast, mapHeader,astr

make_astr,astr_stk,CD=astr.CD,DELT=astr.CDELT,$
  CRPIX=(size(imgStack,/dimension)-1)/2+1,$
  CRVAL=[0.,0.],CTYPE=astr.CTYPE,$
  LATPOLE=astr.LATPOLE,LONGPOLE=astr.LONGPOLE

mkhdr,stackHeader,imgStack
putast,stackHeader,astr_stk

END