from herschel.pacs.spg.phot import PhotHelper
from herschel.ia.dataset import StringParameter
from herschel.ia.dataset import LongParameter

"""
  Load the red map
"""
cameraName = "red"
fitsFile = dataDir + "map-" + str(obsids[0]) +  "-" + str(obsids[1]) + "-" + cameraName + ".fits"
redMap = FitsArchive().load(fitsFile)

"""
  Modifiy the metadata
"""
meta = redMap.meta
meta.remove("obsid")

for i in range(len(obsids)):
  meta.set("obsid00"+str(i+1), LongParameter(obsids[i]))

meta.set("associatedObsids", StringParameter(str(borderObsids)[1:-1], "Other associated obsids that contributed to the map"))

"""
 Make sure the ra and dec metaData information is correct
"""
redMap = centerRaDecMetaData(redMap)

"""
  Save the map
"""
fitsFile = dataDir + str(obsids[0]) +  "-" + str(obsids[1]) + "_JScanamHPDP_" + cameraName + ".fits"
FitsArchive().save(fitsFile, redMap)


"""
  Load the blue map
"""
try:
  cameraName = "blue"
  fitsFile = dataDir + "map-" + str(obsids[0]) +  "-" + str(obsids[1]) + "-" + cameraName + ".fits"
  blueMap = FitsArchive().load(fitsFile)
except:
  cameraName = "green"
  fitsFile = dataDir + "map-" + str(obsids[0]) +  "-" + str(obsids[1]) + "-" + cameraName + ".fits"
  blueMap = FitsArchive().load(fitsFile)

"""
  Modifiy the metadata
"""
meta = blueMap.meta
meta.remove("obsid")

for i in range(len(obsids)):
  meta.set("obsid00"+str(i+1), LongParameter(obsids[i]))

meta.set("associatedObsids", StringParameter(str(borderObsids)[1:-1], "Other associated obsids that contributed to the map"))

"""
 Make sure the ra and dec metaData information is correct
"""
blueMap = centerRaDecMetaData(blueMap)

"""
  Save the map
"""
fitsFile = dataDir + str(obsids[0]) +  "-" + str(obsids[1]) + "_JScanamHPDP_" + cameraName + ".fits"
FitsArchive().save(fitsFile, blueMap)