pro pstat, filename
  sdssproc, filename, im, hdr=hdr

  djs_iterstat, im, mean=mn

  print, filename, mn, hdr['HARTMANN'], hdr['EXPTIME']
end