def draw_otf(cube, flag=None, ii=None, rms=None, suptitle='', show=True, *args, **kwargs): import pylab import analyse if ii is None: ii = analyse.make_2d_map(cube, flag) if rms is None: rms = analyse.make_2d_map(cube, flag, 'rms') fig = pylab.figure(figsize=(10,8)) fig.suptitle(suptitle) plot = custom_draw_map(figure=fig, show=False) plot(ii, subplot=221, title='mom0', *args, **kwargs) plot(rms, subplot=222, title='rms', yaxis_label=False, *args, **kwargs) draw_otf_spectrum(cube, figure=fig, subplot=212, title='spectra', *args, **kwargs) if show: pylab.show() return fig
def easy_analyse(fitspath, output_dir, flag=None, plot=False, save=False): import analyse print('make_cube: %s'%(fitspath.split('/')[-1])) raw_data = analyse.loadfits(fitspath) cw_data = analyse.makespec(raw_data) if flag is None: fitted_data, flag = analyse.basefit(cw_data) else: flag = analyse.loadfits(flag) fitted_data = analyse.basefit_flag(cw_data, None, flag) pass convolved_data = analyse.convolve(fitted_data, 2) fitted_data = analyse.basefit_flag(convolved_data, None, flag) ii = analyse.make_2d_map(fitted_data, flag) rms = analyse.make_2d_map(fitted_data, flag, 'rms') if save: savepath = output_dir + fitspath.split('/')[-1].split('.fits')[0] analyse.savefits(fitted_data, savepath+'.qlook.data.fits', clobber=True) analyse.savefits(flag, savepath+'.qlook.flag.fits', clobber=True) analyse.savefits(ii, savepath+'.qlook.ii.fits', clobber=True) analyse.savefits(rms, savepath+'.qlook.rms.fits', clobber=True) pass if plot: savepath = output_dir + fitspath.split('/')[-1].split('.fits')[0] isotope = fitspath.split('_')[-2] plot = analyse.custom_draw_map(figure=(8,8), xspacing=0.2, yspacing=0.2, tick_labels_size=12, colorber_font_size=11, show=False) plot(ii, title='mom0: '+isotope).save(savepath+'.qlook.ii.png') plot(rms, title='rms: '+isotope).save(savepath+'.qlook.rms.png') analyse.draw_otf_spectrum(fitted_data, figure=(12,5), title='spectra: '+isotope, show=False).savefig(savepath+'.qlook.spectrum.png') pass raw_data.data = None cw_data.data = None return fitted_data, flag
def generate_rms(self): import analyse self.rms = analyse.make_2d_map(self.cube, self.flag, 'rms') return
def generate_ii(self): import analyse self.ii = analyse.make_2d_map(self.cube, self.flag) return