나는 최근에 python을 사용하여 healpy에서 synfast
함수로 재생했습니다. 특히 나는 최대 모드가 계산되도록 설정하는 매개 변수 lmax (기본값 : lmax = 3 * nside-1)으로 연주했습니다. 또한 당신이 볼 수있는 것처럼 ud_grade
함수를 사용했습니다.healpy의 기능 synfast에서 증가 된 lmax의 편향된 스펙트럼
내가 온도계 TT와 편광 EE와 BB에서 균일 한 평평한 스펙트럼으로부터 synfast
으로 매우 간단한지도를 생성 한 몇 가지 사례를 살펴 보았습니다. 그런 다음 anafast
을 사용하여 해당지도에서 스펙트럼을 계산합니다. 다음 플롯에서 볼 수 있듯이 모든 경우에 초기 입력 스펙트럼을 복구하지 못하는 것으로 나타났습니다.
synfast
의
pixwin 인수 설정 기능을 사용
pixwin(nside=128)
anafast
의 출력 신호를 보정.
- 점선 회색 : 입력 스펙트럼 모델
- 블루 (일반 기본 경우)
synfast(nside=128, lmax=3*128-1)
=>anafast(lmax=3*128-1)
- 녹색 :
synfast(nside=256, lmax=3*256-1)
=>ud_grade(256->128)
=>anafast(lmax=3*128-1)
- 마젠타 :
synfast(nside=256, lmax=3*128-1)
=>ud_grade(256->128)
=>anafast(lmax=3*128-1)
- 빨강 :
synfast(nside=128, lmax=4*128-1)
=>anafast(lmax=3*128-1)
- 시안은 : 당신이 볼 수 있듯이
synfast(nside=128, lmax=10*128-1)
=>anafast(lmax=3*128-1)
Plot : spectra estimations from different kind of map using syfast, ud_grade and anafast
- , 나는 대부분의 경우에 편견 스펙트럼을 얻는다.
- 정상적인 경우 (파란색)는 높은 모드의 경우 약간 아래쪽으로 기울어 져 있습니다. 이것은 healpy가> 2 * nside-1 모드에 대한 스펙트럼 추정을 잘 처리하지 못하기 때문에 명백하게 정상입니다.
- 결과는 온도와 편광 스펙트럼 (TT, EE 및 BB)에서 동일합니다.
- 보시다시피, lmax = 4 * 128-1 및 lmax = 10 * 128-1의 경우는 동일한 결과를 제공합니다. 이것은 lmax을 4 * 128-1보다 높게 설정할 수없고 lmax을 요구하더라도 lmax = 4 * 128-1을 설정하는
synfast
에서 비롯됩니다.
누군가이 동작에 대한 설명이 있습니까?
초기 입력 스펙트럼을 복구하는 방법은 무엇입니까?
미리 ud_grade의 예상되는 영향이다