Question

我适合加尔西姆的星系并收到以下错误：

File "/home/luis/Documents/SRC2014/galsim_work/noiseLibrary.py", line 50, in create_galaxy
gal = gal.shear(g1=e1, g2=e2)

File "/usr/lib/python2.7/dist-packages/galsim/base.py", line 572, in shear
shear = galsim.Shear(**kwargs)

File "/usr/lib/python2.7/dist-packages/galsim/shear.py", line 111, in __init__
raise ValueError("Requested shear exceeds 1: %f"%g)

ValueError: Requested shear exceeds 1: 1.007171

我试图通过检查e1和e2的值以及我的程序中通过以下方式创建星系的幅度来解决这个问题：

if e1 > 1:
     print "The e1 component of ellipticity is greater than 1."
     e1 = 0.99
elif e1 < -1:
     print "The e1 component of ellipticity is less than -1."
     e1 = -0.99
if e2 > 1:
     print "The e2 component of ellipticity is greater than 1."
     e2 + 0.99
elif e2 < -1:
     print "The e2 component of ellipticity is less than -1."
     e2 = -0.99

if np.sqrt(e1**2 + e2**2) > 1:
     return create_galaxy(flux, hlr, e1*0.5, e2*0.5, galtype_gal=galtype_gal, sersic_index=sersic_index,
                          psf_flag=psf_flag, psf_type=psf_type, beta=beta, size_psf=size_psf, flux_psf=flux_psf,
                          x_len=x_len, y_len=y_len, scale=scale, method=method,seed=seed)

虽然此解决方案适用于各个组件，但我怀疑递归调用是否会产生以下错误：

File "/usr/lib/python2.7/dist-packages/galsim/base.py", line 1196, in drawImage
image.added_flux = prof.SBProfile.draw(imview.image, gain, wmult)

MemoryError

有关如何解决这个问题的任何建议？请注意，对于我的拟合参数化，我将e1和e2绑定在介于-1和1之间。我怀疑我应该限制椭圆度的大小（根据单位圆限制），而不是单位平方。

谢谢！

Answer 1

你是正确的，椭圆度必须限制在单位圆内，因此每个组件上的单独边界将无法完成工作。这就是关于请求的剪切超过1的ValueError的原因。

但是，在没有更多细节的情况下，很难诊断出内存错误。这是一直发生还是偶尔发生一次？你能给我们一个简短的脚本来重现这个错误，包括所有galaxy和PSF参数的特定值以及你用来制作正在绘制的对象的命令吗？

解读关于同时适应加尔西姆星系的错误

1 个答案: