在 Galsim 中破译有关同时拟合星系的错误答案

【问题标题】：Deciphering errors regarding simultaneous fitting to galaxies in Galsim在 Galsim 中破译有关同时拟合星系的错误
【发布时间】：2014-08-22 17:52:54
【问题描述】：

我正在拟合 galsim 中的星系并收到以下错误：

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 之间以进行拟合参数化。我怀疑我应该限制我的椭圆度的大小（根据单位圆来限制），而不是单位平方。

谢谢！

【问题讨论】：

标签： python galsim

【解决方案1】：

您是正确的，椭圆度必须在单位圆内有界，因此每个分量上的单独边界将无法完成这项工作。这就是请求剪切超过 1 的 ValueError 的原因。

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

【讨论】：

我收到了一个更新，不根据幅度和极角参数化，我仍然可以使用上面提到的幅度约束进行拟合。我目前无法重现内存问题，但我可以收到一封电子邮件来发送脚本吗？
抱歉回复慢；直到现在我才注意到这个评论。您是说现在配件可以正常工作并且没有内存错误吗？在这种情况下，发送脚本毫无意义，因为我的想法是我可能会使用它来尝试诊断内存错误。