err.print_progress_message("Writing mask cube", t0) writemask.writeMask(mask, dict_Header, Parameters, outputMaskCube, Parameters["writeCat"]["compress"], Parameters["writeCat"]["overwrite"]) if Parameters["pipeline"]["trackMemory"]: print_memory_usage(t0) # ------------------------ # ---- STORE CUBELETS ---- # ------------------------ if Parameters["steps"]["doCubelets"] and object_array_exists: err.print_progress_message("Writing cubelets", t0) objects = np.array(objects) cathead = np.array(catParNames) cubelets.writeSubcube(np_Cube, dict_Header, mask, objects, cathead, outroot, outputCubeletsDir, Parameters["writeCat"]["compress"], Parameters["writeCat"]["overwrite"]) if Parameters["pipeline"]["trackMemory"]: print_memory_usage(t0) # ---------------------------- # ---- MAKE MOM0 and MOM1 ---- # ---------------------------- if (Parameters["steps"]["doMom0"] or Parameters["steps"]["doMom1"]): err.print_progress_message("Writing moment maps", t0) debug = 0 write_mom = [ Parameters["steps"]["doMom0"], Parameters["steps"]["doMom1"], False ] writemoment2.writeMoments(np_Cube, mask, outroot, debug, dict_Header, Parameters["writeCat"]["compress"], write_mom[0],
# debug = 0 # writemoment2.writeMoment1(np_Cube, mask, outroot, debug, dict_Header, mom0_Image,Parameters['writeCat']['compress']) # ------------------------ # ---- STORE CUBELETS ---- # ------------------------ if Parameters['steps']['doCubelets'] and Parameters['steps']['doMerge'] and NRdet: print "\n--- %.3f seconds since start"%(time()-t0) print "\n--- SoFiA: Writing cubelets ---" sys.stdout.flush() objects = np.array(objects) cathead = np.array(catParNames) cubelets.writeSubcube(np_Cube, dict_Header, mask, objects, cathead, outroot,Parameters['writeCat']['compress']) # --------------------------------------------------- # ---- APPEND PARAMETER VALUES IN PHYSICAL UNITS ---- # --------------------------------------------------- if Parameters['steps']['doMerge'] and NRdet: print "\n--- %.3f seconds since start"%(time()-t0) print "\n--- SoFiA: Adding WCS position to catalogue ---" sys.stdout.flush() objects, catParNames, catParFormt, catParUnits = wcs_coordinates.add_wcs_coordinates(objects,catParNames,catParFormt,catParUnits,Parameters)
if Parameters["steps"]["doWriteMask"]: err.print_progress_message("Writing mask cube", t0) writemask.writeMask(mask, dict_Header, Parameters, outputMaskCube, Parameters["writeCat"]["compress"], Parameters["writeCat"]["overwrite"]) if Parameters["pipeline"]["trackMemory"]: print_memory_usage(t0) # ------------------------ # ---- STORE CUBELETS ---- # ------------------------ if Parameters["steps"]["doCubelets"] and object_array_exists: err.print_progress_message("Writing cubelets", t0) objects = np.array(objects) cathead = np.array(catParNames) cubelets.writeSubcube(np_Cube, dict_Header, mask, objects, cathead, outroot, outputCubeletsDir, Parameters["writeCat"]["compress"], Parameters["writeCat"]["overwrite"]) if Parameters["pipeline"]["trackMemory"]: print_memory_usage(t0) # ---------------------------- # ---- MAKE MOM0 and MOM1 ---- # ---------------------------- if (Parameters["steps"]["doMom0"] or Parameters["steps"]["doMom1"]): err.print_progress_message("Writing moment maps", t0) debug = 0 write_mom = [Parameters["steps"]["doMom0"], Parameters["steps"]["doMom1"], False] writemoment2.writeMoments(np_Cube, mask, outroot, debug, dict_Header, Parameters["writeCat"]["compress"], write_mom[0], write_mom[1], Parameters["writeCat"]["overwrite"]) if Parameters["pipeline"]["trackMemory"]: print_memory_usage(t0) #writemoment.writeMoments(np_Cube, mask, outroot, debug, dict_Header, Parameters["writeCat"]["compress"], write_mom, Parameters["writeCat"]["overwrite"])
print "\n--- SoFiA: Writing moment-1 map ---" sys.stdout.flush() debug = 0 writemoment.writeMoment1(np_Cube, mask, outroot, debug, dict_Header, mom0_Image) # ------------------------ # ---- STORE CUBELETS ---- # ------------------------ if Parameters["steps"]["doCubelets"] and Parameters["steps"]["doMerge"] and NRdet: print "\n--- SoFiA: Writing cubelets ---" sys.stdout.flush() objects = np.array(objects) cathead = np.array(catParNames) cubelets.writeSubcube(np_Cube, dict_Header, mask, objects, cathead, outroot) # --------------------------------------------------- # ---- APPEND PARAMETER VALUES IN PHYSICAL UNITS ---- # --------------------------------------------------- if Parameters["steps"]["doMerge"] and NRdet: print "\n--- SoFiA: Adding WCS position to catalog ---" sys.stdout.flush() objects, catParNames, catParFormt, catParUnits = wcs_coordinates.add_wcs_coordinates( objects, catParNames, catParFormt, catParUnits, Parameters ) # --------------------