def blank_uvcal_from_uvdata(uvdata): """initialize UVCal object with same times, antennas, and frequencies as uvdata. Parameters ---------- uvdata: UVData object UVData object that you wish to generate a blanck UVCal object from. Returns ------- uvcal: UVCal object UVCal object with all flags set to False and all gains set to unity with same antennas, freqs, jones, and times as uvdata. """ uvcal = UVCal() uvcal.Nfreqs = uvdata.Nfreqs uvcal.Njones = uvdata.Npols uvcal.Ntimes = uvdata.Ntimes uvcal.Nspws = uvdata.Nspws uvcal.history = "" uvcal.Nspws = uvdata.Nspws uvcal.telescope_name = uvdata.telescope_name uvcal.telescope_location = uvdata.telescope_location uvcal.Nants_data = uvdata.Nants_data uvcal.Nants_telescope = uvdata.Nants_telescope uvcal.ant_array = np.asarray(list(set(uvdata.ant_1_array).union(set(uvdata.ant_2_array)))) uvcal.antenna_names = uvdata.antenna_names uvcal.antenna_numbers = uvdata.antenna_numbers uvcal.antenna_positions = uvdata.antenna_positions uvcal.spw_array = uvdata.spw_array uvcal.freq_array = uvdata.freq_array uvcal.jones_array = uvdata.polarization_array uvcal.time_array = np.unique(uvdata.time_array) uvcal.integration_time = np.mean(uvdata.integration_time) uvcal.lst_array = np.unique(uvdata.lst_array) uvcal.gain_convention = "divide" # always use divide for this package. uvcal.flag_array = np.zeros( (uvcal.Nants_data, uvcal.Nspws, uvcal.Nfreqs, uvcal.Ntimes, uvcal.Njones), dtype=np.bool, ) uvcal.quality_array = np.zeros_like(uvcal.flag_array, dtype=np.float64) uvcal.x_orientation = uvdata.x_orientation uvcal.gain_array = np.ones_like(uvcal.flag_array, dtype=np.complex128) uvcal.cal_style = "redundant" uvcal.cal_type = "gain" uvcal.time_range = ( uvcal.time_array.min() - uvcal.integration_time / 2.0, uvcal.time_array.max() + uvcal.integration_time / 2.0, ) uvcal.channel_width = np.median(np.diff(uvcal.freq_array)) return uvcal
pass else: setattr(uvc, var, attr) # Change what's left manually uvc.cal_style = 'redundant' uvc.cal_type = 'gain' uvc.gain_convention = 'multiply' command_given = sys.argv[1:] command_given = ["'" + x + "'" if ' ' in x else x for x in command_given] command_given = ' '.join(command_given) script_path = os.path.dirname(os.path.abspath(__file__)) git_hash = sub.check_output(['git', '-C', script_path, 'rev-parse', 'HEAD']).strip().decode('UTF-8') uvc.history = 'Created using mkgains.py\nCommand run: mkgains.py %s\nmkgains.py Git Hash: %s' % ( command_given, git_hash) uvc.Njones = uvd.Npols uvc.jones_array = uvd.polarization_array uvc.ant_array = np.arange(uvc.Nants_data) uvc.time_range = [ uvc.time_array[0] - (uvc.integration_time / 172800.), uvc.time_array[-1] + (uvc.integration_time / 172800.) ] uvc.quality_array = np.zeros( (uvc.Nants_data, uvc.Nspws, uvc.Nfreqs, uvc.Ntimes, uvc.Njones), dtype='float64') uvc.flag_array = np.zeros( (uvc.Nants_data, uvc.Nspws, uvc.Nfreqs, uvc.Ntimes, uvc.Njones), dtype='bool') uvc.time_array = np.unique(uvd.time_array)
gains[index524, :, poi_i, 1] = y524 datafile.close() ### cal = UVCal() cal.cal_type = 'gain' cal.set_gain() cal.Nfreqs = Nfreqs cal.Njones = Njones cal.Ntimes = Ntimes # # Change the history comment to list field, freq range name, instrument, averaging sample set, pointing JD reference, # calibration catalogs, and whatever else is important. # cal.history = 'EXAMPLE HISTORY, PLEASE CHANGE: EoR0 highband per frequency, per pointing, per polarization bandpass for MWA, averaged per cable over Season 1 using an early version of KGS. Pointing JD is referenced from Aug 23,2013.' # cal.Nspws = 1 cal.freq_array = freq_array.reshape(cal.Nspws, -1) cal.freq_range = [freq_array[0], freq_array[-1]] # valid frequencies for solutions. cal.channel_width = np.diff(freq_array)[0] cal.jones_array = jones_array cal.time_array = time_array # # Pointing integration time # cal.integration_time = 1800. # cal.gain_convention = 'divide' # Use this operation to apply gain solution. cal.x_orientation = 'east' # orientation of 1st jones parameter. #