this_gal = this_gal.dilate(hlr)
this_gal = this_gal.shear(e1=xz[p] * e10[l],
e2=xz[p] * e20[l]).shear(g1=g10[i],
g2=g20[i])
final = galsim.Convolve([this_gal, psf])
print(rng())
image = galsim.ImageF(nx, ny, scale=pixel_scale)
fft_image = image[galsim.BoundsI(1, nx, 1, ny)]
final.drawImage(fft_image, method='fft')
sky_level_pixel = sky_level * pixel_scale**2
fft_image.addNoise(
galsim.PoissonNoise(rng,
sky_level=sky_level_pixel)) #add noise
#image
rst = hsm.EstimateShear(image,
final_epsf_image,
shear_est='BJ')
BJe1array[j].append(rst.corrected_e1)
BJe2array[j].append(rst.corrected_e2)
rst = hsm.EstimateShear(image,
final_epsf_image,
shear_est='REGAUSS')
RGe1array[j].append(rst.corrected_e1)
RGe2array[j].append(rst.corrected_e2)
e1, e2 = elli(image.array)
Me1array[j].append(e1)
Me2array[j].append(e2)
BJe1 = []
BJe2 = []
RGe1 = []
RGe2 = []
e1,e2 = elli(galarray1[i]) gal1e1arr.append(e1) gal1e2arr.append(e2) e1,e2 = elli(galarray2[i]) gal2e1arr.append(e1) gal2e2arr.append(e2) e1,e2 = elli(galarray3[i]) gal3e1arr.append(e1) gal3e2arr.append(e2) e1,e2 = elli(galarray4[i]) gal4e1arr.append(e1) gal4e2arr.append(e2) rst = hsm.EstimateShear(galarray1Sd[i],final_epsf_image,shear_est='REGAUSS') rggal1e1.append(rst.corrected_e1) rggal1e2.append(rst.corrected_e2) rst = hsm.EstimateShear(galarray1Sd[i],final_epsf_image,shear_est='BJ') bjgal1e1.append(rst.corrected_e1) bjgal1e2.append(rst.corrected_e2) rst = hsm.EstimateShear(galarray2Sd[i],final_epsf_image,shear_est='REGAUSS') rggal2e1.append(rst.corrected_e1) rggal2e2.append(rst.corrected_e2) rst = hsm.EstimateShear(galarray2Sd[i],final_epsf_image,shear_est='BJ') bjgal2e1.append(rst.corrected_e1) bjgal2e2.append(rst.corrected_e2)
