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)