如何在模块化星系图像模拟工具包Galsim中设置随机种子

Question

如何为如下函数设置随机种子:

image = galsim.imageD(x_lim,y_lim,scale)
image_gal = gal_obj.drawShoot(image=image)

我需要一个确定性的种子，但我不知道如何设置它。

Answer 1

在绘图时使用"rng"关键字参数。例如，你可以输入

image = galsim.ImageD(x_lim, y_lim, scale=scale)
image_gal = gal_obj.drawShoot(image=image, rng=galsim.BaseDeviate(1234567))

如果你想使用种子1234567进行确定性光子拍摄。

在一个更小的注意事项，使用drawShoot仍然是支持的，但不鼓励GalSim v1.1。相反，我们鼓励使用统一的绘图方法，允许光子拍摄和DFT渲染。在这种情况下，代码片段的最后一行将变成

image_gal = gal_obj.drawImage(image=image, rng=galsim.BaseDeviate(1234567), method='phot')

(也就是说，假设您已经安装了GalSim v1.1)