from herschel.pacs.spg.phot import PhotHelper

for field in fields:
  # Load the maps
  redMapCombined = FitsArchive().load(dataDir + "field" + str(field) + "-red-corr.fits")
  try:
    blueMapCombined = FitsArchive().load(dataDir + "field" + str(field) + "-blue-corr.fits")
  except:
    blueMapCombined = FitsArchive().load(dataDir + "field" + str(field) + "-green-corr.fits")
  # Create the rgb maps
  rgbMap = PhotHelper.createRgbMap(redMapCombined, blueMapCombined, maxSignal, exponent, multiplicativeFactor)
  d = Display(rgbMap)
  d.saveAsJPG(dataDir + "field" + str(field) +".jpg")
  d.saveAsPNG(dataDir + "field" + str(field) +".png")