САС и ПСФ в Галсим

Question

Настроить

• MacOS 10.12.6
• Python 2.7
• Галсим 1.4.4

---

Цель

Я хочу вставить два SED (Balge и Disk) и использовать их для создания WFIRST PSF.

Пока что я могу импортировать SED с помощью galsim.SED () и создать PSF с помощью wfirst.getPSF()

---

проблема

Применяя wfirst.getPSF(), я не могу принять во внимание желаемый SED.

---

Я даже попробовал:

Я также пытался использовать galsim.Bandpass()

---

Больше деталей:

Основываясь на рецепте, приведенном в примере (#13), можно создать PSF с помощью wfirst.getPSF(), а затем свернуть его в SED.

Я следовал этой рутине:

PSFs = wfirst.getPSF(SCAs=use_SCA, approximate_struts=True,_waves=10,logger=logger)
point = galsim.Gaussian(sigma=1.e-8, flux=1.)
star_sed = galsim.SED(lambda x:1, 'nm', 'flambda').withFlux(1.,filter_)  
star = galsim.Convolve(point*star_sed, PSF)

Мне было интересно, есть ли вариант, в котором мы можем принять во внимание SED, когда мы хотим сделать PSF.

-Спасибо

Answer 1

Ключевым моментом путаницы является то, что у PSF нет SED; только астрономические объекты, такие как звезды и галактики, имеют САС. Процесс, на который вы указали в demo13.py, является правильным способом включения SED: вы присоединяете его к рассматриваемому астрономическому объекту (в данном случае к звезде, но вы также можете назначить SED галактике или назначить различные САС для разделения компонентов галактики).

Так что, если у вас был ахроматический galsim.GSObjects называется bulge а также disk для двух компонентов галактики и отдельных SED для каждого (bulge_sed а также disk_sed), хроматический WFIRST PSF называется psfи галсим. filterтогда вы просто делаете:

galaxy = bulge*bulge_sed + disk*disk_sed
object = galsim.Convolve(galaxy, psf)
object.drawImage(bandpass=filter, scale=wfirst.pixel_scale)

PSF для выпуклости и диска будет отличаться, потому что PSF является хроматическим, и вы дали разным SED для выпуклости и диска разные значения, поэтому это должно охватывать описанный вами вариант использования. См. Demo12.py для большего количества примеров того, как использовать хроматическую функциональность, о которой вы спрашиваете (особенно пример C в этой демонстрации относится к вашему вопросу).