for i in range(len(band_run_list)) ] QSO_im_list, err_map_list, PSF_list = [], [], [] zp_list = [] qso_center_list, frm_c_RA_DEC_list = [], [] #%% for i in range(len(band_seq)): # The pixel scale is all 0.168 if len(glob.glob(image_folder + filename_list[i])) == 0: print(filename_list[i] + " DOES NOT EXIST!!!") image_not_exit.write(filename_list[i] + ' DOES NOT EXIST!!!' + '\n') QSO_im, err_map, PSF, pix_scale, zp, qso_center, fr_c_RA_DEC = [], [], [], [], [], [], [] run_list.remove(i) else: QSO_im, err_map, PSF, pix_scale, zp, qso_center, fr_c_RA_DEC = gen_fit_id( image_folder, image_RA, image_DEC, filename_list[i], cut_frame=120) QSO_im_list.append(QSO_im) err_map_list.append(err_map) qso_center_list.append(qso_center) frm_c_RA_DEC_list.append(fr_c_RA_DEC) if len(PSF) != 0 and PSF.shape[0] != PSF.shape[1]: cut = int((PSF.shape[0] - PSF.shape[1]) / 2) if cut > 0: PSF = PSF[cut:-cut, :] elif cut < 0: PSF = PSF[:, -cut:cut] PSF /= PSF.sum() if PSF.shape[0] != PSF.shape[1]: raise ValueError("PSF shape is not a square.") PSF_list.append(PSF) zp_list.append(zp)
filename_list = [QSO_id+'_HSC-{0}.fits'.format(band_seq[i]) for i in range(5)] QSO_im_list, err_map_list, QSO_bkg_list, PSF_list=[], [], [], [] zp_list = [] qso_center_list, frm_c_RA_DEC_list = [], [] for i in range(len(band_seq)): # The pixel scale is all 0.168 if len(glob.glob(image_folder+filename_list[i])) == 0: image_not_exit = open(filename_image_not_exit,"r+") image_not_exit.read() print filename_list[i] + " DOES NOT EXIST!!!" image_not_exit.write(filename_list[i] + ' DOES NOT EXIST!!!'+ '\n') QSO_im, err_map, QSO_bkg, PSF, pix_scale, zp, qso_fr_center, fr_c_RA_DEC = [], [], [], [], [], [], [], [] run_list.remove(i) else: QSO_im, err_map, QSO_bkg, PSF, pix_scale, zp, qso_fr_center, fr_c_RA_DEC = gen_fit_id(image_folder, QSO_RA, QSO_DEC, filename_list[i],cut_frame=120, subbkl=sub_bkg) QSO_im_list.append(QSO_im) err_map_list.append(err_map) QSO_bkg_list.append(QSO_bkg) qso_center_list.append(qso_fr_center) frm_c_RA_DEC_list.append(fr_c_RA_DEC) if len(PSF) != 0 and PSF.shape[0] != PSF.shape[1]: cut = ((PSF.shape[0] - PSF.shape[1])/2) if cut>0: PSF = PSF[cut:-cut,:] elif cut<0: PSF = PSF[:,-cut:cut] PSF /= PSF.sum() if PSF.shape[0] != PSF.shape[1]: raise ValueError("PSF shape is not a square.") PSF_list.append(PSF)
band_run_list = [0] #run I band first run_list = copy.deepcopy(band_run_list) band_seq = ['I'] filename_list = [ galaxy_ID + '_HSC-{0}.fits'.format(band_seq[i]) for i in range(1) ] galaxy_im_list, err_map_list, PSF_list = [], [], [] zp_list = [] galaxy_center_list, frm_c_RA_DEC_list = [], [] for i in range(len(band_seq)): # The pixel scale is all 0.168 galaxy_im, err_map, PSF, pix_scale, zp, galaxy_fr_center, fr_c_RA_DEC = gen_fit_id( image_folder, galaxy_RA, galaxy_DEC, filename_list[i], cut_frame=120, pltshow=0, subbkl='y') galaxy_im_list.append(galaxy_im) err_map_list.append(err_map) galaxy_center_list.append(galaxy_fr_center) frm_c_RA_DEC_list.append(fr_c_RA_DEC) if len(PSF) != 0 and PSF.shape[0] != PSF.shape[1]: cut = ((PSF.shape[0] - PSF.shape[1]) / 2) if cut > 0: PSF = PSF[cut:-cut, :] elif cut < 0: PSF = PSF[:, -cut:cut] PSF /= PSF.sum() if PSF.shape[0] != PSF.shape[1]:
# #import lenstronomy.Plots.output_plots as plot_mcmc_behaviour fig = plt.figure(figsize=(20, 15)) ax = fig.add_subplot(111) out_plot.plot_mcmc_behaviour(ax, samples_mcmc, param_mcmc, dist_mcmc) plt.show() #%% Plot the image again: band_seq = ['G', 'R', 'I', 'Z', 'Y'] filename = galaxy_ID + '_HSC-{0}.fits'.format(band) sub_bkg = True from flux_profile import galaxy_total_compare galaxy_im, err_map, galaxy_bkg, PSF, pix_scale, zp, qso_fr_center, fr_c_RA_DEC = gen_fit_id( image_folder, galaxy_RA, galaxy_DEC, filename, cut_frame=120, subbkl=sub_bkg) ct = (len(galaxy_im) - len(image_host[0])) / 2 # If want to cut to 61, QSO_im[ct:-ct,ct:-ct] if len(image_host) == 1: host = image_host[0] label = ['data', 'QSO', 'host', 'model', 'normalized residual'] elif len(image_host) > 1: host = np.zeros_like(image_host[0]) for i in range(len(image_host)): host += image_host[i] label = [ 'data', 'QSO', 'host as {0} components'.format(i + 1), 'model',