outfits = ufile.replace('.fits', '_old.fits') outfitsname = os.path.split(outfits)[-1] # log that we are saving rotated image WLOG(p, '', 'Saving Rotated Image in ' + outfitsname) # add keys from original header file hdict = spirouImage.CopyOriginalKeys(hdr) # write to file p = spirouImage.WriteImage(p, outfits, newimage, hdict) # ---------------------------------------------------------------------- # End Message # ---------------------------------------------------------------------- wmsg = 'Recipe {0} has been successfully completed' WLOG(p, 'info', wmsg.format(p['PROGRAM'])) # return a copy of locally defined variables in the memory return dict(locals()) # ============================================================================= # Start of code # ============================================================================= if __name__ == "__main__": # run main with no arguments (get from command line - sys.argv) ll = main() # exit message if in debug mode spirouStartup.Exit(ll, has_plots=False) # ============================================================================= # End of code # =============================================================================
# add berv values hdict = spirouImage.AddKey(p, hdict, p['KW_BERV'], value=loc['BERV']) hdict = spirouImage.AddKey(p, hdict, p['KW_BJD'], value=loc['BJD']) hdict = spirouImage.AddKey(p, hdict, p['KW_BERV_MAX'], value=loc['BERV_MAX']) # write image and add header keys (via hdict) p = spirouImage.WriteImage(p, corfile, data, hdict) # ---------------------------------------------------------------------- # End Message # ---------------------------------------------------------------------- p = spirouStartup.End(p) # return a copy of locally defined variables in the memory return dict(locals()) # ============================================================================= # Start of code # ============================================================================= if __name__ == "__main__": # run main with no arguments (get from command line - sys.argv) ll = main() # exit message if in debug mode spirouStartup.Exit(ll) # ============================================================================= # End of code # =============================================================================
centpart = data2[:, col - slice0:col + slice0] # centpart = data2[col - slice:col + slice,:] # weights = np.where((centpart < satseuil) & (centpart > 0), 1, 0.0001) # y = np.average(centpart, axis=1, weights=weights) ## weighted average y = np.nanmedian(centpart, axis=1) # y=average(centpart,axis=1,weights=where((centpart>0),1,0.0001)) # ## weighted average plt.plot(np.arange(ny), y) # ---------------------------------------------------------------------- # End Message # ---------------------------------------------------------------------- p = spirouStartup.End(p, outputs=None) # return a copy of locally defined variables in the memory return dict(locals()) # ============================================================================= # Start of code # ============================================================================= if __name__ == "__main__": # run main with no arguments (get from command line - sys.argv) ll = main() # exit message if in debug mode spirouStartup.Exit(ll, has_plots=True) # ============================================================================= # End of code # =============================================================================