fig = img.axs[1][0].imshow(shear_field_pdf, vmin=g_min, vmax=g_max,cmap="jet") img.figure.colorbar(fig, ax=img.axs[1][0]) fig = img.axs[1][1].imshow(shear_field-shear_field_pdf, cmap="jet") img.figure.colorbar(fig, ax=img.axs[1][1]) yfmt = ["%.1f" % (8 - 4 * k) for k in range(5)] xfmt = ["%.1f" % (-8 + 4 * k) for k in range(5)] # img.set_ticklabel(1, 0, 0, len(xfmt), xfmt) # img.set_ticklabel(1, 0, 1, len(xfmt), xfmt) # # img.set_ticklabel(1, 1, 0, len(xfmt), xfmt) # img.set_ticklabel(1, 1, 1, len(xfmt), xfmt) titles = ["Input", "MCMC Fit", "True surface", "True - Fitted"] for i in range(2): for j in range(2): img.set_ticklabel(i, j, 0, len(yfmt), yfmt) img.set_ticklabel(i, j, 1, len(xfmt), xfmt) img.set_label(i, j, 0, "DEC.[arcmin]",size=img.xy_lb_size-5) img.set_label(i, j, 1, "R.A.[arcmin]",size=img.xy_lb_size-5) img.axs[i][j].set_title(titles[i*2+j]) pic_name = parent_path + "pic/%s_mcmc_fit_expo_%d_slope.png"%(shear_cmd, expo_num) pic_name = tool_box.file_name(pic_name) img.save_img(pic_name) img.close_img() print("Fit. %s"%pic_name)
# else: # fig = img.axs[0][i].imshow(show_data[i], cmap="jet") # img.figure.colorbar(fig, ax=img.axs[0][i]) # figs.append(fig) fig = img.axs[0][i].imshow(show_data[i], cmap="jet") img.figure.colorbar(fig, ax=img.axs[0][i]) img.set_label(0, i, 0, "DEC. [arcmin]") img.set_label(0, i, 1, "R.A. [arcmin]") if i < 2: pass fmt = ["%.1f"%(8 - 4*k) for k in range(5)] img.set_ticklabel(0, i, 0, len(fmt), fmt) fmt = ["%.1f"%(-8 + 4*k) for k in range(5)] img.set_ticklabel(0, i, 1, len(fmt), fmt) # img.axs[0][i].set_yticks([int(nx - 1 - k * nx / 5) for k in range(6)]) # img.axs[0][i].set_yticklabels(["%.1f" % (k * nx / 5. * pixel_scale) for k in range(6)]) # img.axs[0][i].set_xticks([int(k * ny / 5) for k in range(6)]) # img.axs[0][i].set_xticklabels(["%.1f" % (k * ny / 5. * pixel_scale) for k in range(6)]) img.axs[0][i].set_title(titles[i]) img.save_img(parent_path + "pic/shear_field_slope_test.png") img.show_img() #