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()
#
