count_uniq_bll = count_uniq_bll[sortind]

total_baselines = bl_length.size
baseline_bin_indices = range(0, int(NP.ceil(1.0*total_baselines/baseline_chunk_size)+1)*baseline_chunk_size, baseline_chunk_size)
if n_bl_chunks is None:
    n_bl_chunks = int(NP.ceil(1.0*total_baselines/baseline_chunk_size))

bl_chunk = range(len(baseline_bin_indices)-1)
bl_chunk = bl_chunk[:n_bl_chunks]
bl = bl[:min(baseline_bin_indices[n_bl_chunks], total_baselines),:]
bl_length = bl_length[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
bl_id = bl_id[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
    
n_channels = nchan

frac_width = DSP.window_N2width(n_window=None, shape=bpass_shape, fftpow=fftpow, area_normalize=True, power_normalize=False)
# window = NP.sqrt(frac_width * n_channels) * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=False, power_normalize=True)
PDB.set_trace()
# window = n_channels * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
window = n_channels * DSP.window_fftpow(n_channels, shape=bpass_shape, fftpow=fftpow, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
bw = n_channels * freq_resolution
bandpass_str = '{0:0d}x{1:.1f}_kHz'.format(nchan, freq_resolution/1e3)

pfb_instr = ''
pfb_outstr = ''
if pfb_method is None:
    pfb_instr = '_no_pfb'
    pfb_outstr = 'no_pfb_'

snapshot_type_str = ''
if avg_drifts:
コード例 #2
0
    count_uniq_bll = count_uniq_bll[sortind]

total_baselines = bl_length.size
baseline_bin_indices = range(0, int(NP.ceil(1.0*total_baselines/baseline_chunk_size)+1)*baseline_chunk_size, baseline_chunk_size)
if n_bl_chunks is None:
    n_bl_chunks = int(NP.ceil(1.0*total_baselines/baseline_chunk_size))

bl_chunk = range(len(baseline_bin_indices)-1)
bl_chunk = bl_chunk[:n_bl_chunks]
bl = bl[:min(baseline_bin_indices[n_bl_chunks], total_baselines),:]
bl_length = bl_length[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
bl_id = bl_id[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
    
n_channels = nchan

frac_width = DSP.window_N2width(shape=bpass_shape, area_normalize=False, power_normalize=True)
# window = NP.sqrt(frac_width * n_channels) * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=False, power_normalize=True)
window = n_channels * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=True, power_normalize=False) 
bw = n_channels * freq_resolution
bandpass_str = '{0:0d}x{1:.1f}_kHz'.format(nchan, freq_resolution/1e3)

pfb_instr = ''
pfb_outstr = ''
if pfb_method is None:
    pfb_instr = '_no_pfb'
    pfb_outstr = 'no_pfb_'

snapshot_type_str = ''
if avg_drifts:
    snapshot_type_str = 'drift_averaged_'
if beam_switch:
コード例 #3
0
    int(NP.ceil(1.0 * total_baselines / baseline_chunk_size) + 1) *
    baseline_chunk_size, baseline_chunk_size)
if n_bl_chunks is None:
    n_bl_chunks = int(NP.ceil(1.0 * total_baselines / baseline_chunk_size))

bl_chunk = range(len(baseline_bin_indices) - 1)
bl_chunk = bl_chunk[:n_bl_chunks]
bl = bl[:min(baseline_bin_indices[n_bl_chunks], total_baselines), :]
bl_length = bl_length[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
bl_id = bl_id[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]

n_channels = nchan

frac_width = DSP.window_N2width(n_window=None,
                                shape=bpass_shape,
                                fftpow=fftpow,
                                area_normalize=True,
                                power_normalize=False)
# window = NP.sqrt(frac_width * n_channels) * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=False, power_normalize=True)
PDB.set_trace()
# window = n_channels * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
window = n_channels * DSP.window_fftpow(n_channels,
                                        shape=bpass_shape,
                                        fftpow=fftpow,
                                        pad_width=0,
                                        centering=True,
                                        area_normalize=True,
                                        power_normalize=False)
bw = n_channels * freq_resolution
bandpass_str = '{0:0d}x{1:.1f}_kHz'.format(nchan, freq_resolution / 1e3)